Changeset 2400


Ignore:
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:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/README

    r2398 r2400  
    31213121== 06/07/2020 == EM
    31223122Some code cleanup: use "call abort_physic()" instead of "stop" or "call abort"
     3123
     3124== 07/07/2020 == DB
     3125Add 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

    r2399 r2400  
    3939      use assert_m, only : assert
    4040      use vertical_layers_mod, only : presnivs
     41      use geometry_mod, only: cell_area
    4142      implicit none
    4243
     
    8889    REAL, parameter:: kmax = 7.e-4 ! Max horizontal wavenumber
    8990    REAL, parameter:: kmin = 2.e-5 ! Min horizontal wavenumber
     91    REAL kstar                     ! Min horizontal wavenumber constrained by horizontal resolution
    9092    REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity
    9193    REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity   
    92     REAL CPHA                   ! absolute PHASE VELOCITY frequency
     94    REAL CPHA                    ! absolute PHASE VELOCITY frequency
    9395    REAL ZK(NW, ngrid)           ! Horizontal wavenumber amplitude
    9496    REAL ZP(NW, ngrid)           ! Horizontal wavenumber angle       
     
    263265                     * PI / 2.
    264266                ! Horizontal wavenumber amplitude
    265                 ZK(JW, II) = KMIN + (KMAX - KMIN) *RAN_NUM_2
     267                ! From Venus model: TN+GG April/June 2020 (rev2381)
     268                ! "Individual waves are not supposed to occupy the
     269                ! entire domain, but only a fraction of it" (Lott+2012)
     270               
     271                !ZK(JW, II) = KMIN + (KMAX - KMIN) *RAN_NUM_2
     272                KSTAR = PI/SQRT(cell_area(II))
     273                ZK(JW, II) = MAX(KMIN,KSTAR) + (KMAX - MAX(KMIN,KSTAR)) *RAN_NUM_2
    266274                ! Horizontal phase speed
    267275                CPHA = 0.
Note: See TracChangeset for help on using the changeset viewer.