Changeset 2299 for trunk/LMDZ.GENERIC/libf/phystd/phyredem.F90

Timestamp:

Apr 27, 2020, 4:39:37 PM (5 years ago)

Author:

dbardet

Message:

27/04/2017 (r2299) == DB

Add non-orographic gravity waves drag parameterization (Flott scheme adpated from the Earth GCM)

It can be called using 'calllott_nonoro=.true.' in callphys.def, and set the maximum value of the Eliassen-Plam flux 'epflux_max'.
Cumulated output fields are du_nonoro, dv_nonoro (winds tendency due to GW drag), east_gwstress and west_gwstress (stress profile in esatward and westward direction due to GW drag)
These variables are added in Xhisitns and start files.

File:

• trunk/LMDZ.GENERIC/libf/phystd/phyredem.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/phyredem.F90

@@ -142 +142 @@
   use tracer_h, only: noms
   use slab_ice_h, only: noceanmx
-  use callkeys_mod, only: ok_slab_ocean
+  use callkeys_mod, only: ok_slab_ocean, calllott_nonoro
+  use nonoro_gwd_ran_mod, only: du_nonoro_gwd, dv_nonoro_gwd, &
+                                east_gwstress, west_gwstress
 
   implicit none
@@ -208 +210 @@
     enddo
   endif ! of if (nq>0)
+
+! Non-orographic gavity waves
+if (calllott_nonoro) then
+   call put_field("du_nonoro_gwd","Zonal wind tendency due to GW",du_nonoro_gwd)
+   call put_field("dv_nonoro_gwd","Meridional wind tendency due to GW",dv_nonoro_gwd)
+   call put_field("east_gwstress","Eastward stress profile due to GW",east_gwstress)
+   call put_field("west_gwstress","Westward stress profile due to GW",west_gwstress)
+endif
+
 ! close file
       CALL close_restartphy