Changeset 2220 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Jan 22, 2020, 11:41:33 AM (5 years ago)

Author:

dbardet

Message:

Update the nonorographic gravity waves drag parametrization (from the r3599 of Earth s model LMDZ6): 1) add east_ and west_gwdstress variables, 2) delete aleas function: random waves are producted by the MOD function, 3) reproductibility concerning the number of procs is now validated, 4) changing name of some variables like RUW, RVW, ZOP, ZOM,..., 5) tendency of winds due to GW drag are now module variables and written in restart files.

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -17 +17 @@
                      inquire_dimension, inquire_dimension_length
   use ioipsl_getincom, only : getin
-
+  use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd
   implicit none
   
@@ -415 +415 @@
 endif
 
+! Non-orographic gravity waves 
+call get_field("du_nonoro_gwd",du_nonoro_gwd,found,indextime)
+if (.not.found) then
+   write(*,*) "phyetat0: <du_nonoro_gwd> not in file"
+   du_nonoro_gwd(:,:)=0.
+else
+   write(*,*) "phyetat0: Memory of zonal wind tendency due to non-orographic GW"
+   write(*,*) " <du_nonoro_gwd> range:", &
+          minval(du_nonoro_gwd), maxval(du_nonoro_gwd)
+endif
+
+call get_field("dv_nonoro_gwd",dv_nonoro_gwd,found,indextime)
+if (.not.found) then
+   write(*,*) "phyetat0: <dv_nonoro_gwd> not in file"
+   dv_nonoro_gwd(:,:)=0.
+else
+   write(*,*) "phyetat0: Memory of meridional wind tendency due to non-orographic GW"
+   write(*,*) " <dv_nonoro_gwd> range:", &
+          minval(dv_nonoro_gwd), maxval(dv_nonoro_gwd)
+endif
 
 ! tracer on surface