Changeset 2437 for trunk

Timestamp:

Dec 1, 2020, 3:55:37 PM (4 years ago)

Author:

mvals

Message:

Mars GCM:

• improvedclouds_mod.F: update of the nucleation equation with its analytical resolution
• writediagmicrofi.F: makes it possible to get outputs from the microphysics (call example in watercloud_mod.F)

MV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/improvedclouds_mod.F (modified) (2 diffs)
• libf/phymars/watercloud_mod.F (modified) (1 diff)
• libf/phymars/writediagmicrofi.F (added)

trunk/LMDZ.MARS/README

@@ -3239 +3239 @@
 Add stormdust in the aerosol opacities computed by aeroptical.F90
 + add a description of aeroptical in util/README
+
+== 01/12/2020 == MV
+- improvedclouds_mod.F: update of the nucleation equation with its analytical resolution
+- writediagmicrofi.F: makes it possible to get outputs from the microphysics (call example in watercloud_mod.F)

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

@@ -326 +326 @@
         dM = 0.
         do i = 1, nbin_cld
-          n_aer(i) = n_aer(i)/( 1. + rate(i)*microtimestep)
-          m_aer(i) = m_aer(i)/( 1. + rate(i)*microtimestep)
-          dN       = dN + n_aer(i) * rate(i) * microtimestep
-          dM       = dM + m_aer(i) * rate(i) * microtimestep
+          dN       = dN + n_aer(i)*(exp(-rate(i)*microtimestep)-1.)
+          dM       = dM + m_aer(i)*(exp(-rate(i)*microtimestep)-1.)
         enddo
 
@@ -335 +333 @@
 c       Update Dust particles
         zq(ig,l,igcm_dust_mass)   = 
-     &  zq(ig,l,igcm_dust_mass)   - dM/ tauscaling(ig) !max(tauscaling(ig),1.e-10) 
+     &  zq(ig,l,igcm_dust_mass)   + dM/ tauscaling(ig) !max(tauscaling(ig),1.e-10) 
         zq(ig,l,igcm_dust_number) = 
-     &  zq(ig,l,igcm_dust_number) - dN/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
+     &  zq(ig,l,igcm_dust_number) + dN/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
 c       Update CCNs
         zq(ig,l,igcm_ccn_mass)   = 
-     &  zq(ig,l,igcm_ccn_mass)   + dM/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
+     &  zq(ig,l,igcm_ccn_mass)   - dM/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
         zq(ig,l,igcm_ccn_number) = 
-     &  zq(ig,l,igcm_ccn_number) + dN/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
+     &  zq(ig,l,igcm_ccn_number) - dN/ tauscaling(ig) !max(tauscaling(ig),1.e-10)
         
         ENDIF ! of is satu >1

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

@@ -392 +392 @@
           ENDDO
         ENDIF
-     
+
+!      !! Example of how to use writediagmicrofi useful to 
+!      !! get outputs at each microphysical sub-timestep (better to be used in 1D)
+!            CALL WRITEDIAGMICROFI(ngrid,imicro,microstep,
+!     &       microtimestep,'subpdtcloud',
+!     &      'subpdtcloud','K/s',1,subpdtcloud(:,:))     
  
       ENDDO ! of DO microstep=1,imicro