Changeset 2072 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Jun 23, 2014, 4:53:24 PM (10 years ago)
Author:
lguez
Message:

Chose more specific variable names.

Location:
LMDZ5/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/YOEGWD.h

    r1938 r2072  
    1414      real GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC
    1515
    16       REAL RUWMAX
     16      REAL GWD_RANDO_RUWMAX
    1717!     Maximum Eliassen-Palm flux at launch level, in "FLOTT_GWD_rando"
    1818
    19       REAL SAT ! saturation parameter in "FLOTT_GWD_rando"
     19      REAL GWD_RANDO_SAT ! saturation parameter in "FLOTT_GWD_rando"
    2020!     S_c in equation (12) of Lott (JGR, vol 118, page 8897, 2013)
    2121
     
    2323      COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT,        &
    2424     &     GHMAX,GRAHILO,GSIGCR,NKTOPG,NSTRA,GSSEC,GTSEC,GVSEC,         &
    25      &     RUWMAX, sat
     25     &     GWD_RANDO_RUWMAX, gwd_rando_sat
    2626
    2727      save /YOEGWD/
  • LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

    r2057 r2072  
    186186    LOGICAL,SAVE :: ok_strato_omp
    187187    LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp
    188     real, SAVE:: RUWMAX_omp, sat_omp
     188    real, SAVE:: GWD_RANDO_RUWMAX_omp, gwd_rando_sat_omp
    189189    LOGICAL,SAVE      :: carbon_cycle_tr_omp
    190190    LOGICAL,SAVE      :: carbon_cycle_cpl_omp
     
    16871687    CALL getin('ok_gwd_rando', ok_gwd_rando_omp)
    16881688
    1689     RUWMAX_omp = 0.01
    1690     CALL getin('ruwmax', RUWMAX_omp)
    1691 
    1692     sat_omp = 0.25
    1693     CALL getin('sat', sat_omp)
     1689    GWD_RANDO_RUWMAX_omp = 0.01
     1690    CALL getin('gwd_rando_ruwmax', GWD_RANDO_RUWMAX_omp)
     1691
     1692    gwd_rando_sat_omp = 0.25
     1693    CALL getin('gwd_rando_sat', gwd_rando_sat_omp)
    16941694
    16951695    !Config Key  = OK_LES                                               
     
    19171917    ok_hines = ok_hines_omp
    19181918    ok_gwd_rando = ok_gwd_rando_omp
    1919     RUWMAX = RUWMAX_omp
    1920     sat = sat_omp
     1919    GWD_RANDO_RUWMAX = GWD_RANDO_RUWMAX_omp
     1920    gwd_rando_sat = gwd_rando_sat_omp
    19211921    ok_LES = ok_LES_omp
    19221922    callstats = callstats_omp
     
    21092109    write(lunout,*) 'ok_hines = ',  ok_hines
    21102110    write(lunout,*) 'ok_gwd_rando = ',  ok_gwd_rando
    2111     write(lunout,*) 'RUWMAX = ', RUWMAX
    2112     write(lunout,*) 'sat = ', sat
     2111    write(lunout,*) 'GWD_RANDO_RUWMAX = ', GWD_RANDO_RUWMAX
     2112    write(lunout,*) 'gwd_rando_sat = ', gwd_rando_sat
    21132113    write(lunout,*) 'read_climoz = ', read_climoz
    21142114    write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr
  • LMDZ5/trunk/libf/phylmd/flott_gwd_rando_m.F90

    r1938 r2072  
    238238                     + ZK(JW, II) * SIN(ZP(JW, II)) * VH(II, LAUNCH)
    239239                ! Momentum flux at launch lev
    240                 RUW0(JW, II) = RUWMAX
     240                RUW0(JW, II) = GWD_RANDO_RUWMAX
    241241             ENDDO
    242242          end DO
     
    318318          WWP(JW, :) = min(WWP(JW, :), MAX(0., &
    319319               SIGN(1., ZOP(JW, :) * ZOM(JW, :))) * ABS(ZOP(JW, :))**3 &
    320                / BV(:, LL + 1) * EXP(- ZH(:, LL + 1) / H0) * SAT**2 * KMIN**2 &
    321                / ZK(JW, :)**4)
     320               / BV(:, LL + 1) * EXP(- ZH(:, LL + 1) / H0) * GWD_RANDO_SAT**2 &
     321               * KMIN**2 / ZK(JW, :)**4)
    322322       end DO
    323323
Note: See TracChangeset for help on using the changeset viewer.