Changeset 2381


Ignore:
Timestamp:
Jul 1, 2020, 6:11:48 PM (4 years ago)
Author:
slebonnois
Message:

GG+TN+SL : includes a limitation in the horizontal wavelength of non-orographic gravity waves, related to GCM resolution

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90

    r2198 r2381  
    1515
    1616      use dimphy
     17      USE geometry_mod, only: cell_area, latitude_deg
    1718      implicit none
    1819
     
    5758    INTEGER JK, JP, JO, JW
    5859    REAL KMIN, KMAX ! Min and Max horizontal wavenumbers
     60    REAL KSTAR ! Min horizontal wavenumber constrained by resolution
    5961    REAL CMIN, CMAX ! Min and Max absolute ph. vel.
    6062    REAL CPHA ! absolute PHASE VELOCITY frequency
     
    129131
    130132    !  PARAMETERS CORRESPONDING TO V3:
    131     RUWMAX = 0.005      ! Max EP-Flux at Launch altitude
     133    RUWMAX = 0.005      ! Max EP-Flux at Launch altitude (Nominal as Gilli2017)
     134!    RUWMAX = 0.0025    ! Half-nominal
     135!    RUWMAX = 0.01       ! Double-nominal
    132136    SAT    = 0.85       ! Saturation parameter: Sc in (12)
    133137    RDISS  = 0.1        ! Diffusion parameter
     
    135139    DELTAT=24.*3600.    ! Time scale of the waves (first introduced in 9b)
    136140
    137 !!!! TEST GG Values corresponding to min/max horizontal wavel 50-500 km
    138 !(similar to observations)
    139     KMIN = 1.E-5        ! Min horizontal wavenumber
    140 !    KMIN = 6.3E-6       ! Min horizontal wavenumber
    141     KMAX = 1.E-4        ! Max horizontal wavenumber
     141!!!! TESTS
     142!    KMIN = 1.E-5         ! Min horizontal wavenumber    (lambda max 628)
     143!    KMIN = 6.3E-6       ! Min horizontal wavenumber    (lambda max 1000 km)
     144!    KMIN = 1.2E-6       ! Min horizontal wavenumber     (lambda max 5000 km)
     145!    KMAX = 1.E-4        ! Max horizontal wavenumber     (lambda min 50 km)
     146
     147!!! TEST Kobs  and BestFit Diogo
     148    KMIN = 1.E-5         ! Min horizontal wavenumber    (lambda max 628)
     149    KMAX = 1.E-4        ! Max horizontal wavenumber     (lambda min 50 km)
     150
     151!!!  TEST GW8
     152
     153!    KMIN = 1.25E-5       ! Min horizontal wavenumber     (lambda max 500 km)
     154!    KMAX = 1.E-4        ! Max horizontal wavenumber     (lambda min 50 km)
     155
     156!!!  TEST GW9
     157
     158!    KMIN = 1.E-6       ! Min horizontal wavenumber     (lambda max 6000 km)
     159!    KMAX = 2.E-5       ! Max horizontal wavenumber     (lambda min 300 km)
     160
     161!!!  TEST BestFit GG
     162!    KMIN = 6.3E-6       ! Min horizontal wavenumber    (lambda max 1000 km)
     163!    KMAX = 1.E-4        ! Max horizontal wavenumber     (lambda min 50 km)
     164
    142165    !Online output: one value only
    143166    if (output) then
     
    145168      KMAX = 3E-5
    146169    endif
     170
     171!!
    147172    CMIN = 1.           ! Min phase velocity
    148     CMAX = 61.          ! Max phase speed velocity
    149 !    XLAUNCH=0.6         ! Parameter that control launching altitude
    150     XLAUNCH=5e-3        ! Value for top of cloud convective region
     173    CMAX = 61.          ! Max phase speed velocity  TestGW6
     174!     CMAX= 111.         ! TestGW2
     175!    XLAUNCH=0.6        ! Parameter that control launching altitude
     176     XLAUNCH=5e-3        ! Value for top of cloud convective region
     177!    XLAUNCH=2.e-4         ! Value for ~ 70 km? Review this
    151178
    152179!    PR = 9.2e6          ! Reference pressure    ! VENUS!!
     
    246273                ! Horizontal wavenumber amplitude
    247274!                ZK(JW, II) = KMIN + (KMAX - KMIN) * MOD(TT(II, JW) * 100., 1.)
    248                 ZK(JW, II) = KMIN + (KMAX - KMIN) * ALEAS(0.)
     275                ! TN+GG April/June 2020 - "Individual waves are not supposed to occupy
     276                ! the entire domain, but only a fraction of it" [Lott+2012]:
     277                KSTAR = RPI/SQRT(cell_area(II))     ! KSTAR should be < KMAX...
     278                ZK(JW, II) = MAX(KMIN,KSTAR) + (KMAX - MAX(KMIN,KSTAR)) * ALEAS(0.)
     279!               ZK(JW, II) = KMIN + (KMAX - KMIN) * ALEAS(0.)
    249280                ! Horizontal phase speed
    250281!                CPHA = CMIN + (CMAX - CMIN) * MOD(TT(II, JW)**2, 1.)
Note: See TracChangeset for help on using the changeset viewer.