Changeset 2403 for trunk


Ignore:
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)

Location:
trunk/LMDZ.GENERIC
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/README

    r2354 r2403  
    15871587Major cleanup: remove obsolete compilation scripts (makegcm*) and old dynamical
    15881588core, as it is obsolete with respect to the one provide in LMDZ.COMMON.
     1589
     1590
     1591== 09/07/2020 == DB
     1592Add 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)
     1593

  • trunk/LMDZ.GENERIC/deftank/def_generic_giants/PROTO_Saturn_Article/z2sig.def

    r1164 r2403  
    1919179.0
    2020191.0
    21 203.0
    22 215.0
    23 227.0
    24 239.0
     21201.0
     22211.0
     23221.0
     24231.0
     25241.0
    2526251.0
    26 263.0
     27259.0
     28267.0
    2729275.0
    28 287.0
     30283.0
     31291.0
    2932299.0
    30 311.0
    31 323.0
    32 335.0
    33 347.0
    34 359.0
    35 371.0
    36 383.0
    37 395.0
    38 407.0
    39 419.0
    40 431.0
    41 443.0
    42 455.0
    43 467.0
    44 479.0
    45 491.0
    46 503.0
    47 515.0
    48 527.0
    49 539.0
    50 551.0
    51 563.0
    52 575.0
    53 587.0
    54 599.0
    55 611.0
    56 623.0
    57 635.0
    58 647.0
    59 659.0
    60 671.0
    61 683.0
     33307.0
     34313.0
     35319.0
     36325.0
     37331.0
     38337.0
     39343.0
     40349.0
     41355.0
     42361.0
     43367.0
     44373.0
     45379.0
     46385.0
     47391.0
     48397.0
     49403.0
     50409.0
     51415.0
     52421.0
     53427.0
     54433.0
     55439.0
     56445.0
     57451.0
     58457.0
     59463.0
     60469.0
     61475.0
     62481.0
     63487.0
     64493.0
     65499.0
     66505.0
     67511.0
     68517.0
     69523.0
     70529.0
     71535.0
     72541.0
     73547.0
     74553.0
     75559.0
     76571.0
     77577.0
     78583.0
     79589.0
     80595.0
     81601.0
     82607.0
     83613.0
     84619.0
     85625.0
     86631.0
     87637.0
     88643.0
     89649.0
     90655.0
     91661.0
     92667.0
     93673.0
     94679.0
     95685.0
     96691.0
    6297695.0
    6398707.0

  • trunk/LMDZ.GENERIC/deftank/saturn1d/profile

    r1150 r2403  
    6464130
    6565130
     66130
     67130
     68130
     69130
     70130
     71130
     72130
     73130
     74130
     75130
     76130
     77130
     78130
     79130
     80130
     81130
     82130
     83130
     84130
     85130
     86130
     87130
     88130
     89130
     90130
     91130
     92130
     93130
     94130
     95130
     96130
     97130
     98130
     99130
     100130
     101130

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

    r2305 r2403  
    4949                              gwd_convective_source
    5050     use inifis_mod
     51     use geometry_mod, only: cell_area
    5152#ifdef CPP_XIOS
    5253     use xios_output_mod, only: send_xios_field
    … …  
    9798     INTEGER, parameter:: nw = nk * np *no ! Total numbers of gravity waves
    9899     INTEGER jk, jp, jo, jw
    99      REAL, parameter:: kmax = 7.e-4 ! Max horizontal wavenumber
    100      REAL, parameter:: kmin = 2.e-5 ! Min horizontal wavenumber
    101      !REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity
    102      REAL, parameter:: cmax = 100.  ! Test for Saturn: Max horizontal absolute phase velocity
     100     REAL kstar                  ! Control value to constrain the min horizontal wavenumber by the horizontal resolution
     101     REAL, parameter:: kmax = 1.e-4 ! Max horizontal wavenumber=N/u, u(=30) zonal wind velocity at launch
     102     !REAL, parameter:: kmax = 4.e-5 ! Max horizontal wavenumber=N/u, u(=70) zonal wind velocity at launch
     103     REAL, parameter:: kmin = 2.e-6 ! Min horizontal wavenumber=1/sqrt(DxDy) Dx and Dy horizontal grid
     104     REAL, parameter:: cmax = 30.   ! Max horizontal absolute phase velocity=zonal wind at launch
     105     !REAL, parameter:: cmax = 100.  ! Test for Saturn: Max horizontal absolute phase velocity
     106     !REAL, parameter:: cmax = 70.   ! Test for Saturn: Max horizontal absolute phase velocity
    103107     REAL, parameter:: cmin = 1.    ! Min horizontal absolute phase velocity
    104108     REAL cpha                      ! absolute phase velocity frequency
    … …  
    117121     REAL epflux_0(nw, ngrid)             ! Fluxes at launching level (previous name: RUW0)
    118122     INTEGER launch                       ! Launching altitude
    119      REAL, parameter:: xlaunch = 0.4      ! Control the launching altitude
    120      REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer (approx.)
     123     REAL, parameter:: xlaunch = 0.2      ! Control the launching altitude
     124     REAL, parameter:: zmaxth_top = 8000. ! Top of convective layer in m (approx.)
    121125
    122126     REAL prec(ngrid)     ! precipitations
    … …  
    131135     REAL H0bis(ngrid, nlayer)       ! characteristic height of the atmosphere
    132136     REAL, save::  H0                ! characteristic height of the atmosphere
    133      REAL, parameter:: pr = 250      ! Reference pressure [Pa]
    134      REAL, parameter:: tr = 220.     ! Reference temperature [K]
     137     !REAL, parameter:: pr = 250      ! Reference pressure [Pa]
     138     !REAL, parameter:: tr = 220.     ! Reference temperature [K]
     139     ! TEST FOR SATURN STRATOSPHERE
     140     REAL, parameter:: pr = 110      ! Reference pressure [Pa]
     141     REAL, parameter:: tr = 130.     ! Reference temperature [K]
    135142     REAL zh(ngrid, nlayer + 1)      ! Log-pressure altitude (constant H0)
    136143     REAL zhbis(ngrid, nlayer + 1)   ! Log-pressure altitude (varying H)
    … …  
    237244                        + 1.) * pi / 2.
    238245           ! horizontal wavenumber amplitude
    239            zk(jw, ii) = kmin + (kmax - kmin) * ran_num_2
     246           ! TN+GG April/June 2020
     247           ! "Individual waves are not supposed
     248           ! to occupy the entire domain, but
     249           ! only a faction of it" Lott+2012
     250
     251           !zk(jw, ii) = kmin + (kmax - kmin) * ran_num_2
     252           kstar = pi / sqrt(cell_area(ii))
     253           zk(jw, ii) = max(kmin,kstar) + (kmax - max(kmin,kstar)) * ran_num_2
    240254           ! horizontal phase speed
    241255           cpha = 0.
Note: See TracChangeset for help on using the changeset viewer.