Changeset 1989 for LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

Timestamp:

Mar 3, 2014, 10:57:40 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Inclusion du code RRTM

---

Adding RRTM code

MPL

File:

• LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -93 +93 @@
       REAL, SAVE, ALLOCATABLE :: topswcf_aero(:,:),  solswcf_aero(:,:)  ! diag
       !$OMP THREADPRIVATE(topswcf_aero,solswcf_aero)
+! Special RRTM
+      REAL, SAVE, ALLOCATABLE :: ZLWFT0_i(:,:),  ZSWFT0_i(:,:)      ! diag
+      !$OMP THREADPRIVATE(ZLWFT0_i,ZSWFT0_i)
+      REAL, SAVE, ALLOCATABLE :: ZFLDN0(:,:),  ZFLUP0(:,:)      ! diag
+      !$OMP THREADPRIVATE(ZFLDN0,ZFLUP0)
+      REAL, SAVE, ALLOCATABLE :: ZFSDN0(:,:),  ZFSUP0(:,:)      ! diag
+      !$OMP THREADPRIVATE(ZFSDN0,ZFSUP0)
+!
       REAL, SAVE, ALLOCATABLE :: tausum_aero(:,:,:) 
       !$OMP THREADPRIVATE(tausum_aero) 
@@ -243 +251 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq, ref_ice, theta, zphi
 !$OMP THREADPRIVATE(ref_liq, ref_ice, theta, zphi)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: ref_liq_pi, ref_ice_pi
+!$OMP THREADPRIVATE(ref_liq_pi, ref_ice_pi)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: wake_omg, zx_rh
 !$OMP THREADPRIVATE(wake_omg, zx_rh)
@@ -313 +323 @@
       allocate(d_u_lif(klon,klev),d_v_lif(klon,klev))
       allocate(d_ts(klon,nbsrf), d_tr(klon,klev,nbtr))
+! Special RRTM
+      allocate(ZLWFT0_i(klon,klev+1),ZSWFT0_i(klon,klev+1),ZFLDN0(klon,klev+1))
+      allocate(ZFLUP0(klon,klev+1),ZFSDN0(klon,klev+1),ZFSUP0(klon,klev+1))
+!
       allocate(topswad_aero(klon), solswad_aero(klon))
       allocate(topswai_aero(klon), solswai_aero(klon))
@@ -400 +414 @@
       ALLOCATE(fl(klon, klev), re(klon, klev), flwc(klon, klev))
       ALLOCATE(ref_liq(klon, klev), ref_ice(klon, klev), theta(klon, klev))
+      ALLOCATE(ref_liq_pi(klon, klev), ref_ice_pi(klon, klev))
       ALLOCATE(zphi(klon, klev), wake_omg(klon, klev), zx_rh(klon, klev))
       ALLOCATE(pmfd(klon, klev), pmfu(klon, klev))
@@ -550 +565 @@
       DEALLOCATE(fl, re, flwc)
       DEALLOCATE(ref_liq, ref_ice, theta)
+      DEALLOCATE(ref_liq_pi, ref_ice_pi)
       DEALLOCATE(zphi, wake_omg, zx_rh)
       DEALLOCATE(pmfd, pmfu)