Changeset 2403 for trunk/LMDZ.GENERIC/libf

Timestamp:

Jul 9, 2020, 11:10:48 AM (4 years ago)

Author:

dbardet

Message:

Add kstar parameter to control kmin value (set value in the parameterization). kmin=1/lambda_max, to ensure the subgrid scale characteristic, we have to constrain the maximum GW's wavelength (lambda_max) by the size of the mesh. So, kmin=max(kstar, kmin prescribed value)

File:

• trunk/LMDZ.GENERIC/libf/phystd/nonoro_gwd_ran_mod.F90 (modified) (5 diffs)

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

@@ -49 +49 @@
                               gwd_convective_source
      use inifis_mod
+     use geometry_mod, only: cell_area
 #ifdef CPP_XIOS
      use xios_output_mod, only: send_xios_field
@@ -97 +98 @@
      INTEGER, parameter:: nw = nk * np *no ! Total numbers of gravity waves
      INTEGER jk, jp, jo, jw
-     REAL, parameter:: kmax = 7.e-4 ! Max horizontal wavenumber
-     REAL, parameter:: kmin = 2.e-5 ! Min horizontal wavenumber
-     !REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity
-     REAL, parameter:: cmax = 100.  ! Test for Saturn: Max horizontal absolute phase velocity
+     REAL kstar                  ! Control value to constrain the min horizontal wavenumber by the horizontal resolution
+     REAL, parameter:: kmax = 1.e-4 ! Max horizontal wavenumber=N/u, u(=30) zonal wind velocity at launch
+     !REAL, parameter:: kmax = 4.e-5 ! Max horizontal wavenumber=N/u, u(=70) zonal wind velocity at launch
+     REAL, parameter:: kmin = 2.e-6 ! Min horizontal wavenumber=1/sqrt(DxDy) Dx and Dy horizontal grid 
+     REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity=zonal wind at launch
+     !REAL, parameter:: cmax = 100.  ! Test for Saturn: Max horizontal absolute phase velocity
+     !REAL, parameter:: cmax = 70.   ! Test for Saturn: Max horizontal absolute phase velocity
      REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity
      REAL cpha                      ! absolute phase velocity frequency 
@@ -117 +121 @@
      REAL epflux_0(nw, ngrid)             ! Fluxes at launching level (previous name: RUW0)
      INTEGER launch                       ! Launching altitude
-     REAL, parameter:: xlaunch = 0.4      ! Control the launching altitude
-     REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer (approx.)
+     REAL, parameter:: xlaunch = 0.2      ! Control the launching altitude
+     REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer in m (approx.)
 
      REAL prec(ngrid)     ! precipitations
@@ -131 +135 @@
      REAL H0bis(ngrid, nlayer)       ! characteristic height of the atmosphere
      REAL, save::  H0                ! characteristic height of the atmosphere
-     REAL, parameter:: pr = 250      ! Reference pressure [Pa]
-     REAL, parameter:: tr = 220.     ! Reference temperature [K]
+     !REAL, parameter:: pr = 250      ! Reference pressure [Pa]
+     !REAL, parameter:: tr = 220.     ! Reference temperature [K]
+     ! TEST FOR SATURN STRATOSPHERE
+     REAL, parameter:: pr = 110      ! Reference pressure [Pa]
+     REAL, parameter:: tr = 130.     ! Reference temperature [K]
      REAL zh(ngrid, nlayer + 1)      ! Log-pressure altitude (constant H0) 
      REAL zhbis(ngrid, nlayer + 1)   ! Log-pressure altitude (varying H)
@@ -237 +244 @@
                         + 1.) * pi / 2.
            ! horizontal wavenumber amplitude
-           zk(jw, ii) = kmin + (kmax - kmin) * ran_num_2
+           ! TN+GG April/June 2020
+           ! "Individual waves are not supposed
+           ! to occupy the entire domain, but
+           ! only a faction of it" Lott+2012
+
+           !zk(jw, ii) = kmin + (kmax - kmin) * ran_num_2
+           kstar = pi / sqrt(cell_area(ii))
+           zk(jw, ii) = max(kmin,kstar) + (kmax - max(kmin,kstar)) * ran_num_2
            ! horizontal phase speed
            cpha = 0.