Changeset 1675 for trunk/LMDZ.VENUS
- Timestamp:
- Mar 16, 2017, 9:51:39 AM (8 years ago)
- Location:
- trunk/LMDZ.VENUS
- Files:
-
- 5 edited
-
deftank/context_lmdz_physics.xml (modified) (6 diffs)
-
libf/phyvenus/flott_gwd_ran.F90 (modified) (1 diff)
-
libf/phyvenus/hgardfou.F (modified) (3 diffs)
-
libf/phyvenus/load_ksi.F (modified) (2 diffs)
-
libf/phyvenus/moldiff_red.F90 (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/deftank/context_lmdz_physics.xml
r1642 r1675 49 49 <!-- 2D variables --> 50 50 <field_group id="fields_2D" domain_ref="dom_glo"> 51 <field id="phis" 52 long_name="Surface geopotential (gz)" 53 unit="m2/s2" /> 54 <field id="aire" 55 long_name="Mesh area" 56 unit="m2" /> 57 <field id="tsol" 58 long_name="Surface Temperature" 59 unit="K" /> 51 60 <field id="psol" 52 61 long_name="Surface Pressure" 53 62 unit="Pa" /> 54 <field id="aire" 55 long_name="Mesh area" 56 unit="m2" /> 63 <field id="cdragh" 64 long_name="Drag coef on T" 65 unit="-" /> 66 <field id="cdragm" 67 long_name="Drag coef on U" 68 unit="-" /> 69 <field id="tops" 70 long_name="Solar rad. at TOA" 71 unit="W/m2" /> 72 <field id="topl" 73 long_name="IR rad. at TOA" 74 unit="W/m2" /> 75 <field id="sols" 76 long_name="Solar rad. at surface" 77 unit="W/m2" /> 78 <field id="soll" 79 long_name="IR rad. at surface" 80 unit="W/m2" /> 57 81 </field_group> 58 82 … … 62 86 long_name="Atmospheric temperature" 63 87 unit="K" /> 88 <field id="pres" 89 long_name="Atmospheric pressure" 90 unit="Pa" /> 91 <field id="geop" 92 long_name="Geopotential (gz)" 93 unit="m2/s2" /> 64 94 <field id="vitu" 65 95 long_name="Eastward Zonal Wind" … … 68 98 long_name="Northward Meridional Wind" 69 99 unit="m/s" /> 70 </field_group> 100 <field id="vitw" 101 long_name="Vertical wind" 102 unit="Pa/s" /> 103 <field id="Kz" 104 long_name="vertical diffusion coef" 105 unit="m2/s" /> 106 <field id="mmean" 107 long_name="Mean molecular mass" 108 unit="g/mol" /> 109 <field id="rho" 110 long_name="Atmospheric density" 111 unit="kg/m3" /> 112 <field id="dudyn" 113 long_name="Dynamics du" 114 unit="m/s2" /> 115 <field id="duvdf" 116 long_name="Boundary layer du" 117 unit="m/s2" /> 118 <field id="dvvdf" 119 long_name="Boundary layer dv" 120 unit="m/s2" /> 121 <field id="duajs" 122 long_name="Dry adj. du" 123 unit="m/s2" /> 124 <field id="dugwo" 125 long_name="GW oro du" 126 unit="m/s2" /> 127 <field id="dugwno" 128 long_name="GW non-oro du" 129 unit="m/s2" /> 130 <field id="dumolvis" 131 long_name="molec viscosity (u)" 132 unit="m/s2" /> 133 <field id="dvmolvis" 134 long_name="molec viscosity (v)" 135 unit="m/s2" /> 136 <field id="dtdyn" 137 long_name="Dynamics dT" 138 unit="K/s" /> 139 <field id="dtphy" 140 long_name="Physics dT" 141 unit="K/s" /> 142 <field id="dtvdf" 143 long_name="Boundary layer dT" 144 unit="K/s" /> 145 <field id="dtajs" 146 long_name="Dry adj. dT" 147 unit="K/s" /> 148 <field id="dtswr" 149 long_name="SW radiation dT" 150 unit="K/s" /> 151 <field id="dtswrNLTE" 152 long_name="SW radiation (NLTE) dT" 153 unit="K/s" /> 154 <field id="dtswrLTE" 155 long_name="SW radiation (LTE) dT" 156 unit="K/s" /> 157 <field id="dtlwr" 158 long_name="LW radiation dT" 159 unit="K/s" /> 160 <field id="dtlwrNLTE" 161 long_name="LW radiation (NLTE) dT" 162 unit="K/s" /> 163 <field id="dtlwrLTE" 164 long_name="LW radiation (LTE) dT" 165 unit="K/s" /> 166 <field id="dteuv" 167 long_name="EUV radiation dT" 168 unit="K/s" /> 169 <field id="dtcond" 170 long_name="Thermal conduction dT" 171 unit="K/s" /> 172 <field id="dtec" 173 long_name="Kinetic dissip dT" 174 unit="K/s" /> 175 <field id="SWnet" 176 long_name="Net SW flux" 177 unit="W/m2" /> 178 <field id="LWnet" 179 long_name="Net LW flux" 180 unit="W/m2" /> 181 <field id="fluxvdf" 182 long_name="PBL net flux" 183 unit="W/m2" /> 184 <field id="fluxdyn" 185 long_name="Dyn net flux" 186 unit="W/m2" /> 187 <field id="fluxajs" 188 long_name="Dry adj. net flux" 189 unit="W/m2" /> 190 <field id="fluxec" 191 long_name="Kinetic Energy net flux" 192 unit="W/m2" /> 193 <field id="co2" 194 long_name="CO2" 195 unit="mol/mol" /> 196 <field id="co" 197 long_name="CO" 198 unit="mol/mol" /> 199 <field id="h2" 200 long_name="H2" 201 unit="mol/mol" /> 202 <field id="h2o" 203 long_name="H2O" 204 unit="mol/mol" /> 205 <field id="o1d" 206 long_name="O(1D)" 207 unit="mol/mol" /> 208 <field id="o" 209 long_name="O" 210 unit="mol/mol" /> 211 <field id="o2" 212 long_name="O2" 213 unit="mol/mol" /> 214 <field id="o2dg" 215 long_name="O2(DG)" 216 unit="mol/mol" /> 217 <field id="o3" 218 long_name="O3" 219 unit="mol/mol" /> 220 <field id="h" 221 long_name="H" 222 unit="mol/mol" /> 223 <field id="oh" 224 long_name="OH" 225 unit="mol/mol" /> 226 <field id="ho2" 227 long_name="HO2" 228 unit="mol/mol" /> 229 <field id="h2o2" 230 long_name="H2O2" 231 unit="mol/mol" /> 232 <field id="cl" 233 long_name="Cl" 234 unit="mol/mol" /> 235 <field id="clo" 236 long_name="ClO" 237 unit="mol/mol" /> 238 <field id="cl2" 239 long_name="Cl2" 240 unit="mol/mol" /> 241 <field id="hcl" 242 long_name="HCl" 243 unit="mol/mol" /> 244 <field id="hocl" 245 long_name="HOCl" 246 unit="mol/mol" /> 247 <field id="clco" 248 long_name="ClCO" 249 unit="mol/mol" /> 250 <field id="clco3" 251 long_name="ClCO3" 252 unit="mol/mol" /> 253 <field id="cocl2" 254 long_name="COCl2" 255 unit="mol/mol" /> 256 <field id="s" 257 long_name="S" 258 unit="mol/mol" /> 259 <field id="so" 260 long_name="SO" 261 unit="mol/mol" /> 262 <field id="so2" 263 long_name="SO2" 264 unit="mol/mol" /> 265 <field id="so3" 266 long_name="SO3" 267 unit="mol/mol" /> 268 <field id="s2o2" 269 long_name="S2O2" 270 unit="mol/mol" /> 271 <field id="ocs" 272 long_name="OCS" 273 unit="mol/mol" /> 274 <field id="hso3" 275 long_name="HSO3" 276 unit="mol/mol" /> 277 <field id="h2so4" 278 long_name="H2SO4" 279 unit="mol/mol" /> 280 <field id="s2" 281 long_name="S2" 282 unit="mol/mol" /> 283 <field id="clso2" 284 long_name="ClSO2" 285 unit="mol/mol" /> 286 <field id="oscl" 287 long_name="OSCl" 288 unit="mol/mol" /> 289 <field id="n2" 290 long_name="N2" 291 unit="mol/mol" /> 292 <field id="h2oliq" 293 long_name="H2Oliq" 294 unit="mol/mol" /> 295 <field id="h2so4liq" 296 long_name="H2SO4liq" 297 unit="mol/mol" /> 298 <field id="d_qmoldifCO2" 299 long_name="Diff molec CO2" 300 unit="kg/kg" /> 301 <field id="d_qmoldifO3p" 302 long_name="Diff molec 0" 303 unit="kg/kg" /> 304 <field id="d_qmoldifN2" 305 long_name="Diff molec N2" 306 unit="kg/kg" /> 307 </field_group> 71 308 72 309 </field_definition> 73 310 74 <!-- Definition of output files --> 311 <!-- Definition of output files 312 output_freq: ts is physics time steps, d is day 313 enabled=".true." to output a file 314 --> 315 75 316 <file_definition> 76 <!-- Instantaneous outputs; Here output_freq inphysics time steps -->77 <file id=" histins"78 name="X histins"317 <!-- Instantaneous outputs; every physics time steps --> 318 <file id="ins" 319 name="Xins" 79 320 output_freq="1ts" 80 321 type="one_file" 81 enabled=". true.">322 enabled=".false."> 82 323 83 324 <!-- VARS 0D … … 90 331 <field_group operation="instant" 91 332 freq_op="1ts"> 333 <field field_ref="phis" operation="once" /> 92 334 <field field_ref="aire" operation="once" /> 335 <field field_ref="tops" /> 336 <field field_ref="tsol" /> 93 337 <field field_ref="psol" /> 94 338 </field_group> … … 98 342 freq_op="1ts"> 99 343 <field field_ref="temp" /> 344 <field field_ref="pres" /> 345 <field field_ref="geop" /> 100 346 <field field_ref="vitu" /> 101 347 <field field_ref="vitv" /> 348 <field field_ref="vitw" /> 102 349 </field_group> 103 350 </file> 104 351 105 <!-- diurnal averages outputs; output_freq is every day --> 352 <!-- Instantaneous outputs; 2000/Vd (if day_step=240000) --> 353 <file id="ins-2kperVd" 354 name="Xins-2kperVd" 355 output_freq="24ts" 356 type="one_file" 357 enabled=".false."> 358 359 <!-- VARS 2D --> 360 <field_group operation="instant" 361 freq_op="1ts"> 362 <field field_ref="phis" operation="once" /> 363 <field field_ref="aire" operation="once" /> 364 <field field_ref="tops" /> 365 <field field_ref="tsol" /> 366 <field field_ref="psol" /> 367 </field_group> 368 369 <!-- VARS 3D --> 370 <field_group operation="instant" 371 freq_op="1ts"> 372 <field field_ref="temp" /> 373 <field field_ref="pres" /> 374 <field field_ref="geop" /> 375 <field field_ref="vitu" /> 376 <field field_ref="vitv" /> 377 <field field_ref="vitw" /> 378 </field_group> 379 </file> 380 381 <!-- averages outputs; every Vday --> 106 382 <file id="diurnalave" 107 383 name="Xdiurnalave" 108 384 output_freq="1d" 109 385 type="one_file" 386 output_level="1" 387 enabled=".false."> 388 389 <!-- VARS 2D --> 390 <field_group operation="average" 391 freq_op="1ts"> 392 <field field_ref="phis" level="1" operation="once" /> 393 <field field_ref="aire" level="1" operation="once" /> 394 <field field_ref="tops" level="1" /> 395 <field field_ref="tsol" level="1" /> 396 <field field_ref="psol" level="1" /> 397 </field_group> 398 399 <!-- VARS 3D --> 400 <field_group operation="average" 401 freq_op="1ts"> 402 <field field_ref="temp" level="1" /> 403 <field field_ref="pres" level="1" /> 404 <field field_ref="geop" level="1" /> 405 <field field_ref="vitu" level="1" /> 406 <field field_ref="vitv" level="1" /> 407 <field field_ref="vitw" level="3" /> 408 <field field_ref="rho" level="3" /> 409 </field_group> 410 </file> 411 412 <!-- averages outputs; 1000/Vd (if day_step=240000) --> 413 <file id="ave-1kperVd" 414 name="Xave-1kperVd" 415 output_freq="48ts" 416 type="one_file" 417 output_level="1" 418 enabled=".false."> 419 420 <!-- VARS 2D --> 421 <field_group operation="average" 422 freq_op="1ts"> 423 <field field_ref="phis" level="1" operation="once" /> 424 <field field_ref="aire" level="1" operation="once" /> 425 <field field_ref="tops" level="1" /> 426 <field field_ref="tsol" level="1" /> 427 <field field_ref="psol" level="1" /> 428 </field_group> 429 430 <!-- VARS 3D --> 431 <field_group operation="average" 432 freq_op="1ts"> 433 <field field_ref="temp" level="1" /> 434 <field field_ref="pres" level="1" /> 435 <field field_ref="geop" level="1" /> 436 <field field_ref="vitu" level="1" /> 437 <field field_ref="vitv" level="1" /> 438 <field field_ref="vitw" level="3" /> 439 <field field_ref="rho" level="3" /> 440 </field_group> 441 </file> 442 443 <!-- averages outputs; 100/Vd (if day_step=240000) --> 444 <file id="ave-100perVd" 445 name="Xave-100perVd" 446 output_freq="480ts" 447 type="one_file" 448 output_level="1" 449 enabled=".false."> 450 451 <!-- VARS 2D --> 452 <field_group operation="average" 453 freq_op="1ts"> 454 <field field_ref="phis" level="1" operation="once" /> 455 <field field_ref="aire" level="1" operation="once" /> 456 <field field_ref="tops" level="1" /> 457 <field field_ref="tsol" level="1" /> 458 <field field_ref="psol" level="1" /> 459 </field_group> 460 461 <!-- VARS 3D --> 462 <field_group operation="average" 463 freq_op="1ts"> 464 <field field_ref="temp" level="1" /> 465 <field field_ref="pres" level="1" /> 466 <field field_ref="geop" level="1" /> 467 <field field_ref="vitu" level="1" /> 468 <field field_ref="vitv" level="1" /> 469 <field field_ref="dudyn" level="2" /> 470 <field field_ref="duvdf" level="2" /> 471 <field field_ref="duajs" level="2" /> 472 <field field_ref="dtdyn" level="2" /> 473 <field field_ref="dtvdf" level="2" /> 474 <field field_ref="dtajs" level="2" /> 475 <field field_ref="dtswr" level="2" /> 476 <field field_ref="dtlwr" level="2" /> 477 <field field_ref="vitw" level="3" /> 478 <field field_ref="rho" level="3" /> 479 <field field_ref="mmean" level="3" /> 480 </field_group> 481 </file> 482 483 <!-- averages outputs; 24/Vd (if day_step=240000) --> 484 <file id="ave-24perVd" 485 name="Xave-24perVd" 486 output_freq="2000ts" 487 type="one_file" 110 488 output_level="2" 111 489 enabled=".true."> 112 490 113 <!-- VARS 0D114 <field_group operation="average"115 freq_op="1ts">116 <field field_ref="ls" />117 </field_group> -->118 119 491 <!-- VARS 2D --> 120 492 <field_group operation="average" 121 493 freq_op="1ts"> 494 <field field_ref="phis" level="1" operation="once" /> 122 495 <field field_ref="aire" level="1" operation="once" /> 496 <field field_ref="tops" level="1" /> 497 <field field_ref="tsol" level="1" /> 123 498 <field field_ref="psol" level="1" /> 124 499 </field_group> … … 127 502 <field_group operation="average" 128 503 freq_op="1ts"> 129 <field field_ref="temp" level="1" /> 130 <field field_ref="vitu" level="1" /> 131 <field field_ref="vitv" level="1" /> 504 <field field_ref="temp" level="1" /> 505 <field field_ref="pres" level="1" /> 506 <field field_ref="geop" level="1" /> 507 <field field_ref="vitu" level="1" /> 508 <field field_ref="vitv" level="1" /> 509 <field field_ref="dudyn" level="2" /> 510 <field field_ref="duvdf" level="2" /> 511 <field field_ref="duajs" level="2" /> 512 <field field_ref="dtdyn" level="2" /> 513 <field field_ref="dtvdf" level="2" /> 514 <field field_ref="dtajs" level="2" /> 515 <field field_ref="dtswr" level="2" /> 516 <field field_ref="dtlwr" level="2" /> 517 <field field_ref="vitw" level="3" /> 518 <field field_ref="rho" level="3" /> 519 <field field_ref="mmean" level="3" /> 520 <field field_ref="co" level="4" /> 521 <field field_ref="h2o" level="4" /> 522 <field field_ref="o3" level="4" /> 523 <field field_ref="hcl" level="4" /> 524 <field field_ref="so" level="4" /> 525 <field field_ref="so2" level="4" /> 526 <field field_ref="ocs" level="4" /> 527 <field field_ref="h2so4" level="4" /> 528 <field field_ref="so3" level="5" /> 529 <field field_ref="s2" level="5" /> 530 <field field_ref="h2oliq" level="5" /> 531 <field field_ref="h2so4liq" level="5" /> 132 532 </field_group> 133 533 </file> 534 134 535 </file_definition> 135 536 -
trunk/LMDZ.VENUS/libf/phyvenus/flott_gwd_ran.F90
r1621 r1675 138 138 !(similar to observations) 139 139 KMIN = 1.E-5 ! Min horizontal wavenumber 140 ! KMIN = 6.3E-6 ! Min horizontal wavenumber 140 141 KMAX = 1.E-4 ! Max horizontal wavenumber 141 142 !Online output: one value only -
trunk/LMDZ.VENUS/libf/phyvenus/hgardfou.F
r1530 r1675 22 22 DATA firstcall /.TRUE./ 23 23 IF (firstcall) THEN 24 PRINT*, 'hgardfou garantit la temperature dans [ 20,1200] K'24 PRINT*, 'hgardfou garantit la temperature dans [15,1200] K' 25 25 firstcall = .FALSE. 26 26 ENDIF … … 29 29 DO k = 1, klev 30 30 DO i = 1, klon 31 !!!! MODIF GG to avoid crash after 78--> 95 extension!! 32 !!! WARNING: it has to be review/removed when the extension to the 33 !!! thermosphere is completed (physical processes and ionosphere added) 34 IF (k.LT.85) THEN 31 35 zt(i) = t(i,k) 36 ENDIF 32 37 ENDDO 33 38 #ifdef CRAY … … 49 54 ENDIF 50 55 #ifdef CRAY 51 CALL WHENFLT(klon, zt, 1, 20.0, jadrs, jbad)56 CALL WHENFLT(klon, zt, 1, 15.0, jadrs, jbad) 52 57 #else 53 58 jbad = 0 54 59 DO i = 1, klon 55 IF (zt(i).LT. 20.0) THEN60 IF (zt(i).LT.15.0) THEN 56 61 jbad = jbad + 1 57 62 jadrs(jbad) = i -
trunk/LMDZ.VENUS/libf/phyvenus/load_ksi.F
r1530 r1675 45 45 nlve = klev 46 46 47 cc GG modif below 48 c---------------------------------- 49 c Initialisation of values to 0 50 c (for all vertical levels) 51 c---------------------------------- 52 53 ksive(0:klev+1,0:klev+1,nnuve,nbmat) = 0.0 54 47 55 c ------------------------ 48 56 c Loading the ksi file … … 77 85 read(10,*) 78 86 read(10,*) m,Nb 79 if (m.ne.nlve) then 87 cc GG modif below 88 if (nlve.le.78.and.m.ne.nlve) then 80 89 write(*,*) 'This is subroutine load_ksi' 81 90 print*,'Dimension problem between ksi.txt and nlve' -
trunk/LMDZ.VENUS/libf/phyvenus/moldiff_red.F90
r1621 r1675 1462 1462 1463 1463 compteur=0 1464 do while ( pnew2 .ge. pp(il))1464 do while (Pnew2 .ge. pp(il)) 1465 1465 compteur=compteur+1 1466 1466 do nn=1,nq … … 1474 1474 if (compteur .ge. 100000) then 1475 1475 print*,'error moldiff_red infinite loop' 1476 print*,ig,il,pp(il),tt(nl), pnew2,qnew(il,:),Znew21476 print*,ig,il,pp(il),tt(nl),Pnew2,qnew(il,:),Znew2 1477 1477 stop 1478 1478 endif … … 1564 1564 1565 1565 compteur=0 1566 do while ( pnew2 .ge. pp(il))1566 do while (Pnew2 .ge. pp(il)) 1567 1567 compteur=compteur+1 1568 1568 do nn=1,nq … … 1576 1576 if (compteur .ge. 100000) then 1577 1577 print*,'pb moldiff_red infinite loop' 1578 print*,ig,nl,T(nl), pnew2,qnew(il,:),Znew21578 print*,ig,nl,T(nl),Pnew2,qnew(il,:),Znew2 1579 1579 stop 1580 1580 endif
Note: See TracChangeset
for help on using the changeset viewer.
