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

Timestamp:

Jul 22, 2015, 4:14:51 PM (9 years ago)

Author:

lguez

Message:

New parameterization of gravity wave drag due to front/jet systems, by

1. de la Camara and F. Lott. The new Camara-Lott parameterization

replaces the Hines parameterization so it is activated if not ok_hines
and ok_gwd_rando.

Also changed distribution of phase speeds in FLOTT_GWD_rando, from
uniform to Gaussian. Bug fix in sugwd_strato. Bug fix in the arguments
of the call to add_phys_tend for methane oxydation.

For the new Camara-Lott parameterization, we need to compute relative
vorticity in calfis and pass it as a new argument "rot" to
physiq. Interpolation of relative vorticity to the physics grid is not
optimal for now: it is not weighted by cell areas.

Alvaro de la Camara, Fran\c{}cois Lott

File:

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

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -81 +81 @@
       !$OMP THREADPRIVATE(d_u_lif, d_v_lif)
 ! Tendances Ondes de G non oro (runs strato).
-      REAL, SAVE, ALLOCATABLE :: d_u_hin(:,:)
-      !$OMP THREADPRIVATE(d_u_hin)
-      REAL, SAVE, ALLOCATABLE :: d_v_hin(:,:)
-      !$OMP THREADPRIVATE(d_v_hin)
+      REAL, SAVE, ALLOCATABLE :: du_gwd_hines(:,:)
+      !$OMP THREADPRIVATE(du_gwd_hines)
+      REAL, SAVE, ALLOCATABLE :: dv_gwd_hines(:,:)
+      !$OMP THREADPRIVATE(dv_gwd_hines)
+      REAL, SAVE, ALLOCATABLE :: dv_gwd_rando(:,:)
+      !$OMP THREADPRIVATE(dv_gwd_rando)
+      REAL, SAVE, ALLOCATABLE :: dv_gwd_front(:,:)
+      !$OMP THREADPRIVATE(dv_gwd_front)
+      REAL, SAVE, ALLOCATABLE :: east_gwstress(:,:)
+      !$OMP THREADPRIVATE(east_gwstress)
+      REAL, SAVE, ALLOCATABLE :: west_gwstress(:,:)
+      !$OMP THREADPRIVATE(west_gwstress)
       REAL, SAVE, ALLOCATABLE :: d_t_hin(:,:)
       !$OMP THREADPRIVATE(d_t_hin)
@@ -434 +442 @@
       allocate(topsw0_aero(klon,naero_grp), solsw0_aero(klon,naero_grp))
       allocate(topswcf_aero(klon,3), solswcf_aero(klon,3))
-      allocate(d_u_hin(klon,klev),d_v_hin(klon,klev),d_t_hin(klon,klev))
+      allocate(du_gwd_hines(klon,klev),dv_gwd_hines(klon,klev))
+      allocate(dv_gwd_rando(klon,klev),dv_gwd_front(klon,klev))
+      allocate(east_gwstress(klon,klev),west_gwstress(klon,klev))
+      allocate(d_t_hin(klon,klev))
       allocate(d_q_ch4(klon,klev))
 !      allocate(tausum_aero(klon,nwave,naero_spc))
@@ -679 +690 @@
       deallocate(load_tmp6)
       deallocate(load_tmp7)
-      deallocate(d_u_hin,d_v_hin,d_t_hin)
+      deallocate(du_gwd_hines,dv_gwd_hines,d_t_hin)
       deallocate(d_q_ch4)
+      deallocate(dv_gwd_rando,dv_gwd_front)
+      deallocate(east_gwstress,west_gwstress)
 
 !IM ajout variables CFMIP2/CMIP5