Changeset 726 for LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h
- Timestamp:
- Nov 6, 2006, 4:51:16 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h
r704 r726 2 2 ! $Header$ 3 3 ! 4 #ifndef INCA 4 5 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) 5 6 c … … 18 19 . klev, presnivs, nvert) 19 20 20 #ifdef INCA_AER21 cym CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,22 cym . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,23 cym . nhori, nid_tra2)24 CALL histbeg_phy("histrac_aer", itau_phy, zjulian, pdtphys,25 . nhori, nid_tra2)26 27 cym CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,28 cym . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,29 cym . nhori, nid_tra3)30 CALL histbeg_phy("histrac_inst",itau_phy, zjulian, pdtphys,31 . nhori, nid_tra3)32 33 call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",34 . klev, presnivs, nvert)35 call histvert(nid_tra3, "presnivs", "presnivs", "mb",36 . klev, presnivs, nvert)37 #endif38 21 39 22 … … 48 31 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 49 32 . "once", zsto,zout) 50 #ifdef INCA51 CALL histdef(nid_tra, "ps", "Surface pressure", "Pa",52 . iim,jjphy_nb,nhori, 1,1,1,-99, 32,53 . "ave(X)", zsto,zout)54 55 CALL histdef(nid_tra, "ptrop", "Tropopause pressure", "Pa",56 . iim,jjphy_nb,nhori, 1,1,1,-99, 32,57 . "ave(X)", zsto,zout)58 59 C 3d FIELDS60 CALL histdef(nid_tra, "temp", "Air temperature", "K",61 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,62 . "ave(X)", zsto,zout)63 64 CALL histdef(nid_tra, "u", "zonal wind component", "m/s",65 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,66 . "ave(X)", zsto,zout)67 68 CALL histdef(nid_tra, "v", "zonal wind component", "m/s",69 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,70 . "ave(X)", zsto,zout)71 72 CALL histdef(nid_tra, "h2o", "Specific Humidity", "MMR",73 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,74 . "ave(X)", zsto,zout)75 76 CALL histdef(nid_tra, "pmid", "Pressure", "Pa",77 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,78 . "ave(X)", zsto,zout)79 80 CALL histdef(nid_tra, "pdel", "Delta Pressure", "Pa",81 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,82 . "ave(X)", zsto,zout)83 84 ! MS info85 #ifdef INCA_AER86 CALL histdef(nid_tra, "airm", "Air mass", "kg",87 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,88 . "ave(X)", zsto,zout)89 #endif90 #ifdef INCA91 #ifdef INCAINFO92 93 #ifdef INCA_CH494 DO it=1, phtcnt95 WRITE(str2,'(i2.2)') it96 CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",97 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,98 . "ave(X)", zsto,zout)99 ENDDO100 DO it=1, hetcnt101 WRITE(str2,'(i2.2)') it102 CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",103 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,104 . "ave(X)", zsto,zout)105 ENDDO106 DO it=1, extcnt107 WRITE(str2,'(i2.2)') it108 CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",109 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,110 . "ave(X)", zsto,zout)111 ENDDO112 #endif113 #ifdef INCA_NMHC114 DO it=1, phtcnt115 WRITE(str2,'(i2.2)') it116 CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",117 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,118 . "ave(X)", zsto,zout)119 ENDDO120 DO it=1, hetcnt121 WRITE(str2,'(i2.2)') it122 CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",123 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,124 . "ave(X)", zsto,zout)125 ENDDO126 127 DO it=1, extcnt128 WRITE(str2,'(i2.2)') it129 CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",130 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,131 . "ave(X)", zsto,zout)132 ENDDO133 #endif134 DO it=1, nfs135 WRITE(str2,'(i2.2)') it136 CALL histdef(nid_tra, "INV"//str2, "INV"//str2, "CM-3",137 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,138 . "ave(X)", zsto,zout)139 ENDDO140 141 #else142 143 #ifdef INCA_NMHC144 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",145 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,146 . "ave(X)", zsto,zout)147 CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",148 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,149 . "ave(X)", zsto,zout)150 CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",151 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,152 . "ave(X)", zsto,zout)153 CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",154 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,155 . "ave(X)", zsto,zout)156 CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",157 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,158 . "ave(X)", zsto,zout)159 CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",160 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,161 . "ave(X)", zsto,zout)162 #endif163 #ifdef INCA_CH4164 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",165 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,166 . "ave(X)", zsto,zout)167 CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",168 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,169 . "ave(X)", zsto,zout)170 CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",171 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,172 . "ave(X)", zsto,zout)173 CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",174 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,175 . "ave(X)", zsto,zout)176 CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",177 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,178 . "ave(X)", zsto,zout)179 CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",180 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,181 . "ave(X)", zsto,zout)182 #endif183 #endif184 185 ! DO it=1, grpcnt186 ! CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",187 ! . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,188 ! . "ave(X)", zsto,zout)189 ! ENDDO190 #endif191 192 #ifdef INCA_AER193 194 CALL histdef(nid_tra2, "scavcoef_st","scavcoef_st", "S-1",195 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,196 . "ave(X)", zsto,zout)197 CALL histdef(nid_tra2, "scavcoef_cv","scavcoef_cv", "S-1",198 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,199 . "ave(X)", zsto,zout)200 CALL histdef(nid_tra2, "AngstroemComp","AngstroemComp",201 . "angs comp", iim,jjphy_nb,nhori, 1,1,1, -99, 32,202 . "ave(X)", zsto,zout)203 204 CALL histdef(nid_tra2, "TOTAERH2O", "Total aerosol water",205 . "kg/m2", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,206 . "ave(X)", zsto,zout)207 #endif208 #endif209 33 DO it=1,nqmax 210 34 C champ 2D 211 #ifdef INCA212 IF ( prt_flag_ts(it) == 0 ) CYCLE213 CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),214 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,215 . "ave(X)", zsto,zout)216 CALL histdef(nid_tra, "Dep_"//solsym(it), "Dep_"//solsym(it),217 . "cm/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,218 . "ave(X)", zsto,zout)219 #ifdef INCA_AER220 221 CALL histdef(nid_tra2, "Dry_"//solsym(it), "Dry_"//solsym(it),222 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,223 . "ave(X)", zsto,zout)224 225 IF ((it .ge. trmx) .and. (it .le. trnx)) then226 CALL histdef(nid_tra2, "Sed_"//solsym(it), "Sed_"//solsym(it),227 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,228 . "ave(X)", zsto,zout)229 CALL histdef(nid_tra2, "Wet_"//solsym(it), "Wet_"//solsym(it),230 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,231 . "ave(X)", zsto,zout)232 CALL histdef(nid_tra2, "WetST_"//solsym(it), "WetST_"//solsym(it),233 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,234 . "ave(X)", zsto,zout)235 CALL histdef(nid_tra2, "WetCV_"//solsym(it), "WetCV_"//solsym(it),236 . "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,237 . "ave(X)", zsto,zout)238 CALL histdef(nid_tra2, "Emi_alt_"//solsym(it), "Emi_alt_"//solsym(it),239 . "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,240 . "ave(X)", zsto,zout)241 CALL histdef(nid_tra2, "Load_"//solsym(it), "Load_"//solsym(it),242 . "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,243 . "ave(X)", zsto,zout)244 CALL histdef(nid_tra2, "SConc_"//solsym(it), "SConc_"//solsym(it),245 . "kg/m3", iim,jjphy_nb,nhori, 1,1,1, -99, 32,246 . "ave(X)", zsto,zout)247 do la=1,las248 CALL histdef(nid_tra2, "OD"//cla(la)//"_"//solsym(it), "OD"//cla(la)//"_"//solsym(it),249 . "opt. depth", iim,jjphy_nb,nhori, 1,1,1, -99, 32,250 . "ave(X)", zsto,zout)251 252 CALL histdef(nid_tra2, "OD3D"//cla(la)//"_"//solsym(it), "OD3d_"//cla(la)//"_"//solsym(it),253 . "opt. depth 3D", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,254 . "ave(X)", zsto,zout)255 enddo256 257 CALL histdef(nid_tra2, "CLOAD_"//solsym(it), "3D LOAD_"//solsym(it),258 . "kg/m2 resp. #/m2", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,259 . "ave(X)", zsto,zout)260 261 CALL histdef(nid_tra2, "MD_"//solsym(it), "MD_"//solsym(it),262 . "median diameter [m]", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,263 . "ave(X)", zsto,zout)264 265 CALL histdef(nid_tra2, "MDW_"//solsym(it), "MDW_"//solsym(it),266 . "wet median diameter [m]", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,267 . "ave(X)", zsto,zout)268 269 270 271 CALL histdef(nid_tra3, "Inst_Load_"//solsym(it), "Inst_Load_"//solsym(it),272 . "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,273 . "inst(X)", zout,zout)274 275 endif276 277 #endif278 CALL histdef(nid_tra, solsym(it), solsym(it), "VMR",279 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,280 . "ave(X)", zsto,zout)281 #else282 35 iq=it+2 283 36 iiq=niadv(iq) … … 306 59 c---fin Olivia 307 60 308 #endif309 61 ENDDO 310 62 311 #ifdef INCA312 CALL histdef(nid_tra, "O3_column", "O3_column",313 . "DU", iim,jjphy_nb,nhori, 1,1,1, -99, 32,314 . "ave(X)", zsto,zout)315 CALL histdef(nid_tra, "CO_column", "CO_column",316 . "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,317 . "ave(X)", zsto,zout)318 CALL histdef(nid_tra, "CH4_column", "CH4_column",319 . "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,320 . "ave(X)", zsto,zout)321 CALL histdef(nid_tra, "NO2_column", "NO2_column",322 . "10^15 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,323 . "ave(X)", zsto,zout)324 CALL histdef(nid_tra, "O3_ste", "O3_ste",325 . "CM-2 S-1", iim,jjphy_nb,nhori, 1,1,1, -99, 32,326 . "ave(X)", zsto,zout)327 CALL histdef(nid_tra, "O3_prod", "O3_prod", "CM-3 S-1",328 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,329 . "ave(X)", zsto,zout)330 CALL histdef(nid_tra, "O3_loss", "O3_loss", "CM-3 S-1",331 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,332 . "ave(X)", zsto,zout)333 334 ! Special variables for daytime averaging335 ! CALL histdef(nid_tra, "day_cnt", "day_cnt", "-",336 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,337 ! . "t_sum(X)", zsto,zout)338 ! CALL histdef(nid_tra, "NO_day", "NO_day", "VMR",339 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,340 ! . "t_sum(X)", zsto,zout)341 342 343 344 !! Ajout Anne345 #ifdef INCA_AER346 ! for sulfur cycle347 CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",348 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,349 . "ave(X)", zsto,zout)350 CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",351 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,352 . "ave(X)", zsto,zout)353 CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",354 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,355 . "ave(X)", zsto,zout)356 CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",357 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,358 . "ave(X)", zsto,zout)359 CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",360 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,361 . "ave(X)", zsto,zout)362 CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",363 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,364 . "ave(X)", zsto,zout)365 CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",366 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,367 . "ave(X)", zsto,zout)368 CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",369 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,370 . "ave(X)", zsto,zout)371 CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",372 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,373 . "ave(X)", zsto,zout)374 375 CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",376 . "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,377 . "ave(X)", zsto,zout)378 CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",379 . "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,380 . "ave(X)", zsto,zout)381 CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",382 . "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,383 . "ave(X)", zsto,zout)384 CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",385 . "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,386 . "ave(X)", zsto,zout)387 388 CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",389 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,390 . "ave(X)", zsto,zout)391 #endif392 393 #ifdef INCA_CH4394 DO it=1, grpcnt395 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",396 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,397 . "ave(X)", zsto,zout)398 ENDDO399 400 #endif401 402 403 #ifdef INCA_NMHC404 CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",405 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,406 . "ave(X)", zsto,zout)407 CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",408 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,409 . "ave(X)", zsto,zout)410 CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",411 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,412 . "ave(X)", zsto,zout)413 CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",414 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,415 . "ave(X)", zsto,zout)416 CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",417 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,418 . "ave(X)", zsto,zout)419 CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",420 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,421 . "ave(X)", zsto,zout)422 CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",423 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,424 . "ave(X)", zsto,zout)425 426 DO it=1, grpcnt427 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",428 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,429 . "ave(X)", zsto,zout)430 ENDDO431 432 433 #endif434 ! Fin ajout Anne435 #else436 63 CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", 437 64 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, … … 493 120 . "ave(X)", zsto,zout) 494 121 495 #endif496 122 c 497 123 CALL histend(nid_tra) 498 #ifdef INCA_AER499 CALL histend(nid_tra2)500 CALL histend(nid_tra3)501 #endif502 124 ndex2d = 0 503 125 ndex3d = 0 504 126 ndex = 0 127 #endif
Note: See TracChangeset
for help on using the changeset viewer.