Changeset 2392 for trunk

Timestamp:

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

Author:

cmathe

Message:

improvedco2clouds_mod: 4.2 => if zq(ccnco2_number) + threshold > 1

move the end if of 4.2

File:

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

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

@@ -491 +491 @@
 !----------------------------------------------------------------------------------------------------------------------!
       ! we trigger crystal growth
-      if (zq(ig,l,igcm_ccnco2_number) * tauscaling(ig) >= threshold) then
+      if (zq(ig,l,igcm_ccnco2_number) * tauscaling(ig) +threshold >= 1) then
 
         call updaterice_microco2(zq(ig,l,igcm_co2_ice), zq(ig,l,igcm_ccnco2_mass), zq(ig,l,igcm_ccnco2_number), &
@@ -528 +528 @@
           zq(ig,l,igcm_co2) = zq(ig,l,igcm_co2) - dMice
         end if
+      end if ! nouveau end if
 !----------------------------------------------------------------------------------------------------------------------!
 ! 4.3 Dust cores releasing if no more co2 ice
 !----------------------------------------------------------------------------------------------------------------------!
         ! On sublime tout
-        if (zq(ig,l,igcm_co2_ice) <= threshold) then
+        if ((zq(ig,l,igcm_co2_ice) <= threshold).or.(zq(ig,l,igcm_ccnco2_number) < 1.)) then
 
           if (co2useh2o) then
@@ -565 +566 @@
           riceco2(ig,l) = 0.
         end if !of if co2_ice < threshold
-      end if ! of if NCCN > 1
+!      end if ! of if NCCN > 1
     end do ! of ig loop
   end do ! of nlayer loop