Changeset 2403
- Timestamp:
- Jul 9, 2020, 11:10:48 AM (4 years ago)
- Location:
- trunk/LMDZ.GENERIC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/README
r2354 r2403 1587 1587 Major cleanup: remove obsolete compilation scripts (makegcm*) and old dynamical 1588 1588 core, as it is obsolete with respect to the one provide in LMDZ.COMMON. 1589 1590 1591 == 09/07/2020 == DB 1592 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) 1593 -
trunk/LMDZ.GENERIC/deftank/def_generic_giants/PROTO_Saturn_Article/z2sig.def
r1164 r2403 19 19 179.0 20 20 191.0 21 203.0 22 215.0 23 227.0 24 239.0 21 201.0 22 211.0 23 221.0 24 231.0 25 241.0 25 26 251.0 26 263.0 27 259.0 28 267.0 27 29 275.0 28 287.0 30 283.0 31 291.0 29 32 299.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 33 307.0 34 313.0 35 319.0 36 325.0 37 331.0 38 337.0 39 343.0 40 349.0 41 355.0 42 361.0 43 367.0 44 373.0 45 379.0 46 385.0 47 391.0 48 397.0 49 403.0 50 409.0 51 415.0 52 421.0 53 427.0 54 433.0 55 439.0 56 445.0 57 451.0 58 457.0 59 463.0 60 469.0 61 475.0 62 481.0 63 487.0 64 493.0 65 499.0 66 505.0 67 511.0 68 517.0 69 523.0 70 529.0 71 535.0 72 541.0 73 547.0 74 553.0 75 559.0 76 571.0 77 577.0 78 583.0 79 589.0 80 595.0 81 601.0 82 607.0 83 613.0 84 619.0 85 625.0 86 631.0 87 637.0 88 643.0 89 649.0 90 655.0 91 661.0 92 667.0 93 673.0 94 679.0 95 685.0 96 691.0 62 97 695.0 63 98 707.0 -
trunk/LMDZ.GENERIC/deftank/saturn1d/profile
r1150 r2403 64 64 130 65 65 130 66 130 67 130 68 130 69 130 70 130 71 130 72 130 73 130 74 130 75 130 76 130 77 130 78 130 79 130 80 130 81 130 82 130 83 130 84 130 85 130 86 130 87 130 88 130 89 130 90 130 91 130 92 130 93 130 94 130 95 130 96 130 97 130 98 130 99 130 100 130 101 130 -
trunk/LMDZ.GENERIC/libf/phystd/nonoro_gwd_ran_mod.F90
r2305 r2403 49 49 gwd_convective_source 50 50 use inifis_mod 51 use geometry_mod, only: cell_area 51 52 #ifdef CPP_XIOS 52 53 use xios_output_mod, only: send_xios_field … … 97 98 INTEGER, parameter:: nw = nk * np *no ! Total numbers of gravity waves 98 99 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 103 107 REAL, parameter:: cmin = 1. ! Min horizontal absolute phase velocity 104 108 REAL cpha ! absolute phase velocity frequency … … 117 121 REAL epflux_0(nw, ngrid) ! Fluxes at launching level (previous name: RUW0) 118 122 INTEGER launch ! Launching altitude 119 REAL, parameter:: xlaunch = 0. 4! Control the launching altitude120 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.) 121 125 122 126 REAL prec(ngrid) ! precipitations … … 131 135 REAL H0bis(ngrid, nlayer) ! characteristic height of the atmosphere 132 136 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] 135 142 REAL zh(ngrid, nlayer + 1) ! Log-pressure altitude (constant H0) 136 143 REAL zhbis(ngrid, nlayer + 1) ! Log-pressure altitude (varying H) … … 237 244 + 1.) * pi / 2. 238 245 ! 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 240 254 ! horizontal phase speed 241 255 cpha = 0.
Note: See TracChangeset
for help on using the changeset viewer.