Changeset 2333 for LMDZ5/trunk/libf/phylmd/phys_state_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_state_var_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -395 +395 @@
 !!! fin nrlmd le 10/04/2012
 
-      REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), dv_gwd_rando(:, :)
-      !$OMP THREADPRIVATE(du_gwd_rando, dv_gwd_rando)
+      REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), du_gwd_front(:, :)
+      !$OMP THREADPRIVATE(du_gwd_rando, du_gwd_front)
       ! tendencies on wind due to gravity waves
 
@@ -581 +581 @@
       ALLOCATE(ale_bl_trig(klon))
 !!! fin nrlmd le 10/04/2012
-      if (ok_gwd_rando) allocate(du_gwd_rando(klon, klev), &
-           dv_gwd_rando(klon, klev))
+      if (ok_gwd_rando) allocate(du_gwd_rando(klon, klev))
+      if (.not. ok_hines .and. ok_gwd_rando) allocate(du_gwd_front(klon, klev))
 
 END SUBROUTINE phys_state_var_init
@@ -705 +705 @@
       deallocate(tau_aero_lw_rrtm,piz_aero_lw_rrtm,cg_aero_lw_rrtm)
       deallocate(ccm)
-      if (ok_gwd_rando) deallocate(du_gwd_rando, dv_gwd_rando)
+      if (ok_gwd_rando) deallocate(du_gwd_rando)
+      if (.not. ok_hines .and. ok_gwd_rando) deallocate(du_gwd_front)
        
 !!! nrlmd le 10/04/2012