Changeset 1938 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Jan 20, 2014, 5:20:05 PM (10 years ago)

Author:

lguez

Message:

New run-time parameters: ok_gwd_rando, ruwmax and sat. Logical
variable ok_gwd_rando activates the new stochastic parameterization of
the drag of non-orographic gravity waves, by F. Lott. Default value of
ok_gwd_rando is false. A run with ok_gwd_rando creates new variables
in restartphy.nc: du_gwd_rando and dv_gwd_rando. Possible output in
history files of new variables du_gwd_rando, dv_gwd_rando and
vstr_gwd_rando.

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -146 +146 @@
          o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
          o_dtr_ls, o_dtr_trsp, o_dtr_sscav, &
-         o_dtr_sat, o_dtr_uscav, o_trac_cum
+         o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
+         o_vstr_gwd_rando
 
     USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &
@@ -167 +168 @@
          vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
          vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
-         T2sumSTD, nlevSTD
+         T2sumSTD, nlevSTD, du_gwd_rando, dv_gwd_rando
 
     USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &
@@ -216 +217 @@
 
     USE phys_output_var_mod, only: vars_defined, snow_o, zfra_o, bils_diss, &
-         bils_ec, bils_tke, bils_kinetic, bils_latent, &
-         bils_enthalp, itau_con, nfiles, clef_files, &
-         nid_files
+         bils_ec, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
+         itau_con, nfiles, clef_files, nid_files, zvstr_gwd_rando
     USE indice_sol_mod, only: nbsrf
     USE infotrac, only: nqtot
@@ -892 +892 @@
        CALL histwrite_phy(o_rnebcon, rnebcon)
        CALL histwrite_phy(o_rnebls, rneb)
-IF (vars_defined)  THEN
-      DO k=1, klev
-         DO i=1, klon
-            zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
-         ENDDO
-      ENDDO
-ENDIF
+       IF (vars_defined)  THEN
+          DO k=1, klev
+             DO i=1, klon
+                zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
+             ENDDO
+          ENDDO
+       ENDIF
        CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d)
        CALL histwrite_phy(o_rhum, zx_rh)
@@ -1119 +1119 @@
           CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
        ENDIF
+
+       IF (ok_gwd_rando) then
+          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_vstr_gwd_rando, zvstr_gwd_rando)
+       end IF
+
        CALL histwrite_phy(o_rsu, swup)
        CALL histwrite_phy(o_rsd, swdn)