Changeset 624 for LMDZ4/trunk/libf


Ignore:
Timestamp:
Apr 29, 2005, 9:58:12 AM (20 years ago)
Author:
lmdzadmin
Message:

Rajout de sorties pour INCA -Anne Cozic
MAF

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

Legend:

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

    r616 r624  
    277277!#ENDIF
    278278
    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
     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
    352352#else
    353353         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
  • LMDZ4/trunk/libf/phylmd/write_histrac.h

    r619 r624  
    5454! Ajout Anne
    5555!#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)
     56      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,airm, zx_tmp_3d)
     57      CALL histwrite(nid_tra,"airm",itau_w,zx_tmp_3d,
     58     .                                   iim*(jjm+1)*klev,ndex)
    5959!#endif
    6060
     
    290290      CALL histwrite(nid_tra,"O3_loss",itau_w,zx_tmp_3d,
    291291     .                                   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
     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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,zx_tmp_3d,
     338     .                                   iim*(jjm+1)*klev,ndex)
     339
     340c 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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,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",itau_w,zx_tmp_3d,
     395     .                                   iim*(jjm+1)*klev,ndex)
     396
     397#endif
     398
     399! Fin ajout Anne
    400400
    401401!     ... Special section for daytime averaging
Note: See TracChangeset for help on using the changeset viewer.