Changeset 506 for LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h
- Timestamp:
- Mar 24, 2004, 9:43:54 AM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histmth.h
r486 r506 1 1 IF (ok_mensuel) THEN 2 c 3 IF(region.EQ.'3d') THEN 4 zsto = dtime 5 zout = dtime * ecrit_mth 6 c 7 c zsto1: pour des valeurs "instantannees" mensuelles 8 zsto1 = dtime * ecrit_mth 9 c zsto2: pour des flux radiatifs calcules tous les 2 heures 10 zsto2 = dtime * radpas 11 PRINT*,' zsto,zsto1,zsto2,zout=',zsto, zsto1, zsto2,zout 12 ELSEIF(region.EQ.'1d') THEN 13 zsto = dtime 14 c zout = dtime 15 zout = dtime *4. 16 ENDIF 2 17 c 3 18 idayref = day_ref … … 14 29 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 15 30 CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 16 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 31 c071103 32 c . 1,iim,1,jjmp1, 33 . imin_debut,nbpti, 34 . jmin_debut,nbptj, 35 c071103 36 . itau_phy, zjulian, dtime, 17 37 . nhori, nid_mth) 18 38 write(*,*)'Mensuel ', itau_phy, zjulian … … 22 42 c . klev, znivsig, nvert) 23 43 c 24 zsto = dtime 25 zout = dtime * ecrit_mth 26 c 44 c 45 IF(lev_histmth.GE.1) THEN 27 46 CALL histdef(nid_mth, "phis", "Surface geop. height", "-", 28 47 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 33 52 . "once", zsto,zout) 34 53 c 54 call histdef(nid_mth, "pourc_"//clnsurf(is_ter), 55 $ "Fraction "//clnsurf(is_ter), "W/m2", 56 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 57 $ "once", zsto,zout) 58 c 35 59 c Champs 2D: 36 60 c … … 46 70 . iim,jjmp1,nhori, 1,1,1, -99, 32, 47 71 . "ave(X)", zsto,zout) 48 c 72 c ENSEMBLES BEG 73 CALL histdef(nid_mth, "t2m_min", "Temp. 2m min.", 74 . "K", 75 . iim,jjmp1,nhori, 1,1,1, -99, 32, 76 . t2mincels, zsto,zout) 77 c 78 CALL histdef(nid_mth, "t2m_max", "Temp. 2m max.", 79 . "K", 80 . iim,jjmp1,nhori, 1,1,1, -99, 32, 81 . t2maxcels, zsto,zout) 82 c 83 c CALL histdef(nid_mth, "tsoil", "Sfce soil Temperature", 84 c . "K", 85 c . iim,jjmp1,nhori, 1,1,1, -99, 32, 86 c . "ave(X)", zsto,zout) 87 c 88 CALL histdef(nid_mth, "wind10m","10-m wind speed","m/s", 89 . iim,jjmp1,nhori, 1,1,1, -99, 32, 90 . "ave(X)", zsto,zout) 91 c 92 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-", 93 . iim,jjmp1,nhori, 1,1,1, -99, 32, 94 . "ave(X)", zsto,zout) 95 96 c 97 c ENSEMBLES END 49 98 CALL histdef(nid_mth, "q2m", "Specific humidity 2m", "Kg/Kg", 50 99 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 58 107 . iim,jjmp1,nhori, 1,1,1, -99, 32, 59 108 . "ave(X)", zsto,zout) 60 c61 109 c 62 110 CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", … … 74 122 endif 75 123 c 124 CALL histdef(nid_mth, "ndayrain", 125 . "Number of day with rain (liq+sol)", "-", 126 . iim,jjmp1,nhori, 1,1,1, -99, 32, 127 . "inst(X)", zsto1,zout) 128 c 76 129 CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 77 . " kg/(s*m2)",130 . "Kg/(s*m2)", 78 131 . iim,jjmp1,nhori, 1,1,1, -99, 32, 79 132 . "ave(X)", zsto,zout) 80 133 c 81 134 CALL histdef(nid_mth, "plul", "Large-scale Precip.", 82 . " kg/(s*m2)",135 . "Kg/(s*m2)", 83 136 . iim,jjmp1,nhori, 1,1,1, -99, 32, 84 137 . "ave(X)", zsto,zout) 85 138 c 86 139 CALL histdef(nid_mth, "pluc", "Convective Precip.", 87 . "kg/(s*m2)", 88 . iim,jjmp1,nhori, 1,1,1, -99, 32, 89 . "ave(X)", zsto,zout) 90 c 91 CALL histdef(nid_mth, "snow", "Snow fall", "kg/(s*m2)", 92 . iim,jjmp1,nhori, 1,1,1, -99, 32, 93 . "ave(X)", zsto,zout) 94 c 95 CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2", 96 . iim,jjmp1,nhori, 1,1,1, -99, 32, 97 . "ave(X)", zsto,zout) 98 c 99 CALL histdef(nid_mth, "evap", "Evaporation", "kg/(s*m2)", 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", 100 155 . iim,jjmp1,nhori, 1,1,1, -99, 32, 101 156 . "ave(X)", zsto,zout) … … 105 160 . "ave(X)", zsto,zout) 106 161 c 162 CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2", 163 . iim,jjmp1,nhori, 1,1,1, -99, 32, 164 . "ave(X)", zsto,zout) 165 c 107 166 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 108 167 . iim,jjmp1,nhori, 1,1,1, -99, 32, 109 168 . "ave(X)", zsto,zout) 110 169 c 170 CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2", 171 . iim,jjmp1,nhori, 1,1,1, -99, 32, 172 . "ave(X)", zsto,zout) 173 c 174 CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2", 175 . iim,jjmp1,nhori, 1,1,1,-99, 176 . 32, "ave(X)", zsto2,zout) 177 c 178 CALL histdef(nid_mth, "SWupTOAclr", 179 . "SWup clear sky at TOA","W/m2", 180 . iim,jjmp1,nhori, 1,1,1,-99, 181 . 32, "ave(X)", zsto2,zout) 182 c 183 CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2", 184 . iim,jjmp1,nhori, 1,1,1,-99, 185 . 32, "ave(X)", zsto2,zout) 186 c 187 CALL histdef(nid_mth, "SWdnTOAclr", 188 . "SWdn clear sky at TOA","W/m2", 189 . iim,jjmp1,nhori, 1,1,1,-99, 190 . 32, "ave(X)", zsto2,zout) 191 c 192 CALL histdef(nid_mth, "SWup200", "SWup at 200hPa","W/m2", 193 . iim,jjmp1,nhori, 1,1,1,-99, 194 . 32, "ave(X)", zsto2,zout) 195 c 196 CALL histdef(nid_mth, "SWup200clr", 197 . "SWup clear sky at 200hPa","W/m2", 198 . iim,jjmp1,nhori, 1,1,1,-99, 199 . 32, "ave(X)", zsto2,zout) 200 c 201 CALL histdef(nid_mth, "SWdn200", "SWdn at 200hPa","W/m2", 202 . iim,jjmp1,nhori, 1,1,1,-99, 203 . 32, "ave(X)", zsto2,zout) 204 c 205 CALL histdef(nid_mth, "SWdn200clr", 206 . "SWdn clear sky at 200hPa","W/m2", 207 . iim,jjmp1,nhori, 1,1,1,-99, 208 . 32, "ave(X)", zsto2,zout) 209 c 210 CALL histdef(nid_mth, "LWup200", "LWup at 200hPa","W/m2", 211 . iim,jjmp1,nhori, 1,1,1,-99, 212 . 32, "ave(X)", zsto2,zout) 213 c 214 CALL histdef(nid_mth, "LWup200clr", 215 . "LWup clear sky at 200hPa","W/m2", 216 . iim,jjmp1,nhori, 1,1,1,-99, 217 . 32, "ave(X)", zsto2,zout) 218 c 219 CALL histdef(nid_mth, "LWdn200", "LWdn at 200hPa","W/m2", 220 . iim,jjmp1,nhori, 1,1,1,-99, 221 . 32, "ave(X)", zsto2,zout) 222 c 223 CALL histdef(nid_mth, "LWdn200clr", 224 . "LWdn clear sky at 200hPa","W/m2", 225 . iim,jjmp1,nhori, 1,1,1,-99, 226 . 32, "ave(X)", zsto2,zout) 227 c 111 228 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 112 229 . iim,jjmp1,nhori, 1,1,1, -99, 32, 113 230 . "ave(X)", zsto,zout) 114 231 c 232 CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2", 233 . iim,jjmp1,nhori, 1,1,1, -99, 32, 234 . "ave(X)", zsto,zout) 235 c 115 236 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 116 237 . iim,jjmp1,nhori, 1,1,1, -99, 32, 117 238 . "ave(X)", zsto,zout) 118 cIM: 071003 239 c 240 CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2", 241 . iim,jjmp1,nhori, 1,1,1, -99, 32, 242 . "ave(X)", zsto,zout) 243 c 244 CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2", 245 . iim,jjmp1,nhori, 1,1,1,-99, 246 . 32, "ave(X)", zsto2,zout) 247 c 248 CALL histdef(nid_mth, "SWupSFCclr", 249 . "SWup clear sky at surface","W/m2", 250 . iim,jjmp1,nhori, 1,1,1,-99, 251 . 32, "ave(X)", zsto2,zout) 252 c 253 CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2", 254 . iim,jjmp1,nhori, 1,1,1,-99, 255 . 32, "ave(X)", zsto2,zout) 256 c 257 CALL histdef(nid_mth, "SWdnSFCclr", 258 . "SWdn clear sky at surface","W/m2", 259 . iim,jjmp1,nhori, 1,1,1,-99, 260 . 32, "ave(X)", zsto2,zout) 261 c 262 CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 263 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 264 . "ave(X)", zsto2,zout) 265 c 119 266 CALL histdef(nid_mth,"LWdnSFC","Down. IR rad. at surface", 120 267 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 121 . "ave(X)", zsto,zout) 122 cIM: 071003 123 CALL histdef(nid_mth,"LWupSFC","Upwd. IR rad. at surface", 268 . "ave(X)", zsto2,zout) 269 c 270 CALL histdef(nid_mth,"LWupSFCclr", 271 . "CS Upwd. IR rad. at surface", 124 272 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 125 . "ave(X)", zsto,zout) 126 c 127 CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2", 128 . iim,jjmp1,nhori, 1,1,1, -99, 32, 129 . "ave(X)", zsto,zout) 130 c 131 CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2", 132 . iim,jjmp1,nhori, 1,1,1, -99, 32, 133 . "ave(X)", zsto,zout) 134 c 135 CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2", 136 . iim,jjmp1,nhori, 1,1,1, -99, 32, 137 . "ave(X)", zsto,zout) 138 c 139 CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2", 140 . iim,jjmp1,nhori, 1,1,1, -99, 32, 141 . "ave(X)", zsto,zout) 273 . "ave(X)", zsto2,zout) 274 c 275 CALL histdef(nid_mth,"LWdnSFCclr", 276 . "Down. CS IR rad. at surface", 277 . "W/m2", iim,jjmp1,nhori, 1,1,1, -99, 32, 278 . "ave(X)", zsto2,zout) 142 279 c 143 280 CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2", … … 166 303 c 167 304 CALL histdef(nid_mth, "fqcalving","Ice Calving", 168 . "kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32, 169 . "ave(X)", zsto,zout) 170 c 305 . "Kg/m2/s",iim,jjmp1,nhori, 1,1,1, -99, 32, 306 . "ave(X)", zsto,zout) 307 cIM: 171003 308 DO nsrf = 1, nbsrf 309 call histdef(nid_mth, "taux_"//clnsurf(nsrf), 310 $ "Zonal wind stress"//clnsurf(nsrf), "Pa", 311 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 312 $ "ave(X)", zsto,zout) 313 314 call histdef(nid_mth, "tauy_"//clnsurf(nsrf), 315 $ "Meridional xind stress "//clnsurf(nsrf), "Pa", 316 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 317 $ "ave(X)", zsto,zout) 318 ENDDO 319 cIM: 171003 320 c 321 cIM if ( lev_histday.gt.1 ) then 171 322 DO nsrf = 1, nbsrf 172 323 C 324 IF(nsrf.GT.1) THEN 173 325 call histdef(nid_mth, "pourc_"//clnsurf(nsrf), 174 326 $ "Fraction "//clnsurf(nsrf), "W/m2", 175 327 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 176 328 $ "ave(X)", zsto,zout) 329 ENDIF !nsrf.GT.1 177 330 C 178 331 call histdef(nid_mth, "tsol_"//clnsurf(nsrf), … … 194 347 $ "LW "//clnsurf(nsrf), "W/m2", 195 348 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 196 $ "ave(X)", zsto ,zout)349 $ "ave(X)", zsto2,zout) 197 350 c 198 351 call histdef(nid_mth, "fsw_"//clnsurf(nsrf), 199 352 $ "SW "//clnsurf(nsrf), "W/m2", 200 353 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 201 $ "ave(X)", zsto ,zout)354 $ "ave(X)", zsto2,zout) 202 355 C 203 356 call histdef(nid_mth, "wbils_"//clnsurf(nsrf), … … 206 359 $ "ave(X)", zsto,zout) 207 360 C 208 call histdef(nid_mth, "taux_"//clnsurf(nsrf), 209 $ "Zonal wind stress"//clnsurf(nsrf), "Pa", 210 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 211 $ "ave(X)", zsto,zout) 212 213 call histdef(nid_mth, "tauy_"//clnsurf(nsrf), 214 $ "Meridional xind stress "//clnsurf(nsrf), "Pa", 215 $ iim,jjmp1,nhori, 1,1,1, -99, 32, 216 $ "ave(X)", zsto,zout) 361 END DO 362 cIM endif !lev_histday 363 c 364 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 365 . iim,jjmp1,nhori, 1,1,1, -99, 32, 366 . "ave(X)", zsto,zout) 367 c 368 CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 369 . iim,jjmp1,nhori, 1,1,1, -99, 32, 370 . "ave(X)", zsto,zout) 371 c 372 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-", 373 . iim,jjmp1,nhori, 1,1,1, -99, 32, 374 . "ave(X)", zsto,zout) 375 c 376 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-", 377 . iim,jjmp1,nhori, 1,1,1, -99, 32, 378 . "ave(X)", zsto,zout) 379 c 380 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-", 381 . iim,jjmp1,nhori, 1,1,1, -99, 32, 382 . "ave(X)", zsto,zout) 383 c 384 CALL histdef(nid_mth, "cldt", "Total cloudiness", "-", 385 . iim,jjmp1,nhori, 1,1,1, -99, 32, 386 . "ave(X)", zsto,zout) 387 c 388 CALL histdef(nid_mth,"cldq","Cloud liquid water path","Kg/m2", 389 . iim,jjmp1,nhori, 1,1,1, -99, 32, 390 . "ave(X)", zsto,zout) 391 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", 397 . iim,jjmp1,nhori, 1,1,1, -99, 32, 398 . "ave(X)", zsto,zout) 399 cIM: 071003 400 c 401 CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 402 . iim,jjmp1,nhori, 1,1,1, -99, 32, 403 . "ave(X)", zsto,zout) 404 c 405 CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 406 . iim,jjmp1,nhori, 1,1,1, -99, 32, 407 . "ave(X)", zsto,zout) 408 c 409 CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-", 410 . iim,jjmp1,nhori, 1,1,1, -99, 32, 411 . "ave(X)", zsto,zout) 412 c 413 CALL histdef(nid_mth, "vq", "Merid humidity transport", "-", 414 . iim,jjmp1,nhori, 1,1,1, -99, 32, 415 . "ave(X)", zsto,zout) 416 cKE43 417 IF(iflag_con.GE.3) THEN ! sb 418 c 419 CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg", 420 . iim,jjmp1,nhori, 1,1,1, -99, 32, 421 . "ave(X)", zsto,zout) 422 c 423 CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa", 424 . iim,jjmp1,nhori, 1,1,1, -99, 32, 425 . "ave(X)", zsto,zout) 426 c 427 CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa", 428 . iim,jjmp1,nhori, 1,1,1, -99, 32, 429 . "ave(X)", zsto,zout) 430 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", 436 . iim,jjmp1,nhori, 1,1,1, -99, 32, 437 . "ave(X)", zsto,zout) 438 c 439 ENDIF !iflag_con .GE. 3 440 c34EK 441 c 442 c Champs interpolles sur des niveaux de pression 443 DO k=1, nlevENS 444 CALL histdef(nid_mth, "u"//clev(k), 445 . "Zonal wind"//clev(k)//"mb","m/s", 446 . iim,jjmp1,nhori, 1,1,1, -99, 32, 447 . "ave(X)", zsto,zout) 448 c 449 CALL histdef(nid_mth, "v"//clev(k), 450 . "Meridional wind"//clev(k)//"mb","m/s", 451 . iim,jjmp1,nhori, 1,1,1, -99, 32, 452 . "ave(X)", zsto,zout) 453 c 454 CALL histdef(nid_mth, "w"//clev(k), 455 . "Vertical wind"//clev(k)//"mb","m/s", 456 . iim,jjmp1,nhori, 1,1,1, -99, 32, 457 . "ave(X)", zsto,zout) 458 c 459 CALL histdef(nid_mth, "phi"//clev(k), 460 . "Geopotential"//clev(k)//"mb","m2/s2", 461 . iim,jjmp1,nhori, 1,1,1, -99, 32, 462 . "ave(X)", zsto,zout) 463 c 464 ENDDO 465 ENDIF !lev_histmth.GE.1 466 c 467 IF(lev_histmth.GE.2) THEN 468 c 469 c Champs 3D: 470 c 471 c 472 CALL histdef(nid_mth,"lwcon","Cloud water content","Kg/Kg", 473 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 474 . "ave(X)", zsto,zout) 475 c 476 CALL histdef(nid_mth,"iwcon","Cloud ice water content","Kg/Kg", 477 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 478 . "ave(X)", zsto,zout) 479 c 480 CALL histdef(nid_mth, "temp", "Air temperature", "K", 481 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 482 . "ave(X)", zsto,zout) 483 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", 489 c . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 490 c . "ave(X)", zsto,zout) 491 c 492 CALL histdef(nid_mth, "geop", "Geopotential height", "m", 493 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 494 . "ave(X)", zsto,zout) 495 c 496 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", 497 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 498 . "ave(X)", zsto,zout) 499 c 500 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", 501 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 502 . "ave(X)", zsto,zout) 503 c 504 CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s", 505 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 506 . "ave(X)", zsto,zout) 507 c 508 CALL histdef(nid_mth, "pres", "Air pressure", "Pa", 509 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 510 . "ave(X)", zsto,zout) 511 c 512 CALL histdef(nid_mth, "rneb", "Cloud fraction", "-", 513 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 514 . "ave(X)", zsto,zout) 515 c 516 CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction" 517 . , "-", 518 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 519 . "ave(X)", zsto,zout) 520 c 521 CALL histdef(nid_mth, "rhum", "Relative humidity", "-", 522 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 523 . "ave(X)", zsto,zout) 524 c 525 CALL histdef(nid_mth, "ozone", "Ozone concentration", "-", 526 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 527 . "ave(X)", zsto,zout) 528 c 529 CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s", 530 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 531 . "ave(X)", zsto,zout) 532 c 533 CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s", 534 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 535 . "ave(X)", zsto,zout) 536 c 537 CALL histdef(nid_mth, "dqphy", "Physics dQ", "Kg/Kg/s", 538 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 539 . "ave(X)", zsto,zout) 540 c 541 c#define histmthNMC 542 c#ifdef histmthNMC 543 c#include "ini_histmthNMC.h" 544 c#endif 545 c 546 ENDIF !lev_histmth.GE.2 547 c 548 IF(lev_histmth.GE.3) THEN 549 c 550 DO nsrf=1, nbsrf 217 551 c 218 552 call histdef(nid_mth, "albe_"//clnsurf(nsrf), … … 229 563 . iim,jjmp1,nhori, 1,1,1, -99, 32, 230 564 . "ave(X)", zsto,zout) 231 232 END DO 233 C 234 c CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-", 235 c . iim,jjmp1,nhori, 1,1,1, -99, 32, 236 c . "ave(X)", zsto,zout) 565 c 566 ENDDO !nsrf=1, nbsrf 237 567 c 238 568 CALL histdef(nid_mth, "albs", "Surface albedo", "-", … … 243 573 . "ave(X)", zsto,zout) 244 574 c 245 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 246 . iim,jjmp1,nhori, 1,1,1, -99, 32, 247 . "ave(X)", zsto,zout) 248 c 249 CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 250 . iim,jjmp1,nhori, 1,1,1, -99, 32, 251 . "ave(X)", zsto,zout) 252 c 253 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-", 254 . iim,jjmp1,nhori, 1,1,1, -99, 32, 255 . "ave(X)", zsto,zout) 256 c 257 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-", 258 . iim,jjmp1,nhori, 1,1,1, -99, 32, 259 . "ave(X)", zsto,zout) 260 c 261 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-", 262 . iim,jjmp1,nhori, 1,1,1, -99, 32, 263 . "ave(X)", zsto,zout) 264 c 265 CALL histdef(nid_mth, "cldt", "Total cloudiness", "-", 266 . iim,jjmp1,nhori, 1,1,1, -99, 32, 267 . "ave(X)", zsto,zout) 268 c 269 CALL histdef(nid_mth,"cldq","Cloud liquid water path","Kg/m2", 270 . iim,jjmp1,nhori, 1,1,1, -99, 32, 271 . "ave(X)", zsto,zout) 272 cIM: 071003 273 CALL histdef(nid_mth,"lwp","Cloud water path","Kg/m2", 274 . iim,jjmp1,nhori, 1,1,1, -99, 32, 275 . "ave(X)", zsto,zout) 276 c 277 CALL histdef(nid_mth,"iwp","Cloud ice water path","Kg/m2", 278 . iim,jjmp1,nhori, 1,1,1, -99, 32, 279 . "ave(X)", zsto,zout) 280 cIM: 071003 281 CALL histdef(nid_mth,"lwcon","Cloud water content","Kg/Kg", 282 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 283 . "ave(X)", zsto,zout) 284 c 285 CALL histdef(nid_mth,"iwcon","Cloud ice water content","Kg/Kg", 286 . iim,jjmp1,nhori, klev,1,klev, nvert, 32, 287 . "ave(X)", zsto,zout) 288 c 289 CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 290 . iim,jjmp1,nhori, 1,1,1, -99, 32, 291 . "ave(X)", zsto,zout) 292 c 293 CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 294 . iim,jjmp1,nhori, 1,1,1, -99, 32, 295 . "ave(X)", zsto,zout) 296 c 297 CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-", 298 . iim,jjmp1,nhori, 1,1,1, -99, 32, 299 . "ave(X)", zsto,zout) 300 c 301 CALL histdef(nid_mth, "vq", "Merid humidity transport", "-", 302 . iim,jjmp1,nhori, 1,1,1, -99, 32, 303 . "ave(X)", zsto,zout) 304 cKE43 305 IF (iflag_con .GE. 3) THEN ! sb 306 c 307 CALL histdef(nid_mth, "cape", "Conv avlbl pot ener", "J/Kg", 308 . iim,jjmp1,nhori, 1,1,1, -99, 32, 309 . "ave(X)", zsto,zout) 310 c 311 CALL histdef(nid_mth, "pbase", "Cld base pressure", "hPa", 312 . iim,jjmp1,nhori, 1,1,1, -99, 32, 313 . "ave(X)", zsto,zout) 314 c 315 CALL histdef(nid_mth, "ptop", "Cld top pressure", "hPa", 316 . iim,jjmp1,nhori, 1,1,1, -99, 32, 317 . "ave(X)", zsto,zout) 318 c 319 CALL histdef(nid_mth, "fbase", "Cld base mass flux", "Kg/m2/s", 320 . iim,jjmp1,nhori, 1,1,1, -99, 32, 321 . "ave(X)", zsto,zout) 322 c 323 c 324 CALL histdef(nid_mth, "upwd", "saturated updraft", "Kg/m2/s", 575 ENDIF !lev_histmth.GE.3 576 c 577 IF(lev_histmth.GE.4) THEN 578 c 579 CALL histdef(nid_mth, "clwcon", 580 . "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", 325 586 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 326 587 . "ave(X)", zsto,zout) … … 331 592 c 332 593 CALL histdef(nid_mth, "dnwd0", "unsat. downdraft", "Kg/m2/s", 333 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,334 . "ave(X)", zsto,zout)335 c336 CALL histdef(nid_mth,"Ma","undilute adiab updraft","Kg/m2/s",337 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,338 . "ave(X)", zsto,zout)339 c340 CALL histdef(nid_mth, "prw", "Precipitable water", "kg/m2",341 . iim,jjmp1,nhori, 1,1,1, -99, 32,342 . "ave(X)", zsto,zout)343 c344 ENDIF345 c34EK346 c347 c Champs 3D:348 c349 CALL histdef(nid_mth, "temp", "Air temperature", "K",350 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,351 . "ave(X)", zsto,zout)352 c353 CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg",354 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,355 . "ave(X)", zsto,zout)356 cIM: 071003357 CALL histdef(nid_mth,"wvap","Water vapor mixing ratio","Kg/Kg",358 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,359 . "ave(X)", zsto,zout)360 c361 CALL histdef(nid_mth, "geop", "Geopotential height", "m",362 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,363 . "ave(X)", zsto,zout)364 c365 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s",366 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,367 . "ave(X)", zsto,zout)368 c369 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s",370 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,371 . "ave(X)", zsto,zout)372 c373 CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s",374 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,375 . "ave(X)", zsto,zout)376 c377 CALL histdef(nid_mth, "pres", "Air pressure", "Pa",378 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,379 . "ave(X)", zsto,zout)380 c381 CALL histdef(nid_mth, "rneb", "Cloud fraction", "-",382 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,383 . "ave(X)", zsto,zout)384 c385 CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction"386 . , "-",387 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,388 . "ave(X)", zsto,zout)389 c390 CALL histdef(nid_mth, "rhum", "Relative humidity", "-",391 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,392 . "ave(X)", zsto,zout)393 c394 CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content"395 . , "kg/kg",396 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,397 . "ave(X)", zsto,zout)398 c399 CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg",400 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,401 . "ave(X)", zsto,zout)402 c403 CALL histdef(nid_mth, "dtphy", "Physics dT", "K/s",404 594 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 405 595 . "ave(X)", zsto,zout) … … 519 709 ENDIF 520 710 C 521 CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",522 . iim,jjmp1,nhori, klev,1,klev,nvert, 32,523 . "ave(X)", zsto,zout)524 c525 711 if (nqmax.GE.3) THEN 526 712 DO iq=1,nqmax-2 … … 537 723 ENDIF 538 724 c 539 CALL histdef(nid_mth, "SWupTOA", "SWup at TOA","W/m2", 540 . iim,jjmp1,nhori, 1,1,1,-99, 541 . 32, "ave(X)", zsto,zout) 542 c 543 CALL histdef(nid_mth, "SWupSFC", "SWup at surface","W/m2", 544 . iim,jjmp1,nhori, 1,1,1,-99, 545 . 32, "ave(X)", zsto,zout) 546 c 547 CALL histdef(nid_mth, "SWdnTOA", "SWdn at TOA","W/m2", 548 . iim,jjmp1,nhori, 1,1,1,-99, 549 . 32, "ave(X)", zsto,zout) 550 c 551 CALL histdef(nid_mth, "SWdnSFC", "SWdn at surface","W/m2", 552 . iim,jjmp1,nhori, 1,1,1,-99, 553 . 32, "ave(X)", zsto,zout) 554 c 555 CALL histdef(nid_mth, "SWupTOAclr", 556 . "SWup clear sky at TOA","W/m2", 557 . iim,jjmp1,nhori, 1,1,1,-99, 558 . 32, "ave(X)", zsto,zout) 559 c 560 CALL histdef(nid_mth, "SWupSFCclr", 561 . "SWup clear sky at surface","W/m2", 562 . iim,jjmp1,nhori, 1,1,1,-99, 563 . 32, "ave(X)", zsto,zout) 564 c 565 CALL histdef(nid_mth, "SWdnTOAclr", 566 . "SWdn clear sky at TOA","W/m2", 567 . iim,jjmp1,nhori, 1,1,1,-99, 568 . 32, "ave(X)", zsto,zout) 569 c 570 CALL histdef(nid_mth, "SWdnSFCclr", 571 . "SWdn clear sky at surface","W/m2", 572 . iim,jjmp1,nhori, 1,1,1,-99, 573 . 32, "ave(X)", zsto,zout) 574 575 c Champs interpolles sur des niveaux de pression 576 577 CALL histdef(nid_mth, "u850", "Zonal wind 850mb", "m/s", 578 . iim,jjmp1,nhori, 1,1,1, -99, 32, 579 . "ave(X)", zsto,zout) 580 581 CALL histdef(nid_mth, "v850", "Meridional wind 850mb", "m/s", 582 . iim,jjmp1,nhori, 1,1,1, -99, 32, 583 . "ave(X)", zsto,zout) 584 c 585 CALL histdef(nid_mth, "u500", "Zonal wind 500mb", "m/s", 586 . iim,jjmp1,nhori, 1,1,1, -99, 32, 587 . "ave(X)", zsto,zout) 588 589 CALL histdef(nid_mth, "v500", "Meridional wind 500mb", "m/s", 590 . iim,jjmp1,nhori, 1,1,1, -99, 32, 591 . "ave(X)", zsto,zout) 592 593 CALL histdef(nid_mth, "u200", "Zonal wind 200mb", "m/s", 594 . iim,jjmp1,nhori, 1,1,1, -99, 32, 595 . "ave(X)", zsto,zout) 596 597 CALL histdef(nid_mth, "v200", "Meridional wind 200mb", "m/s", 598 . iim,jjmp1,nhori, 1,1,1, -99, 32, 599 . "ave(X)", zsto,zout) 600 601 CALL histdef(nid_mth, "phi500","Geopotentiel à 500mb","m2/s2", 602 . iim,jjmp1,nhori, 1,1,1, -99, 32, 603 . "ave(X)", zsto,zout) 604 605 CALL histdef(nid_mth, "w500", "Vertical wind 500mb", "Pa/s", 606 . iim,jjmp1,nhori, 1,1,1, -99, 32, 607 . "ave(X)", zsto,zout) 725 ENDIF !lev_histmth.GE.4 608 726 c 609 727 CALL histend(nid_mth)
Note: See TracChangeset
for help on using the changeset viewer.