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

Timestamp:

Jul 22, 2015, 4:14:51 PM (10 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_output_var_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_var_mod.F90

@@ -59 +59 @@
   LOGICAL, SAVE :: vars_defined = .FALSE. ! ug PAS THREADPRIVATE ET C'EST NORMAL
 
+  REAL, allocatable:: zustr_gwd_hines(:), zvstr_gwd_hines(:) ! (klon)
+  REAL, allocatable:: zustr_gwd_front(:), zvstr_gwd_front(:) ! (klon)
   REAL, allocatable:: zustr_gwd_rando(:), zvstr_gwd_rando(:) ! (klon)
   !$OMP THREADPRIVATE(zustr_gwd_rando, zvstr_gwd_rando)
@@ -86 +88 @@
     allocate (bils_ec(klon),bils_ech(klon),bils_tke(klon),bils_diss(klon),bils_kinetic(klon),bils_enthalp(klon),bils_latent(klon))
 
+    IF (ok_hines) allocate(zustr_gwd_hines(klon), zvstr_gwd_hines(klon))
+    IF (.not.ok_hines.and.ok_gwd_rando) &
+                  allocate(zustr_gwd_front(klon), zvstr_gwd_front(klon))
     IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))