Changeset 511 for LMDZ.3.3/branches
- Timestamp:
- Mar 24, 2004, 11:41:29 AM (21 years ago)
- Location:
- LMDZ.3.3/branches/rel-LF/libf/phylmd
- Files:
-
- 7 edited
-
ini_histISCCP.h (modified) (3 diffs)
-
ini_histday.h (modified) (2 diffs)
-
ini_histhf.h (modified) (3 diffs)
-
ini_histhf3d.h (modified) (3 diffs)
-
ini_histins.h (modified) (2 diffs)
-
ini_histmth.h (modified) (23 diffs)
-
ini_histmthNMC.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histISCCP.h
r500 r511 1 1 IF (ok_isccp) THEN 2 2 c 3 IF (region.EQ.'3d') THEN !region4 3 c pour les champs instantannes, il faut mettre la meme valeur pour 5 4 c zout et zsto. … … 18 17 c ecriture mensuelle 19 18 zout = dtime * ecrit_mth 20 c21 ELSE IF (region.EQ.'1d') THEN !region22 zsto = dtime23 zout = dtime *4.24 c zout = dtime25 ENDIF26 19 c 27 20 print*,'ISCCP zout,zsto=',zout,zsto … … 45 38 c 46 39 CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:), 47 . imin_debut,nbpti, 48 . jmin_debut,nbptj, 40 . 1, iim, 1, jjmp1, 49 41 . itau_phy, zjulian, dtime, 50 42 . nhori, nid_isccp) -
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histday.h
r506 r511 1 1 IF (ok_journe) THEN 2 2 c 3 IF(region.EQ.'3d') THEN4 3 zsto = dtime 5 c zout = dtime * FLOAT(ecrit_day) 6 c zsto1= dtime * FLOAT(ecrit_day) 7 zout = dtime * REAL(ecrit_day) 8 zsto1= dtime * REAL(ecrit_day) 9 ELSEIF(region.EQ.'1d') THEN 10 zsto = dtime 11 c zout = dtime 12 zout = dtime *4. 13 ENDIF 4 zout = dtime * FLOAT(ecrit_day) 5 zsto1= dtime * FLOAT(ecrit_day) 6 c zout = dtime * REAL(ecrit_day) 7 c zsto1= dtime * REAL(ecrit_day) 14 8 c 15 9 idayref = day_ref … … 25 19 ENDDO 26 20 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 27 write(*,*)'zx_lon = ',zx_lon(:,1)28 write(*,*)'zx_lat = ',zx_lat(1,:)29 21 CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 30 c071103 31 c . 1,iim,1,jjmp1, 32 . imin_debut,nbpti, 33 . jmin_debut,nbptj, 34 c071103 35 . itau_phy, zjulian, dtime, 22 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 36 23 . nhori, nid_day) 37 24 write(*,*)'Journee ', itau_phy, zjulian -
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf.h
r506 r511 1 1 2 2 IF (ok_hf) THEN 3 c4 IF(region.EQ.'3d') THEN5 3 c 6 4 zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) … … 9 7 PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf 10 8 c 11 ELSEIF(region.EQ.'1d') THEN12 c13 c pour les champs instantannes, il faut mettre la meme valeur pour14 c zout et tsto.15 c dtime est passe par ailleurs a histbeg16 zsto = dtime17 c zout = dtime18 zout = dtime *4.19 ENDIF20 print*,'zout,zsto=',zout,zsto21 cIM cf LF22 9 idayref = day_ref 23 10 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 24 25 cccIM CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)26 c CALL ymds2ju(annee_ref, 1, 1, 0.0, zjulian)27 c zjulian = zjulian + day_ini28 cIM cf LF29 11 30 12 c … … 39 21 cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat, 40 22 CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 41 cIM cf LF . 1,iim,1,jjmp1, 0, zjulian, dtime, 42 c071103 43 c . 1,iim,1,jjmp1, 44 . imin_debut,nbpti, 45 . jmin_debut,nbptj, 46 c071103 47 . itau_phy, zjulian, dtime, 23 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 48 24 . nhori, nid_hf) 49 25 -
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf3d.h
r510 r511 1 1 2 2 c IF (ok_hf) THEN 3 c4 IF(region.EQ.'3d') THEN5 3 c 6 4 zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) … … 9 7 PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf 10 8 c 11 ELSEIF(region.EQ.'1d') THEN12 c13 c pour les champs instantannes, il faut mettre la meme valeur pour14 c zout et tsto.15 c dtime est passe par ailleurs a histbeg16 zsto = dtime17 c zout = dtime18 zout = dtime *4.19 ENDIF20 print*,'zout,zsto=',zout,zsto21 9 cIM cf LF 22 10 idayref = day_ref … … 39 27 cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat, 40 28 CALL histbeg("histhf3d", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 41 cIM cf LF . 1,iim,1,jjmp1, 0, zjulian, dtime, 42 c071103 43 c . 1,iim,1,jjmp1, 44 . imin_debut,nbpti, 45 . jmin_debut,nbptj, 46 c071103 47 . itau_phy, zjulian, dtime, 29 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 48 30 . nhori, nid_hf3d) 49 31 -
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histins.h
r506 r511 1 1 IF (ok_instan) THEN 2 2 c 3 IF(region.EQ.'3d') THEN4 3 zsto = dtime * ecrit_ins 5 4 zout = dtime * ecrit_ins 6 ELSEIF(REGION.eq.'1d') THEN7 zsto = dtime8 c zout = dtime9 zout = dtime *4.10 ENDIF11 5 c 12 6 idayref = day_ref … … 23 17 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 24 18 CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 25 c071103 26 c . 1,iim,1,jjmp1, 27 . imin_debut,nbpti, 28 . jmin_debut,nbptj, 29 c071103 30 . itau_phy, zjulian, dtime, 19 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 31 20 . nhori, nid_ins) 32 21 write(*,*)'Inst ', itau_phy, zjulian -
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) -
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmthNMC.h
r508 r511 1 c IF (ok_mensuelNMC) THEN2 1 IF (ok_mensuel) THEN 3 2 c 4 IF(region.EQ.'3d') THEN5 3 zsto = dtime 6 cIM 150304 BEG7 4 zout = dtime * ecrit_mth 8 c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))9 cIM 150304 END10 cIM 120304 BEG11 c zout = dtime * REAL(NINT(86400./dtime*ecrit_hf))12 c zsto = zout13 c zsto = dtime14 cIM 120304 END15 ELSEIF(region.EQ.'1d') THEN16 zsto = dtime17 c zout = dtime18 zout = dtime *4.19 ENDIF20 5 c 21 6 idayref = day_ref … … 32 17 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 33 18 CALL histbeg("histNMC.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 34 c071103 35 c . 1,iim,1,jjmp1, 36 . imin_debut,nbpti, 37 . jmin_debut,nbptj, 38 c071103 39 . itau_phy, zjulian, dtime, 19 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 40 20 . nhori, nid_nmc) 41 21 write(*,*)'Mensuel NMC ', itau_phy, zjulian
Note: See TracChangeset
for help on using the changeset viewer.
