Changeset 759 for LMDZ4/branches/LMDZ4_V2_patch
- Timestamp:
- May 22, 2007, 3:55:22 PM (18 years ago)
- Location:
- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd
- Files:
-
- 3 edited
-
ini_histrac.h (modified) (9 diffs)
-
phytrac.F (modified) (4 diffs)
-
write_histrac.h (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h
r740 r759 17 17 18 18 #ifdef INCA_AER 19 CALL histbeg("histrac_aer", iim,zx_lon , jjm+1,zx_lat,19 CALL histbeg("histrac_aer", iim,zx_lon(:,1), jjm+1,zx_lat(1,:), 20 20 . 1,iim,1,jjm+1, itau_phy, zjulian, pdtphys, 21 21 . nhori, nid_tra2) … … 25 25 . nhori, nid_tra3) 26 26 27 call histvert(nid_tra2, "presnivs", " presnivs", "mb",27 call histvert(nid_tra2, "presnivs", "Vertical levels", "mb", 28 28 . klev, presnivs, nvert) 29 29 call histvert(nid_tra3, "presnivs", "presnivs", "mb", … … 92 92 #ifdef INCA 93 93 #ifdef INCAINFO 94 95 #ifdef INCA_CH4 94 96 DO it=1, phtcnt 95 97 WRITE(str2,'(i2.2)') it … … 111 113 . "ave(X)", zsto,zout) 112 114 ENDDO 113 115 #endif 116 #ifdef INCA_NMHC 117 DO it=1, phtcnt 118 WRITE(str2,'(i2.2)') it 119 CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1", 120 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 121 . "ave(X)", zsto,zout) 122 ENDDO 123 DO it=1, hetcnt 124 WRITE(str2,'(i2.2)') it 125 CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1", 126 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 127 . "ave(X)", zsto,zout) 128 ENDDO 129 DO it=1, extcnt 130 WRITE(str2,'(i2.2)') it 131 CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1", 132 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 133 . "ave(X)", zsto,zout) 134 ENDDO 135 #endif 114 136 DO it=1, nfs 115 137 WRITE(str2,'(i2.2)') it … … 120 142 121 143 #else 144 145 #ifdef INCA_NMHC 122 146 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1", 123 147 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, … … 139 163 . "ave(X)", zsto,zout) 140 164 #endif 141 142 DO it=1, grpcnt 143 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR", 144 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 145 . "ave(X)", zsto,zout) 146 ENDDO 165 #ifdef INCA_CH4 166 CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1", 167 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 168 . "ave(X)", zsto,zout) 169 CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1", 170 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 171 . "ave(X)", zsto,zout) 172 CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1", 173 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 174 . "ave(X)", zsto,zout) 175 CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1", 176 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 177 . "ave(X)", zsto,zout) 178 CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1", 179 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 180 . "ave(X)", zsto,zout) 181 CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1", 182 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 183 . "ave(X)", zsto,zout) 184 #endif 185 #endif 186 187 ! DO it=1, grpcnt 188 ! CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR", 189 ! . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 190 ! . "ave(X)", zsto,zout) 191 ! ENDDO 147 192 #endif 148 193 … … 349 394 #endif 350 395 396 #ifdef INCA_CH4 397 DO it=1, grpcnt 398 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR", 399 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 400 . "ave(X)", zsto,zout) 401 ENDDO 402 403 #endif 404 405 351 406 #ifdef INCA_NMHC 352 407 CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1", … … 371 426 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 372 427 . "ave(X)", zsto,zout) 428 429 DO it=1, grpcnt 430 CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR", 431 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 432 . "ave(X)", zsto,zout) 433 ENDDO 434 435 373 436 #endif 374 437 ! Fin ajout Anne … … 376 439 CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", 377 440 . iim,jjm+1,nhori, 1,1,1, -99, 32, 378 . "inst(X)", zsto,zout)441 . "inst(X)", zout,zout) 379 442 380 443 CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", 381 444 . iim,jjm+1,nhori, 1,1,1, -99, 32, 382 . "inst(X)", z sto,zout)445 . "inst(X)", zout,zout) 383 446 CALL histdef(nid_tra, "psrf1", "nature sol", "-", 384 447 . iim,jjm+1,nhori, 1,1,1, -99, 32, 385 . "inst(X)", z sto,zout)448 . "inst(X)", zout,zout) 386 449 CALL histdef(nid_tra, "psrf2", "nature sol", "-", 387 450 . iim,jjm+1,nhori, 1,1,1, -99, 32, 388 . "inst(X)", z sto,zout)451 . "inst(X)", zout,zout) 389 452 CALL histdef(nid_tra, "psrf3", "nature sol", "-", 390 453 . iim,jjm+1,nhori, 1,1,1, -99, 32, 391 . "inst(X)", z sto,zout)454 . "inst(X)", zout,zout) 392 455 CALL histdef(nid_tra, "psrf4", "nature sol", "-", 393 456 . iim,jjm+1,nhori, 1,1,1, -99, 32, 394 . "inst(X)", z sto,zout)457 . "inst(X)", zout,zout) 395 458 CALL histdef(nid_tra, "ftsol1", "temper sol", "-", 396 459 . iim,jjm+1,nhori, 1,1,1, -99, 32, 397 . "inst(X)", z sto,zout)460 . "inst(X)", zout,zout) 398 461 CALL histdef(nid_tra, "ftsol2", "temper sol", "-", 399 462 . iim,jjm+1,nhori, 1,1,1, -99, 32, 400 . "inst(X)", z sto,zout)463 . "inst(X)", zout,zout) 401 464 CALL histdef(nid_tra, "ftsol3", "temper sol", "-", 402 465 . iim,jjm+1,nhori, 1,1,1, -99, 32, 403 . "inst ", zsto,zout)466 . "inst(X)", zout,zout) 404 467 CALL histdef(nid_tra, "ftsol4", "temper sol", "-", 405 468 . iim,jjm+1,nhori, 1,1,1, -99, 32, 406 . "inst(X)", z sto,zout)469 . "inst(X)", zout,zout) 407 470 CALL histdef(nid_tra, "pplay", "flux u mont","-", 408 471 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 409 . "inst(X)", z sto,zout)472 . "inst(X)", zout,zout) 410 473 CALL histdef(nid_tra, "t", "flux u mont","-", 411 474 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, 412 . "inst(X)", z sto,zout)475 . "inst(X)", zout,zout) 413 476 CALL histdef(nid_tra, "mfu", "flux u mont","-", 414 477 . iim,jjm+1,nhori, klev,1,klev,nvert, 32, -
LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/phytrac.F
r740 r759 68 68 USE species_names 69 69 USE chem_mods 70 #ifdef INCA_NMHC 70 71 USE pht_tables, ONLY : jrates 72 USE lightning, ONLY : prod_light 73 #endif 74 #ifdef INCA_CH4 75 USE pht_tables, ONLY : jrates 76 USE lightning, ONLY : prod_light 77 #endif 71 78 USE transport_controls, ONLY : conv_flg, pbl_flg 72 79 USE airplane_src, ONLY : ptrop 73 USE lightning, ONLY : prod_light74 80 #ifdef INCA_AER 75 81 USE AEROSOL_MOD, only : ntr,trmx,trnx 76 82 USE AEROSOL_DIAG,only : cla,las,tausum,angst,aload,cload,totaerh2o,tau, 77 $ emiss20,sconc,scavcoef_st,scavcoef_cv 83 $ emiss20,sconc,scavcoef_st,scavcoef_cv, lambda 78 84 $ ,cload05ss ,cload05bc ,cload05pom ,cload05dust ,cload05so4 79 85 $ ,cload125ss ,cload125bc ,cload125pom ,cload125dust ,cload125so4 … … 108 114 #include "comgeomphy.h" 109 115 #include "advtrac.h" 116 110 117 c====================================================================== 111 118 … … 335 342 INTEGER :: ncsec 336 343 337 INTEGER :: prt_flag_ts(nbtr)=(/ 1,1,1344 INTEGER :: prt_flag_ts(nbtr)=(/ 338 345 #ifdef INCA_CH4 339 . 346 . 1,1,1,0,0,1,1,1,1,1, 340 347 . 0,1,0,0,0,0,0,1,0,0, 341 348 . 0,1,1,1,1,0,1,1,1,0, 342 349 . 1,1,1,1,1,1,1,1,1,1, 343 350 . 1,0,0 344 #ifdef INCA_AER 345 . ,1,1,1,1,0,1,1,1 346 #endif 347 #endif 348 #ifdef INCA_AER 349 c aerosol tracers 350 . ,1,0,1,1,1,1,1,1,0,1, 351 . 0,1,1,1,1,1,0,1,0,1,1,1 351 #ifdef INCA_AER 352 . ,1,1,1,1,0,1,1,1,1,0, 353 . 1,1,1,1,1,1,0,1,0,1, 354 . 1,1,1,1,0,1,0,1,1,1 355 #endif 352 356 #endif 353 357 #ifdef INCA_NMHC 354 . , 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, 1, 1, 1, 1, 360 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 361 . 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 362 . 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,1,1,1,1, 360 . 1,1,1,1,1,1,1,1,1,1, 361 . 1,1,1,1,1,1,1,1,1,1, 362 . 1,1,1,1,1,1,1,1,1,1, 363 . 1,1,1,1,1,1,1,1,1,1, 364 . 1,1,1,1,1,1,1,1,1,1, 365 . 1,1,1,1,1,1,1,1,1,1, 366 . 1,1,1,1,1,1,1 367 #ifdef INCA_AER 368 . ,1,1,1,1,0,1,1,1,1,0, 369 . 1,1,1,1,1,1,0,1,0,1, 370 . 1,1,1,1,0,1,0,1,1,1 371 #endif 372 #endif 373 #if defined(INCA_AER) && !defined(INCA_CH4) && !defined(INCA_NMHC) 374 . 1,1,1,0,1,1,1,1,0,1, 375 . 1,1,1,1,1,0,1,0,1,1, 376 . 1,1,1,0,1,0,1,1,1 377 #endif 378 #if defined(INCA) && !defined(INCA_CH4) && !defined(INCA_NMHC) && !defined(INCA_AER) 379 . 1,1,1,1,1,1,1,1,1,1, 380 . 1 363 381 #endif 364 382 . /) … … 374 392 #ifdef INCA_AER 375 393 integer la 394 character*20 form 376 395 #endif 377 396 c -
LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h
r740 r759 62 62 #ifdef INCA 63 63 #ifdef INCAINFO 64 #ifdef INCA_NMHC 64 65 DO it=1, phtcnt 65 66 WRITE(str2,'(i2.2)') it … … 85 86 . iim*(jjm+1)*klev,ndex3d) 86 87 ENDDO 88 #endif 89 #ifdef INCA_CH4 90 DO it=1, phtcnt 91 WRITE(str2,'(i2.2)') it 92 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it), 93 . zx_tmp_3d) 94 CALL histwrite(nid_tra,"j"//str2,itau_w,zx_tmp_3d, 95 . iim*(jjm+1)*klev,ndex3d) 96 ENDDO 97 98 DO it=1, hetcnt 99 WRITE(str2,'(i2.2)') it 100 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it), 101 . zx_tmp_3d) 102 CALL histwrite(nid_tra,"w"//str2,itau_w,zx_tmp_3d, 103 . iim*(jjm+1)*klev,ndex3d) 104 ENDDO 105 106 DO it=1, extcnt 107 WRITE(str2,'(i2.2)') it 108 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it), 109 . zx_tmp_3d) 110 CALL histwrite(nid_tra,"ext"//str2,itau_w,zx_tmp_3d, 111 . iim*(jjm+1)*klev,ndex3d) 112 ENDDO 113 #endif 87 114 88 115 DO it=1, nfs … … 94 121 ENDDO 95 122 #else 123 #ifdef INCA_NMHC 96 124 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2), 97 125 . zx_tmp_3d) … … 124 152 . iim*(jjm+1)*klev,ndex3d) 125 153 #endif 126 DO it=1, grpcnt 127 128 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 129 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 130 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 131 . iim*(jjm+1)*klev,ndex3d) 132 ENDDO 154 #ifdef INCA_CH4 155 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2), 156 . zx_tmp_3d) 157 CALL histwrite(nid_tra,"jO3",itau_w,zx_tmp_3d, 158 . iim*(jjm+1)*klev,ndex3d) 159 160 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4), 161 . zx_tmp_3d) 162 CALL histwrite(nid_tra,"jNO2",itau_w,zx_tmp_3d, 163 . iim*(jjm+1)*klev,ndex3d) 164 165 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13), 166 . zx_tmp_3d) 167 CALL histwrite(nid_tra,"jH2O2",itau_w,zx_tmp_3d, 168 . iim*(jjm+1)*klev,ndex3d) 169 170 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1), 171 . zx_tmp_3d) 172 CALL histwrite(nid_tra,"wHNO3",itau_w,zx_tmp_3d, 173 . iim*(jjm+1)*klev,ndex3d) 174 175 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1), 176 . zx_tmp_3d) 177 CALL histwrite(nid_tra,"kN2O5",itau_w,zx_tmp_3d, 178 . iim*(jjm+1)*klev,ndex3d) 179 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1), 180 . zx_tmp_3d) 181 CALL histwrite(nid_tra,"LghtNO",itau_w,zx_tmp_3d, 182 . iim*(jjm+1)*klev,ndex3d) 183 #endif 184 #endif 185 ! DO it=1, grpcnt 186 ! 187 ! CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 188 ! zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 189 ! CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 190 ! . iim*(jjm+1)*klev,ndex3d) 191 ! ENDDO 133 192 #endif 134 193 … … 228 287 #endif 229 288 else 230 #ifdef INCA 289 #ifdef INCA_CH4 290 if ( solsym(it) == 'OX' ) then 291 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3) 292 end if 293 #endif 294 #ifdef INCA_NMHC 231 295 if ( solsym(it) == 'OX' ) then 232 296 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3) … … 264 328 265 329 #ifdef INCA 266 !#ifdef INCA_CH4267 330 CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_col(1), zx_tmp_2d) 268 331 CALL histwrite(nid_tra,"O3_column",itau_w,zx_tmp_2d, … … 363 426 #endif 364 427 428 #ifdef INCA_CH4 429 DO it=1, grpcnt 430 431 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 432 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 433 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 434 . iim*(jjm+1)*klev,ndex3d) 435 ENDDO 436 437 #endif 438 439 365 440 #ifdef INCA_NMHC 366 441 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1), … … 398 473 CALL histwrite(nid_tra,"CO_loss",itau_w,zx_tmp_3d, 399 474 . iim*(jjm+1)*klev,ndex) 475 476 DO it=1, grpcnt 477 478 CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d) 479 zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it) 480 CALL histwrite(nid_tra,grpsym(it),itau_w,zx_tmp_3d, 481 . iim*(jjm+1)*klev,ndex3d) 482 ENDDO 483 400 484 401 485 #endif
Note: See TracChangeset
for help on using the changeset viewer.
