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

Rajout de sorties pour INCA -Anne Cozic
MAF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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.