Changeset 2400 for trunk

Timestamp:

Jul 7, 2020, 3:11:53 PM (4 years ago)

Author:

dbardet

Message:

Add kstar parameter to control kmin value. kmi=1/lambda_max, to ensure the subgrid scale characteristic, we have to constrain the GWs wavelength by the size of the mesh

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/nonoro_gwd_ran_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/README

@@ -3121 +3121 @@
 == 06/07/2020 == EM
 Some code cleanup: use "call abort_physic()" instead of "stop" or "call abort"
+
+== 07/07/2020 == DB
+Add kstar parameter to control kmin value. kmi=1/lambda_max, to ensure the "subgrid scale" characteristic, we have to constrain the GW's wavelength by the size of the mesh

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

@@ -39 +39 @@
       use assert_m, only : assert
       use vertical_layers_mod, only : presnivs
+      use geometry_mod, only: cell_area
       implicit none
 
@@ -88 +89 @@
     REAL, parameter:: kmax = 7.e-4 ! Max horizontal wavenumber
     REAL, parameter:: kmin = 2.e-5 ! Min horizontal wavenumber
+    REAL kstar                     ! Min horizontal wavenumber constrained by horizontal resolution
     REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity
     REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity    
-    REAL CPHA                   ! absolute PHASE VELOCITY frequency
+    REAL CPHA                    ! absolute PHASE VELOCITY frequency
     REAL ZK(NW, ngrid)           ! Horizontal wavenumber amplitude
     REAL ZP(NW, ngrid)           ! Horizontal wavenumber angle        
@@ -263 +265 @@
                      * PI / 2.
                 ! Horizontal wavenumber amplitude
-                ZK(JW, II) = KMIN + (KMAX - KMIN) *RAN_NUM_2
+                ! From Venus model: TN+GG April/June 2020 (rev2381)
+                ! "Individual waves are not supposed to occupy the 
+                ! entire domain, but only a fraction 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.