Changeset 4031 for LMDZ6/trunk/libf/phylmd/ecrad/radiation_scheme.F90

Timestamp:

Nov 29, 2021, 3:29:39 PM (3 years ago)

Author:

idelkadi

Message:

Corrections in the LMDZ-ECRAN interface:

• gas concentrations in VMR instead of MMR
• use of a constant value for each gas instead of a table of values
• specify in the 14 SW bands of ECRAD those corresponding to the 6 SW bands for diffuse and direct albedo in LMDZ
• passage of all processing of input fields for Ecrad (conversions, calculations, ...) from radiation.F90 to radlw.F90

File:

• LMDZ6/trunk/libf/phylmd/ecrad/radiation_scheme.F90 (modified) (10 diffs)

LMDZ6/trunk/libf/phylmd/ecrad/radiation_scheme.F90

@@ -22 +22 @@
      &  PPRESSURE, PTEMPERATURE, &
      &  PPRESSURE_H, PTEMPERATURE_H, PQ, PQSAT, &
-     &  PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4, PO3_DP, &
+     &  PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, &
+     &  PCCL4, PO3, PO2, &
      &  PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW, &
      &  ZRE_LIQUID_UM, ZRE_ICE_UM, &
@@ -151 +152 @@
 ! AI
 REAL(KIND=JPRB),   INTENT(IN) :: PQSAT(KLON,KLEV)
-REAL(KIND=JPRB),   INTENT(IN) :: PCO2(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PCH4(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PN2O(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PNO2(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PCFC11(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PCFC12(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PHCFC22(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PCCL4(KLON,KLEV) 
-REAL(KIND=JPRB),   INTENT(IN) :: PO3_DP(KLON,KLEV) ! AI (kg/kg) ATTENTION (Pa*kg/kg) 
+REAL(KIND=JPRB),   INTENT(IN) :: PCO2
+REAL(KIND=JPRB),   INTENT(IN) :: PCH4 
+REAL(KIND=JPRB),   INTENT(IN) :: PN2O 
+REAL(KIND=JPRB),   INTENT(IN) :: PNO2 
+REAL(KIND=JPRB),   INTENT(IN) :: PCFC11
+REAL(KIND=JPRB),   INTENT(IN) :: PCFC12
+REAL(KIND=JPRB),   INTENT(IN) :: PHCFC22
+REAL(KIND=JPRB),   INTENT(IN) :: PCCL4
+REAL(KIND=JPRB),   INTENT(IN) :: PO3(KLON,KLEV) ! AI (kg/kg) ATTENTION (Pa*kg/kg) 
+REAL(KIND=JPRB),   INTENT(IN) :: PO2
 
 ! *** Cloud fraction and hydrometeor mass mixing ratios
@@ -326 +328 @@
   print*,'PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4 =', &
         PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4
-  print*,'PO3_DP =',PO3_DP
+  print*,'PO3 =',PO3
   print*,'PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW =', &
         PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW
@@ -345 +347 @@
 ! AI appel radiation_setup
 call SETUP_RADIATION_SCHEME(loutput)
+!! Les 6 bandes SW pour l'albedo :
+!! 0.185-0.25, 0.25-0.4, 0.4-0.69 , 0.69-1.19, 1.19-2.38, 2.38-4.00 micro-metre
+  call rad_config%define_sw_albedo_intervals(6, &
+             &  (/ 0.25e-6_jprb, 0.44e-6_jprb, 1.19e-6_jprb, &
+             &     2.38e-6_jprb, 4.00e-6_jprb /),  (/ 1,2,3,4,5,6 /))
 
 if (lprint_config) then
@@ -372 +379 @@
   print*,'rad_config%i_emiss_from_band_lw =', rad_config%i_emiss_from_band_lw
 endif
-!stop
-! A EFFACER
-!print*,'n_g_lw, n_g_sw =', rad_config%n_g_lw, rad_config%n_g_sw
-!print*,'use_canopy_full_spectrum_lw = ', rad_config%use_canopy_full_spectrum_lw
-!print*,'rad_config%i_band_from_reordered_g_lw =', &
-!       rad_config%i_band_from_reordered_g_lw
-!print*,'use_canopy_full_spectrum_lw =', rad_config%use_canopy_full_spectrum_lw
-!rad_config%use_canopy_full_spectrum_lw = .TRUE.
-! AI ATTENTION
-!rad_config%i_band_from_reordered_g_lw = 1
-!rad_config%use_spectral_solar_scaling = .true.
-!endif
-! AI ATTENTION test
-!rad_config%i_gas_model = IGasModelMonochromatic
 
 ! AI ATTENTION
 ! Allocate memory in radiation objects
-CALL single_level%allocate(KLON, NSW, 2, &
+CALL single_level%allocate(KLON, NSW, 1, &
      &                     use_sw_albedo_direct=.TRUE.)
 
@@ -420 +413 @@
 ! temperature at the half-level corresponding to the surface as
 ! follows:
-thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) &
-     &  = PTEMPERATURE(KIDIA:KFDIA,KLEV) &
-     &  + 0.5_JPRB * (PTEMPERATURE_H(KIDIA:KFDIA,KLEV+1) &
-     &               -PTEMPERATURE_H(KIDIA:KFDIA,KLEV))
+!thermodynamics%temperature_hl(KIDIA:KFDIA,KLEV+1) &
+!     &  = PTEMPERATURE(KIDIA:KFDIA,KLEV) &
+!     &  + 0.5_JPRB * (PTEMPERATURE_H(KIDIA:KFDIA,KLEV+1) &
+!     &               -PTEMPERATURE_H(KIDIA:KFDIA,KLEV))
 
 ! Alternatively we respect the model's atmospheric temperature in the
@@ -461 +454 @@
 !single_level%lw_emissivity(KIDIA:KFDIA,1)  = 1.0_JPRB
 single_level%lw_emissivity(KIDIA:KFDIA,1)  = PEMIS(KIDIA:KFDIA)
-single_level%lw_emissivity(KIDIA:KFDIA,2)  = PEMIS_WINDOW(KIDIA:KFDIA)
+!single_level%lw_emissivity(KIDIA:KFDIA,2)  = PEMIS_WINDOW(KIDIA:KFDIA)
 
 ! Create the relevant seed from date and time get the starting day
@@ -499 +492 @@
 !     &  PLAND_SEA_MASK, PCCN_LAND, PCCN_SEA, &
 !     &  ZRE_LIQUID_UM)
-cloud%re_liq(KIDIA:KFDIA,:) = ZRE_LIQUID_UM(KIDIA:KFDIA,:) * 1.0e-6_JPRB
+cloud%re_liq(KIDIA:KFDIA,:) = ZRE_LIQUID_UM(KIDIA:KFDIA,:)
 
 !CALL ICE_EFFECTIVE_RADIUS(KIDIA, KFDIA, KLON, KLEV, &
 !     &  PPRESSURE, PTEMPERATURE, PCLOUD_FRAC, PQ_ICE, PQ_SNOW, PGEMU, &
 !     &  ZRE_ICE_UM)
-cloud%re_ice(KIDIA:KFDIA,:) = ZRE_ICE_UM(KIDIA:KFDIA,:) * 1.0e-6_JPRB
+cloud%re_ice(KIDIA:KFDIA,:) = ZRE_ICE_UM(KIDIA:KFDIA,:)
 
 ! Get the cloud overlap decorrelation length (for cloud boundaries),
@@ -539 +532 @@
   ENDDO
 ENDDO
+!AI ATTENTION meme traitement dans le version offline
+!call cloud%create_inv_cloud_effective_size_eta(ncol, nlev, &
+!               &  thermodynamics%pressure_hl, &
+!               &  low_inv_effective_size, &
+!               &  middle_inv_effective_size, &
+!               &  high_inv_effective_size, 0.8_jprb, 0.45_jprb)      
 
 print*,'******** AEROSOLS (allocate + input) **************************************'
@@ -608 +607 @@
 !  ENDDO
 !ENDDO
-ZO3 = PO3_DP
 
 !  Insert gas mixing ratios
 print*,'Insert gas mixing ratios'
 CALL gas%put(IH2O,    IMassMixingRatio, PQ)
-CALL gas%put(ICO2,    IMassMixingRatio, PCO2)
-CALL gas%put(ICH4,    IMassMixingRatio, PCH4)
-CALL gas%put(IN2O,    IMassMixingRatio, PN2O)
-CALL gas%put(ICFC11,  IMassMixingRatio, PCFC11)
-CALL gas%put(ICFC12,  IMassMixingRatio, PCFC12)
-CALL gas%put(IHCFC22, IMassMixingRatio, PHCFC22)
-CALL gas%put(ICCL4,   IMassMixingRatio, PCCL4)
-CALL gas%put(IO3,     IMassMixingRatio, ZO3)
-CALL gas%put_well_mixed(IO2, IVolumeMixingRatio, 0.20944_JPRB)
+CALL gas%put(IO3,     IMassMixingRatio, PO3)
+CALL gas%put_well_mixed(ICO2,    IVolumeMixingRatio, PCO2)
+CALL gas%put_well_mixed(ICH4,    IVolumeMixingRatio, PCH4)
+CALL gas%put_well_mixed(IN2O,    IVolumeMixingRatio, PN2O)
+CALL gas%put_well_mixed(ICFC11,  IVolumeMixingRatio, PCFC11)
+CALL gas%put_well_mixed(ICFC12,  IVolumeMixingRatio, PCFC12)
+CALL gas%put_well_mixed(IHCFC22, IVolumeMixingRatio, PHCFC22)
+CALL gas%put_well_mixed(ICCL4,   IVolumeMixingRatio, PCCL4)
+CALL gas%put_well_mixed(IO2,     IVolumeMixingRatio, PO2)
 ! Ensure the units of the gas mixing ratios are what is required by
 ! the gas absorption model