Changeset 616 for LMDZ4


Ignore:
Timestamp:
Apr 13, 2005, 10:58:34 AM (20 years ago)
Author:
lmdzadmin
Message:

Mise a jour pour INCA.2.0 Anne C
MAFi+LF

Location:
LMDZ4/trunk/libf
Files:
2 added
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/trunk/libf/dyn3d/advtrac.F

    r566 r616  
    44c
    55c
    6 #ifdef INCA_CH4
     6#ifdef INCA
    77      SUBROUTINE advtrac(pbaru,pbarv ,
    88     *                   p,  masse,q,iapptrac,teta,
    … …  
    1616     *                  pk)
    1717#endif
    18 
    1918c     Auteur :  F. Hourdin
    2019c
    … …  
    5150      REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
    5251      REAL pk(ip1jmp1,llm)
    53 #ifdef INCA_CH4
     52#ifdef INCA
    5453      INTEGER            :: hadv_flg(nqmx)
    55       REAL               :: mmt_adj(ip1jmp1,llm)
     54      REAL               :: mmt_adj(ip1jmp1,llm,1)
    5655      REAL               :: flxw(ip1jmp1,llm)
    5756#endif
    … …  
    129128        CALL groupe( massem, pbaruc,pbarvc, pbarug,pbarvg,wg )
    130129
    131 #ifdef INCA_CH4
     130#ifdef INCA
    132131      ! ... Flux de masse diaganostiques traceurs
    133132      flxw = wg / FLOAT(iapp_tracvl)
    … …  
    214213        else if (iadv(iq).eq.20) then
    215214            call pentes_ini (q(1,1,iq),wg,massem,pbarug,pbarvg,0)
    216 #ifdef INCA_CH4
     215#ifdef INCA
    217216       do iiq = iq+1, iq+3
    218217         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
    … …  
    232231           call  prather(q(1,1,iq),wg,massem,pbarug,pbarvg,
    233232     s                     n,dtbon)
    234 #ifdef INCA_CH4
     233#ifdef INCA
    235234       do iiq = iq+1, iq+9
    236235         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)

  • LMDZ4/trunk/libf/dyn3d/caladvtrac.F

    r567 r616  
    44c
    55c
    6 #ifdef INCA_CH4
     6#ifdef INCA
    77            SUBROUTINE caladvtrac(q,pbaru,pbarv ,
    88     *                   p ,masse, dq ,  teta,
    … …  
    4141      REAL p( ip1jmp1,llmp1),q( ip1jmp1,llm,nqmx),dq( ip1jmp1,llm,2 )
    4242      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    43 #ifdef INCA_CH4
     43#ifdef INCA
    4444      INTEGER            :: hadv_flg(nqmx)
    4545      REAL               :: mmt_adj(iip1,jjp1,llm,1)
    … …  
    7373c   advection
    7474
    75 #ifdef INCA_CH4
     75#ifdef INCA
    7676      CALL advtrac( pbaru,pbarv,
    7777     *             p,  masse,q,iapptrac, teta,

  • LMDZ4/trunk/libf/dyn3d/calfis.F

    r524 r616  
    2323     $                  pdq,
    2424     $                  pw,
    25 #ifdef INCA_CH4
     25#ifdef INCA
    2626     $                  flxw,
    2727#endif
    … …  
    163163      REAL unskap, pksurcp
    164164
    165 #ifdef INCA_CH4
     165#ifdef INCA
    166166      REAL flxw(iip1,jjp1,llm)
    167167      REAL flxwfi(ngridmx,llm)
    … …  
    432432
    433433
    434 #ifdef INCA_CH4
     434#ifdef INCA
    435435      CALL gr_dyn_fi(llm,iip1,jjp1,ngridmx,flxw,flxwfi)
    436436#endif
    … …  
    461461     .             zqfi,
    462462     .             pvervel,
    463 #ifdef INCA_CH4
     463#ifdef INCA
    464464     .             flxwfi,
    465465#endif

  • LMDZ4/trunk/libf/dyn3d/iniadvtrac.F

    r543 r616  
    88#ifdef INCA
    99      USE transport_controls, only : hadv_flg, vadv_flg
     10      USE species_names
    1011      USE chemshut
    1112#endif
    … …  
    192193         str1=tnom(iq)
    193194         tname(iiq)=tnom(iq)
    194          ttext(iiq)=str1(1:lnblnk(str1))//descrq(iadv(iiq))
     195         IF (iadv(iiq).eq.0) THEN
     196           ttext(iiq)=str1(1:lnblnk(str1))
     197         ELSE
     198           ttext(iiq)=str1(1:lnblnk(str1))//descrq(iadv(iiq))
     199         endif
    195200         str2=ttext(iiq)
    196201c   schemas tenant compte des moments d'ordre superieur.

  • LMDZ4/trunk/libf/dyn3d/leapfrog.F

    r566 r616  
    132132
    133133      REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm)
    134 #ifdef INCA_CH4
     134#ifdef INCA
    135135      REAL :: flxw(ip1jmp1,llm)
    136136#endif
    … …  
    274274
    275275c
    276 #ifdef INCA_CH4
     276#ifdef INCA
    277277             CALL caladvtrac(q,pbaru,pbarv,
    278278     *                      p, masse, dq,  teta,
    … …  
    358358     $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    359359     $               du,dv,dteta,dq,w,
    360 #ifdef INCA_CH4
     360#ifdef INCA
    361361     $               flxw,
    362362#endif

  • LMDZ4/trunk/libf/phylmd/ini_histrac.h

    r549 r616  
    8787     .                "ave(X)", zsto,zout)
    8888
    89 #ifdef INCA_CH4
     89#ifdef INCA
    9090#ifdef INCAINFO
    9191         DO it=1, phtcnt
    … …  
    244244
    245245#ifdef INCA
    246 #ifdef INCA_CH4
     246!#ifdef INCA_CH4
    247247         CALL histdef(nid_tra, "O3_column", "O3_column",
    248248     .           "DU", iim,jjm+1,nhori, 1,1,1, -99, 32,
    … …  
    275275!    .                "t_sum(X)", zsto,zout)
    276276
    277 #endif
     277!#ENDIF
     278
     279!!! Ajout Anne
     280!#ifdef INCA_AER
     281!! for sulfur cycle
     282!         CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",
     283!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     284!     .                "ave(X)", zsto,zout)
     285!         CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",
     286!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     287!     .                "ave(X)", zsto,zout)
     288!         CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",
     289!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     290!     .                "ave(X)", zsto,zout)
     291!         CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",
     292!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     293!     .                "ave(X)", zsto,zout)
     294!         CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",
     295!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     296!     .                "ave(X)", zsto,zout)
     297!         CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",
     298!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     299!     .                "ave(X)", zsto,zout)
     300!         CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",
     301!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     302!     .                "ave(X)", zsto,zout)
     303!         CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",
     304!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     305!     .                "ave(X)", zsto,zout)
     306!         CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",
     307!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     308!     .                "ave(X)", zsto,zout)
     309!
     310!         CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",
     311!     .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     312!     .           "ave(X)", zsto,zout)
     313!         CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",
     314!     .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     315!     .           "ave(X)", zsto,zout)
     316!         CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",
     317!     .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     318!     .           "ave(X)", zsto,zout)
     319!         CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",
     320!     .           "kg/m2/s", iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     321!     .           "ave(X)", zsto,zout)
     322!
     323!         CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",
     324!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     325!     .                "ave(X)", zsto,zout)
     326!#endif
     327!
     328!#ifdef INCA_NMHC
     329!         CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
     330!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     331!     .                "ave(X)", zsto,zout)
     332!         CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",
     333!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     334!     .                "ave(X)", zsto,zout)
     335!         CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",
     336!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     337!     .                "ave(X)", zsto,zout)
     338!         CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",
     339!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     340!     .                "ave(X)", zsto,zout)
     341!         CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",
     342!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     343!     .                "ave(X)", zsto,zout)
     344!         CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",
     345!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     346!     .                "ave(X)", zsto,zout)
     347!         CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",
     348!     .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
     349!     .                "ave(X)", zsto,zout)
     350!#endif
     351!! Fin ajout Anne
    278352#else
    279353         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",

  • LMDZ4/trunk/libf/phylmd/physiq.F

    r594 r616  
    88     .            u,v,t,qx,
    99     .            omega,
    10 #ifdef INCA_CH4
     10#ifdef INCA
    1111     .            flxmass_w,
    1212#endif
    … …  
    1717#ifdef INCA
    1818      USE chemshut
     19      USE species_names
    1920#ifdef INCA_CH4
    2021!      USE obs_pos
    … …  
    8283#include "clesphys.h"
    8384#include "control.h"
     85#include "logic.h"
    8486#include "temps.h"
    8587#include "comgeomphy.h"
    … …  
    202204      REAL omega(klon,klev)
    203205
    204 #ifdef INCA_CH4
     206#ifdef INCA
    205207      REAL flxmass_w(klon,klev)
    206208#endif
    … …  
    10101012      solswad(:)=0.
    10111013!rv
     1014! anne
     1015      d_u_con(:,:) = 0.0
     1016      d_v_con(:,:) = 0.0
     1017      rnebcon0(:,:) = 0.0
     1018      clwcon0(:,:) = 0.0
     1019      rnebcon(:,:) = 0.0
     1020      clwcon(:,:) = 0.0
     1021      paire_ter(:) = 0.0
     1022      nhistoW(:,:,:,:) = 0.0
     1023      histoW(:,:,:,:) = 0.0
     1024! fin anne
    10121025
    10131026cym
    … …  
    12961309     $                   presnivs,
    12971310     $                   calday,
    1298      $                   tracnam,
    1299      $                   natsnam,
    1300 c     $                   mxoutflds,
    1301 c     $                   outinst,
    1302 c     $                   outtimav,
    13031311     $                   klon,
    13041312     $                   nqmax,
    … …  
    24222430     $                          iim,
    24232431     $                          jjm,
    2424      $                          tr_seri)
     2432#ifdef INCA_AER
     2433     $                          tr_seri,
     2434     $                          ftsol,
     2435     $                          paprs,
     2436     $                          cdragh,
     2437     $                          cdragm,
     2438     $                          pctsrf,
     2439     $                          pdtphys,
     2440     $                          itap)
     2441#else
     2442     $                          tr_seri)     
     2443#endif       
     2444
    24252445#ifdef INCAINFO
    24262446           WRITE(lunout,*)'OK.'
    … …  
    26602680     I                   ibas_con,
    26612681     I                   pmflxr,pmflxs,prfl,psfl,
    2662 #ifdef INCA_CH4
     2682#ifdef INCA
    26632683     I                   flxmass_w,
    26642684#endif

  • LMDZ4/trunk/libf/phylmd/phytrac.F

    r593 r616  
    4848     I                    prfl,
    4949     I                    psfl,
    50 #ifdef INCA_CH4
     50#ifdef INCA
    5151     I                    flxmass_w,
    5252#endif
    … …  
    6666#ifdef INCA_AER
    6767      USE AEROSOL_MOD, only : ntr,trmx,trnx
    68       USE AEROSOL_DIAG, only : cla,las,tausum,angst,aload,scon
    69      .                            ,scavcoef_st,scavcoef_cv
    70       USE AEROSOL_PROGNOS, ONLY : md
     68      USE AEROSOL_DIAG,only : cla,las,tausum,angst,aload,cload,totaerh2o,tau,
     69     $  emiss20,sconc,scavcoef_st,scavcoef_cv
     70     $  ,cload05ss  ,cload05bc  ,cload05pom  ,cload05dust  ,cload05so4
     71     $  ,cload125ss  ,cload125bc  ,cload125pom  ,cload125dust  ,cload125so4
     72      USE AEROSOL_PROGNOS, ONLY : md,mdw
     73      USE AEROSOL_METEO, only : airm
     74#endif
     75#ifdef INCA_NMHC
     76      USE RESISTANCE_DIAGNOSE, ONLY : surf_alb, sol_irrad, surf_temp, surf_wind,
     77     $                                aero_resist, lamin_resist, surf_resist
    7178#endif
    7279#endif
    … …  
    136143      REAL prfl(klon,klev+1),   psfl(klon,klev+1)     !--lessivage large-scale
    137144
    138 #ifdef INCA_CH4
     145#ifdef INCA
    139146      REAL flxmass_w(klon,klev)
    140147#endif
    … …  
    315322
    316323#ifdef INCA
     324      INTEGER           :: lastgas
    317325      INTEGER           :: ncsec
     326
     327      INTEGER           :: prt_flag_ts(nbtr)=(/1,1,1
    318328#ifdef INCA_CH4
    319 #ifdef  INCA_AER
    320       INTEGER           :: prt_flag_ts(51)=(/1,1,1
    321 #else
    322       INTEGER           :: prt_flag_ts(43)=(/1,1,1
    323 #endif
    324 #else
    325 #ifdef INCA_AER
    326       INTEGER           :: prt_flag_ts(11)=(/1,1,1
    327 #endif
    328 #endif
    329 
    330 #ifdef INCA_CH4
    331      .                                              ,1,0,1,1,0,1,0,
    332      .                                         0,0,0,0,0,0,0,1,0,0,
     329     .                                              ,0,0,1,1,1,1,1,
     330     .                                         0,1,0,0,0,0,0,1,0,0,
    333331     .                                         0,1,1,1,1,0,1,1,1,0,
    334332     .                                         1,1,1,1,1,1,1,1,1,1,
    335333     .                                         1,0,0
    336 #endif
    337334#ifdef INCA_AER
    338      .                                              ,1,1,1,1,1,1,1,1
     335     .                                         ,1,1,1,1,0,1,1,1
     336#endif
     337#endif
     338#ifdef INCA_AER
     339c aerosol tracers
     340     .                                        ,1,0,1,1,1,1,1,1,0,1,
     341     .                                         0,1,1,1,1,1,0,1,0,1,1,1
     342#endif
     343#ifdef INCA_NMHC
     344     .                                                 , 1, 1, 1, 1, 1, 1, 1,
     345     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     346     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     347     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     348     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     349     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     350     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     351     .                                          1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
     352     .                                          1, 1, 1, 1, 1, 1, 1
    339353#endif
    340354     .                                         /)
    … …  
    444458
    445459#ifdef INCA_AER
     460
     461! Changement Anne 01/04/2005
    446462        CALL aerosolmain (tr_seri,
    447463     $                 pdtphys,
    448464     $                 pplay,
     465     $                 pdel,
    449466     $                 prfl,
    450467     $                 pmflxr,
    … …  
    455472     $                 ibas_con,
    456473     $                 pphi,
    457      $                 airephy,
    458      $                 nstep)
     474     $                 airephy, ! paire,
     475     $                 nstep,
     476     $                 rneb,         ! for chimiaq
     477     $                 t_seri,       ! for chimiaq
     478     $                 rh)
     479! fin changement anne
     480
    459481#endif
    460482
    … …  
    470492     $                 pplay,      !pmid
    471493     $                 pdel,       !pdel
     494     $                 airephy,
    472495     $                 pctsrf(1,1),!oro
    473496     $                 ftsol,      !tsurf

  • LMDZ4/trunk/libf/phylmd/write_histrac.h

    r541 r616  
    5252     .                                   iim*(jjm+1)*klev,ndex3d)
    5353
    54 #ifdef INCA_CH4
     54! Ajout Anne
     55!#ifdef INCA_AER
     56!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,airm, zx_tmp_3d)
     57!      CALL histwrite(nid_tra,"airm",itra,zx_tmp_3d,
     58!     .                                   iim*(jjm+1)*klev,ndex)
     59!#endif
     60
     61! Fin ajout Anne
     62
     63#ifdef INCA
    5564#ifdef INCAINFO
    5665      DO it=1, phtcnt
    … …  
    130139       CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_st(1,1,it),
    131140     .                  zx_tmp_3d)
    132        CALL histwrite(nid_tra,"scavcoef_st",itau_w,zx_tmp_3d,
     141       CALL histwrite(nid_tra2,"scavcoef_st",itau_w,zx_tmp_3d,
    133142     .                  iim*(jjm+1)*klev,ndex3d)
    134143       CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_cv(1,1,it),
    135144     .                  zx_tmp_3d)
    136        CALL histwrite(nid_tra,"scavcoef_cv",itau_w,zx_tmp_3d,
     145       CALL histwrite(nid_tra2,"scavcoef_cv",itau_w,zx_tmp_3d,
    137146     .                  iim*(jjm+1)*klev,ndex3d)
    138147
    … …  
    215224#endif
    216225       else
    217 #ifdef INCA_CH4
     226#ifdef INCA
    218227       if ( solsym(it) == 'OX' ) then
    219228       zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(id_o3)
    … …  
    251260
    252261#ifdef INCA
    253 #ifdef INCA_CH4
    254       CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_tr_col(1), zx_tmp_2d)
     262!#ifdef INCA_CH4
     263      CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_col(1), zx_tmp_2d)
    255264      CALL histwrite(nid_tra,"O3_column",itau_w,zx_tmp_2d,
    256265     .     iim*(jjm+1),ndex2d)
    257266
    258       CALL gr_fi_ecrit(1, klon,iim,jjm+1, co_tr_col(1), zx_tmp_2d)
     267      CALL gr_fi_ecrit(1, klon,iim,jjm+1, co_col(1), zx_tmp_2d)
    259268      CALL histwrite(nid_tra,"CO_column",itau_w,zx_tmp_2d,
    260269     .     iim*(jjm+1),ndex2d)
    261270
    262       CALL gr_fi_ecrit(1, klon,iim,jjm+1, ch4_tr_col(1), zx_tmp_2d)
     271      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ch4_col(1), zx_tmp_2d)
    263272      CALL histwrite(nid_tra,"CH4_column",itau_w,zx_tmp_2d,
    264273     .     iim*(jjm+1),ndex2d)
    265274
    266       CALL gr_fi_ecrit(1, klon,iim,jjm+1, no2_tr_col(1), zx_tmp_2d)
     275      CALL gr_fi_ecrit(1, klon,iim,jjm+1, no2_col(1), zx_tmp_2d)
    267276      CALL histwrite(nid_tra,"NO2_column",itau_w,zx_tmp_2d,
    268277     .     iim*(jjm+1),ndex2d)
    … …  
    281290      CALL histwrite(nid_tra,"O3_loss",itau_w,zx_tmp_3d,
    282291     .                                   iim*(jjm+1)*klev,ndex3d)
     292!! Ajout Anne
     293!#ifdef INCA_AER
     294!! for sulfur cycle
     295!
     296!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsoh(1,1),
     297!     .     zx_tmp_3d)
     298!      CALL histwrite(nid_tra,"SO2_p_dmsoh",itra,zx_tmp_3d,
     299!     .                                   iim*(jjm+1)*klev,ndex)
     300!
     301!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsno3(1,1),
     302!     .     zx_tmp_3d)
     303!      CALL histwrite(nid_tra,"SO2_p_dmsno3",itra,zx_tmp_3d,
     304!     .                                   iim*(jjm+1)*klev,ndex)
     305!
     306!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_h2soh(1,1),
     307!     .     zx_tmp_3d)
     308!      CALL histwrite(nid_tra,"SO2_p_h2soh",itra,zx_tmp_3d,
     309!     .                                   iim*(jjm+1)*klev,ndex)
     310!
     311!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsooh(1,1),
     312!     .     zx_tmp_3d)
     313!      CALL histwrite(nid_tra,"SO2_p_dmsooh",itra,zx_tmp_3d,
     314!     .                                   iim*(jjm+1)*klev,ndex)
     315!
     316!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,DMSO_p_dmsoh(1,1),
     317!     .     zx_tmp_3d)
     318!      CALL histwrite(nid_tra,"DMSO_p_dmsoh",itra,zx_tmp_3d,
     319!     .                                   iim*(jjm+1)*klev,ndex)
     320!
     321!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASMSAM_p_dmsooh(1,1),
     322!     .     zx_tmp_3d)
     323!      CALL histwrite(nid_tra,"ASMSAM_p_dmsooh",itra,zx_tmp_3d,
     324!     .                                   iim*(jjm+1)*klev,ndex)
     325!
     326!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2oh(1,1),
     327!     .     zx_tmp_3d)
     328!      CALL histwrite(nid_tra,"ASSO4M_p_so2oh",itra,zx_tmp_3d,
     329!     .                                   iim*(jjm+1)*klev,ndex)
     330!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2h2o2(1,1),
     331!     .     zx_tmp_3d)
     332!      CALL histwrite(nid_tra,"ASSO4M_p_so2h2o2",itra,zx_tmp_3d,
     333!     .                                   iim*(jjm+1)*klev,ndex)
     334!
     335!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2o3(1,1),
     336!     .     zx_tmp_3d)
     337!      CALL histwrite(nid_tra,"ASSO4M_p_so2o3",itra,zx_tmp_3d,
     338!     .                                   iim*(jjm+1)*klev,ndex)
     339!
     340!c closing the sulfur budget
     341!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_so2(1,1), zx_tmp_3d)
     342!      CALL histwrite(nid_tra,"Wet3D_SO2",itra,zx_tmp_3d,
     343!     .     iim*(jjm+1)*klev,ndex)
     344!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_dms(1,1), zx_tmp_3d)
     345!      CALL histwrite(nid_tra,"Wet3D_DMS",itra,zx_tmp_3d,
     346!     .     iim*(jjm+1)*klev,ndex)
     347!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_hno3(1,1), zx_tmp_3d)
     348!      CALL histwrite(nid_tra,"Wet3D_HNO3",itra,zx_tmp_3d,
     349!     .     iim*(jjm+1)*klev,ndex)
     350!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_h2o2(1,1), zx_tmp_3d)
     351!      CALL histwrite(nid_tra,"Wet3D_H2O2",itra,zx_tmp_3d,
     352!     .     iim*(jjm+1)*klev,ndex)
     353!
     354!
     355!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,PH_HIST(1,1),
     356!     .     zx_tmp_3d)
     357!      CALL histwrite(nid_tra,"PH_HIST",itra,zx_tmp_3d,
     358!     .                                   iim*(jjm+1)*klev,ndex)
     359!#endif
     360!
     361!#ifdef INCA_NMHC
     362!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1),
     363!     .     zx_tmp_3d)
     364!      CALL histwrite(nid_tra,"CO2_basprod",itra,zx_tmp_3d,
     365!     .                                   iim*(jjm+1)*klev,ndex)
     366!
     367!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_nmhcprod(1,1),
     368!     .     zx_tmp_3d)
     369!      CALL histwrite(nid_tra,"CO2_nmhcprod",itra,zx_tmp_3d,
     370!     .                                   iim*(jjm+1)*klev,ndex)
     371!
     372!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_radicalprod(1,1),
     373!     .     zx_tmp_3d)
     374!      CALL histwrite(nid_tra,"CO2_radicalprod",itra,zx_tmp_3d,
     375!     .                                   iim*(jjm+1)*klev,ndex)
     376!
     377!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_prod(1,1),
     378!     .     zx_tmp_3d)
     379!      CALL histwrite(nid_tra,"HNO3_prod",itra,zx_tmp_3d,
     380!     .                                   iim*(jjm+1)*klev,ndex)
     381!
     382!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_loss(1,1),
     383!     .     zx_tmp_3d)
     384!      CALL histwrite(nid_tra,"HNO3_loss",itra,zx_tmp_3d,
     385!     .                                   iim*(jjm+1)*klev,ndex)
     386!
     387!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_prod(1,1),
     388!     .     zx_tmp_3d)
     389!      CALL histwrite(nid_tra,"CO_prod",itra,zx_tmp_3d,
     390!     .                                   iim*(jjm+1)*klev,ndex)
     391!
     392!      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_loss(1,1),
     393!     .     zx_tmp_3d)
     394!      CALL histwrite(nid_tra,"CO_loss",itra,zx_tmp_3d,
     395!     .                                   iim*(jjm+1)*klev,ndex)
     396!
     397!#endif
     398!
     399!! Fin ajout Anne
    283400
    284401!     ... Special section for daytime averaging
    … …  
    292409!    .                                  iim*(jjm+1)*klev,ndex3d)
    293410
    294 #endif
     411!#endif
    295412#else
    296413
Note: See TracChangeset for help on using the changeset viewer.