Changeset 726 for LMDZ4/branches/V3_test/libf/phylmd/phytrac.F
- Timestamp:
- Nov 6, 2006, 4:51:16 PM (18 years ago)
- File:
-
- 1 edited
-
LMDZ4/branches/V3_test/libf/phylmd/phytrac.F (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/V3_test/libf/phylmd/phytrac.F
r704 r726 67 67 USE iophy 68 68 USE vampir 69 #ifdef INCA 70 USE sflx 71 USE chem_tracnm 72 USE species_names 73 USE chem_mods 74 #ifdef INCA_NMHC 75 USE pht_tables, ONLY : jrates 76 USE lightning, ONLY : prod_light 77 #endif 78 #ifdef INCA_CH4 79 USE pht_tables, ONLY : jrates 80 USE lightning, ONLY : prod_light 81 #endif 82 USE transport_controls, ONLY : conv_flg, pbl_flg 83 USE airplane_src, ONLY : ptrop 84 #ifdef INCA_AER 85 USE AEROSOL_MOD, only : ntr,trmx,trnx 86 USE AEROSOL_DIAG,only : cla,las,tausum,angst,aload,cload,totaerh2o,tau, 87 $ emiss20,sconc,scavcoef_st,scavcoef_cv 88 $ ,cload05ss ,cload05bc ,cload05pom ,cload05dust ,cload05so4 89 $ ,cload125ss ,cload125bc ,cload125pom ,cload125dust ,cload125so4 90 USE AEROSOL_PROGNOS, ONLY : md,mdw 91 USE AEROSOL_METEO, only : airm 92 #endif 93 #ifdef INCA_NMHC 94 USE RESISTANCE_DIAGNOSE, ONLY : surf_alb, sol_irrad, surf_temp, surf_wind, 95 $ aero_resist, lamin_resist, surf_resist 96 #endif 97 #endif 69 98 70 IMPLICIT none 99 71 c====================================================================== … … 165 137 #ifdef INCA 166 138 REAL flxmass_w(klon,klev) 139 CHARACTER(len=8) :: solsym(nqmax) 167 140 #endif 168 141 c integer iflag_con … … 230 203 cAA Pour l'instant seuls les cas du rn et du pb ont ete envisages. 231 204 232 REAL source(klon ) ! a voir lorsque le flux est prescrit205 REAL source(klon,nqmax) ! a voir lorsque le flux est prescrit 233 206 cAA 234 207 cAA Pour la source de radon et son reservoir de sol … … 365 338 INTEGER :: lastgas 366 339 INTEGER :: ncsec 367 368 INTEGER :: prt_flag_ts(nbtr)=(/ 369 #ifdef INCA_CH4 370 . 1,1,1,0,0,1,1,1,1,1, 371 . 0,1,0,0,0,0,0,1,0,0, 372 . 0,1,1,1,1,0,1,1,1,0, 373 . 1,1,1,1,1,1,1,1,1,1, 374 . 1,0,0 375 #ifdef INCA_AER 376 . ,1,1,1,1,0,1,1,1,1,0, 377 . 1,1,1,1,1,1,0,1,0,1, 378 . 1,1,1,1,0,1,0,1,1,1 379 #endif 380 #endif 381 #ifdef INCA_NMHC 382 . 1,1,1,1,1,1,1,1,1,1, 383 . 1,1,1,1,1,1,1,1,1,1, 384 . 1,1,1,1,1,1,1,1,1,1, 385 . 1,1,1,1,1,1,1,1,1,1, 386 . 1,1,1,1,1,1,1,1,1,1, 387 . 1,1,1,1,1,1,1,1,1,1, 388 . 1,1,1,1,1,1,1,1,1,1, 389 . 1,1,1,1,1,1,1,1,1,1, 390 . 1,1,1,1,1,1,1 391 #ifdef INCA_AER 392 . ,1,1,1,1,0,1,1,1,1,0, 393 . 1,1,1,1,1,1,0,1,0,1, 394 . 1,1,1,1,0,1,0,1,1,1 395 #endif 396 #endif 397 #if defined(INCA_AER) && !defined(INCA_CH4) && !defined(INCA_NMHC) 398 . 1,1,1,1,1,1,1,1,1,1, 399 . 1,1,1,1,1,1,1,1,1,1, 400 . 1,1,1,1,1,1,1,1,1 401 #endif 402 #if defined(INCA) && !defined(INCA_CH4) && !defined(INCA_NMHC) && !defined(INCA_AER) 403 . 1,1,1,1,1,1,1,1,1,1, 404 . 1 405 #endif 406 407 . /) 408 340 INTEGER :: prt_flag_ts(nbtr) 409 341 410 342 REAL, PARAMETER :: dry_mass = 28.966 … … 413 345 REAL :: calday 414 346 REAL :: pdel(klon,klev) 415 REAL :: dummy(klon,klev) = 0.347 REAL :: dummy(klon,klev) 416 348 #endif 417 349 #ifdef INCA_AER … … 420 352 c 421 353 c====================================================================== 354 355 #ifdef INCA 356 prt_flag_ts(:)=(/ 357 #ifdef INCA_CH4 358 . 1,1,1,0,0,1,1,1,1,1, 359 . 0,1,0,0,0,0,0,1,0,0, 360 . 0,1,1,1,1,0,1,1,1,0, 361 . 1,1,1,1,1,1,1,1,1,1, 362 . 1,0,0 363 #ifdef INCA_AER 364 . ,1,1,1,1,0,1,1,1,1,0, 365 . 1,1,1,1,1,1,0,1,0,1, 366 . 1,1,1,1,0,1,0,1,1,1 367 #endif 368 #endif 369 #ifdef INCA_NMHC 370 . 1,1,1,1,1,1,1,1,1,1, 371 . 1,1,1,1,1,1,1,1,1,1, 372 . 1,1,1,1,1,1,1,1,1,1, 373 . 1,1,1,1,1,1,1,1,1,1, 374 . 1,1,1,1,1,1,1,1,1,1, 375 . 1,1,1,1,1,1,1,1,1,1, 376 . 1,1,1,1,1,1,1,1,1,1, 377 . 1,1,1,1,1,1,1,1,1,1, 378 . 1,1,1,1,1,1,1 379 #ifdef INCA_AER 380 . ,1,1,1,1,0,1,1,1,1,0, 381 . 1,1,1,1,1,1,0,1,0,1, 382 . 1,1,1,1,0,1,0,1,1,1 383 #endif 384 #endif 385 #if defined(INCA_AER) && !defined(INCA_CH4) && !defined(INCA_NMHC) 386 . 1,1,1,1,1,1,1,1,1,1, 387 . 1,1,1,1,1,1,1,1,1,1, 388 . 1,1,1,1,1,1,1,1,1 389 #endif 390 #if defined(INCA) && !defined(INCA_CH4) && !defined(INCA_NMHC) && !defined(INCA_AER) 391 . 1,1,1,1,1,1,1,1,1,1, 392 . 1 393 #endif 394 395 . /) 396 dummy(:,:) = 0. 397 398 #endif 422 399 modname='phytrac' 423 400 … … 549 526 $ rneb, ! for chimiaq 550 527 $ t_seri, ! for chimiaq 551 $ rh) 528 $ rh, 529 $ lafin) 552 530 ! fin changement anne 553 531 … … 555 533 556 534 CALL chemmain (tr_seri, !mmr 557 $ nas, !nas558 535 $ nstep, !nstep 559 536 $ calday, !calday … … 590 567 $ obuf, !obuf 591 568 $ iip1, !nx 592 $ jjp1) !ny 569 $ jjp1, !ny 570 $ source, 571 $ solsym) 593 572 #ifdef INCAINFO 594 573 #ifdef INCA_AER … … 797 776 C CALL minmaxqfi(tr_seri(1,1,it),0.,1.e33,'cltracrn it='//itn) 798 777 else ! couche limite avec flux prescrit 799 #ifdef INCA 800 DO k = 1, klon 801 source(k) = eflux(k,it)-dflux(k,it) 802 END DO 803 #else 778 #ifndef INCA 804 779 805 780 Cmaf provisoire source / traceur a creer 806 781 DO i=1, klon 807 source(i ) = 0.0 ! pas de source, pour l'instant782 source(i,it) = 0.0 ! pas de source, pour l'instant 808 783 ENDDO 809 784 C 810 785 #endif 811 786 CALL cltrac(pdtphys, coefh,t_seri, 812 s tr_seri(1,1,it), source ,787 s tr_seri(1,1,it), source(:,it), 813 788 e paprs, pplay, delp, 814 789 s d_tr_cl(1,1,it))
Note: See TracChangeset
for help on using the changeset viewer.
