Changeset 2299 for trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.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/inifis_mod.F90 (modified) (2 diffs)

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

@@ -22 +22 @@
   use planetwide_mod, only: planetwide_sumval
   use callkeys_mod
+  use ioipsl_getin_p_mod, only : getin_p
   use mod_phys_lmdz_para, only : is_parallel
 
@@ -290 +291 @@
      write(*,*) " calladj = ",calladj
 
+     write(*,*)"call Lott's non-oro GWs parameterisation scheme ?"
+     calllott_nonoro=.false. ! default value
+     call getin_p("calllott_nonoro",calllott_nonoro)
+     write(*,*)" calllott_nonoro = ",calllott_nonoro
+     
+     write(*,*)"Max Value of EP flux at launching altitude"
+     epflux_max=0.0                ! default value
+     call getin_p("epflux_max",epflux_max)
+     write(*,*)"epflux_max = ",epflux_max
+
      write(*,*) "call thermal plume model ?"
      calltherm=.false. ! default value