Changeset 673 for LMDZ4


Ignore:
Timestamp:
Nov 16, 2005, 3:46:34 PM (19 years ago)
Author:
lmdzadmin
Message:

suppression de prt_flag_ts (INCA) AC
MAF

Location:
LMDZ4/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/ini_histrac.h

    r666 r673  
    2121     .                 nhori, nid_tra2)
    2222
    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,
    2424     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
    2525     .                 nhori, nid_tra3)
     
    2727         call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",
    2828     .                 klev, presnivs, nvert)
    29 
    30          call histvert(nid_tra3, "presnivs", "Vertical levels", "mb",
     29         call histvert(nid_tra3, "presnivs", "presnivs", "mb",
    3130     .                 klev, presnivs, nvert)
    3231#endif
     32
    3333
    3434         zsto = pdtphys
     
    8484#ifdef INCA
    8585#ifdef INCAINFO
     86
     87#ifdef INCA_CH4
    8688         DO it=1, phtcnt
    8789         WRITE(str2,'(i2.2)') it
     
    9698     .                "ave(X)", zsto,zout)
    9799         ENDDO
    98 
    99100         DO it=1, extcnt
    100101         WRITE(str2,'(i2.2)') it
     
    103104     .                "ave(X)", zsto,zout)
    104105         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
    106127         DO it=1, nfs
    107128         WRITE(str2,'(i2.2)') it
     
    112133
    113134#else
     135
     136#ifdef INCA_NMHC
    114137         CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
    115138     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     
    131154     .                "ave(X)", zsto,zout)
    132155#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
    139183#endif
    140184
     
    159203C champ 2D
    160204#ifdef INCA
    161          IF ( prt_flag_ts(it) == 0 ) CYCLE
    162205
    163206         CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),
     
    341384#endif
    342385
     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
    343396#ifdef INCA_NMHC
    344397         CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
     
    363416     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    364417     .                "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
    365426#endif
    366427! Fin ajout Anne
  • LMDZ4/trunk/libf/phylmd/phytrac.F

    r666 r673  
    6969      USE species_names
    7070      USE chem_mods
     71#ifdef INCA_NMHC
    7172      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
    7279      USE transport_controls, ONLY : conv_flg, pbl_flg
    7380      USE airplane_src, ONLY : ptrop
    74       USE lightning, ONLY : prod_light
    7581#ifdef INCA_AER
    7682      USE AEROSOL_MOD, only : ntr,trmx,trnx
     
    340346      INTEGER           :: ncsec
    341347
    342       INTEGER           :: prt_flag_ts(nbtr)=(/1,1,1
    343 #ifdef INCA_CH4
    344      .                                              ,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,0
    349 #endif
    350 #ifdef INCA_NMHC
    351      .                                                 , 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, 1
    360 #endif
    361 #ifdef INCA_AER
    362      .                                         ,1,1,1,1,0,1,1,1
    363 #endif
    364 #ifdef INCA_AER
    365 c aerosol tracers
    366      .                                        ,1,0,1,1,1,1,1,1,0,1,
    367      .                                         0,1,1,1,1,1,0,1,0,1,1,1
    368 #endif
    369      .                                         /)
    370348
    371349
  • LMDZ4/trunk/libf/phylmd/write_histrac.h

    r646 r673  
    6363#ifdef INCA
    6464#ifdef INCAINFO
     65#ifdef INCA_NMHC
    6566      DO it=1, phtcnt
    6667      WRITE(str2,'(i2.2)') it
     
    8687     .                                   iim*(jjm+1)*klev,ndex3d)
    8788      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
    88115
    89116      DO it=1, nfs
     
    95122      ENDDO
    96123#else
     124#ifdef INCA_NMHC
    97125      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
    98126     .     zx_tmp_3d)
     
    125153     .                                   iim*(jjm+1)*klev,ndex3d)
    126154#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
    134193#endif
    135194
     
    160219C champs 2D
    161220#ifdef INCA
    162       IF ( prt_flag_ts(it) == 0 ) CYCLE
     221
    163222      CALL gr_fi_ecrit(1, klon,iim,jjm+1, eflux(1,it),zx_tmp_2d)
    164223      CALL histwrite(nid_tra,"Emi_"//solsym(it),itau_w,zx_tmp_2d,
     
    229288#endif
    230289       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
    232296       if ( solsym(it) == 'OX' ) then
    233297       zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
     
    364428#endif
    365429
     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
    366442#ifdef INCA_NMHC
    367443      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1),
     
    399475      CALL histwrite(nid_tra,"CO_loss",itau_w,zx_tmp_3d,
    400476     .                                   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
    401486
    402487#endif
Note: See TracChangeset for help on using the changeset viewer.