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

Timestamp:

Mar 31, 2023, 8:42:57 PM (18 months ago)

Author:

lguez

Message:

Merge LMDZ_ECRad branch back into trunk!

Location:

LMDZ6/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/ecrad/radiation_scheme.F90 (modified) (14 diffs)

LMDZ6/trunk/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, &
-     &  PMU0, PTEMPERATURE_SKIN, PALBEDO_DIF, PALBEDO_DIR, &
+     &  PMU0, PTEMPERATURE_SKIN, &
+     &  PALBEDO_DIF, PALBEDO_DIR, &
      &  PEMIS, PEMIS_WINDOW, &
-     &  PCCN_LAND, PCCN_SEA, &
-     &  PGELAM, PGEMU, PLAND_SEA_MASK, &
-     &  PPRESSURE, PTEMPERATURE, &
+     &  PGELAM, PGEMU, &
      &  PPRESSURE_H, PTEMPERATURE_H, PQ, PQSAT, &
      &  PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, &
      &  PCCL4, PO3, PO2, &
-     &  PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW, &
+     &  PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_SNOW, &
      &  ZRE_LIQUID_UM, ZRE_ICE_UM, &
      &  PAEROSOL_OLD, PAEROSOL, &
@@ -117 +116 @@
 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
@@ -139 +139 @@
 REAL(KIND=JPRB),   INTENT(IN) :: PGEMU(KLON)
 ! Land-sea mask
-REAL(KIND=JPRB),   INTENT(IN) :: PLAND_SEA_MASK(KLON) 
+!REAL(KIND=JPRB),   INTENT(IN) :: PLAND_SEA_MASK(KLON) 
 
 ! *** Variables on full levels
-REAL(KIND=JPRB),   INTENT(IN) :: PPRESSURE(KLON,KLEV)    ! (Pa)
-REAL(KIND=JPRB),   INTENT(IN) :: PTEMPERATURE(KLON,KLEV) ! (K)
+!REAL(KIND=JPRB),   INTENT(IN) :: PPRESSURE(KLON,KLEV)    ! (Pa)
+!REAL(KIND=JPRB),   INTENT(IN) :: PTEMPERATURE(KLON,KLEV) ! (K)
 ! *** Variables on half levels
 REAL(KIND=JPRB),   INTENT(IN) :: PPRESSURE_H(KLON,KLEV+1)    ! (Pa)
@@ -167 +167 @@
 REAL(KIND=JPRB),   INTENT(IN) :: PQ_LIQUID(KLON,KLEV)
 REAL(KIND=JPRB),   INTENT(IN) :: PQ_ICE(KLON,KLEV)
-REAL(KIND=JPRB),   INTENT(IN) :: PQ_RAIN(KLON,KLEV)
+!REAL(KIND=JPRB),   INTENT(IN) :: PQ_RAIN(KLON,KLEV)
 REAL(KIND=JPRB),   INTENT(IN) :: PQ_SNOW(KLON,KLEV)
 
@@ -174 +174 @@
 REAL(KIND=JPRB),   INTENT(IN) :: PAEROSOL(KLON,KLEV,KAEROSOL)
 
-REAL(KIND=JPRB),   INTENT(IN) :: PCCN_LAND(KLON) 
-REAL(KIND=JPRB),   INTENT(IN) :: PCCN_SEA(KLON) 
+!REAL(KIND=JPRB),   INTENT(IN) :: PCCN_LAND(KLON) 
+!REAL(KIND=JPRB),   INTENT(IN) :: PCCN_SEA(KLON) 
 
 !AI mars 2021
@@ -313 +313 @@
 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
@@ -320 +320 @@
   print*,'PTEMPERATURE_SKIN =',PTEMPERATURE_SKIN
   print*,'PEMIS, PEMIS_WINDOW =', PEMIS, PEMIS_WINDOW
-  print*,'PCCN_LAND, PCCN_SEA =', PCCN_LAND, PCCN_SEA
   print*,'PGELAM, PGEMU =', PGELAM, PGEMU
-  print*,'PPRESSURE =', PPRESSURE
-  print*,'PTEMPERATURE =', PTEMPERATURE
   print*,'PPRESSURE_H =', PPRESSURE_H
   print*,'PTEMPERATURE_H =', PTEMPERATURE_H
@@ -331 +328 @@
         PCO2, PCH4, PN2O, PNO2, PCFC11, PCFC12, PHCFC22, PCCL4
   print*,'PO3 =',PO3
-  print*,'PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW =', &
-        PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_RAIN, PQ_SNOW
+  print*,'PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_SNOW =', &
+        PCLOUD_FRAC, PQ_LIQUID, PQ_ICE, PQ_SNOW
   print*,'ZRE_LIQUID_UM, ZRE_ICE_UM =', &
         ZRE_LIQUID_UM, ZRE_ICE_UM
@@ -407 +404 @@
 ! Alternative approximate version using temperature and pressure from
 ! the thermodynamics structure
-CALL thermodynamics%calc_saturation_wrt_liquid(KIDIA, KFDIA)
+!CALL thermodynamics%calc_saturation_wrt_liquid(KIDIA, KFDIA)
+!AI ATTENTION
+thermodynamics%h2o_sat_liq = PQSAT
 
 print*,'********** SINGLE LEVEL VARS **********************************'
 !AI ATTENTION
-!thermodynamics%h2o_sat_liq = PQSAT
 ! Set single-level fileds
 single_level%solar_irradiance              = PSOLAR_IRRADIANCE
@@ -510 +508 @@
 
 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
@@ -644 +642 @@
 
 ! Compute UV fluxes as weighted sum of appropriate shortwave bands
+!AI ATTENTION
+if (0.eq.1) then
 PFLUX_UV       (KIDIA:KFDIA) = 0.0_JPRB
 DO JBAND = 1,NWEIGHT_UV
@@ -660 +660 @@
        &  * flux%sw_dn_surf_clear_band(IBAND_PAR(JBAND),KIDIA:KFDIA)
 ENDDO
-
+endif
 ! Compute effective broadband emissivity
 ZBLACK_BODY_NET_LW = flux%lw_dn(KIDIA:KFDIA,KLEV+1) &
@@ -679 +679 @@
 !AI ATTENTION
 !IF (YRERAD%LAPPROXSWUPDATE) THEN
+if (0.eq.1) then
 IF (rad_config%do_surface_sw_spectral_flux) THEN
   PSWDIFFUSEBAND(KIDIA:KFDIA,:) = 0.0_JPRB
@@ -693 +694 @@
   ENDDO
 ENDIF
-
+endif
 CALL single_level%deallocate
 CALL thermodynamics%deallocate