Changeset 4950 for LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.F90

Timestamp:

May 22, 2024, 3:16:36 PM (6 weeks ago)

Author:

lebasn

Message:

StratAer?: New model version (microphysic, composition routine, code cleaning, new params...)

File:

• LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.F90

@@ -51 +51 @@
   
   !============= NUCLEATION VARS =============
+  ! MOLECULAR ACCOMODATION OF H2SO4 (Raes and Van Dingen)
+  REAL,SAVE    :: ALPH2SO4               ! H2SO4 accommodation  coefficient [condensation/evaporation]
+  !$OMP THREADPRIVATE(ALPH2SO4)
+  
   ! flag to constraint nucleation rate in a lat/pres box
   LOGICAL,SAVE :: flag_nuc_rate_box      ! Nucleation rate limit or not to a lat/pres
@@ -64 +68 @@
   INTEGER,SAVE :: flh2o  ! ds stratemit : flh2o =0 (tr_seri), flh2o=1 (dq) 
   !$OMP THREADPRIVATE(flh2o)
-!  REAL,ALLOCATABLE,SAVE    :: d_q_emiss(:,:)
-!  !$OMP THREADPRIVATE(d_q_emiss)
   
   REAL,ALLOCATABLE,SAVE    :: budg_emi(:,:)            !DIMENSION(klon,n)
@@ -144 +146 @@
   !$OMP THREADPRIVATE(day_emit_roc)
   
+  REAL,ALLOCATABLE,SAVE    :: RRSI(:) ! radius [cm] for each aerosol size
+  REAL,ALLOCATABLE,SAVE    :: Vbin(:) ! volume [m3] for each aerosol size  
+  !$OMP THREADPRIVATE(RRSI, Vbin)
   REAL,SAVE    :: dlat, dlon             ! delta latitude and d longitude of grid in degree
   !$OMP THREADPRIVATE(dlat, dlon)
@@ -153 +158 @@
     USE print_control_mod, ONLY : lunout
     USE mod_phys_lmdz_para, ONLY : is_master
-    USE infotrac_phy, ONLY: id_OCS_strat,id_SO2_strat,id_H2SO4_strat,nbtr_sulgas
+    USE infotrac_phy, ONLY: id_OCS_strat,id_SO2_strat,id_H2SO4_strat,nbtr_sulgas,nbtr_bin
+    USE phys_local_var_mod, ONLY : mdw
+    USE aerophys, ONLY: mdwmin, V_rat
+    USE YOMCST  , ONLY : RPI
+    
+    INTEGER :: it
     
     WRITE(lunout,*) 'IN STRATAER_LOCAL_VAR INIT WELCOME!'
@@ -185 +195 @@
     
     ! nuc init
+    ALPH2SO4 = 0.1
     flag_nuc_rate_box = .FALSE.
     nuclat_min=0  ; nuclat_max=0
@@ -238 +249 @@
     ENDIF ! if master
     
+    !--initialising dry diameters to geometrically spaced mass/volume (see Jacobson 1994)
+    mdw(1)=mdwmin
+    IF (V_rat.LT.1.62) THEN ! compensate for dip in second bin for lower volume ratio
+       mdw(2)=mdw(1)*2.**(1./3.)
+       DO it=3, nbtr_bin
+          mdw(it)=mdw(it-1)*V_rat**(1./3.)
+       ENDDO
+    ELSE
+       DO it=2, nbtr_bin
+          mdw(it)=mdw(it-1)*V_rat**(1./3.)
+       ENDDO
+    ENDIF
+    IF (is_master) WRITE(lunout,*) 'init mdw=', mdw
+    
+    !   compute particle radius RRSI [cm] and volume Vbin [m3] from diameter mdw [m]
+    ALLOCATE(RRSI(nbtr_bin), Vbin(nbtr_bin))
+    
+    DO it=1,nbtr_bin
+       !     [cm]
+       RRSI(it)=mdw(it)/2.*100.
+       !     [m3]
+       Vbin(it)=4.0*RPI*((mdw(it)/2.)**3)/3.0
+    ENDDO
+    
+    IF (is_master) THEN
+       WRITE(lunout,*) 'init RRSI=', RRSI
+       WRITE(lunout,*) 'init Vbin=', Vbin
+    ENDIF
+    
     WRITE(lunout,*) 'IN STRATAER INIT END'