Changeset 766 for LMDZ4/trunk/libf/phylmd/ini_histrac.h
- Timestamp:
- Jun 4, 2007, 4:34:47 PM (17 years ago)
- File:
-
- 1 edited
-
LMDZ4/trunk/libf/phylmd/ini_histrac.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ini_histrac.h
r684 r766 2 2 ! $Header$ 3 3 ! 4 #ifndef INCA 5 c$OMP MASTER 4 6 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) 5 7 c 6 CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon) 7 DO i = 1, iim 8 zx_lon(i,1) = xlon(i+1) 9 zx_lon(i,jjm+1) = xlon(i+1) 10 ENDDO 11 CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat) 12 CALL histbeg("histrac", iim,zx_lon(:,1), jjm+1,zx_lat(1,:), 13 . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys, 8 cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon) 9 cym DO i = 1, iim 10 cym zx_lon(i,1) = xlon(i+1) 11 cym zx_lon(i,jjm+1) = xlon(i+1) 12 cym ENDDO 13 cym CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat) 14 cym CALL histbeg("histrac", iim,zx_lon(:,1), jjm+1,zx_lat(1,:), 15 cym . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys, 16 cym . nhori, nid_tra) 17 CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys, 14 18 . nhori, nid_tra) 15 19 CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", 16 20 . klev, presnivs, nvert) 17 21 18 #ifdef INCA_AER19 CALL histbeg("histrac_aer", iim,zx_lon(:,1), jjm+1,zx_lat(1,:),20 . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,21 . nhori, nid_tra2)22 23 CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,24 . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys,25 . nhori, nid_tra3)26 27 call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",28 . klev, presnivs, nvert)29 call histvert(nid_tra3, "presnivs", "presnivs", "mb",30 . klev, presnivs, nvert)31 #endif32 22 33 23 34 24 zsto = pdtphys 35 zout = ecrit_tra25 zout = pdtphys * FLOAT(ecrit_tra) 36 26 c 37 27 CALL histdef(nid_tra, "phis", "Surface geop. height", "-", 38 . iim,jj m+1,nhori, 1,1,1, -99, 32,28 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 39 29 . "once", zsto,zout) 40 30 c 41 31 CALL histdef(nid_tra, "aire", "Grid area", "-", 42 . iim,jj m+1,nhori, 1,1,1, -99, 32,32 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 43 33 . "once", zsto,zout) 44 #ifdef INCA45 CALL histdef(nid_tra, "ps", "Surface pressure", "Pa",46 . iim,jjm+1,nhori, 1,1,1,-99, 32,47 . "ave(X)", zsto,zout)48 49 CALL histdef(nid_tra, "ptrop", "Tropopause pressure", "Pa",50 . iim,jjm+1,nhori, 1,1,1,-99, 32,51 . "ave(X)", zsto,zout)52 53 C 3d FIELDS54 CALL histdef(nid_tra, "temp", "Air temperature", "K",55 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,56 . "ave(X)", zsto,zout)57 58 CALL histdef(nid_tra, "u", "zonal wind component", "m/s",59 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,60 . "ave(X)", zsto,zout)61 62 CALL histdef(nid_tra, "v", "zonal wind component", "m/s",63 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,64 . "ave(X)", zsto,zout)65 66 CALL histdef(nid_tra, "h2o", "Specific Humidity", "MMR",67 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,68 . "ave(X)", zsto,zout)69 70 CALL histdef(nid_tra, "pmid", "Pressure", "Pa",71 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,72 . "ave(X)", zsto,zout)73 74 CALL histdef(nid_tra, "pdel", "Delta Pressure", "Pa",75 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,76 . "ave(X)", zsto,zout)77 78 ! MS info79 #ifdef INCA_AER80 CALL histdef(nid_tra, "airm", "Air mass", "kg",81 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,82 . "ave(X)", zsto,zout)83 #endif84 #ifdef INCA85 #ifdef INCAINFO86 87 #ifdef INCA_CH488 DO it=1, phtcnt89 WRITE(str2,'(i2.2)') it90 CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",91 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,92 . "ave(X)", zsto,zout)93 ENDDO94 DO it=1, hetcnt95 WRITE(str2,'(i2.2)') it96 CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",97 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,98 . "ave(X)", zsto,zout)99 ENDDO100 DO it=1, extcnt101 WRITE(str2,'(i2.2)') it102 CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",103 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,104 . "ave(X)", zsto,zout)105 ENDDO106 #endif107 #ifdef INCA_NMHC108 DO it=1, phtcnt109 WRITE(str2,'(i2.2)') it110 CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",111 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,112 . "ave(X)", zsto,zout)113 ENDDO114 DO it=1, hetcnt115 WRITE(str2,'(i2.2)') it116 CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",117 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,118 . "ave(X)", zsto,zout)119 ENDDO120 DO it=1, extcnt121 WRITE(str2,'(i2.2)') it122 CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",123 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,124 . "ave(X)", zsto,zout)125 ENDDO126 #endif127 DO it=1, nfs128 WRITE(str2,'(i2.2)') it129 CALL histdef(nid_tra, "INV"//str2, "INV"//str2, "CM-3",130 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,131 . "ave(X)", zsto,zout)132 ENDDO133 134 #else135 136 #ifdef INCA_NMHC137 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",138 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,139 . "ave(X)", zsto,zout)140 CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",141 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,142 . "ave(X)", zsto,zout)143 CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",144 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,145 . "ave(X)", zsto,zout)146 CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",147 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,148 . "ave(X)", zsto,zout)149 CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",150 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,151 . "ave(X)", zsto,zout)152 CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",153 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,154 . "ave(X)", zsto,zout)155 #endif156 #ifdef INCA_CH4157 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",158 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,159 . "ave(X)", zsto,zout)160 CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",161 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,162 . "ave(X)", zsto,zout)163 CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",164 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,165 . "ave(X)", zsto,zout)166 CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",167 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,168 . "ave(X)", zsto,zout)169 CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",170 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,171 . "ave(X)", zsto,zout)172 CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",173 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,174 . "ave(X)", zsto,zout)175 #endif176 #endif177 178 ! DO it=1, grpcnt179 ! CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",180 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,181 ! . "ave(X)", zsto,zout)182 ! ENDDO183 #endif184 185 #ifdef INCA_AER186 187 CALL histdef(nid_tra2, "scavcoef_st","scavcoef_st", "S-1",188 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,189 . "ave(X)", zsto,zout)190 CALL histdef(nid_tra2, "scavcoef_cv","scavcoef_cv", "S-1",191 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,192 . "ave(X)", zsto,zout)193 CALL histdef(nid_tra2, "AngstroemComp","AngstroemComp",194 . "angs comp", iim,jjm+1,nhori, 1,1,1, -99, 32,195 . "ave(X)", zsto,zout)196 197 CALL histdef(nid_tra2, "TOTAERH2O", "Total aerosol water",198 . "kg/m2", iim,jjm+1,nhori, klev,1,klev, nvert, 32,199 . "ave(X)", zsto,zout)200 #endif201 #endif202 34 DO it=1,nqmax 203 35 C champ 2D 204 #ifdef INCA205 IF ( prt_flag_ts(it) == 0 ) CYCLE206 CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),207 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,208 . "ave(X)", zsto,zout)209 CALL histdef(nid_tra, "Dep_"//solsym(it), "Dep_"//solsym(it),210 . "cm/s", iim,jjm+1,nhori, 1,1,1, -99, 32,211 . "ave(X)", zsto,zout)212 #ifdef INCA_AER213 214 CALL histdef(nid_tra2, "Dry_"//solsym(it), "Dry_"//solsym(it),215 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,216 . "ave(X)", zsto,zout)217 218 IF ((it .ge. trmx) .and. (it .le. trnx)) then219 CALL histdef(nid_tra2, "Sed_"//solsym(it), "Sed_"//solsym(it),220 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,221 . "ave(X)", zsto,zout)222 CALL histdef(nid_tra2, "Wet_"//solsym(it), "Wet_"//solsym(it),223 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,224 . "ave(X)", zsto,zout)225 CALL histdef(nid_tra2, "WetST_"//solsym(it), "WetST_"//solsym(it),226 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,227 . "ave(X)", zsto,zout)228 CALL histdef(nid_tra2, "WetCV_"//solsym(it), "WetCV_"//solsym(it),229 . "kg/m2/s", iim,jjm+1,nhori, 1,1,1, -99, 32,230 . "ave(X)", zsto,zout)231 CALL histdef(nid_tra2, "Emi_alt_"//solsym(it), "Emi_alt_"//solsym(it),232 . "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,233 . "ave(X)", zsto,zout)234 CALL histdef(nid_tra2, "Load_"//solsym(it), "Load_"//solsym(it),235 . "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,236 . "ave(X)", zsto,zout)237 CALL histdef(nid_tra2, "SConc_"//solsym(it), "SConc_"//solsym(it),238 . "kg/m3", iim,jjm+1,nhori, 1,1,1, -99, 32,239 . "ave(X)", zsto,zout)240 do la=1,las241 CALL histdef(nid_tra2, "OD"//cla(la)//"_"//solsym(it), "OD"//cla(la)//"_"//solsym(it),242 . "opt. depth", iim,jjm+1,nhori, 1,1,1, -99, 32,243 . "ave(X)", zsto,zout)244 245 CALL histdef(nid_tra2, "OD3D"//cla(la)//"_"//solsym(it), "OD3d_"//cla(la)//"_"//solsym(it),246 . "opt. depth 3D", iim,jjm+1,nhori, klev,1,klev, nvert, 32,247 . "ave(X)", zsto,zout)248 enddo249 250 CALL histdef(nid_tra2, "CLOAD_"//solsym(it), "3D LOAD_"//solsym(it),251 . "kg/m2 resp. #/m2", iim,jjm+1,nhori, klev,1,klev, nvert, 32,252 . "ave(X)", zsto,zout)253 254 CALL histdef(nid_tra2, "MD_"//solsym(it), "MD_"//solsym(it),255 . "median diameter [m]", iim,jjm+1,nhori, klev,1,klev, nvert, 32,256 . "ave(X)", zsto,zout)257 258 CALL histdef(nid_tra2, "MDW_"//solsym(it), "MDW_"//solsym(it),259 . "wet median diameter [m]", iim,jjm+1,nhori, klev,1,klev, nvert, 32,260 . "ave(X)", zsto,zout)261 262 263 264 CALL histdef(nid_tra3, "Inst_Load_"//solsym(it), "Inst_Load_"//solsym(it),265 . "kg/m2", iim,jjm+1,nhori, 1,1,1, -99, 32,266 . "inst(X)", zout,zout)267 268 endif269 270 #endif271 CALL histdef(nid_tra, solsym(it), solsym(it), "VMR",272 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,273 . "ave(X)", zsto,zout)274 #else275 36 iq=it+2 276 37 iiq=niadv(iq) 277 38 CALL histdef(nid_tra, tnom(iq), ttext(iiq), "U/kga", 278 . iim,jj m+1,nhori, klev,1,klev,nvert, 32,39 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 279 40 . "ave(X)", zsto,zout) 280 41 if (lessivage) THEN 281 42 CALL histdef(nid_tra, "fl"//tnom(iq),"Flux "//ttext(iiq), 282 . "U/m2/s",iim,jj m+1,nhori, klev,1,klev,nvert, 32,43 . "U/m2/s",iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 283 44 . "ave(X)", zsto,zout) 284 45 endif … … 287 48 CALL histdef(nid_tra, "d_tr_th_"//tnom(iq), 288 49 . "tendance thermique"// ttext(iiq), "?", 289 . iim,jj m+1,nhori, klev,1,klev,nvert, 32,50 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 290 51 . "ave(X)", zsto,zout) 291 52 CALL histdef(nid_tra, "d_tr_cv_"//tnom(iq), 292 53 . "tendance convection"// ttext(iiq), "?", 293 . iim,jj m+1,nhori, klev,1,klev,nvert, 32,54 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 294 55 . "ave(X)", zsto,zout) 295 56 CALL histdef(nid_tra, "d_tr_cl_"//tnom(iq), 296 57 . "tendance couche limite"// ttext(iiq), "?", 297 . iim,jj m+1,nhori, klev,1,klev,nvert, 32,58 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 298 59 . "ave(X)", zsto,zout) 299 60 c---fin Olivia 300 61 301 #endif302 62 ENDDO 303 63 304 #ifdef INCA 305 CALL histdef(nid_tra, "O3_column", "O3_column", 306 . "DU", iim,jjm+1,nhori, 1,1,1, -99, 32, 307 . "ave(X)", zsto,zout) 308 CALL histdef(nid_tra, "CO_column", "CO_column", 309 . "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32, 310 . "ave(X)", zsto,zout) 311 CALL histdef(nid_tra, "CH4_column", "CH4_column", 312 . "10^18 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32, 313 . "ave(X)", zsto,zout) 314 CALL histdef(nid_tra, "NO2_column", "NO2_column", 315 . "10^15 CM-2", iim,jjm+1,nhori, 1,1,1, -99, 32, 316 . "ave(X)", zsto,zout) 317 CALL histdef(nid_tra, "O3_ste", "O3_ste", 318 . "CM-2 S-1", iim,jjm+1,nhori, 1,1,1, -99, 32, 319 . "ave(X)", zsto,zout) 320 CALL histdef(nid_tra, "O3_prod", "O3_prod", "CM-3 S-1", 321 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 64 CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", 65 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 66 . "inst(X)", zout,zout) 67 68 CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", 69 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 70 . "inst(X)", zout,zout) 71 CALL histdef(nid_tra, "psrf1", "nature sol", "-", 72 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 73 . "inst(X)", zout,zout) 74 CALL histdef(nid_tra, "psrf2", "nature sol", "-", 75 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 76 . "inst(X)", zout,zout) 77 CALL histdef(nid_tra, "psrf3", "nature sol", "-", 78 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 79 . "inst(X)", zout,zout) 80 CALL histdef(nid_tra, "psrf4", "nature sol", "-", 81 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 82 . "inst(X)", zout,zout) 83 CALL histdef(nid_tra, "ftsol1", "temper sol", "-", 84 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 85 . "inst(X)", zout,zout) 86 CALL histdef(nid_tra, "ftsol2", "temper sol", "-", 87 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 88 . "inst(X)", zout,zout) 89 CALL histdef(nid_tra, "ftsol3", "temper sol", "-", 90 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 91 . "inst", zout,zout) 92 CALL histdef(nid_tra, "ftsol4", "temper sol", "-", 93 . iim,jjphy_nb,nhori, 1,1,1, -99, 32, 94 . "inst(X)", zout,zout) 95 CALL histdef(nid_tra, "pplay", "flux u mont","-", 96 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 97 . "inst(X)", zout,zout) 98 CALL histdef(nid_tra, "t", "flux u mont","-", 99 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 100 . "inst(X)", zout,zout) 101 CALL histdef(nid_tra, "mfu", "flux u mont","-", 102 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 322 103 . "ave(X)", zsto,zout) 323 CALL histdef(nid_tra, "O3_loss", "O3_loss", "CM-3 S-1", 324 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 104 CALL histdef(nid_tra, "mfd", "flux u decen","-", 105 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 106 . "ave(X)", zsto,zout) 107 CALL histdef(nid_tra, "en_u", "flux u mont","-", 108 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 109 . "ave(X)", zsto,zout) 110 CALL histdef(nid_tra, "en_d", "flux u mont","-", 111 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 112 . "ave(X)", zsto,zout) 113 CALL histdef(nid_tra, "de_d", "flux u mont","-", 114 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 115 . "ave(X)", zsto,zout) 116 CALL histdef(nid_tra, "de_u", "flux u decen","-", 117 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 118 . "ave(X)", zsto,zout) 119 CALL histdef(nid_tra, "coefh", "turbulent coef","-", 120 . iim,jjphy_nb,nhori, klev,1,klev,nvert, 32, 325 121 . "ave(X)", zsto,zout) 326 122 327 ! Special variables for daytime averaging328 ! CALL histdef(nid_tra, "day_cnt", "day_cnt", "-",329 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,330 ! . "t_sum(X)", zsto,zout)331 ! CALL histdef(nid_tra, "NO_day", "NO_day", "VMR",332 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32,333 ! . "t_sum(X)", zsto,zout)334 335 336 337 !! Ajout Anne338 #ifdef INCA_AER339 ! for sulfur cycle340 CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",341 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,342 . "ave(X)", zsto,zout)343 CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",344 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,345 . "ave(X)", zsto,zout)346 CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",347 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,348 . "ave(X)", zsto,zout)349 CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",350 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,351 . "ave(X)", zsto,zout)352 CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",353 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,354 . "ave(X)", zsto,zout)355 CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",356 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,357 . "ave(X)", zsto,zout)358 CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",359 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,360 . "ave(X)", zsto,zout)361 CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",362 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,363 . "ave(X)", zsto,zout)364 CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",365 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,366 . "ave(X)", zsto,zout)367 368 CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",369 . "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,370 . "ave(X)", zsto,zout)371 CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",372 . "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,373 . "ave(X)", zsto,zout)374 CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",375 . "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,376 . "ave(X)", zsto,zout)377 CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",378 . "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,379 . "ave(X)", zsto,zout)380 381 CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",382 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,383 . "ave(X)", zsto,zout)384 #endif385 386 #ifdef INCA_CH4387 DO it=1, grpcnt388 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",389 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,390 . "ave(X)", zsto,zout)391 ENDDO392 393 #endif394 395 396 #ifdef INCA_NMHC397 CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",398 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,399 . "ave(X)", zsto,zout)400 CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",401 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,402 . "ave(X)", zsto,zout)403 CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",404 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,405 . "ave(X)", zsto,zout)406 CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",407 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,408 . "ave(X)", zsto,zout)409 CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",410 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,411 . "ave(X)", zsto,zout)412 CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",413 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,414 . "ave(X)", zsto,zout)415 CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",416 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,417 . "ave(X)", zsto,zout)418 419 DO it=1, grpcnt420 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",421 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,422 . "ave(X)", zsto,zout)423 ENDDO424 425 426 #endif427 ! Fin ajout Anne428 #else429 CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",430 . iim,jjm+1,nhori, 1,1,1, -99, 32,431 . "inst(X)", zout,zout)432 433 CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",434 . iim,jjm+1,nhori, 1,1,1, -99, 32,435 . "inst(X)", zout,zout)436 CALL histdef(nid_tra, "psrf1", "nature sol", "-",437 . iim,jjm+1,nhori, 1,1,1, -99, 32,438 . "inst(X)", zout,zout)439 CALL histdef(nid_tra, "psrf2", "nature sol", "-",440 . iim,jjm+1,nhori, 1,1,1, -99, 32,441 . "inst(X)", zout,zout)442 CALL histdef(nid_tra, "psrf3", "nature sol", "-",443 . iim,jjm+1,nhori, 1,1,1, -99, 32,444 . "inst(X)", zout,zout)445 CALL histdef(nid_tra, "psrf4", "nature sol", "-",446 . iim,jjm+1,nhori, 1,1,1, -99, 32,447 . "inst(X)", zout,zout)448 CALL histdef(nid_tra, "ftsol1", "temper sol", "-",449 . iim,jjm+1,nhori, 1,1,1, -99, 32,450 . "inst(X)", zout,zout)451 CALL histdef(nid_tra, "ftsol2", "temper sol", "-",452 . iim,jjm+1,nhori, 1,1,1, -99, 32,453 . "inst(X)", zout,zout)454 CALL histdef(nid_tra, "ftsol3", "temper sol", "-",455 . iim,jjm+1,nhori, 1,1,1, -99, 32,456 . "inst(X)", zout,zout)457 CALL histdef(nid_tra, "ftsol4", "temper sol", "-",458 . iim,jjm+1,nhori, 1,1,1, -99, 32,459 . "inst(X)", zout,zout)460 CALL histdef(nid_tra, "pplay", "flux u mont","-",461 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,462 . "inst(X)", zout,zout)463 CALL histdef(nid_tra, "t", "flux u mont","-",464 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,465 . "inst(X)", zout,zout)466 CALL histdef(nid_tra, "mfu", "flux u mont","-",467 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,468 . "ave(X)", zsto,zout)469 CALL histdef(nid_tra, "mfd", "flux u decen","-",470 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,471 . "ave(X)", zsto,zout)472 CALL histdef(nid_tra, "en_u", "flux u mont","-",473 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,474 . "ave(X)", zsto,zout)475 CALL histdef(nid_tra, "en_d", "flux u mont","-",476 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,477 . "ave(X)", zsto,zout)478 CALL histdef(nid_tra, "de_d", "flux u mont","-",479 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,480 . "ave(X)", zsto,zout)481 CALL histdef(nid_tra, "de_u", "flux u decen","-",482 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,483 . "ave(X)", zsto,zout)484 CALL histdef(nid_tra, "coefh", "turbulent coef","-",485 . iim,jjm+1,nhori, klev,1,klev,nvert, 32,486 . "ave(X)", zsto,zout)487 488 #endif489 123 c 490 124 CALL histend(nid_tra) 491 #ifdef INCA_AER492 CALL histend(nid_tra2)493 CALL histend(nid_tra3)494 #endif495 125 ndex2d = 0 496 126 ndex3d = 0 497 127 ndex = 0 128 c$OMP END MASTER 129 #endif
Note: See TracChangeset
for help on using the changeset viewer.
