3 FUNCTION cooling_casiana (self, temp)
5 class(cooling_t):: self
6 real(kind=dp) cooling_casiana, temp, qq
7 real(kind=dp),
parameter:: &
43 hm=4.98282e-05*4e-28, &
59 real(kind=8),
parameter:: &
62 if (temp .le. t_mc)
then 64 else if (temp .le. te0)
then 66 else if (temp .le. te1)
then 68 else if (temp .le. te2)
then 70 else if (temp .le. te3)
then 72 else if (temp .le. te4)
then 74 else if (temp .le. te5)
then 76 else if (temp .le. te6)
then 78 else if (temp .le. te7)
then 80 else if (temp .le. te8)
then 81 qq = (h7*temp**c7)*2.5e27
82 else if (temp .le. te9)
then 84 else if (temp .le. te10)
then 86 else if (temp .le. te11)
then 88 else if (temp .le. te12)
then 90 else if (temp .le. te13)
then 92 else if (temp .le. te14)
then