Changeset 673 for LMDZ4/trunk/libf
- Timestamp:
- Nov 16, 2005, 3:46:34 PM (19 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 3 edited
-
ini_histrac.h (modified) (10 diffs)
-
phytrac.F (modified) (2 diffs)
-
write_histrac.h (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ini_histrac.h
r666 r673 21 21 . nhori, nid_tra2) 22 22 23 CALL histbeg("histrac_inst", iim,zx_lon (:,1), jjm+1,zx_lat(1,:),23 CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat, 24 24 . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys, 25 25 . nhori, nid_tra3) … … 27 27 call histvert(nid_tra2, "presnivs", "Vertical levels", "mb", 28 28 . klev, presnivs, nvert) 29 30 call histvert(nid_tra3, "presnivs", "Vertical levels", "mb", 29 call histvert(nid_tra3, "presnivs", "presnivs", "mb", 31 30 . klev, presnivs, nvert) 32 31 #endif 32 33 33 34 34 zsto = pdtphys … … 84 84 #ifdef INCA 85 85 #ifdef INCAINFO 86 87 #ifdef INCA_CH4 86 88 DO it=1, phtcnt 87 89 WRITE(str2,'(i2.2)') it … … 96 98 . "ave(X)", zsto,zout) 97 99 ENDDO 98 99 100 DO it=1, extcnt 100 101 WRITE(str2,'(i2.2)') it … … 103 104 . "ave(X)", zsto,zout) 104 105 ENDDO 105 106 #endif 107 #ifdef INCA_NMHC 108 DO it=1, phtcnt 109 WRITE(str2,'(i2.2)') it 110 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 ENDDO 114 DO it=1, hetcnt 115 WRITE(str2,'(i2.2)') it 116 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 ENDDO 120 DO it=1, extcnt 121 WRITE(str2,'(i2.2)') it 122 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 ENDDO 126 #endif 106 127 DO it=1, nfs 107 128 WRITE(str2,'(i2.2)') it … … 112 133 113 134 #else 135 136 #ifdef INCA_NMHC 114 137 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1", 115 138 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, … … 131 154 . "ave(X)", zsto,zout) 132 155 #endif 133 134 DO it=1, grpcnt 135 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR", 136 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 137 . "ave(X)", zsto,zout) 138 ENDDO 156 #ifdef INCA_CH4 157 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 #endif 176 #endif 177 178 ! DO it=1, grpcnt 179 ! 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 ! ENDDO 139 183 #endif 140 184 … … 159 203 C champ 2D 160 204 #ifdef INCA 161 IF ( prt_flag_ts(it) == 0 ) CYCLE162 205 163 206 CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it), … … 341 384 #endif 342 385 386 #ifdef INCA_CH4 387 DO it=1, grpcnt 388 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 ENDDO 392 393 #endif 394 395 343 396 #ifdef INCA_NMHC 344 397 CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1", … … 363 416 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 364 417 . "ave(X)", zsto,zout) 418 419 DO it=1, grpcnt 420 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 ENDDO 424 425 365 426 #endif 366 427 ! Fin ajout Anne -
LMDZ4/trunk/libf/phylmd/phytrac.F
r666 r673 69 69 USE species_names 70 70 USE chem_mods 71 #ifdef INCA_NMHC 71 72 USE pht_tables, ONLY : jrates 73 USE lightning, ONLY : prod_light 74 #endif 75 #ifdef INCA_CH4 76 USE pht_tables, ONLY : jrates 77 USE lightning, ONLY : prod_light 78 #endif 72 79 USE transport_controls, ONLY : conv_flg, pbl_flg 73 80 USE airplane_src, ONLY : ptrop 74 USE lightning, ONLY : prod_light75 81 #ifdef INCA_AER 76 82 USE AEROSOL_MOD, only : ntr,trmx,trnx … … 340 346 INTEGER :: ncsec 341 347 342 INTEGER :: prt_flag_ts(nbtr)=(/1,1,1343 #ifdef INCA_CH4344 . ,0,0,1,1,1,1,1,345 . 0,1,0,0,0,0,0,1,0,0,346 . 0,1,1,1,1,0,1,1,1,0,347 . 1,1,1,1,1,1,1,1,1,1,348 . 1,0,0349 #endif350 #ifdef INCA_NMHC351 . , 1, 1, 1, 1, 1, 1, 1,352 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,353 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,354 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,355 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,356 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,357 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,358 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,359 . 1, 1, 1, 1, 1, 1, 1360 #endif361 #ifdef INCA_AER362 . ,1,1,1,1,0,1,1,1363 #endif364 #ifdef INCA_AER365 c aerosol tracers366 . ,1,0,1,1,1,1,1,1,0,1,367 . 0,1,1,1,1,1,0,1,0,1,1,1368 #endif369 . /)370 348 371 349 -
LMDZ4/trunk/libf/phylmd/write_histrac.h
r646 r673 63 63 #ifdef INCA 64 64 #ifdef INCAINFO 65 #ifdef INCA_NMHC 65 66 DO it=1, phtcnt 66 67 WRITE(str2,'(i2.2)') it … … 86 87 . iim*(jjm+1)*klev,ndex3d) 87 88 ENDDO 89 #endif 90 #ifdef INCA_CH4 91 DO it=1, phtcnt 92 WRITE(str2,'(i2.2)') it 93 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it), 94 . zx_tmp_3d) 95 CALL histwrite(nid_tra,"j"//str2,itau_w,zx_tmp_3d, 96 . iim*(jjm+1)*klev,ndex3d) 97 ENDDO 98 99 DO it=1, hetcnt 100 WRITE(str2,'(i2.2)') it 101 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it), 102 . zx_tmp_3d) 103 CALL histwrite(nid_tra,"w"//str2,itau_w,zx_tmp_3d, 104 . iim*(jjm+1)*klev,ndex3d) 105 ENDDO 106 107 DO it=1, extcnt 108 WRITE(str2,'(i2.2)') it 109 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it), 110 . zx_tmp_3d) 111 CALL histwrite(nid_tra,"ext"//str2,itau_w,zx_tmp_3d, 112 . iim*(jjm+1)*klev,ndex3d) 113 ENDDO 114 #endif 88 115 89 116 DO it=1, nfs … … 95 122 ENDDO 96 123 #else 124 #ifdef INCA_NMHC 97 125 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2), 98 126 . zx_tmp_3d) … … 125 153 . iim*(jjm+1)*klev,ndex3d) 126 154 #endif 127 DO it=1, grpcnt 128 129 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 130 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 131 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 132 . iim*(jjm+1)*klev,ndex3d) 133 ENDDO 155 #ifdef INCA_CH4 156 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2), 157 . zx_tmp_3d) 158 CALL histwrite(nid_tra,"jO3",itau_w,zx_tmp_3d, 159 . iim*(jjm+1)*klev,ndex3d) 160 161 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4), 162 . zx_tmp_3d) 163 CALL histwrite(nid_tra,"jNO2",itau_w,zx_tmp_3d, 164 . iim*(jjm+1)*klev,ndex3d) 165 166 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13), 167 . zx_tmp_3d) 168 CALL histwrite(nid_tra,"jH2O2",itau_w,zx_tmp_3d, 169 . iim*(jjm+1)*klev,ndex3d) 170 171 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1), 172 . zx_tmp_3d) 173 CALL histwrite(nid_tra,"wHNO3",itau_w,zx_tmp_3d, 174 . iim*(jjm+1)*klev,ndex3d) 175 176 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1), 177 . zx_tmp_3d) 178 CALL histwrite(nid_tra,"kN2O5",itau_w,zx_tmp_3d, 179 . iim*(jjm+1)*klev,ndex3d) 180 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1), 181 . zx_tmp_3d) 182 CALL histwrite(nid_tra,"LghtNO",itau_w,zx_tmp_3d, 183 . iim*(jjm+1)*klev,ndex3d) 184 #endif 185 #endif 186 ! DO it=1, grpcnt 187 ! 188 ! CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 189 ! zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 190 ! CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 191 ! . iim*(jjm+1)*klev,ndex3d) 192 ! ENDDO 134 193 #endif 135 194 … … 160 219 C champs 2D 161 220 #ifdef INCA 162 IF ( prt_flag_ts(it) == 0 ) CYCLE 221 163 222 CALL gr_fi_ecrit(1, klon,iim,jjm+1, eflux(1,it),zx_tmp_2d) 164 223 CALL histwrite(nid_tra,"Emi_"//solsym(it),itau_w,zx_tmp_2d, … … 229 288 #endif 230 289 else 231 #ifdef INCA 290 #ifdef INCA_CH4 291 if ( solsym(it) == 'OX' ) then 292 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3) 293 end if 294 #endif 295 #ifdef INCA_NMHC 232 296 if ( solsym(it) == 'OX' ) then 233 297 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3) … … 364 428 #endif 365 429 430 #ifdef INCA_CH4 431 DO it=1, grpcnt 432 433 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 434 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 435 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 436 . iim*(jjm+1)*klev,ndex3d) 437 ENDDO 438 439 #endif 440 441 366 442 #ifdef INCA_NMHC 367 443 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1), … … 399 475 CALL histwrite(nid_tra,"CO_loss",itau_w,zx_tmp_3d, 400 476 . iim*(jjm+1)*klev,ndex) 477 478 DO it=1, grpcnt 479 480 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 481 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 482 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 483 . iim*(jjm+1)*klev,ndex3d) 484 ENDDO 485 401 486 402 487 #endif
Note: See TracChangeset
for help on using the changeset viewer.
