Changeset 766 for LMDZ4/trunk/libf/phylmd/ini_histhf.h
- Timestamp:
- Jun 4, 2007, 4:34:47 PM (17 years ago)
- File:
-
- 1 edited
-
LMDZ4/trunk/libf/phylmd/ini_histhf.h (modified) (17 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ini_histhf.h
r684 r766 4 4 5 5 IF (ok_hf) THEN 6 c$OMP MASTER 6 7 c 7 8 zstohf = ecrit_hf … … 14 15 15 16 c 16 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 17 DO i = 1, iim 18 zx_lon(i,1) = rlon(i+1) 19 zx_lon(i,jjmp1) = rlon(i+1) 20 ENDDO 21 22 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 23 CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 24 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 17 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 18 cym DO i = 1, iim 19 cym zx_lon(i,1) = rlon(i+1) 20 cym zx_lon(i,jjmp1) = rlon(i+1) 21 cym ENDDO 22 23 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 24 25 cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat, 26 cym CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 27 cym . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 28 cym . nhori, nid_hf) 29 30 CALL histbeg_phy("histhf", itau_phy, zjulian, dtime, 25 31 . nhori, nid_hf) 26 32 … … 34 40 c 35 41 CALL histdef(nid_hf, "aireTER","Grid area CONT","-", 36 . iim,jj mp1,nhori, 1,1,1, -99, 32,42 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 37 43 . "once", zstohf,zout) 38 44 c 39 45 CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-", 40 . iim,jj mp1,nhori, 1,1,1, -99, 32,46 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 41 47 . "once", zstohf,zout) 42 48 c 43 49 CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-", 44 . iim,jj mp1,nhori, 1,1,1, -99, 32,50 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 45 51 . "once", zstohf,zout) 46 52 c … … 48 54 c 49 55 CALL histdef(nid_hf, "t2m", "Temperature 2m", "K", 50 . iim,jj mp1,nhori, 1,1,1, -99, 32,56 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 51 57 . "inst(X)", zstohf,zout) 52 58 c 53 59 CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg", 54 . iim,jj mp1,nhori, 1,1,1, -99, 32,60 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 55 61 . "inst(X)", zstohf,zout) 56 62 c 57 63 CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa", 58 . iim,jj mp1,nhori, 1,1,1, -99, 32,64 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 59 65 . "inst(X)", zstohf,zout) 60 66 c 61 67 CALL histdef(nid_hf, "precip", "Precipitation", "kg/(s*m2)", 62 . iim,jj mp1,nhori, 1,1,1, -99, 32,68 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 63 69 . "ave(X)", zstophy,zout) 64 70 c 65 71 cIM ENSEMBLES BEG 66 72 CALL histdef(nid_hf, "tsol", "Surface Temperature", "K", 67 . iim,jj mp1,nhori, 1,1,1, -99, 32,73 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 68 74 . "inst(X)", zstohf,zout) 69 75 c 70 76 CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa", 71 . iim,jj mp1,nhori, 1,1,1, -99, 32,77 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 72 78 . "inst(X)", zstohf,zout) 73 79 c 74 80 CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s", 75 . iim,jj mp1,nhori, 1,1,1, -99, 32,81 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 76 82 . "inst(X)", zstohf,zout) 77 83 c 78 84 CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s", 79 . iim,jj mp1,nhori, 1,1,1, -99, 32,85 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 80 86 . "inst(X)", zstohf,zout) 81 87 c 82 88 CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s", 83 . iim,jj mp1,nhori, 1,1,1, -99, 32,89 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 84 90 . "inst(X)", zstohf,zout) 85 91 c … … 99 105 CALL histdef(nid_hf, "phi"//bb2, 100 106 . "Geopotential "//bb2//"mb", "m", 101 . iim,jj mp1,nhori, 1,1,1, -99, 32,107 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 102 108 . "inst(X)", zstohf,zout) 103 109 ENDIF !bb2.EQ."500" … … 109 115 c 110 116 CALL histdef(nid_hf, "cldt", "Total cloudiness", "%", 111 . iim,jj mp1,nhori, 1,1,1, -99, 32,117 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 112 118 . "ave(X)", zstophy,zout) 113 119 c 114 120 CALL histdef(nid_hf, "SWnetOR", 115 121 . "Sfce net SW radiation OR", "W/m2", 116 . iim,jj mp1,nhori, 1,1,1, -99, 32,122 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 117 123 . "ave(X)", zstophy,zout) 118 124 c 119 125 CALL histdef(nid_hf, "SWdownOR", 120 126 . "Sfce incident SW radiation OR", "W/m2", 121 . iim,jj mp1,nhori, 1,1,1, -99, 32,127 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 122 128 . "ave(X)", zstophy,zout) 123 129 c 124 130 CALL histdef(nid_hf, "LWdownOR", 125 131 . "Sfce incident LW radiation OR", "W/m2", 126 . iim,jj mp1,nhori, 1,1,1, -99, 32,132 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 127 133 . "ave(X)", zstophy,zout) 128 134 ENDIF !lev_histhf.GE.2 … … 148 154 CALL histdef(nid_hf, "t"//bb2, 149 155 . "Temperature "//bb2//"mb","K", 150 . iim,jj mp1,nhori, 1,1,1, -99, 32,156 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 151 157 . "inst(X)", zstohf,zout) 152 158 c … … 154 160 CALL histdef(nid_hf, "phi"//bb2, 155 161 . "Geopotential "//bb2//"mb", "m", 156 . iim,jj mp1,nhori, 1,1,1, -99, 32,162 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 157 163 . "inst(X)", zstohf,zout) 158 164 ENDIF … … 160 166 CALL histdef(nid_hf, "q"//bb2, 161 167 . "Specific humidity "//bb2//"mb","kg/kg", 162 . iim,jj mp1,nhori, 1,1,1, -99, 32,168 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 163 169 . "inst(X)", zstohf,zout) 164 170 c 165 171 CALL histdef(nid_hf, "u"//bb2, 166 172 . "Zonal wind "//bb2//"mb","m/s", 167 . iim,jj mp1,nhori, 1,1,1, -99, 32,173 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 168 174 . "inst(X)", zstohf,zout) 169 175 c 170 176 CALL histdef(nid_hf, "v"//bb2, 171 177 . "Meridional wind "//bb2//"mb","m/s", 172 . iim,jj mp1,nhori, 1,1,1, -99, 32,178 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 173 179 . "inst(X)", zstohf,zout) 174 180 c … … 210 216 CALL histdef(nid_hf, "u"//bb2, 211 217 . "Zonal wind "//bb2//"mb","m/s", 212 . iim,jj mp1,nhori, 1,1,1, -99, 32,218 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 213 219 . "inst(X)", zstohf,zstohf) 214 220 c 215 221 CALL histdef(nid_hf, "v"//bb2, 216 222 . "Meridional wind "//bb2//"mb","m/s", 217 . iim,jj mp1,nhori, 1,1,1, -99, 32,223 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 218 224 . "inst(X)", zstohf,zstohf) 219 225 c … … 223 229 CALL histdef(nid_hf, "topl", 224 230 . "Outgoing longwave radiation", "W/m2", 225 . iim,jj mp1,nhori, 1,1,1, -99, 32,231 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 226 232 . "ave(X)", zstophy,zout) 227 233 c 228 234 CALL histdef(nid_hf, "precip","Total precipitation rate", 229 . "kg/(s*m2)",iim,jj mp1,nhori, 1,1,1, -99, 32,235 . "kg/(s*m2)",iim,jjphy_nb,nhori, 1,1,1, -99, 32, 230 236 . "ave(X)", zstophy,zout) 231 237 c 232 238 CALL histdef(nid_hf, "slp", 233 239 . "Mean sea-level pressure", "Pa", 234 . iim,jj mp1,nhori, 1,1,1, -99, 32,240 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 235 241 . "inst(X)", zstohf,zout) 236 242 c … … 255 261 CALL histdef(nid_hf, "t"//bb2, 256 262 . "Temperature "//bb2//"mb","K", 257 . iim,jj mp1,nhori, 1,1,1, -99, 32,263 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 258 264 . "inst(X)", zstohf,zstohf) 259 265 c 260 266 CALL histdef(nid_hf, "q"//bb2, 261 267 . "Specific humidity "//bb2//"mb","kg/kg", 262 . iim,jj mp1,nhori, 1,1,1, -99, 32,268 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 263 269 . "inst(X)", zstohf,zstohf) 264 270 c … … 269 275 CALL histdef(nid_hf, "phi"//bb2, 270 276 . "Geopotential "//bb2//"mb", "m", 271 . iim,jj mp1,nhori, 1,1,1, -99, 32,277 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 272 278 . "inst(X)", zstohf,zstohf) 273 279 c 274 280 CALL histdef(nid_hf, "w"//bb2, 275 281 . "Vertical motion "//bb2//"mb","Pa/s", 276 . iim,jj mp1,nhori, 1,1,1, -99, 32,282 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 277 283 . "inst(X)", zstohf,zstohf) 278 284 ENDIF !bb2.EQ."500" … … 282 288 CALL histdef(nid_hf, "t"//bb3, 283 289 . "Temperature "//bb3//"mb","K", 284 . iim,jj mp1,nhori, 1,1,1, -99, 32,290 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 285 291 . "inst(X)", zstohf,zstohf) 286 292 c 287 293 CALL histdef(nid_hf, "u"//bb3, 288 294 . "Zonal wind "//bb3//"mb","m/s", 289 . iim,jj mp1,nhori, 1,1,1, -99, 32,295 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 290 296 . "inst(X)", zstohf,zstohf) 291 297 c 292 298 CALL histdef(nid_hf, "v"//bb3, 293 299 . "Meridional wind "//bb3//"mb","m/s", 294 . iim,jj mp1,nhori, 1,1,1, -99, 32,300 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 295 301 . "inst(X)", zstohf,zstohf) 296 302 c … … 300 306 c 301 307 CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s", 302 . iim,jj mp1,nhori, 1,1,1, -99, 32,308 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 303 309 . "inst(X)", zstohf,zout) 304 310 c 305 311 CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s", 306 . iim,jj mp1,nhori, 1,1,1, -99, 32,312 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 307 313 . "inst(X)", zstohf,zout) 308 314 c 309 315 CALL histdef(nid_hf, "t2m", "Temperature 2m", "K", 310 . iim,jj mp1,nhori, 1,1,1, -99, 32,316 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 311 317 . "inst(X)", zstohf,zout) 312 318 c 313 319 CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg", 314 . iim,jj mp1,nhori, 1,1,1, -99, 32,320 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 315 321 . "inst(X)", zstohf,zout) 316 322 c 317 323 CALL histdef(nid_hf, "psol", "Surface Pressure", "Pa", 318 . iim,jj mp1,nhori, 1,1,1, -99, 32,324 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 319 325 . "inst(X)", zstohf,zout) 320 326 c 321 327 CALL histdef(nid_hf, "cldt", "Total cloudiness", "%", 322 . iim,jj mp1,nhori, 1,1,1, -99, 32,328 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 323 329 . "ave(X)", zstohf,zout) 324 330 c 325 331 CALL histdef(nid_hf, "taux", 326 332 $ "Zonal wind stress","Pa", 327 $ iim,jj mp1,nhori, 1,1,1, -99, 32,333 $ iim,jjphy_nb,nhori, 1,1,1, -99, 32, 328 334 $ "ave(X)", zstohf,zout) 329 335 c 330 336 CALL histdef(nid_hf, "tauy", 331 337 $ "Meridional wind stress ", "Pa", 332 $ iim,jj mp1,nhori, 1,1,1, -99, 32,338 $ iim,jjphy_nb,nhori, 1,1,1, -99, 32, 333 339 $ "ave(X)", zstohf,zout) 334 340 c 335 341 CALL histdef(nid_hf, "prw", "Precipitable water", "kg/m2", 336 . iim,jj mp1,nhori, 1,1,1, -99, 32,342 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 337 343 cIM 090904 . "ave(X)", zstohf,zout) 338 344 . "ave(X)", zstophy,zout) 339 345 c 340 346 CALL histdef(nid_hf, "sens", "Sensible heat flux", "W/m2", 341 . iim,jj mp1,nhori, 1,1,1, -99, 32,347 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 342 348 cIM 090904 . "ave(X)", zstohf,zout) 343 349 . "ave(X)", zstophy,zout) 344 350 c 345 351 CALL histdef(nid_hf, "flat", "Latent heat flux", "W/m2", 346 . iim,jj mp1,nhori, 1,1,1, -99, 32,352 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 347 353 . "ave(X)", zstohf,zout) 348 354 c 349 355 c CALL histdef(nid_hf, "snow_mass", "Snow Mass", "kg/m2", 350 c . iim,jj mp1,nhori, 1,1,1, -99, 32,356 c . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 351 357 c . "ave(X)", zstohf,zout) 352 358 c 353 359 CALL histdef(nid_hf, "SWdnSFC", "SWdn at surface","W/m2", 354 . iim,jj mp1,nhori, 1,1,1,-99,360 . iim,jjphy_nb,nhori, 1,1,1,-99, 355 361 . 32, "ave(X)", zstohf,zout) 356 362 c 357 363 CALL histdef(nid_hf, "SWupSFC", "SWup at surface","W/m2", 358 . iim,jj mp1,nhori, 1,1,1,-99,364 . iim,jjphy_nb,nhori, 1,1,1,-99, 359 365 . 32, "ave(X)", zstohf,zout) 360 366 c 361 367 CALL histdef(nid_hf,"LWdnSFC","Down. IR rad. at surface", 362 . "W/m2", iim,jj mp1,nhori, 1,1,1, -99, 32,368 . "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 363 369 . "ave(X)", zstohf,zout) 364 370 c 365 371 CALL histdef(nid_hf,"LWupSFC","Upwd. IR rad. at surface", 366 . "W/m2", iim,jj mp1,nhori, 1,1,1, -99, 32,372 . "W/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32, 367 373 . "ave(X)", zstohf,zout) 368 374 c 369 375 CALL histdef(nid_hf, "SWdnTOA", "SWdn at TOA","W/m2", 370 . iim,jj mp1,nhori, 1,1,1,-99,376 . iim,jjphy_nb,nhori, 1,1,1,-99, 371 377 . 32, "ave(X)", zstohf,zout) 372 378 c 373 379 CALL histdef(nid_hf, "SWupTOA", "SWup at TOA","W/m2", 374 . iim,jj mp1,nhori, 1,1,1,-99,380 . iim,jjphy_nb,nhori, 1,1,1,-99, 375 381 . 32, "ave(X)", zstohf,zout) 376 382 c … … 378 384 CALL histdef(nid_hf, "PV"//ctetaSTD(k), 379 385 $ "Vorticite potentielle "//ctetaSTD(k)//"K", 380 $ "1/(Pa*s)",iim,jj mp1,nhori, 1,1,1, -99,386 $ "1/(Pa*s)",iim,jjphy_nb,nhori, 1,1,1, -99, 381 387 $ 32, "inst(X)", zstohf,zout) 382 388 ENDDO !k=1, nbteta … … 393 399 CALL histend(nid_hf) 394 400 c 401 c$OMP END MASTER 395 402 ENDIF ! ok_hf
Note: See TracChangeset
for help on using the changeset viewer.
