Changeset 2957

Timestamp:

May 5, 2023, 9:13:28 AM (18 months ago)

Author:

emillour

Message:

Generic PCM:
Fix a buggy behavior concerning H2O aerosol variance; aeroptproperties is not
designed to handle aerosol variance which is not constant, whereas h2o_reffrad
returns a variance which varies (between 0.09 and 0.13) with location and time.
Revert to a simpler setup where the H2O aerosol variance is uniform and set by
the user (nueff_iaero_h2o flag in callphys.def; default value 0.1)
Also added some "intent()" in optci arguments and increased length of string
to store varaible name in writediagfi.
EM

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/optci.F90 (modified) (3 diffs)
• libf/phystd/radii_mod.F90 (modified) (7 diffs)
• libf/phystd/writediagfi.F (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1792 +1792 @@
 While at it also moved "Nmix_co2" and "radfixed" from callkeys_mod to radii_mod
 since they are only used in that module.
+
+== 05/05/2023 == EM
+Fix a buggy behavior concerning H2O aerosol variance; aeroptproperties is not
+designed to handle aerosol variance which is not constant, whereas h2o_reffrad
+returns a variance which varies (between 0.09 and 0.13) with location and time.
+Revert to a simpler setup where the H2O aerosol variance is uniform and set by
+the user (nueff_iaero_h2o flag in callphys.def; default value 0.1)
+Also added some "intent()" in optci arguments and increased length of string
+to store varaible name in writediagfi.

trunk/LMDZ.GENERIC/libf/phystd/optci.F90

@@ -42 +42 @@
 
 
-  real*8 DTAUI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
+  real*8,intent(out) :: DTAUI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
   real*8 DTAUKI(L_LEVELS,L_NSPECTI,L_NGAUSS)
   real*8 TAUI(L_NLEVRAD,L_NSPECTI,L_NGAUSS)
-  real*8 TAUCUMI(L_LEVELS,L_NSPECTI,L_NGAUSS)
-  real*8 PLEV(L_LEVELS)
-  real*8 TLEV(L_LEVELS)
-  real*8 TMID(L_LEVELS), PMID(L_LEVELS)
-  real*8 COSBI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
-  real*8 WBARI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
+  real*8,intent(out) :: TAUCUMI(L_LEVELS,L_NSPECTI,L_NGAUSS)
+  real*8,intent(in) :: PLEV(L_LEVELS)
+  real*8,intent(in) :: TLEV(L_LEVELS) ! not used
+  real*8,intent(in) :: TMID(L_LEVELS)
+  real*8,intent(in) :: PMID(L_LEVELS)
+  real*8,intent(out) :: COSBI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
+  real*8,intent(out) :: WBARI(L_NLAYRAD,L_NSPECTI,L_NGAUSS)
 
   ! for aerosols
-  real*8  QXIAER(L_LEVELS,L_NSPECTI,NAERKIND)
-  real*8  QSIAER(L_LEVELS,L_NSPECTI,NAERKIND)
-  real*8  GIAER(L_LEVELS,L_NSPECTI,NAERKIND)
-  real*8  TAUAERO(L_LEVELS,NAERKIND)
+  real*8,intent(in) ::  QXIAER(L_LEVELS,L_NSPECTI,NAERKIND)
+  real*8,intent(in) ::  QSIAER(L_LEVELS,L_NSPECTI,NAERKIND)
+  real*8,intent(in) ::  GIAER(L_LEVELS,L_NSPECTI,NAERKIND)
+  real*8,intent(in) ::  TAUAERO(L_LEVELS,NAERKIND)
   real*8  TAUAEROLK(L_LEVELS,L_NSPECTI,NAERKIND)
   real*8  TAEROS(L_LEVELS,L_NSPECTI,NAERKIND)
@@ -66 +67 @@
   real*8  LCOEF(4), LKCOEF(L_LEVELS,4)
 
-  real*8 taugsurf(L_NSPECTI,L_NGAUSS-1)
+  real*8,intent(out) :: taugsurf(L_NSPECTI,L_NGAUSS-1)
   real*8 DCONT,DAERO
   double precision wn_cont, p_cont, p_air, T_cont, dtemp, dtempc
@@ -72 +73 @@
 
   ! variable species mixing ratio variables
-  real*8  QVAR(L_LEVELS), WRATIO(L_LEVELS), MUVAR(L_LEVELS)
+  real*8,intent(in) :: QVAR(L_LEVELS)
+  real*8,intent(in) :: MUVAR(L_LEVELS)
+  real*8  WRATIO(L_LEVELS)
   real*8  KCOEF(4)
   integer NVAR(L_LEVELS)

trunk/LMDZ.GENERIC/libf/phystd/radii_mod.F90

@@ -7 +7 @@
 !     CO2 cloud properties (initialized in inifis)
       real,save :: Nmix_co2 ! Number mixing ratio of CO2 ice particles
+!$OMP THREADPRIVATE(Nmix_co2)
 
       ! flag to specify if we assume a constant fixed radius for particles
       logical,save :: radfixed ! initialized in inifis
+!$OMP THREADPRIVATE(radfixed)
 
 !     water cloud optical properties (initialized in su_aer_radii below)
@@ -17 +19 @@
       real, save ::  Nmix_h2o_ice
 !$OMP THREADPRIVATE(rad_h2o,rad_h2o_ice,Nmix_h2o,Nmix_h2o_ice)
-      real, parameter ::  coef_chaud=0.13
-      real, parameter ::  coef_froid=0.09
+
+      real,save :: nueff_iaero_h2o ! effective variance of H2O aerosol
+                                   ! (initialized in su_aer_radii below)
+!$OMP THREADPRIVATE(nueff_iaero_h2o)
+! coefficients for a variable nueff() for h2o aerosol; disabled for now 
+      real, parameter ::  coef_hot=0.13
+      real, parameter ::  coef_cold=0.09
 
 
@@ -76 +83 @@
 
          if(iaer.eq.iaero_h2o)then ! H2O ice
-            reffrad(1:ngrid,1:nlayer,iaer) = 1.e-5
-            nueffrad(1:ngrid,1:nlayer,iaer) = 0.1 
+           nueff_iaero_h2o=0.1 ! default value for variance of h2o aerosols
+           call getin_p("nueff_iaero_h2o",nueff_iaero_h2o)
+           if (is_master) write(*,*)" nueff_iaero_h2o = ",nueff_iaero_h2o
+           reffrad(1:ngrid,1:nlayer,iaer) = 1.e-5
+           nueffrad(1:ngrid,1:nlayer,iaer) = nueff_iaero_h2o 
          endif
 
@@ -85 +95 @@
          endif
  
-         if(iaer.eq.iaero_h2so4)then ! H2O ice
+         if(iaer.eq.iaero_h2so4)then ! H2SO4 ice
             reffrad(1:ngrid,1:nlayer,iaer) = 1.e-6
             nueffrad(1:ngrid,1:nlayer,iaer) = 0.1 
@@ -214 +224 @@
                zfice = MIN(MAX(zfice,0.0),1.0)
                reffrad(ig,l)= rad_h2o * (1.-zfice) + rad_h2o_ice * zfice
-               nueffrad(ig,l) = coef_chaud * (1.-zfice) + coef_froid * zfice
+!               nueffrad(ig,l) = coef_hot * (1.-zfice) + coef_cold * zfice
             enddo
          enddo
@@ -224 +234 @@
                zrad_liq  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o*pi*rhowater) )
                zrad_ice  = CBRT( 3*pq(ig,l)/(4*Nmix_h2o_ice*pi*rhowaterice) )
-               nueffrad(ig,l) = coef_chaud * (1.-zfice) + coef_froid * zfice
+!               nueffrad(ig,l) = coef_hot * (1.-zfice) + coef_cold * zfice
                zrad = zrad_liq * (1.-zfice) + zrad_ice * zfice
 
@@ -231 +241 @@
             enddo      
       end if
+
+! For now only constant nueff is enabled (otherwise some specific handling
+! of variable nueff is required in aeroptproperties)
+      nueffrad(1:ngrid,1:nlayer)=nueff_iaero_h2o
 
    end subroutine h2o_reffrad

trunk/LMDZ.GENERIC/libf/phystd/writediagfi.F

@@ -78 +78 @@
 
       integer,save :: zitau=0
-      character(len=20),save :: firstnom='1234567890'
+      character(len=40),save :: firstnom='1234567890'
 !$OMP THREADPRIVATE(zitau,firstnom)