Changeset 4677 for LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

Timestamp:

Sep 7, 2023, 1:07:27 PM (9 months ago)

Author:

idelkadi

Message:

Implementation in the LMDZ code of the double call of the ECRAD radiative transfer code to estimate the 3D radiative effect of clouds.

• This double call of Ecrad is controlled by the ok_3Deffect logic key.
• If this key is enabled, 2 files of parameter configuration "namelists" for ECRAD are required at runtime: namelist_ecrad and namelist_ecrad_s2.
• If this key is deactivated, the configuration and initialization part (reading namelist and netcdf files) is performed only once during simulation (1st call to ECRAD). Otherwise, configuration and initialization are performed each time Ecrad is called.

File:

• LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -399 +399 @@
 !$OMP THREADPRIVATE(swdnc0p, swdn0p, swdnp, swupc0p, swup0p, swupp)
 
+!AI ajout variables double appel Ecrad (3Deffect)
+      REAL,ALLOCATABLE,SAVE :: heat_s2(:,:), cool_s2(:,:)
+!$OMP THREADPRIVATE(heat_s2, cool_s2)
+      REAL,ALLOCATABLE,SAVE :: heat0_s2(:,:), cool0_s2(:,:)
+!$OMP THREADPRIVATE(heat0_s2, cool0_s2)
+      REAL,ALLOCATABLE,SAVE :: radsol_s2(:), topsw_s2(:), toplw_s2(:)
+!$OMP THREADPRIVATE(radsol_s2, topsw_s2, toplw_s2)
+      REAL,ALLOCATABLE,SAVE :: albpla_s2(:)
+!$OMP THREADPRIVATE(albpla_s2)
+      REAL,ALLOCATABLE,SAVE :: solsw_s2(:), solswfdiff_s2(:), sollw_s2(:)
+!$OMP THREADPRIVATE(solsw_s2, solswfdiff_s2, sollw_s2)
+      REAL,ALLOCATABLE,SAVE :: sollwdown_s2(:)
+!$OMP THREADPRIVATE(sollwdown_s2)
+      REAL,ALLOCATABLE,SAVE :: topsw0_s2(:),toplw0_s2(:)
+      REAL,ALLOCATABLE,SAVE :: solsw0_s2(:),sollw0_s2(:)
+!$OMP THREADPRIVATE(topsw0_s2,toplw0_s2,solsw0_s2,sollw0_s2)
+      REAL,ALLOCATABLE,SAVE :: lwdnc0_s2(:,:), lwdn0_s2(:,:), lwdn_s2(:,:)
+      REAL,ALLOCATABLE,SAVE :: lwupc0_s2(:,:), lwup0_s2(:,:), lwup_s2(:,:)
+!$OMP THREADPRIVATE(lwdnc0_s2,lwdn0_s2,lwdn_s2,lwupc0_s2,lwup0_s2,lwup_s2)       
+      REAL,ALLOCATABLE,SAVE :: swdnc0_s2(:,:), swdn0_s2(:,:), swdn_s2(:,:)
+      REAL,ALLOCATABLE,SAVE :: swupc0_s2(:,:), swup0_s2(:,:), swup_s2(:,:)
+!$OMP THREADPRIVATE(swdnc0_s2, swdn0_s2, swdn_s2, swupc0_s2, swup0_s2, swup_s2)
+
 ! pbase : cloud base pressure
 ! bbase : cloud base buoyancy
@@ -697 +720 @@
       ALLOCATE(swdnc0p(klon,klevp1), swdn0p(klon,klevp1), swdnp(klon,klevp1))
       ALLOCATE(swupc0p(klon,klevp1), swup0p(klon,klevp1), swupp(klon,klevp1))
+
+!AI Ajout pour Ecrad (3Deffect)       
+      ALLOCATE(heat_s2(klon,klev), cool_s2(klon,klev))
+      ALLOCATE(heat0_s2(klon,klev), cool0_s2(klon,klev))
+      ALLOCATE(radsol_s2(klon), topsw_s2(klon), toplw_s2(klon))
+      ALLOCATE(albpla_s2(klon))
+      ALLOCATE(solsw_s2(klon), solswfdiff_s2(klon), sollw_s2(klon))
+      ALLOCATE(sollwdown_s2(klon))
+      ALLOCATE(topsw0_s2(klon),toplw0_s2(klon))
+      ALLOCATE(solsw0_s2(klon),sollw0_s2(klon))
+      ALLOCATE(lwdnc0_s2(klon,klevp1), lwdn0_s2(klon,klevp1), lwdn_s2(klon,klevp1))
+      ALLOCATE(lwupc0_s2(klon,klevp1), lwup0_s2(klon,klevp1), lwup_s2(klon,klevp1))
+      ALLOCATE(swdnc0_s2(klon,klevp1), swdn0_s2(klon,klevp1), swdn_s2(klon,klevp1))
+      ALLOCATE(swupc0_s2(klon,klevp1), swup0_s2(klon,klevp1), swup_s2(klon,klevp1))
 
       ALLOCATE(cape(klon))
@@ -856 +893 @@
       DEALLOCATE(topsw0,toplw0,solsw0,sollw0)
       DEALLOCATE(albpla)
+
+!AI Ajout pour Ecrad (3Deffect)
+      DEALLOCATE(heat_s2, cool_s2)
+      DEALLOCATE(heat0_s2, cool0_s2)
+      DEALLOCATE(radsol_s2, topsw_s2, toplw_s2)
+      DEALLOCATE(albpla_s2)
+      DEALLOCATE(solsw_s2, solswfdiff_s2, sollw_s2)
+      DEALLOCATE(sollwdown_s2)
+      DEALLOCATE(topsw0_s2,toplw0_s2)
+      DEALLOCATE(solsw0_s2,sollw0_s2)
+      DEALLOCATE(lwdnc0_s2, lwdn0_s2, lwdn_s2)
+      DEALLOCATE(lwupc0_s2, lwup0_s2, lwup_s2)
+      DEALLOCATE(swdnc0_s2, swdn0_s2, swdn_s2)
+      DEALLOCATE(swupc0_s2, swup0_s2, swup_s2)
+
 !IM ajout variables CFMIP2/CMIP5
       DEALLOCATE(heatp, coolp)