Changeset 531

Timestamp:

Feb 15, 2012, 5:22:46 PM (13 years ago)

Author:

tnavarro

Message:

corrected bug

File:

• trunk/LMDZ.MARS/libf/phymars/improvedclouds.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/improvedclouds.F

@@ -170 +170 @@
 
         do i=1,nbin_cld-1
-          rad_cld(i+1)  = rad_cld(i) * CBRT(vrat_cld) !!**(1./3.)
+          rad_cld(i+1)  = rad_cld(i) * vrat_cld**(1./3.)
           vol_cld(i+1)  = vol_cld(i) * vrat_cld
         enddo
         
         do i=1,nbin_cld
-          rb_cld(i+1)= CBRT( (2.*vrat_cld) / (vrat_cld+1.) ) *
+          rb_cld(i+1)= ( (2.*vrat_cld) / (vrat_cld+1.) )**(1./3.) *
      &      rad_cld(i)
           dr_cld(i)  = rb_cld(i+1) - rb_cld(i)
@@ -364 +364 @@
           rhocloud(ig,l) = min(max(rhocloud(ig,l),rho_ice),rho_dust)
           rice(ig,l) =
-     &      CBRT( Mo / No * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
+     &      CBRT( real(Mo)/real(No) * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
 c          nuice(ig,l)=nuice_ref ! used for rad. transfer calculations
           if ((Mo.lt.1.e-20) .or. (No.le.1)) rice(ig,l) = 1.e-8
@@ -484 +484 @@
           rice_out(ig,l)=rice(ig,l)
           rice(ig,l) =
-     &      CBRT( Mo / No * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
+     &      CBRT( real(Mo)/real(No) * 0.75 / pi / rhocloud(ig,l) ) !**(1./3.)
           if ((Mo.lt.1.e-20) .or. (No.le.1)) rice(ig,l) = 1.e-8
           rice_out(ig,l)=rice(ig,l)-rice_out(ig,l)