Changeset 2340 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jun 9, 2020, 11:48:45 AM (4 years ago)

Author:

jvatant

Message:

In addition to r2297, for the n-layer aerosol scheme, enables to set the particle size effective variance with aeronlay_nueff in callphys.def.
--JVO

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callkeys_mod.F90 (modified) (1 diff)
• inifis_mod.F90 (modified) (2 diffs)
• radii_mod.F90 (modified) (2 diffs)

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

@@ -117 +117 @@
       real,dimension(:),allocatable,save :: aeronlay_sclhght
       real,dimension(:),allocatable,save :: aeronlay_size
-!$OMP THREADPRIVATE(aeronlay_tauref,aeronlay_lamref,aeronlay_ptop,aeronlay_pbot,aeronlay_sclhght,aeronlay_size)
+      real,dimension(:),allocatable,save :: aeronlay_nueff
+!$OMP THREADPRIVATE(aeronlay_tauref,aeronlay_lamref,aeronlay_ptop,aeronlay_pbot,aeronlay_sclhght,aeronlay_size,aeronlay_nueff)
       real,save :: co2supsat
       real,save :: pceil

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

@@ -704 +704 @@
      IF(.NOT.ALLOCATED(aeronlay_sclhght))     ALLOCATE(aeronlay_sclhght(nlayaero))     
      IF(.NOT.ALLOCATED(aeronlay_size))        ALLOCATE(aeronlay_size(nlayaero))     
+     IF(.NOT.ALLOCATED(aeronlay_nueff))       ALLOCATE(aeronlay_nueff(nlayaero))     
      IF(.NOT.ALLOCATED(optprop_aeronlay_ir))  ALLOCATE(optprop_aeronlay_ir(nlayaero))     
      IF(.NOT.ALLOCATED(optprop_aeronlay_vis)) ALLOCATE(optprop_aeronlay_vis(nlayaero))     
@@ -739 +740 @@
      write(*,*)" aeronlay_sclhght = ",aeronlay_sclhght
 
-     write(*,*)"Generic n-layer aerosols: particles sizes (m)"
+     write(*,*)"Generic n-layer aerosols: particles effective radii (m)"
      aeronlay_size=1.e-6
      call getin_p("aeronlay_size",aeronlay_size)
      write(*,*)" aeronlay_size = ",aeronlay_size
+
+     write(*,*)"Generic n-layer aerosols: particles radii effective variance"
+     aeronlay_nueff=0.1
+     call getin_p("aeronlay_nueff",aeronlay_nueff)
+     write(*,*)" aeronlay_nueff = ",aeronlay_nueff
 
      write(*,*)"Generic n-layer aerosols: VIS optical properties file"

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

@@ -38 +38 @@
       use aerosol_mod, only: iaero_back2lay, iaero_co2, iaero_dust, &
                              iaero_h2o, iaero_h2so4, iaero_nh3, iaero_nlay, iaero_aurora
-      use callkeys_mod, only: size_nh3_cloud, nlayaero, aeronlay_size
+      use callkeys_mod, only: size_nh3_cloud, nlayaero, aeronlay_size, aeronlay_nueff
 
       Implicit none
@@ -94 +94 @@
               if(iaer.eq.iaero_nlay(ia))then ! N-layer aerosols
                  reffrad(1:ngrid,1:nlayer,iaer) = aeronlay_size(ia)
-                 nueffrad(1:ngrid,1:nlayer,iaer) = 0.1 
+                 nueffrad(1:ngrid,1:nlayer,iaer) = aeronlay_nueff(ia) 
               endif
             enddo