Changeset 3017

Timestamp:

Jul 27, 2023, 12:23:20 AM (16 months ago)

Author:

jnaar

Message:

Change the way zdq is estimated to fasten (by a lot) the adaptative subtimestep of water ice clouds. JN

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/improvedclouds_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4142 +4142 @@
 into module data there. Also convert lteparams.h into module nlteparams_h.F90.
 And while at it also turn nlthermeq.F and blendrad.F into modules.
+
+== 27/07/2023 == JN
+Optimization for adaptative subtimestep of water-ice clouds : should now be
+much faster. + small cleanup

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

@@ -397 +397 @@
        IF ( zq(ig,l,igcm_ccn_number)*tauscaling(ig).ge. 1.) THEN ! we trigger crystal growth 
 
-     
         call updaterice_micro(zq(ig,l,igcm_h2o_ice),
      &    zq(ig,l,igcm_ccn_mass),zq(ig,l,igcm_ccn_number),
@@ -417 +416 @@
 ccccccc  implicit scheme of mass growth
 c         cste here must be computed at each step
-          cste = 4*pi*rho_ice*microtimestep
-c          print*, 'ig',ig
-c          print*, 'l',l
-c          print*, 'cste',cste
-c          print*, 'seq',seq
-c          print*, 'zqsat',zqsat(ig,l)
-c          print*, 'zq vap',zq(ig,l,igcm_h2o_vap)
-c          print*, 'res',res
-c          print*, 'No',No
-c          print*, 'rice',rice(ig,l)
+        cste = 4*pi*rho_ice*microtimestep
 
         dMice =
@@ -578 +568 @@
           ENDDO !i=1,nq
 
+          IF (cloud_adapt_ts) then
+c           Estimation of how much is actually condensing/subliming
+                zdq=(dMice/spenttime)*(ptimestep-spenttime)
+c           Threshold with powerlaw (sanity check)
+                zdq=min(abs(zdq),
+     &            abs(zpotcond(ig,l)*((ptimestep-spenttime)/ptimestep)))
+                call adapt_imicro(ptimestep,zdq,
+     &             zimicro(ig,l))
+          ENDIF! (cloud_adapt_ts) then
 c         Increment time spent in here
           spenttime=spenttime+microtimestep
           count_micro(ig,l)=count_micro(ig,l)+1
-          IF (cloud_adapt_ts) then
-c           Estimation of how much is actually condensing/subliming
-c                print*, 'zq',zq(ig,l,igcm_h2o_vap)
-c                print*, 'zq0',zq0(ig,l,igcm_h2o_vap)
-c                print*, 'ptimestep',ptimestep
-                zdq=(zq(ig,l,igcm_h2o_vap)-zq0(ig,l,igcm_h2o_vap))
-     &                   *ptimestep/microtimestep
-c           Estimation of how much is actually condensing/subliming
-c                zdq=min(abs(dMice)*ptimestep,zpotcond(ig,l))
-                call adapt_imicro(ptimestep,zdq,
-     &             zimicro(ig,l))
-          ENDIF! (cloud_adapt_ts) then
           ENDDO ! while (.not. ending_ts)
         ENDDO ! of ig loop