Changeset 511 for LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h
- Timestamp:
- Mar 24, 2004, 11:41:29 AM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h
r506 r511 1 1 IF (ok_mensuel) THEN 2 2 c 3 IF(region.EQ.'3d') THEN4 3 zsto = dtime 5 4 zout = dtime * ecrit_mth … … 10 9 zsto2 = dtime * radpas 11 10 PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout 12 ELSEIF(region.EQ.'1d') THEN13 zsto = dtime14 c zout = dtime15 zout = dtime *4.16 ENDIF17 11 c 18 12 idayref = day_ref … … 29 23 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 30 24 CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 31 c071103 32 c . 1,iim,1,jjmp1, 33 . imin_debut,nbpti, 34 . jmin_debut,nbptj, 35 c071103 36 . itau_phy, zjulian, dtime, 25 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 37 26 . nhori, nid_mth) 38 27 write(*,*)'Mensuel ', itau_phy, zjulian … … 96 85 c 97 86 c ENSEMBLES END 98 CALL histdef(nid_mth, "q2m", "Specific humidity 2m", " Kg/Kg",87 CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "kg/kg", 99 88 . iim,jjmp1,nhori, 1,1,1, -99, 32, 100 89 . "ave(X)", zsto,zout) … … 112 101 . "ave(X)", zsto,zout) 113 102 c 114 CALL histdef(nid_mth, "qsurf", "Surface Air humidity", " Kg/Kg",103 CALL histdef(nid_mth, "qsurf", "Surface Air humidity", "kg/kg", 115 104 . iim,jjmp1,nhori, 1,1,1, -99, 32, 116 105 . "ave(X)", zsto,zout) … … 128 117 c 129 118 CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 130 . " Kg/(s*m2)",119 . "kg/(s*m2)", 131 120 . iim,jjmp1,nhori, 1,1,1, -99, 32, 132 121 . "ave(X)", zsto,zout) 133 122 c 134 123 CALL histdef(nid_mth, "plul", "Large-scale Precip.", 135 . " Kg/(s*m2)",124 . "kg/(s*m2)", 136 125 . iim,jjmp1,nhori, 1,1,1, -99, 32, 137 126 . "ave(X)", zsto,zout) 138 127 c 139 128 CALL histdef(nid_mth, "pluc", "Convective Precip.", 140 . "Kg/(s*m2)", 141 . iim,jjmp1,nhori, 1,1,1, -99, 32, 142 . "ave(X)", zsto,zout) 143 c 144 CALL histdef(nid_mth, "snow", "Snow fall", "Kg/(s*m2)", 145 . iim,jjmp1,nhori, 1,1,1, -99, 32, 146 . "ave(X)", zsto,zout) 147 c 148 CALL histdef(nid_mth, "snow_mass", "Snow Mass", "Kg/m2", 149 . iim,jjmp1,nhori, 1,1,1, -99, 32, 150 . "ave(X)", zsto,zout) 151 c 152 cIM: 071003 evap(Kg/m2/s) ==> evap(mm/day) 153 CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)", 154 c CALL histdef(nid_mth, "evap", "Evaporation", "mm/day", 129 . "kg/(s*m2)", 130 . iim,jjmp1,nhori, 1,1,1, -99, 32, 131 . "ave(X)", zsto,zout) 132 c 133 CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)", 134 . iim,jjmp1,nhori, 1,1,1, -99, 32, 135 . "ave(X)", zsto,zout) 136 c 137 CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2", 138 . iim,jjmp1,nhori, 1,1,1, -99, 32, 139 . "ave(X)", zsto,zout) 140 c 141 CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)", 155 142 . iim,jjmp1,nhori, 1,1,1, -99, 32, 156 143 . "ave(X)", zsto,zout) … … 303 290 c 304 291 CALL histdef(nid_mth, "fqcalving","Ice Calving", 305 . " Kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32,292 . "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32, 306 293 . "ave(X)", zsto,zout) 307 294 cIM: 171003 … … 319 306 cIM: 171003 320 307 c 321 cIM if ( lev_histday.gt.1 ) then322 308 DO nsrf = 1, nbsrf 323 309 C … … 360 346 C 361 347 END DO 362 cIM endif !lev_histday363 348 c 364 349 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", … … 386 371 . "ave(X)", zsto,zout) 387 372 c 388 CALL histdef(nid_mth,"cldq","Cloud liquid water path"," Kg/m2",373 CALL histdef(nid_mth,"cldq","Cloud liquid water path","kg/m2", 389 374 . iim,jjmp1,nhori, 1,1,1, -99, 32, 390 375 . "ave(X)", zsto,zout) 391 376 cIM: 071003 392 CALL histdef(nid_mth,"lwp","Cloud water path"," Kg/m2",393 . iim,jjmp1,nhori, 1,1,1, -99, 32, 394 . "ave(X)", zsto,zout) 395 c 396 CALL histdef(nid_mth,"iwp","Cloud ice water path"," Kg/m2",377 CALL histdef(nid_mth,"lwp","Cloud water path","kg/m2", 378 . iim,jjmp1,nhori, 1,1,1, -99, 32, 379 . "ave(X)", zsto,zout) 380 c 381 CALL histdef(nid_mth,"iwp","Cloud ice water path","kg/m2", 397 382 . iim,jjmp1,nhori, 1,1,1, -99, 32, 398 383 . "ave(X)", zsto,zout) … … 417 402 IF(iflag_con.GE.3) THEN ! sb 418 403 c 419 CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/ Kg",404 CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/kg", 420 405 . iim,jjmp1,nhori, 1,1,1, -99, 32, 421 406 . "ave(X)", zsto,zout) … … 429 414 . "ave(X)", zsto,zout) 430 415 c 431 CALL histdef(nid_mth, "fbase", "Cld base mass flux", " Kg/m2/s",432 . iim,jjmp1,nhori, 1,1,1, -99, 32, 433 . "ave(X)", zsto,zout) 434 c 435 CALL histdef(nid_mth, "prw", "Precipitable water", " Kg/m2",416 CALL histdef(nid_mth, "fbase", "Cld base mass flux", "kg/m2/s", 417 . iim,jjmp1,nhori, 1,1,1, -99, 32, 418 . "ave(X)", zsto,zout) 419 c 420 CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2", 436 421 . iim,jjmp1,nhori, 1,1,1, -99, 32, 437 422 . "ave(X)", zsto,zout) … … 470 455 c 471 456 c 472 CALL histdef(nid_mth,"lwcon","Cloud water content"," Kg/Kg",457 CALL histdef(nid_mth,"lwcon","Cloud water content","kg/kg", 473 458 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 474 459 . "ave(X)", zsto,zout) 475 460 c 476 CALL histdef(nid_mth,"iwcon","Cloud ice water content"," Kg/Kg",461 CALL histdef(nid_mth,"iwcon","Cloud ice water content","kg/kg", 477 462 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 478 463 . "ave(X)", zsto,zout) … … 482 467 . "ave(X)", zsto,zout) 483 468 c 484 CALL histdef(nid_mth, "ovap", "Specific humidity", " Kg/Kg",485 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 486 . "ave(X)", zsto,zout) 487 c 488 c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio"," Kg/Kg",469 CALL histdef(nid_mth, "ovap", "Specific humidity", "kg/kg", 470 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 471 . "ave(X)", zsto,zout) 472 c 473 c CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","kg/kg", 489 474 c . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 490 475 c . "ave(X)", zsto,zout) … … 527 512 . "ave(X)", zsto,zout) 528 513 c 529 CALL histdef(nid_mth, "upwd", "saturated updraft", " Kg/m2/s",514 CALL histdef(nid_mth, "upwd", "saturated updraft", "kg/m2/s", 530 515 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 531 516 . "ave(X)", zsto,zout) … … 535 520 . "ave(X)", zsto,zout) 536 521 c 537 CALL histdef(nid_mth, "dqphy", "Physics dQ", " Kg/Kg/s",522 CALL histdef(nid_mth, "dqphy", "Physics dQ", "kg/kg/s", 538 523 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 539 524 . "ave(X)", zsto,zout) … … 579 564 CALL histdef(nid_mth, "clwcon", 580 565 . "Convective Cloud Liquid water content" 581 . , " Kg/Kg",582 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 583 . "ave(X)", zsto,zout) 584 c 585 CALL histdef(nid_mth,"Ma","undilute adiab updraft"," Kg/m2/s",586 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 587 . "ave(X)", zsto,zout) 588 c 589 CALL histdef(nid_mth, "dnwd", "saturated downdraft"," Kg/m2/s",590 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 591 . "ave(X)", zsto,zout) 592 c 593 CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", " Kg/m2/s",566 . , "kg/kg", 567 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 568 . "ave(X)", zsto,zout) 569 c 570 CALL histdef(nid_mth,"Ma","undilute adiab updraft","kg/m2/s", 571 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 572 . "ave(X)", zsto,zout) 573 c 574 CALL histdef(nid_mth, "dnwd", "saturated downdraft","kg/m2/s", 575 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 576 . "ave(X)", zsto,zout) 577 c 578 CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "kg/m2/s", 594 579 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 595 580 . "ave(X)", zsto,zout) … … 599 584 . "ave(X)", zsto,zout) 600 585 c 601 CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", " Kg/Kg/s",586 CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "kg/kg/s", 602 587 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 603 588 . "ave(X)", zsto,zout) … … 611 596 c . "ave(X)", zsto,zout) 612 597 c 613 CALL histdef(nid_mth, "dqcon", "Convection dQ", " Kg/Kg/s",598 CALL histdef(nid_mth, "dqcon", "Convection dQ", "kg/kg/s", 614 599 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 615 600 . "ave(X)", zsto,zout) … … 623 608 $ ,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout) 624 609 c 625 CALL histdef(nid_mth, "dqlsc", "Condensation dQ", " Kg/Kg/s",610 CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "kg/kg/s", 626 611 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 627 612 . "ave(X)", zsto,zout) … … 631 616 . "ave(X)", zsto,zout) 632 617 c 633 CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", " Kg/Kg/s",618 CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "kg/kg/s", 634 619 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 635 620 . "ave(X)", zsto,zout) … … 639 624 . "ave(X)", zsto,zout) 640 625 c 641 CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", " Kg/Kg/s",626 CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "kg/kg/s", 642 627 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 643 628 . "ave(X)", zsto,zout) … … 656 641 . "ave(X)", zsto,zout) 657 642 658 CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", " Kg/Kg/s",643 CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "kg/kg/s", 659 644 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 660 645 . "ave(X)", zsto,zout)
Note: See TracChangeset
for help on using the changeset viewer.
