Changeset 2072 for LMDZ5/trunk/libf/phylmd
- Timestamp:
- Jun 23, 2014, 4:53:24 PM (10 years ago)
- Location:
- LMDZ5/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/phylmd/YOEGWD.h
r1938 r2072 14 14 real GHMAX,GRAHILO,GSIGCR,GSSEC,GTSEC,GVSEC 15 15 16 REAL RUWMAX16 REAL GWD_RANDO_RUWMAX 17 17 ! Maximum Eliassen-Palm flux at launch level, in "FLOTT_GWD_rando" 18 18 19 REAL SAT ! saturation parameter in "FLOTT_GWD_rando"19 REAL GWD_RANDO_SAT ! saturation parameter in "FLOTT_GWD_rando" 20 20 ! S_c in equation (12) of Lott (JGR, vol 118, page 8897, 2013) 21 21 … … 23 23 COMMON/YOEGWD/ GFRCRIT,GKWAKE,GRCRIT,GVCRIT,GKDRAG,GKLIFT, & 24 24 & GHMAX,GRAHILO,GSIGCR,NKTOPG,NSTRA,GSSEC,GTSEC,GVSEC, & 25 & RUWMAX,sat25 & GWD_RANDO_RUWMAX, gwd_rando_sat 26 26 27 27 save /YOEGWD/ -
LMDZ5/trunk/libf/phylmd/conf_phys_m.F90
r2057 r2072 186 186 LOGICAL,SAVE :: ok_strato_omp 187 187 LOGICAL,SAVE :: ok_hines_omp, ok_gwd_rando_omp 188 real, SAVE:: RUWMAX_omp,sat_omp188 real, SAVE:: GWD_RANDO_RUWMAX_omp, gwd_rando_sat_omp 189 189 LOGICAL,SAVE :: carbon_cycle_tr_omp 190 190 LOGICAL,SAVE :: carbon_cycle_cpl_omp … … 1687 1687 CALL getin('ok_gwd_rando', ok_gwd_rando_omp) 1688 1688 1689 RUWMAX_omp = 0.011690 CALL getin(' ruwmax',RUWMAX_omp)1691 1692 sat_omp = 0.251693 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) 1694 1694 1695 1695 !Config Key = OK_LES … … 1917 1917 ok_hines = ok_hines_omp 1918 1918 ok_gwd_rando = ok_gwd_rando_omp 1919 RUWMAX =RUWMAX_omp1920 sat =sat_omp1919 GWD_RANDO_RUWMAX = GWD_RANDO_RUWMAX_omp 1920 gwd_rando_sat = gwd_rando_sat_omp 1921 1921 ok_LES = ok_LES_omp 1922 1922 callstats = callstats_omp … … 2109 2109 write(lunout,*) 'ok_hines = ', ok_hines 2110 2110 write(lunout,*) 'ok_gwd_rando = ', ok_gwd_rando 2111 write(lunout,*) ' RUWMAX = ',RUWMAX2112 write(lunout,*) ' sat = ',sat2111 write(lunout,*) 'GWD_RANDO_RUWMAX = ', GWD_RANDO_RUWMAX 2112 write(lunout,*) 'gwd_rando_sat = ', gwd_rando_sat 2113 2113 write(lunout,*) 'read_climoz = ', read_climoz 2114 2114 write(lunout,*) 'carbon_cycle_tr = ', carbon_cycle_tr -
LMDZ5/trunk/libf/phylmd/flott_gwd_rando_m.F90
r1938 r2072 238 238 + ZK(JW, II) * SIN(ZP(JW, II)) * VH(II, LAUNCH) 239 239 ! Momentum flux at launch lev 240 RUW0(JW, II) = RUWMAX240 RUW0(JW, II) = GWD_RANDO_RUWMAX 241 241 ENDDO 242 242 end DO … … 318 318 WWP(JW, :) = min(WWP(JW, :), MAX(0., & 319 319 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) 322 322 end DO 323 323
Note: See TracChangeset
for help on using the changeset viewer.