Changeset 4853 for LMDZ6/trunk/libf/phylmd/radlwsw_m.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/radlwsw_m.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

@@ -1354 +1354 @@
       endif
 ! lldebug_for_offline
- 
+
+  if (namelist_ecrad_file.eq.'namelist_ecrad') then 
+      print*,' 1er apell Ecrad : ok_3Deffect, namelist_ecrad_file = ', &
+          ok_3Deffect, namelist_ecrad_file    
       CALL RADIATION_SCHEME &
       & (ist, iend, klon, klev, naero_grp, NSW, &
@@ -1396 +1399 @@
       & PSFSWDIF, PSFSWDIR, &
       & cloud_cover_sw)
+  else
+   print*,' 2e apell Ecrad : ok_3Deffect, namelist_ecrad_file = ', &
+          ok_3Deffect, namelist_ecrad_file       
+   CALL RADIATION_SCHEME_S2 &
+      & (ist, iend, klon, klev, naero_grp, NSW, &
+      & namelist_ecrad_file, ok_3Deffect, &
+      & debut, ok_volcan, flag_aerosol_strat, &
+      & day_cur, current_time, &
+!       Cste solaire/(d_Terre-Soleil)**2
+      & SOLARIRAD, &
+!       Cos(angle zin), temp sol              
+      & rmu0, tsol, &
+!       Albedo diffuse et directe
+      & PALBD_NEW,PALBP_NEW, &
+!       Emessivite : PEMIS_WINDOW (???), &
+      & ZEMIS, ZEMISW, &
+!       longitude(rad), sin(latitude), PMASQ_ ???
+      & ZGELAM, ZGEMU, &
+!       Temp et pres aux interf, vapeur eau, Satur spec humid 
+      & paprs_i, ZTH_i, q_i, qsat_i, & 
+!       Gas
+       & ZCO2, ZCH4, ZN2O, ZNO2, ZCFC11, ZCFC12, ZHCFC22, &
+       & ZCCL4, POZON_i(:,:,1), ZO2, &
+!       nuages :
+      & cldfra_i, flwc_i, fiwc_i, ZQ_SNOW, &
+!       rayons effectifs des gouttelettes              
+      & ref_liq_i, ref_ice_i, &
+!       aerosols
+     & ZAEROSOL_OLD, ZAEROSOL, &
+! Outputs
+!       Net flux :
+      & ZSWFT_i, ZLWFT_i, ZSWFT0_ii, ZLWFT0_ii, &
+!       DWN flux :
+      & ZFSDWN_i, ZFLUX_i(:,2,:), ZFCDWN_i, ZFLUC_i(:,2,:), &
+!       UP flux :
+      & ZFSUP_i, ZFLUX_i(:,1,:), ZFCUP_i, ZFLUC_i(:,1,:), &
+!       Surf Direct flux : ATTENTION
+      & ZFLUX_DIR, ZFLUX_DIR_CLEAR, ZFLUX_DIR_INTO_SUN, &
+!       UV and para flux
+      & ZFLUX_UV, ZFLUX_PAR, ZFLUX_PAR_CLEAR, &
+!      & ZFLUX_SW_DN_TOA,
+      & ZEMIS_OUT, ZLWDERIVATIVE, &
+      & PSFSWDIF, PSFSWDIR, &
+      & cloud_cover_sw)
+  endif
+
 
       print *,'========= RADLWSW: apres RADIATION_SCHEME ==================== '