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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -1 +1 @@
 !
-! $Id$
+! $Header$
 !
 MODULE phys_output_write_mod
@@ -132 +132 @@
          o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
          o_dtoro, o_dulif, o_dvlif, o_dtlif, &
-         o_duhin, o_dvhin, o_dthin, o_dqch4, o_rsu, &
+         o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &
+         o_du_gwd_front, o_dv_gwd_front, &
+         o_east_gwstress, o_west_gwstress, &
          o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, &
          o_rlucs, o_rldcs, o_tnt, o_tntr, &
@@ -151 +153 @@
          o_dtr_ls, o_dtr_trsp, o_dtr_sscav, &
          o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
-         o_vstr_gwd_rando
+         o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
+         o_ustr_gwd_front,o_vstr_gwd_front
 
     USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &
@@ -174 +177 @@
          vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
          vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
-         T2sumSTD, nlevSTD, du_gwd_rando, dv_gwd_rando, &
+         T2sumSTD, nlevSTD, du_gwd_rando, du_gwd_front, &
          ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
          rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, &
@@ -230 +233 @@
          zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &
          d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
-         d_v_lif, d_t_lif, d_u_hin, d_v_hin, d_t_hin, &
+         d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, d_t_hin, &
+         dv_gwd_rando, dv_gwd_front, &
+         east_gwstress, west_gwstress, &
          d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD
 
     USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &
          bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
-         itau_con, nfiles, clef_files, nid_files, zvstr_gwd_rando
+         itau_con, nfiles, clef_files, nid_files, &
+         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
+         zustr_gwd_front, zvstr_gwd_front                                   
     USE ocean_slab_mod, only: tslab, slab_bils, slab_bilg, tice, seaice
     USE pbl_surface_mod, only: snow
@@ -1217 +1224 @@
 
        IF (ok_hines) THEN
-          IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys
-          CALL histwrite_phy(o_duhin, zx_tmp_fi3d)
-          IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys
-          CALL histwrite_phy(o_dvhin, zx_tmp_fi3d)
+          CALL histwrite_phy(o_du_gwd_hines, du_gwd_hines/pdtphys)
+          CALL histwrite_phy(o_dv_gwd_hines, dv_gwd_hines/pdtphys)
           IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
           CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
+          CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
+          CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
+       end IF
+
+       if (.not. ok_hines .and. ok_gwd_rando) then
+          CALL histwrite_phy(o_du_gwd_front, du_gwd_front / pdtphys)
+          CALL histwrite_phy(o_dv_gwd_front, dv_gwd_front / pdtphys)
+          CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front)
+          CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front)
        ENDIF
 
@@ -1228 +1242 @@
           CALL histwrite_phy(o_du_gwd_rando, du_gwd_rando / pdtphys)
           CALL histwrite_phy(o_dv_gwd_rando, dv_gwd_rando / pdtphys)
+          CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando)
           CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando)
+          CALL histwrite_phy(o_east_gwstress, east_gwstress )
+          CALL histwrite_phy(o_west_gwstress, west_gwstress )
        end IF