Changeset 2388 for trunk

Timestamp:

Jul 2, 2020, 3:42:33 PM (4 years ago)

Author:

cmathe

Message:

mprovedco2clouds_mod: case satu < 1 removed, debug in 4.2 (write)

File:

• trunk/LMDZ.MARS/libf/phymars/improvedco2clouds_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/improvedco2clouds_mod.F90

@@ -482 +482 @@
           end if ! of if co2useh2o
         end if ! of if No > 1e-30
-
-        ! WIP: cas où satu < 1 et donc si on a de la glace de co2 et que T > Tcond, alors on relache notre glace en vap
-      else
-          if (tcond(ig,l) <= zt(ig,l) .and. zq(ig,l,igcm_co2_ice) >= threshold) then
-            zq(ig,l,igcm_co2) = zq(ig,l,igcm_co2) + zq(ig,l,igcm_co2_ice)
-            zq(ig,l,igcm_co2_ice) = 0.
-          end if
       end if ! of is satu > 1
 !----------------------------------------------------------------------------------------------------------------------!
@@ -499 +492 @@
       ! we trigger crystal growth
       if (zq(ig,l,igcm_ccnco2_number) * tauscaling(ig) + threshold >= 1.) then
+        write(*,*)'yes I am in', zq(ig, l, igcm_ccnco2_number), zq(ig, l, igcm_co2_ice), satu
 
         call updaterice_microco2(zq(ig,l,igcm_co2_ice), zq(ig,l,igcm_ccnco2_mass), zq(ig,l,igcm_ccnco2_number), &
@@ -510 +504 @@
         facteurmax = abs(tcond(ig,l)-zt(ig,l)) * (cpp/lw)
 
-        ! call scheme of microphys. mass growth for CO2
+        ! call scheme of microphys. mass growth for CO2 (evaporation/condensation)
         call massflowrateCO2(pplay(ig,l), zt(ig,l), satu, riceco2(ig,l), mmean(ig,l), Ic_rice)