Changeset 4188 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad

Timestamp:

Jun 29, 2022, 9:23:14 AM (2 years ago)

Author:

idelkadi

Message:

Implementation of the Ecrad radiative transfer code in the LMD model (continued) :
Integration of aerosols (direct effect)

Location:

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad

Files:

• radiation_aerosol_optics.F90 (modified) (3 diffs)
• radiation_scheme.F90 (modified) (4 diffs)
• radiation_setup.F90 (modified) (3 diffs)
• readaerosol_optic_ecrad.F90 (added)
• setup_aerosol_optics_lmdz_m.F90 (added)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_aerosol_optics.F90

@@ -89 +89 @@
          &  IAerosolClassUndefined,   IAerosolClassIgnored, &
          &  IAerosolClassHydrophobic, IAerosolClassHydrophilic
+    USE phys_local_var_mod, ONLY: rhcl 
 
     integer, intent(in) :: nlev               ! number of model levels
@@ -187 +188 @@
       end if
 
-      call gas%get(IH2O, IMassMixingRatio, h2o_mmr, istartcol=istartcol)
+!AI juin 2022      
+      !call gas%get(IH2O, IMassMixingRatio, h2o_mmr, istartcol=istartcol)
 
       ! Loop over position
@@ -195 +197 @@
           ! saturation and the index to the relative-humidity index of
           ! hydrophilic-aerosol data
-          rh  = h2o_mmr(jcol,jlev) / thermodynamics%h2o_sat_liq(jcol,jlev)
-          irh = ao%calc_rh_index(rh)
+! AI juin 2022          
+!          rh  = h2o_mmr(jcol,jlev) / thermodynamics%h2o_sat_liq(jcol,jlev)
+!          irh = ao%calc_rh_index(rh)
+          irh = ao%calc_rh_index(rhcl(jcol,jlev))
+!          print*,'irh=',irh
 
           factor = ( thermodynamics%pressure_hl(jcol,jlev+1) &

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_scheme.F90

@@ -13 +13 @@
 SUBROUTINE RADIATION_SCHEME &
 ! Inputs
-     & (KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW, &
+     & (KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW, &
      &  IDAY, TIME, &
      &  PSOLAR_IRRADIANCE, &
@@ -117 +117 @@
 INTEGER(KIND=JPIM),INTENT(IN) :: KLEV     ! Number of levels
 !INTEGER, INTENT(IN) :: KLON, KLEV
-INTEGER(KIND=JPIM),INTENT(IN) :: KAEROLMDZ ! Number of aerosol types
+!INTEGER(KIND=JPIM),INTENT(IN) :: KAEROLMDZ ! Number of aerosol types
+INTEGER(KIND=JPIM),INTENT(IN) :: KAEROSOL
 INTEGER(KIND=JPIM),INTENT(IN) :: NSW ! Numbe of bands
 
 ! AI ATTENTION
-INTEGER, PARAMETER :: KAEROSOL = 12
+!INTEGER, PARAMETER :: KAEROSOL = 12
 
 ! *** Single-level fields
@@ -313 +314 @@
 if (lprint_input) then
   print*,'********** Verification des entrees *************'
-  print*,'KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW =', &
-        KIDIA, KFDIA, KLON, KLEV, KAEROLMDZ, NSW
+  print*,'KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW =', &
+        KIDIA, KFDIA, KLON, KLEV, KAEROSOL, NSW
   print*,'IDAY, TIME =', IDAY, TIME
   print*,'PSOLAR_IRRADIANCE =', PSOLAR_IRRADIANCE
@@ -510 +511 @@
 
 print*,'******** AEROSOLS (allocate + input) **************************************'
-IF (NAERMACC > 0) THEN
+!IF (NAERMACC > 0) THEN
   CALL aerosol%allocate(KLON, 1, KLEV, KAEROSOL) ! MACC climatology
-ELSE
-  CALL aerosol%allocate(KLON, 1, KLEV, 6) ! Tegen climatology
-ENDIF
+!ELSE
+!  CALL aerosol%allocate(KLON, 1, KLEV, 6) ! Tegen climatology
+!ENDIF
 ! Compute the dry mass of each layer neglecting humidity effects, in
 ! kg m-2, needed to scale some of the aerosol inputs

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/ecrad/radiation_setup.F90

@@ -238 +238 @@
       ! hydrophobic black carbon, 11 ammonium sulphate, 12 inactive
       ! SO2
-      rad_config%n_aerosol_types = 12
-      if (lprint_setp) then
-        print*,'rad_config%n_aerosol_types =', rad_config%n_aerosol_types
-      endif
+!      rad_config%n_aerosol_types = 12
+!      if (lprint_setp) then
+!        print*,'rad_config%n_aerosol_types =', rad_config%n_aerosol_types
+!      endif
 
       ! Indices to the aerosol optical properties in
@@ -247 +247 @@
       ! index hydrophilic aerosol types and positive numbers index
       ! hydrophobic aerosol types
-      rad_config%i_aerosol_type_map = 0 ! There can be up to 256 types
+!      rad_config%i_aerosol_type_map = 0 ! There can be up to 256 types
 !      if (lprint_setp) then
 !        print*,'rad_config%i_aerosol_type_map =', rad_config%i_aerosol_type_map
 !      endif
 
-      rad_config%i_aerosol_type_map(1:12) = (/ &
-           &  -1, &  ! Sea salt, size bin 1 (OPAC)
-           &  -2, &  ! Sea salt, size bin 2 (OPAC)
-           &  -3, &  ! Sea salt, size bin 3 (OPAC)
-           &   7, &  ! Desert dust, size bin 1 (Woodward 2001)
-           &   8, &  ! Desert dust, size bin 2 (Woodward 2001)
-           &   9, &  ! Desert dust, size bin 3 (Woodward 2001)
-           &  -4, &  ! Hydrophilic organic matter (OPAC)
-           &  10, &  ! Hydrophobic organic matter (OPAC)
-           &  11, &  ! Black carbon (Boucher)
-           &  11, &  ! Black carbon (Boucher)
-           &  -5, &  ! Ammonium sulphate (OPAC)
-           &  14 /)  ! Stratospheric sulphate (hand edited from OPAC)
+!      rad_config%i_aerosol_type_map(1:12) = (/ &
+!           &  -1, &  ! Sea salt, size bin 1 (OPAC)
+!           &  -2, &  ! Sea salt, size bin 2 (OPAC)
+!           &  -3, &  ! Sea salt, size bin 3 (OPAC)
+!           &   7, &  ! Desert dust, size bin 1 (Woodward 2001)
+!           &   8, &  ! Desert dust, size bin 2 (Woodward 2001)
+!           &   9, &  ! Desert dust, size bin 3 (Woodward 2001)
+!           &  -4, &  ! Hydrophilic organic matter (OPAC)
+!           &  10, &  ! Hydrophobic organic matter (OPAC)
+!           &  11, &  ! Black carbon (Boucher)
+!           &  11, &  ! Black carbon (Boucher)
+!           &  -5, &  ! Ammonium sulphate (OPAC)
+!           &  14 /)  ! Stratospheric sulphate (hand edited from OPAC)
 !     if (lprint_setp) then
 !       print*,'rad_config%i_aerosol_type_map =', rad_config%i_aerosol_type_map
@@ -433 +433 @@
          &  'photosynthetically active radiation, PAR')
 
+     rad_config%i_aerosol_type_map(1:13) = (/ &
+           &  -1, &  ! Sea salt, size bin 1 (OPAC)
+           &  -2, &  ! Sea salt, size bin 2 (OPAC)
+           &  -3, &  ! Sea salt, size bin 3 (OPAC)
+           &  -4, &  ! Hydrophilic organic matter (OPAC)
+           &  -5, &  ! Ammonium sulphate (OPAC)
+           &  -6, &
+           &  -7, &
+           &   1, &
+           &   2, &
+           &   3, &
+           &  -8, &
+           &  -9, &
+           &   4 /)  ! Stratospheric sulphate (hand edited from OPAC)
+   rad_config%aerosol_optics_override_file_name = 'aerosol_optics_lmdz.nc'
+
 !    IF (YRERAD%NAERMACC > 0) THEN
       ! With the MACC aerosol climatology we need to add in the