Changeset 4853 for LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_interface.F90

Timestamp:

Mar 19, 2024, 3:34:21 PM (2 months ago)

Author:

idelkadi

Message:

Ecrad update in LMDZ / Implementation of Ecrad double call in LMDZ

• version 1.6.1 (commit 210d7911380f53a788c3cad73b3cf9b4e022ef87)
• interface routines between lmdz and ecrad grouped in a new "lmdz" directory
• double call method redesigned so as to go through the Ecrad initialization and configuration part only once for the entire simulation
• clean-up in the read.F routine: delete unitules arguments
• modification of default gas model in namelist (default: ECCKD)

File:

• LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_interface.F90 (modified) (12 diffs)

LMDZ6/trunk/libf/phylmd/ecrad/radiation/radiation_interface.F90

@@ -68 +68 @@
 
     ! Load the look-up tables from files in the specified directory
-    if (config%i_gas_model == IGasModelMonochromatic) then
+    if (config%i_gas_model_sw == IGasModelMonochromatic) then
       call setup_gas_optics_mono(config, trim(config%directory_name))
-    else if (config%i_gas_model == IGasModelIFSRRTMG) then
-      call setup_gas_optics_rrtmg(config, trim(config%directory_name))
-    else if (config%i_gas_model == IGasModelECCKD) then
-      call setup_gas_optics_ecckd(config)
+    else
+      ! Note that we can run RRTMG and ECCKD for different parts of
+      ! the spectrum: the setup routines only configure the relevant
+      ! part.
+      if (config%i_gas_model_sw == IGasModelIFSRRTMG .or. config%i_gas_model_lw == IGasModelIFSRRTMG) then
+        call setup_gas_optics_rrtmg(config, trim(config%directory_name))
+      end if
+      if (config%i_gas_model_sw == IGasModelECCKD .or. config%i_gas_model_lw == IGasModelECCKD) then
+        call setup_gas_optics_ecckd(config)
+      end if
     end if
 
@@ -122 +128 @@
 
     if (config%do_clouds) then
-      if (config%i_gas_model == IGasModelMonochromatic) then
+      if (config%i_gas_model_sw == IGasModelMonochromatic) then
         !      call setup_cloud_optics_mono(config)
       elseif (config%use_general_cloud_optics) then
@@ -132 +138 @@
 
     if (config%use_aerosols) then
-      if (config%i_gas_model == IGasModelMonochromatic) then
+      if (config%i_gas_model_sw == IGasModelMonochromatic) then
 !        call setup_aerosol_optics_mono(config)
       else 
@@ -167 +173 @@
     type(gas_type),    intent(inout) :: gas
 
-    if (config%i_gas_model == IGasModelMonochromatic) then
+    if (config%i_gas_model_sw == IGasModelMonochromatic) then
       call set_gas_units_mono(gas)
-    elseif (config%i_gas_model == IGasModelECCKD) then
+    elseif (config%i_gas_model_sw == IGasModelIFSRRTMG &
+         &  .or. config%i_gas_model_lw == IGasModelIFSRRTMG) then
+      ! Convert to mass-mixing ratio for RRTMG; note that ecCKD can
+      ! work with this but performs an internal scaling
+      call set_gas_units_ifs(gas)
+    else
+      ! Use volume mixing ratio preferred by ecCKD
       call set_gas_units_ecckd(gas)
-    else
-      call set_gas_units_ifs(gas)
     end if
 
@@ -196 +206 @@
     use radiation_io,             only : nulout
     use radiation_config,         only : config_type, &
-         &   IGasModelMonochromatic, IGasModelIFSRRTMG, &
+         &   IGasModelMonochromatic, IGasModelIFSRRTMG, IGasModelECCKD, &
          &   ISolverMcICA, ISolverSpartacus, ISolverHomogeneous, &
          &   ISolverTripleclouds
@@ -227 +237 @@
     use radiation_general_cloud_optics, only : general_cloud_optics
     use radiation_aerosol_optics, only : add_aerosol_optics
-    USE mod_phys_lmdz_para
-
 
     ! Inputs
@@ -238 +246 @@
     type(thermodynamics_type),intent(in)   :: thermodynamics
     type(gas_type),           intent(in)   :: gas
-    type(cloud_type),        intent(inout):: cloud
+    type(cloud_type),         intent(inout):: cloud
     type(aerosol_type),       intent(in)   :: aerosol
     ! Output
@@ -297 +305 @@
 
     real(jphook) :: hook_handle
-    integer :: jcol, jlev
 
     if (lhook) call dr_hook('radiation_interface:radiation',0,hook_handle)
-
-    if (config%i_solver_sw == ISolverSPARTACUS) then
-         print*,'Dans radiation, mpi_rank, omp_rank, size, chape inv_cloud = ',&
-              mpi_rank, omp_rank, &
-              shape(cloud%inv_cloud_effective_size), &
-              size(cloud%inv_cloud_effective_size)   
-!         do jcol=istartcol, iendcol
-!              do jlev=1,nlev
-!              print*,'Entree radiation_interf, mpi_rank, omp_rank, jcol, jlev &
-!           &  cloud%inv_cloud_effective_size =',mpi_rank, omp_rank, jcol, jlev, &
-!           &  cloud%inv_cloud_effective_size(jcol,jlev)
-!              enddo
-!          enddo
-    endif
-!    cloud%inv_cloud_effective_size=0.05_jprb
 
     if (thermodynamics%pressure_hl(istartcol,2) &
@@ -339 +331 @@
       ! incoming shortwave flux at each g-point, for the specified
       ! range of atmospheric columns
-      if (config%i_gas_model == IGasModelMonochromatic) then
+      if (config%i_gas_model_sw == IGasModelMonochromatic) then
         call gas_optics_mono(ncol,nlev,istartcol,iendcol, config, &
              &  single_level, thermodynamics, gas, lw_albedo, &
              &  od_lw, od_sw, ssa_sw, &
              &  planck_hl, lw_emission, incoming_sw)
-      else if (config%i_gas_model == IGasModelIFSRRTMG) then
-        call gas_optics_rrtmg(ncol,nlev,istartcol,iendcol, config, &
-             &  single_level, thermodynamics, gas, &
-             &  od_lw, od_sw, ssa_sw, lw_albedo=lw_albedo, &
-             &  planck_hl=planck_hl, lw_emission=lw_emission, &
-             &  incoming_sw=incoming_sw)
       else
-        call gas_optics_ecckd(ncol,nlev,istartcol,iendcol, config, &
-             &  single_level, thermodynamics, gas, &
-             &  od_lw, od_sw, ssa_sw, lw_albedo=lw_albedo, &
-             &  planck_hl=planck_hl, lw_emission=lw_emission, &
-             &  incoming_sw=incoming_sw)
+        if (config%i_gas_model_sw == IGasModelIFSRRTMG &
+             &   .or. config%i_gas_model_lw == IGasModelIFSRRTMG) then
+          call gas_optics_rrtmg(ncol,nlev,istartcol,iendcol, config, &
+               &  single_level, thermodynamics, gas, &
+               &  od_lw, od_sw, ssa_sw, lw_albedo=lw_albedo, &
+               &  planck_hl=planck_hl, lw_emission=lw_emission, &
+               &  incoming_sw=incoming_sw)
+        end if
+        if (config%i_gas_model_sw == IGasModelECCKD &
+             &   .or. config%i_gas_model_lw == IGasModelECCKD) then
+          call gas_optics_ecckd(ncol,nlev,istartcol,iendcol, config, &
+               &  single_level, thermodynamics, gas, &
+               &  od_lw, od_sw, ssa_sw, lw_albedo=lw_albedo, &
+               &  planck_hl=planck_hl, lw_emission=lw_emission, &
+               &  incoming_sw=incoming_sw)
+        end if
       end if
 
@@ -368 +365 @@
         ! Compute hydrometeor absorption/scattering properties in each
         ! shortwave and longwave band
-        if (config%i_gas_model == IGasModelMonochromatic) then
+        if (config%i_gas_model_sw == IGasModelMonochromatic) then
           call cloud_optics_mono(nlev, istartcol, iendcol, &
                &  config, thermodynamics, cloud, &
@@ -387 +384 @@
 
       if (config%use_aerosols) then
-        if (config%i_gas_model == IGasModelMonochromatic) then
+        if (config%i_gas_model_sw == IGasModelMonochromatic) then
 !          call add_aerosol_optics_mono(nlev,istartcol,iendcol, &
 !               &  config, thermodynamics, gas, aerosol, & 
@@ -474 +471 @@
         else if (config%i_solver_sw == ISolverSPARTACUS) then
           ! Compute fluxes using the SPARTACUS shortwave solver
-!             cloud%inv_cloud_effective_size=0.05_jprb
-!             do jcol=istartcol, iendcol
-!              do jlev=1,nlev
-!              print*,'jcol, jlev, dans radiation_interf i &
-!               &       cloud%inv_cloud_effective_size =',jcol, jlev, &
-!               cloud%inv_cloud_effective_size(jcol,jlev)
-!              enddo
-!             enddo
           call solver_spartacus_sw(nlev,istartcol,iendcol, &
                &  config, single_level, thermodynamics, cloud, &