Changeset 897
- Timestamp:
- Feb 6, 2008, 6:48:54 PM (17 years ago)
- Location:
- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd
- Files:
-
- 6 edited
-
ini_histday.h (modified) (8 diffs)
-
ini_histhf.h (modified) (5 diffs)
-
ini_histmth.h (modified) (2 diffs)
-
write_histday.h (modified) (5 diffs)
-
write_histhf.h (modified) (3 diffs)
-
write_histmth.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histday.h
r892 r897 192 192 . "ave(X)", zstophy,zout) 193 193 c 194 CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa", 195 . iim,jj_nb,nhori, 1,1,1, -99, 32, 196 . "ave(X)", zstophy,zout) 197 c 198 CALL histdef(nid_day, "cape_max", "CAPE max.", 199 . "J/kg", 200 . iim,jj_nb,nhori, 1,1,1, -99, 32, 201 . capemaxcels, zstophy,zout) 202 c 203 CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", 204 . "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 205 . "ave(X)", zstophy,zout) 206 c 207 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", 208 . iim,jj_nb,nhori, 1,1,1, -99, 32, 209 . "ave(X)", zstophy,zout) 210 c 211 CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2", 212 . iim,jj_nb,nhori, 1,1,1,-99, 213 . 32, "ave(X)", zstophy,zout) 214 c 194 215 c Champs dynamiques sur niveaux de pression 195 216 c … … 236 257 ENDDO !nlevSTD 237 258 c 238 CALL histdef(nid_day, "slp", "Sea Level Pressure", "Pa", 239 . iim,jj_nb,nhori, 1,1,1, -99, 32, 240 . "ave(X)", zstophy,zout) 241 c 242 CALL histdef(nid_day, "cape_max", "CAPE max.", 243 . "J/kg", 244 . iim,jj_nb,nhori, 1,1,1, -99, 32, 245 . capemaxcels, zstophy,zout) 246 c 247 CALL histdef(nid_day, "solldown", "Down. IR rad. at surface", 248 . "W/m2", iim,jj_nb,nhori, 1,1,1, -99, 32, 249 . "ave(X)", zstophy,zout) 250 c 251 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", 252 . iim,jj_nb,nhori, 1,1,1, -99, 32, 253 . "ave(X)", zstophy,zout) 254 c 255 CALL histdef(nid_day, "SWdnSFC", "SWdn at surface","W/m2", 256 . iim,jj_nb,nhori, 1,1,1,-99, 257 . 32, "ave(X)", zstophy,zout) 259 258 260 c 259 261 cIM diagnostiques flux ocean-atm ou ocean-glace de mer … … 310 312 . "ave(X)", zstophy,zout) 311 313 c 312 ENDIF !lev_histday.GE.2 313 c================================================================= 314 IF(lev_histday.GE.3) THEN 315 c================================================================= 316 c INITIALISATION DES CHAMPS 3D 317 c================================================================= 318 c================================================================= 319 c================================================================= 320 c Champs 3D: 321 c 322 CALL histdef(nid_day, "temp", "Air temperature", "K", 323 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 324 . "ave(X)", zstophy,zout) 325 c 326 CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg", 327 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 328 . "ave(X)", zstophy,zout) 329 c 330 CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2", 331 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 332 . "ave(X)", zstophy,zout) 333 c 334 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 335 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 336 . "ave(X)", zstophy,zout) 337 c 338 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 339 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 340 . "ave(X)", zstophy,zout) 341 c 342 CALL histdef(nid_day, "vitw", "Vertical wind", "Pa/s", 343 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 344 . "ave(X)", zstophy,zout) 345 c 346 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 347 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 348 . "ave(X)", zstophy,zout) 349 c 350 c================================================================= 351 c FIN INITIALISATION DES CHAMPS 3D 352 c================================================================= 353 ENDIF ! lev_histday.GE.3 354 c 355 IF (lev_histday.GE.4) THEN 356 c================================================================= 357 c 358 c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES 359 c 360 c================================================================= 314 CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%", 315 . iim,jj_nb,nhori, 1,1,1, -99, 32, 316 . "ave(X)", zstophy,zout) 361 317 c 362 318 CALL histdef(nid_day, "SWupTOA", "SWup at TOA","W/m2", … … 376 332 . iim,jj_nb,nhori, 1,1,1,-99, 377 333 . 32, "ave(X)", zstophy,zout) 378 334 c 379 335 CALL histdef(nid_day, "SWupSFCclr", 380 336 . "SWup clear sky at surface","W/m2", 381 337 . iim,jj_nb,nhori, 1,1,1,-99, 382 338 . 32, "ave(X)", zstophy,zout) 383 339 c 384 340 CALL histdef(nid_day, "SWdnTOAclr", 385 341 . "SWdn clear sky at TOA","W/m2", 386 342 . iim,jj_nb,nhori, 1,1,1,-99, 387 343 . 32, "ave(X)", zstophy,zout) 388 344 c 389 345 CALL histdef(nid_day, "SWdnSFCclr", 390 346 . "SWdn clear sky at surface","W/m2", … … 399 355 c 400 356 CALL histdef(nid_day, "LWupSFC", 401 . "LW down at surface","W/m2", 402 . iim,jj_nb,nhori, 1,1,1,-99, 403 . 32, "ave(X)", zstophy,zout) 357 . "LW up at surface","W/m2", 358 . iim,jj_nb,nhori, 1,1,1,-99, 359 . 32, "ave(X)", zstophy,zout) 360 c 361 CALL histdef(nid_day,"lwp","Cloud water path","kg/m2", 362 . iim,jj_nb,nhori, 1,1,1, -99, 32, 363 . "ave(X)", zstophy,zout) 364 c 365 CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2", 366 . iim,jj_nb,nhori, 1,1,1, -99, 32, 367 . "ave(X)", zstophy,zout) 368 c 369 CALL histdef(nid_day, "taux","Zonal wind stress","Pa", 370 . iim,jj_nb,nhori, 1,1,1, -99, 32, 371 . "ave(X)", zstophy,zout) 372 c 373 CALL histdef(nid_day, "tauy","Meridional wind stress ", "Pa", 374 . iim,jj_nb,nhori, 1,1,1, -99, 32, 375 . "ave(X)", zstophy,zout) 404 376 c 405 377 cIM IF (iflag_con.GE.3) THEN … … 413 385 cIM ENDIF !(iflag_con.GE.3) THEN 414 386 c 415 c================================================================= 416 c================================================================= 387 ENDIF !lev_histday.GE.2 388 c================================================================= 389 IF(lev_histday.GE.3) THEN 390 c================================================================= 391 c INITIALISATION DES CHAMPS 3D 392 c================================================================= 393 c================================================================= 394 c================================================================= 395 c Champs 3D: 396 c 397 CALL histdef(nid_day, "temp", "Air temperature", "K", 398 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 399 . "ave(X)", zstophy,zout) 400 c 401 CALL histdef(nid_day, "ovap", "Specific humidity", "kg/kg", 402 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 403 . "ave(X)", zstophy,zout) 404 c 405 CALL histdef(nid_day, "geop", "Geopotential height", "m2/s2", 406 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 407 . "ave(X)", zstophy,zout) 408 c 409 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 410 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 411 . "ave(X)", zstophy,zout) 412 c 413 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 414 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 415 . "ave(X)", zstophy,zout) 416 c 417 CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", 418 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 419 . "ave(X)", zstophy,zout) 420 c 421 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 422 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 423 . "ave(X)", zstophy,zout) 424 c 425 c================================================================= 426 c FIN INITIALISATION DES CHAMPS 3D 427 c================================================================= 428 ENDIF ! lev_histday.GE.3 429 c 430 IF (lev_histday.GE.4) THEN 417 431 c================================================================= 418 432 c INITIALISATION DES CHAMPS SUR LES SOUS SURFACES … … 590 604 . "ave(X)", zstophy,zout) 591 605 c 592 CALL histdef(nid_day,"lwp","Cloud water path","kg/m2", 593 . iim,jj_nb,nhori, 1,1,1, -99, 32, 594 . "ave(X)", zstophy,zout) 595 c 596 CALL histdef(nid_day,"iwp","Cloud ice water path","kg/m2", 597 . iim,jj_nb,nhori, 1,1,1, -99, 32, 598 . "ave(X)", zstophy,zout) 599 c 606 IF (ok_isccp) THEN 600 607 CALL histdef(nid_day,"meantaucld", 601 608 . "ISCCP mean cloud optical thickness","1", 602 609 . iim,jj_nb,nhori, 1,1,1, -99, 32, 603 610 . "ave(X)", zstophy,zout) 611 ENDIF 604 612 c 605 613 CALL histdef(nid_day,"cldtau","Cloud optical thickness","1", … … 621 629 . "ave(X)", zstophy,zout) 622 630 c 623 CALL histdef(nid_day, "rh2m", "Relative humidity at 2m", "%",624 . iim,jj_nb,nhori, 1,1,1, -99, 32,625 . "ave(X)", zstophy,zout)626 c627 631 CALL histdef(nid_day, "qsat2m", "Saturant humidity at 2m", "%", 628 632 . iim,jj_nb,nhori, 1,1,1, -99, 32, -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histhf.h
r776 r897 92 92 DO k=1, nlevSTD 93 93 c 94 c bb=clevSTD(k) 95 c 96 IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) 97 IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) 94 bb2=clevSTD(k) 98 95 c 99 96 c IF(k.GE.2) THEN … … 102 99 c ENDIF 103 100 c 104 IF(bb2.EQ."500") THEN 105 CALL histdef(nid_hf, "phi"//bb2, 101 bb2=clevSTD(k) 102 IF(bb2.EQ."500") THEN 103 CALL histdef(nid_hf, "phi"//bb2, 106 104 . "Geopotential "//bb2//"mb", "m", 107 105 . iim,jj_nb,nhori, 1,1,1, -99, 32, 108 106 . "inst(X)", zstohf,zout) 109 107 ENDIF !bb2.EQ."500" 110 108 ENDDO !k=1, nlevSTD 111 109 c … … 139 137 c 140 138 DO k=1, nlevSTD 139 c 140 c bb=clevSTD(k) 141 c 142 c IF(k.GE.2) THEN 143 c IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k) 144 c IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k) 145 c aa=clevSTD(k) 146 c bb=aa(1:lnblnk1(aa)) 147 c ENDIF 148 c 149 bb2=clevSTD(k) 150 IF(bb2.EQ."1000".OR.bb2.EQ."925".OR. 151 $ bb2.EQ."850".OR.bb2.EQ."700".OR. 152 $ bb2.EQ."500".OR.bb2.EQ."300".OR. 153 $ bb2.EQ."250".OR.bb2.EQ."200".OR. 154 $ bb2.EQ."150".OR.bb2.EQ."100".OR. 155 $ bb2.EQ."50".OR.bb2.EQ."30".OR. 156 $ bb2.EQ."10")THEN 157 c 158 CALL histdef(nid_hf, "t"//bb2, 159 . "Temperature "//bb2//"mb","K", 160 . iim,jj_nb,nhori, 1,1,1, -99, 32, 161 . "inst(X)", zstohf,zout) 162 c 163 IF(bb2.NE."500") THEN 164 CALL histdef(nid_hf, "phi"//bb2, 165 . "Geopotential "//bb2//"mb", "m", 166 . iim,jj_nb,nhori, 1,1,1, -99, 32, 167 . "inst(X)", zstohf,zout) 168 ENDIF 169 c 170 IF(bb2.GE."200") THEN 171 CALL histdef(nid_hf, "q"//bb2, 172 . "Specific humidity "//bb2//"mb","kg/kg", 173 . iim,jj_nb,nhori, 1,1,1, -99, 32, 174 . "inst(X)", zstohf,zout) 175 ENDIF 176 c 177 CALL histdef(nid_hf, "u"//bb2, 178 . "Zonal wind "//bb2//"mb","m/s", 179 . iim,jj_nb,nhori, 1,1,1, -99, 32, 180 . "inst(X)", zstohf,zout) 181 c 182 CALL histdef(nid_hf, "v"//bb2, 183 . "Meridional wind "//bb2//"mb","m/s", 184 . iim,jj_nb,nhori, 1,1,1, -99, 32, 185 . "inst(X)", zstohf,zout) 186 c 187 CALL histdef(nid_hf, "w"//bb2, 188 . "Vertical motion "//bb2//"mb","Pa/s", 189 . iim,jj_nb,nhori, 1,1,1, -99, 32, 190 . "inst(X)", zstohf,zstohf) 191 ENDIF 192 c 193 ENDDO !nlevSTD 194 c 195 ENDIF !lev_histhf.GE.3 196 c 197 IF(lev_histhf.GE.4) THEN 198 c 199 cIM sorties hf 3D 200 c 201 #undef histhf3d 202 #define histhf3d 203 #ifdef histhf3d 204 #include "ini_histhf3d.h" 205 #endif 206 ENDIF !lev_histhf.GE.4 207 c 208 ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN 209 c 210 IF(lev_histhf.GE.1) THEN 211 c 212 DO k=1, nlevSTD 213 c 141 214 c 142 215 c bb=clevSTD(k) … … 149 222 c ENDIF 150 223 c 151 IF(bb2.EQ."850".or.bb2.EQ."700".OR. 152 $ bb2.EQ."500".or.bb2.EQ."200") THEN 153 c 154 CALL histdef(nid_hf, "t"//bb2, 155 . "Temperature "//bb2//"mb","K", 156 . iim,jj_nb,nhori, 1,1,1, -99, 32, 157 . "inst(X)", zstohf,zout) 158 c 159 IF(bb2.NE."500") THEN 160 CALL histdef(nid_hf, "phi"//bb2, 161 . "Geopotential "//bb2//"mb", "m", 162 . iim,jj_nb,nhori, 1,1,1, -99, 32, 163 . "inst(X)", zstohf,zout) 164 ENDIF 165 c 166 CALL histdef(nid_hf, "q"//bb2, 167 . "Specific humidity "//bb2//"mb","kg/kg", 168 . iim,jj_nb,nhori, 1,1,1, -99, 32, 169 . "inst(X)", zstohf,zout) 224 IF(bb2.EQ."850".or.bb2.EQ."200") THEN 170 225 c 171 226 CALL histdef(nid_hf, "u"//bb2, 172 227 . "Zonal wind "//bb2//"mb","m/s", 173 . 174 . "inst(X)", zstohf,zout)228 . iim,jj_nb,nhori, 1,1,1, -99, 32, 229 . "inst(X)", zstohf,zstohf) 175 230 c 176 231 CALL histdef(nid_hf, "v"//bb2, 177 232 . "Meridional wind "//bb2//"mb","m/s", 178 . 179 . "inst(X)", zstohf,zout)180 c 181 ENDIF ! IF(bb2.EQ."850".or.bb2.EQ."700".OR.182 c 183 ENDDO !nlevSTD 184 c 185 ENDIF !lev_histhf.GE.3186 c 187 IF(lev_histhf.GE.4) THEN188 c 189 cIM sorties hf 3D 190 c 191 #undef histhf3d 192 #define histhf3d 193 #ifdef histhf3d 194 #include "ini_histhf3d.h" 195 #endif 196 ENDIF !lev_histhf.GE.4197 c 198 ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN 199 c 200 IF(lev_histhf.GE.1) THEN 201 c 202 DO k=1, nlevSTD 203 c 233 . iim,jj_nb,nhori, 1,1,1, -99, 32, 234 . "inst(X)", zstohf,zstohf) 235 c 236 ENDIF !bb2.EQ."850".or.bb2.EQ."200" 237 ENDDO !k=1, nlevSTD 238 c 239 CALL histdef(nid_hf, "topl", 240 . "Outgoing longwave radiation", "W/m2", 241 . iim,jj_nb,nhori, 1,1,1, -99, 32, 242 . "ave(X)", zstophy,zout) 243 c 244 CALL histdef(nid_hf, "precip","Total precipitation rate", 245 . "kg/(s*m2)",iim,jj_nb,nhori, 1,1,1, -99, 32, 246 . "ave(X)", zstophy,zout) 247 c 248 CALL histdef(nid_hf, "slp", 249 . "Mean sea-level pressure", "Pa", 250 . iim,jj_nb,nhori, 1,1,1, -99, 32, 251 . "inst(X)", zstohf,zout) 252 c 253 cIM 280504 BEG 254 ENDIF !(lev_histhf.GE.1) THEN 255 c 256 IF(lev_histhf.GE.2) THEN 257 c 258 DO k=1, nlevSTD 204 259 c 205 260 c bb=clevSTD(k) … … 212 267 c ENDIF 213 268 c 214 IF(bb2.EQ."850".or.bb2.EQ."200") THEN215 c216 CALL histdef(nid_hf, "u"//bb2,217 . "Zonal wind "//bb2//"mb","m/s",218 . iim,jj_nb,nhori, 1,1,1, -99, 32,219 . "inst(X)", zstohf,zstohf)220 c221 CALL histdef(nid_hf, "v"//bb2,222 . "Meridional wind "//bb2//"mb","m/s",223 . iim,jj_nb,nhori, 1,1,1, -99, 32,224 . "inst(X)", zstohf,zstohf)225 c226 ENDIF !bb2.EQ."850".or.bb2.EQ."200"227 ENDDO !k=1, nlevSTD228 c229 CALL histdef(nid_hf, "topl",230 . "Outgoing longwave radiation", "W/m2",231 . iim,jj_nb,nhori, 1,1,1, -99, 32,232 . "ave(X)", zstophy,zout)233 c234 CALL histdef(nid_hf, "precip","Total precipitation rate",235 . "kg/(s*m2)",iim,jj_nb,nhori, 1,1,1, -99, 32,236 . "ave(X)", zstophy,zout)237 c238 CALL histdef(nid_hf, "slp",239 . "Mean sea-level pressure", "Pa",240 . iim,jj_nb,nhori, 1,1,1, -99, 32,241 . "inst(X)", zstohf,zout)242 c243 cIM 280504 BEG244 ENDIF !(lev_histhf.GE.1) THEN245 c246 IF(lev_histhf.GE.2) THEN247 c248 DO k=1, nlevSTD249 c250 c bb=clevSTD(k)251 c252 c IF(k.GE.2) THEN253 IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)254 IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)255 c aa=clevSTD(k)256 c bb=aa(1:lnblnk1(aa))257 c ENDIF258 c259 269 IF(bb2.EQ."850".OR.bb2.EQ."500") THEN 260 270 c -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h
r896 r897 292 292 . "ave(X)", zstophy,zout) 293 293 c 294 c295 294 c Effets des aerosols 296 295 c … … 1286 1285 . "ave(X)", zstophy,zout) 1287 1286 c 1288 c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", 1287 c Effets des aerosols 1288 c 1289 IF (ok_ade.OR.ok_aie) THEN 1290 CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2", 1291 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1292 . "ave(X)", zsto,zout) 1293 c 1294 CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2", 1295 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1296 . "ave(X)", zsto,zout) 1297 c 1298 CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2", 1299 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1300 . "ave(X)", zsto,zout) 1301 c 1302 CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2", 1303 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1304 . "ave(X)", zsto,zout) 1305 ENDIF 1306 c 1307 c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", 1289 1308 c . iim,jj_nb,nhori, 1,1,1, -99, 32, 1290 1309 c . "ave(X)", zstophy,zout) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histday.h
r892 r897 284 284 CALL histwrite_phy(nid_day,"fder",itau_w,fder) 285 285 c 286 ENDIF !lev_histday.GE.2 287 c 288 IF(lev_histday.GE.3) THEN 289 c================================================================= 290 c ECRITURE DES CHAMPS 3D 291 c================================================================= 292 293 c 294 c Champs 3D: 295 c 296 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 297 CALL histwrite_phy(nid_day,"temp",itau_w,t_seri) 298 c 299 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 300 CALL histwrite_phy(nid_day,"ovap",itau_w,qx(:,:,ivap)) 301 c 302 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 303 CALL histwrite_phy(nid_day,"geop",itau_w,zphi) 304 c 305 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 306 CALL histwrite_phy(nid_day,"vitu",itau_w,u_seri) 307 c 308 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 309 CALL histwrite_phy(nid_day,"vitv",itau_w,v_seri) 310 c 311 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 312 CALL histwrite_phy(nid_day,"vitw",itau_w,omega) 313 c 314 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 315 CALL histwrite_phy(nid_day,"pres",itau_w,pplay) 316 c 317 ENDIF !lev_histday.GE.3 318 c================================================================= 319 c FIN ECRITURE DES CHAMPS 3D 320 c================================================================= 321 IF(lev_histday.ge.4) THEN 322 c================================================================= 323 c 324 c ECRITURE DES CHAMPS SUR LES SOUS SURFACES 325 c 326 c================================================================= 286 DO i=1, klon 287 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 288 ENDDO 289 c 290 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1, zx_tmp_2d,zx_tmp_2d) 291 CALL histwrite_phy(nid_day,"rh2m",itau_w,zx_tmp_fi2d) 327 292 zx_tmp_fi2d(1 : klon) = swup( 1 : klon, klevp1 ) 328 293 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) … … 360 325 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) 361 326 CALL histwrite_phy(nid_day,"LWupSFC",itau_w,zx_tmp_fi2d) 327 c 328 zx_tmp_fi2d(1:klon) = flwp(1:klon) 329 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 330 CALL histwrite_phy(nid_day,"lwp",itau_w,zx_tmp_fi2d) 331 c 332 zx_tmp_fi2d(1:klon) = fiwp(1:klon) 333 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 334 CALL histwrite_phy(nid_day,"iwp",itau_w,zx_tmp_fi2d) 335 c 336 zx_tmp_fi2d(1 : klon) = zxfluxu(1 : klon, 1) 337 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 338 CALL histwrite_phy(nid_day,"taux",itau_w,zx_tmp_fi2d) 339 c 340 zx_tmp_fi2d(1 : klon) = zxfluxv(1 : klon, 1) 341 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 342 CALL histwrite_phy(nid_day,"tauy",itau_w,zx_tmp_fi2d) 362 343 c 363 344 cIM IF (iflag_con.GE.3) THEN … … 369 350 cIM ENDIF 370 351 c 352 ENDIF !lev_histday.GE.2 353 c 354 IF(lev_histday.GE.3) THEN 355 c================================================================= 356 c ECRITURE DES CHAMPS 3D 357 c================================================================= 358 c 359 c Champs 3D: 360 c 361 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 362 CALL histwrite_phy(nid_day,"temp",itau_w,t_seri) 363 c 364 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 365 CALL histwrite_phy(nid_day,"ovap",itau_w,qx(:,:,ivap)) 366 c 367 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 368 CALL histwrite_phy(nid_day,"geop",itau_w,zphi) 369 c 370 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 371 CALL histwrite_phy(nid_day,"vitu",itau_w,u_seri) 372 c 373 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 374 CALL histwrite_phy(nid_day,"vitv",itau_w,v_seri) 375 c 376 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 377 CALL histwrite_phy(nid_day,"vitw",itau_w,omega) 378 c 379 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 380 CALL histwrite_phy(nid_day,"pres",itau_w,pplay) 381 c 382 ENDIF !lev_histday.GE.3 383 c================================================================= 384 c FIN ECRITURE DES CHAMPS 3D 385 c================================================================= 386 IF(lev_histday.ge.4) THEN 387 c================================================================= 388 c 389 c ECRITURE DES CHAMPS SUR LES SOUS SURFACES 390 c 391 c================================================================= 371 392 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter) 372 393 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d) … … 522 543 CALL histwrite_phy(nid_day,"iwcon",itau_w,fiwc) 523 544 c 524 zx_tmp_fi2d(1:klon) = flwp(1:klon)525 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)526 CALL histwrite_phy(nid_day,"lwp",itau_w,zx_tmp_fi2d)527 c528 zx_tmp_fi2d(1:klon) = fiwp(1:klon)529 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)530 CALL histwrite_phy(nid_day,"iwp",itau_w,zx_tmp_fi2d)531 c532 545 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1, meantaucld,zx_tmp_2d) 546 IF (ok_isccp) THEN 533 547 CALL histwrite_phy(nid_day,"meantaucld",itau_w,meantaucld) 548 ENDIF 534 549 c 535 550 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau,zx_tmp_3d) … … 548 563 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 549 564 CALL histwrite_phy(nid_day,"dqcon",itau_w,zx_tmp_fi3d) 550 c551 DO i=1, klon552 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)553 ENDDO554 c555 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1, zx_tmp_2d,zx_tmp_2d)556 CALL histwrite_phy(nid_day,"rh2m",itau_w,zx_tmp_fi2d)557 565 c 558 566 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1, qsat2m,zx_tmp_2d) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histhf.h
r782 r897 117 117 c 118 118 c IF(k.GE.2) THEN 119 IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k)120 IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k)119 c IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k) 120 c IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k) 121 121 c aa=clevSTD(k) 122 122 c bb=aa(1:lnblnk1(aa)) 123 123 c ENDIF 124 124 c 125 IF(bb2.EQ."850".OR.bb2.EQ."700".OR. 126 $ bb2.EQ."500".OR.bb2.EQ."200") THEN 125 bb2=clevSTD(k) 126 IF(bb2.EQ."1000".OR.bb2.EQ."925".OR. 127 $ bb2.EQ."850".OR.bb2.EQ."700".OR. 128 $ bb2.EQ."500".OR.bb2.EQ."300".OR. 129 $ bb2.EQ."250".OR.bb2.EQ."200".OR. 130 $ bb2.EQ."150".OR.bb2.EQ."100".OR. 131 $ bb2.EQ."50".OR.bb2.EQ."30".OR. 132 $ bb2.EQ."10")THEN 127 133 c 128 134 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d) … … 135 141 ENDIF !bb2.NE."500" 136 142 c 143 IF(bb2.GE."200") THEN 144 c 137 145 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d) 138 146 CALL histwrite_phy(nid_hf,"q"//bb2,itau_w, qlevSTD(:,k)) 147 ENDIF 139 148 c 140 149 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d) … … 144 153 CALL histwrite_phy(nid_hf,"v"//bb2,itau_w, vlevSTD(:,k)) 145 154 c 146 ENDIF !bb2.EQ."850".OR.bb2.EQ."700" 155 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, wlevSTD(:,k),zx_tmp_2d) 156 CALL histwrite_phy(nid_hf,"w"//bb2,itau_w, wlevSTD(:,k)) 157 c 158 ENDIF 147 159 ENDDO !nlevENS 148 160 c -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histmth.h
r896 r897 321 321 CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w, 322 322 $ zx_tmp_fi2d) 323 ENDIF !nsrf.G T.2323 ENDIF !nsrf.GE.2 324 324 C 325 325 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) … … 1286 1286 CALL histwrite_phy(nid_mth,"fder",itau_w,fder) 1287 1287 c 1288 c Effets des aerosols 1289 c 1290 IF (ok_ade.OR.ok_aie) THEN 1291 zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon) 1292 CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d) 1293 c 1294 zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon) 1295 CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d) 1296 c 1297 zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon) 1298 CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d) 1299 c 1300 zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon) 1301 CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d) 1302 ENDIF 1303 c 1288 1304 c 1289 1305 c DO i = 1, klon
Note: See TracChangeset
for help on using the changeset viewer.
