Changeset 2333 for LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_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_output_ctrlout_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1190 +1190 @@
   TYPE(ctrl_out), SAVE :: o_dvlif = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     'dvlif', 'Orography dV', 'm/s2', (/ ('', i=1, 9) /))
-  TYPE(ctrl_out), SAVE :: o_duhin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
-    'duhin', 'Hines GWD dU', 'm/s2', (/ ('', i=1, 9) /))
-  TYPE(ctrl_out), SAVE :: o_dvhin = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
-    'dvhin', 'Hines GWD dV', 'm/s2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_du_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'du_gwd_hines', 'Hines GWD dU', 'm/s2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_dv_gwd_hines = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'dv_gwd_hines', 'Hines GWD dV', 'm/s2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_du_gwd_front = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'du_gwd_front', 'Fronts GWD dU', 'm/s2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_dv_gwd_front = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'dv_gwd_front', 'Fronts GWD dV', 'm/s2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_east_gwstress = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'east_gwstress', 'Eastward GW Stress', 'Pa', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_west_gwstress = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'west_gwstress', 'Westward GW Stress', 'Pa', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_dtoro = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     'dtoro', 'Orography dT', 'K/s', (/ ('', i=1, 9) /))
@@ -1209 +1217 @@
        = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'dv_gwd_rando', &
        "Random gravity waves dV/dt", "m/s2", (/ ('', i=1, 9) /))
+  type(ctrl_out), save:: o_ustr_gwd_hines &
+       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_hines', &
+       "zonal wind stress Hines gravity waves", "Pa", (/ ('', i=1, 9) /))
+  type(ctrl_out), save:: o_vstr_gwd_hines &
+       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_hines', &
+       "meridional wind stress Hines gravity waves", "Pa", (/ ('', i=1, 9) /))
+  type(ctrl_out), save:: o_ustr_gwd_front &
+       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_front', &
+       "zonal wind stress fronts gravity waves", "Pa", (/ ('', i=1, 9) /))
+  type(ctrl_out), save:: o_vstr_gwd_front &
+       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_front', &
+       "meridional wind stress fronts gravity waves", "Pa", (/ ('', i=1, 9) /))
+  type(ctrl_out), save:: o_ustr_gwd_rando &
+       = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'ustr_gwd_rando', &
+       "zonal wind stress random gravity waves", "Pa", (/ ('', i=1, 9) /))
   type(ctrl_out), save:: o_vstr_gwd_rando &
        = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), 'vstr_gwd_rando', &