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

Timestamp:

Jan 20, 2014, 5:20:05 PM (11 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/physiq.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -55 +55 @@
   !IM stations CFMIP
   USE CFMIP_point_locations
+  use FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
+
   IMPLICIT none
   !>======================================================================
@@ -614 +616 @@
 
   ! tendance nulles
-  REAL du0(klon,klev),dv0(klon,klev),dq0(klon,klev),dql0(klon,klev)
+  REAL, dimension(klon,klev):: du0, dv0, dt0, dq0, dql0
 
   !
@@ -1390 +1392 @@
   du0(:,:)=0.
   dv0(:,:)=0.
+  dt0 = 0.
   dq0(:,:)=0.
   dql0(:,:)=0.
@@ -3251 +3254 @@
   ENDIF
 
-  !
-  !IM cf. FLott BEG
+  if (ok_gwd_rando) then
+     call FLOTT_GWD_rando(DTIME, pplay, t_seri, u_seri, v_seri, &
+          rain_fall + snow_fall, zustr_gwd_rando, zvstr_gwd_rando, &
+          du_gwd_rando, dv_gwd_rando)
+     CALL add_phys_tend(du_gwd_rando, dv_gwd_rando, dt0, dq0, dql0, &
+          'flott_gwd_rando')
+  end if
+
   ! STRESS NECESSAIRES: TOUTE LA PHYSIQUE