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

Timestamp:

Jan 22, 2014, 10:51:36 AM (10 years ago)

Author:

idelkadi

Message:

Changes concerinng the Hourdin et al. 2002 version of the thermals and
Ayotte cases

calltherm.F90

---

call to thermcell_2002 modified :
1) iflag_thermals changed from 1 to >= 1000 ; iflag_thermals-1000
controls sub options.
2) thermals w and fraction in output files.

hbtm.F

---

Singularity in the 1/heat dependency of the Monin Obukov length L when
heat=0. Since 1/L is used rather than L, it is preferable to compute
directly L. There is a dependency in 1/u* then which is treated with a
threshold value.
(+ some cleaning in the syntax).

lmdz1d.F

---

If nday<0, -nday is used as the total number of time steps of the
simulation.
The option with imposed wtsurf and wqsurf read in lmdz1d.def was not
active anymore.
< IF(.NOT.ok_flux_surf) THEN
changed to

IF(.NOT.ok_flux_surf.or.max(abs(wtsurf),abs(wqsurf))>0.) THEN

before

fsens=-wtsurf*rcpd*rho(1)
flat=-wqsurf*rlvtt*rho(1)

phys_output_write.F90 et phys_local_var_mod.F90

---

Removing the d_u_ajsb contribution to duthe (same for dv).
Those tendencies are not computed by the model ...
< zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys - &
< d_u_ajsb(1:klon,1:klev)/pdtphys
---

zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys

! d_u_ajsb(1:klon,1:klev)/pdtphys

thermcell_dq.F90, thermcell_main.F90

---

Some cleaning

thermcell_old.F

---

Control by iflag_thermals >= 1000
wa_moy and fraca in outputs
+ cleaning

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -206 +206 @@
          rneb, rnebjn, zx_rh, frugs, agesno, d_t_dyn, d_q_dyn, &
          d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
-         d_u_ajsb, d_u_ajs, d_v_ajsb, d_v_ajs, &
+         d_u_ajs, d_v_ajs, &
          d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
          d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
@@ -1046 +1046 @@
        CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
        IF (vars_defined) THEN
-          zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys - &
-               d_u_ajsb(1:klon,1:klev)/pdtphys
+          zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
        ENDIF
        CALL histwrite_phy(o_duthe, zx_tmp_fi3d) 
        IF (vars_defined) THEN
-          zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys - &
-               d_v_ajsb(1:klon,1:klev)/pdtphys
+          zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys
        ENDIF
        CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)