Ignore:
Timestamp:
Nov 6, 2006, 4:51:16 PM (18 years ago)
Author:
Laurent Fairhead
Message:

Modifications pour rendre INCA plus independant de LMDZ ACo
LF

Location:
LMDZ4/branches/V3_test/libf
Files:
17 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/V3_test/libf/dyn3d/advtrac.F

    r703 r726  
    88     *                   p,  masse,q,iapptrac,teta,
    99     *                  flxw,
    10      *                  pk,
    11      *                  mmt_adj,
    12      *                  hadv_flg)
     10     *                  pk)
    1311#else
    1412      SUBROUTINE advtrac(pbaru,pbarv ,
    … …  
    5149      REAL pk(ip1jmp1,llm)
    5250#ifdef INCA
    53       INTEGER            :: hadv_flg(nqmx)
    54       REAL               :: mmt_adj(ip1jmp1,llm,1)
    5551      REAL               :: flxw(ip1jmp1,llm)
    5652#endif
    … …  
    215211#ifdef INCA
    216212       do iiq = iq+1, iq+3
    217          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     213c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     214          q(:,:,iiq)=q(:,:,iiq)*1
    218215       enddo
    219216#endif
    … …  
    233230#ifdef INCA
    234231       do iiq = iq+1, iq+9
    235          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     232c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
     233         q(:,:,iiq)=q(:,:,iiq)*1
    236234       enddo
    237235#endif

  • LMDZ4/branches/V3_test/libf/dyn3d/advtrac.h

    r524 r726  
    55c INCLUDE 'advtrac.h'
    66
    7       COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv
     7      COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv,
     8     &     nbtrac, nprath, mmt_adj, hadv_flg, vadv_flg, conv_flg,
     9     &     pbl_flg, tracnam
    810      INTEGER iadv(nqmx) ! indice schema de transport
    911      INTEGER hadv(nqmx) ! indice schema transport horizontal
    … …  
    1315      character*10 tname(nqmx) ! nom du traceur pour restart
    1416      character*13 ttext(nqmx) ! nom long du traceur pour sorties
     17
     18
     19      integer nbtrac
     20      integer nprath
     21      real    mmt_adj(iim+1,jjm+1,llm, 1)
     22      integer hadv_flg(nqmx)
     23      integer vadv_flg(nqmx)
     24      integer conv_flg(nqmx-2)
     25      integer pbl_flg(nqmx-2)
     26      character*8 tracnam(nqmx-2)
    1527c-----------------------------------------------------------------------

  • LMDZ4/branches/V3_test/libf/dyn3d/caladvtrac.F

    r703 r726  
    88     *                   p ,masse, dq ,  teta,
    99     *                   flxw,
    10      *                   pk,
    11      *                   mmt_adj,
    12      *                   hadv_flg)
     10     *                   pk)
    1311#else
    1412            SUBROUTINE caladvtrac(q,pbaru,pbarv ,
    … …  
    4240      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    4341#ifdef INCA
    44       INTEGER            :: hadv_flg(nqmx)
    45       REAL               :: mmt_adj(iip1,jjp1,llm,1)
    4642      REAL               :: flxw(ip1jmp1,llm)
    4743#endif
    … …  
    7773     *             p,  masse,q,iapptrac, teta,
    7874     .             flxw,
    79      .             pk,
    80      .             mmt_adj,
    81      .             hadv_flg)
     75     .             pk)
    8276#else
    8377      CALL advtrac( pbaru,pbarv,

  • LMDZ4/branches/V3_test/libf/dyn3d/gcm.F

    r704 r726  
    5858#include "iniprint.h"
    5959#include "tracstoke.h"
    60 
     60#include "advtrac.h"
    6161
    6262      INTEGER         longcles
    … …  
    145145      dynhistave_file = 'dyn_hist_ave.nc'
    146146
     147c initialisation Anne
     148      hadv_flg(:) = 0.
     149      vadv_flg(:) = 0.
     150      conv_flg(:) = 0.
     151      pbl_flg(:)  = 0.
     152      tracnam(:)  = '        '
     153      nprath = 1
     154      nbtrac = 0
     155      mmt_adj(:,:,:,:) = 1
     156
     157
    147158c--------------------------------------------------------------------------
    148159c   Iflag_phys controle l'appel a la physique :
    … …  
    188199      call init_phys_openmp
    189200      call InitComgeomphy
     201
     202#ifdef INCA
     203      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
     204      call init_inca_para(iim,jjm+1,klon2,phy_size,klon_para_nb)
     205#endif
     206
    190207c
    191208c
    … …  
    222239      endif
    223240
     241#ifdef INCA
     242      call init_inca_dim(klon,llm,iim,jjm,
     243     $     rlonu,rlatu,rlonv,rlatv)
     244#endif
    224245
    225246

  • LMDZ4/branches/V3_test/libf/dyn3d/guide.F

    r703 r726  
    367367                if (first.and.ini_anal) vcov(ij,l)=a
    368368            enddo
    369             if (first.and.ini_anal) vcov(ij,l)=a
    370369         enddo
    371370      endif

  • LMDZ4/branches/V3_test/libf/dyn3d/iniadvtrac.F

    r703 r726  
    66      subroutine iniadvtrac(nq)
    77      USE ioipsl
    8 #ifdef INCA
    9       USE transport_controls, only : hadv_flg, vadv_flg
    10       USE species_names
    11       USE chemshut
    12 #endif
    138      IMPLICIT NONE
    149c=======================================================================
    … …  
    5954      descrq(20)='SLP'
    6055      descrq(30)='PRA'
     56
     57#ifdef INCA
     58
     59      CALL init_transport(
     60     $     hadv_flg,
     61     $     vadv_flg,
     62     $     conv_flg,
     63     $     pbl_flg,
     64     $     tracnam)
     65#endif
    6166
    6267c-----------------------------------------------------------------------

  • LMDZ4/branches/V3_test/libf/dyn3d/leapfrog.F

    r703 r726  
    77     &                    time_0)
    88
    9 #ifdef INCA
    10       USE transport_controls, ONLY : hadv_flg, mmt_adj
    11 #endif
    129
    1310cIM : pour sortir les param. du modele dans un fis. netcdf 110106
    … …  
    6158#include "com_io_dyn.h"
    6259#include "iniprint.h"
    63 
     60#include "advtrac.h"
    6461c#include "tracstoke.h"
    6562
    … …  
    297294     *                      p, masse, dq,  teta,
    298295     .             flxw,
    299      .             pk,
    300      .             mmt_adj,
    301      .             hadv_flg)
     296     .             pk)
    302297#else
    303298             CALL caladvtrac(q,pbaru,pbarv,

  • LMDZ4/branches/V3_test/libf/dyn3dpar/advtrac.h

    r709 r726  
    55c INCLUDE 'advtrac.h'
    66
    7       COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv
     7      COMMON/advtr/iadv,hadv,vadv,tnom,tname,ttext,niadv, 
     8     &     nbtrac, nprath, mmt_adj, hadv_flg, vadv_flg, conv_flg,
     9     &     pbl_flg, tracnam
    810      INTEGER iadv(nqmx) ! indice schema de transport
    911      INTEGER hadv(nqmx) ! indice schema transport horizontal
    … …  
    1315      character*10 tname(nqmx) ! nom du traceur pour restart
    1416      character*13 ttext(nqmx) ! nom long du traceur pour sorties
     17
     18      integer nbtrac
     19      integer nprath
     20      real    mmt_adj(iim+1,jjm+1,llm, 1)
     21      integer hadv_flg(nqmx)
     22      integer vadv_flg(nqmx)
     23      integer conv_flg(nqmx-2)
     24      integer pbl_flg(nqmx-2)
     25      character*8 tracnam(nqmx-2)
    1526c-----------------------------------------------------------------------

  • LMDZ4/branches/V3_test/libf/dyn3dpar/caladvtrac_p.F

    r709 r726  
    99     *                   flxw,
    1010     *                   pk,
    11      *                   mmt_adj,
    12      *                   hadv_flg,iapptrac)
     11     *                   iapptrac)
    1312#else
    1413            SUBROUTINE caladvtrac_p(q,pbaru,pbarv ,
    … …  
    4140      REAL p( ip1jmp1,llmp1),q( ip1jmp1,llm,nqmx),dq( ip1jmp1,llm,2 )
    4241      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    43 #ifdef INCA_CH4
    44 cym      INTEGER            :: hadv_flg(nq)
    45       INTEGER            :: hadv_flg(nqmx)
    46       REAL               :: mmt_adj(iip1,jjp1,llm)
     42#ifdef INCA
    4743      REAL               :: flxw(ip1jmp1,llm)
    4844#endif
    … …  
    7874c      print *,'appel a advtrac'
    7975
    80 #ifdef INCA_CH4
     76#ifdef INCA
    8177      CALL advtrac_p( pbaru,pbarv,
    8278     *             p,  masse,q,iapptrac, teta,
    8379     .             flxw,
    84      .             pk,
    85      .             mmt_adj,
    86      .             hadv_flg)
     80     .             pk)
    8781#else
    8882      CALL advtrac_p( pbaru,pbarv,

  • LMDZ4/branches/V3_test/libf/dyn3dpar/gcm.F

    r709 r726  
    1414      USE mod_hallo
    1515      USE Bands
    16 #ifdef INCA
    17       USE inca_dim
    18 #endif
    1916      IMPLICIT NONE
    2017
    … …  
    6461#include "iniprint.h"
    6562#include "tracstoke.h"
    66 
     63#include "advtrac.h"
    6764
    6865      INTEGER         longcles
    … …  
    157154      dynhist_file = 'dyn_hist'
    158155      dynhistave_file = 'dyn_hist_ave'
     156
     157
     158c initialisation Anne
     159      hadv_flg(:) = 0.
     160      vadv_flg(:) = 0.
     161      conv_flg(:) = 0.
     162      pbl_flg(:)  = 0.
     163      tracnam(:)  = '        '
     164      nprath = 1
     165      nbtrac = 0
     166      mmt_adj(:,:,:,:) = 1
     167
    159168
    160169c--------------------------------------------------------------------------
    … …  
    221230c$OMP END PARALLEL
    222231#ifdef INCA
    223       call init_inca_dim
     232      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
     233      call init_inca_para(iim,jjm+1,klon2,phy_size,klon_para_nb)
    224234#endif
    225235
    … …  
    250260      endif
    251261
     262#ifdef INCA
     263      call init_inca_dim(klon,llm,iim,jjm,
     264     $     rlonu,rlatu,rlonv,rlatv)
     265#endif
    252266
    253267

  • LMDZ4/branches/V3_test/libf/dyn3dpar/iniadvtrac.F

    r709 r726  
    66      subroutine iniadvtrac(nq)
    77      USE ioipsl
    8 #ifdef INCA
    9       USE transport_controls, only : hadv_flg, vadv_flg
    10 cym      USE chemshut
    11       USE species_names
    12 #endif
    138      IMPLICIT NONE
    149c=======================================================================
    … …  
    6055      descrq(30)='PRA'
    6156
     57#ifdef INCA
     58
     59      CALL init_transport(
     60     $     hadv_flg,
     61     $     vadv_flg,
     62     $     conv_flg,
     63     $     pbl_flg,
     64     $     tracnam)
     65#endif
     66
    6267c-----------------------------------------------------------------------
    6368c        Choix  des schemas d'advection pour l'eau et les traceurs
    … …  
    110115      tnom(2)='H2Ol'
    111116      nq=nbtrac+2
     117     
    112118       if (nq.gt.nqmx) then
    113        print*,'nombre de traceurs incompatible INCA/LMDZT'
     119       print*,'nombre de traceurs incompatible INCA/LMDZT', nq, nbtrac
    114120       stop
    115121       endif

  • LMDZ4/branches/V3_test/libf/dyn3dpar/leapfrog_p.F

    r709 r726  
    1919       USE Write_Field_p
    2020       USE vampir
    21        
    22 #ifdef INCA
    23       USE transport_controls, ONLY : hadv_flg, mmt_adj
    24 #endif
    2521
    2622      IMPLICIT NONE
    … …  
    7672#include "academic.h"
    7773#include "clesphys.h"
    78 
     74#include "advtrac.h"
    7975     
    8076      include 'mpif.h'
    … …  
    148144
    149145      REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm)
    150 #ifdef INCA_CH4
     146#ifdef INCA
    151147      REAL :: flxw(ip1jmp1,llm)
    152148#endif
    … …  
    528524c$OMP PARALLEL DEFAULT(SHARED)
    529525c
    530 #ifdef INCA_CH4
     526#ifdef INCA
    531527             CALL caladvtrac_p(q,pbaru,pbarv,
    532528     *                      p, masse, dq,  teta,
    533529     .             flxw,
    534      .             pk,
    535      .             mmt_adj,
    536      .             hadv_flg,iapptrac)
     530     .             pk,
     531     .             iapptrac)
    537532#else
    538533             CALL caladvtrac_p(q,pbaru,pbarv,
    … …  
    680675     *                               jj_Nb_physic,2,2,Request_physic)
    681676        enddo
    682 #ifdef INCA_CH4
     677#ifdef INCA
    683678        call Register_SwapFieldHallo(flxw,flxw,ip1jmp1,llm,
    684679     *                               jj_Nb_physic,2,2,Request_physic)
    … …  
    698693     $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    699694     $               du,dv,dteta,dq,w,
    700 #ifdef INCA_CH4
     695#ifdef INCA
    701696     $               flxw,
    702697#endif

  • LMDZ4/branches/V3_test/libf/dyn3dpar/mod_hallo.F90

    r718 r726  
    33implicit none
    44!  include 'mpif.h'
    5   integer, parameter :: MaxRequest=80
     5  integer, parameter :: MaxRequest=200
    66  integer, parameter :: MaxProc=80
    77  integer, parameter :: MaxBufferSize=1024*1024*16
    … …  
    430430      type(request_SR),pointer :: Req
    431431      type(Hallo),pointer :: PtrHallo
    432       integer, dimension(2*mpi_size) :: TabRequest
    433       integer, dimension(MPI_STATUS_SIZE,2*mpi_size) :: TabStatus
     432      integer, dimension(4) :: TabRequest
     433      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    434434      integer :: NbRequest
    435435      integer :: i,rank,pos,ij,l,ierr
    … …  
    512512      type(request_SR),pointer :: Req
    513513      type(Hallo),pointer :: PtrHallo
    514       integer, dimension(mpi_size) :: TabRequest
    515       integer, dimension(MPI_STATUS_SIZE,mpi_size) :: TabStatus
     514      integer, dimension(4) :: TabRequest
     515      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    516516      integer :: NbRequest
    517517      integer :: i,rank,pos,ij,l,ierr
    … …  
    553553      type(request_SR),pointer :: Req
    554554      type(Hallo),pointer :: PtrHallo
    555       integer, dimension(mpi_size) :: TabRequest
    556       integer, dimension(MPI_STATUS_SIZE,mpi_size) :: TabStatus
     555      integer, dimension(4) :: TabRequest
     556      integer, dimension(MPI_STATUS_SIZE,4) :: TabStatus
    557557      integer :: NbRequest
    558558      integer :: i,rank,pos,ij,l,ierr

  • LMDZ4/branches/V3_test/libf/phylmd/ini_histrac.h

    r704 r726  
    22! $Header$
    33!
     4#ifndef INCA
    45         CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
    56c
    … …  
    1819     .                 klev, presnivs, nvert)
    1920
    20 #ifdef INCA_AER
    21 cym         CALL histbeg("histrac_aer", iim,zx_lon, jjm+1,zx_lat,
    22 cym     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
    23 cym     .                 nhori, nid_tra2)
    24          CALL histbeg_phy("histrac_aer", itau_phy, zjulian, pdtphys,
    25      .                 nhori, nid_tra2)
    26 
    27 cym         CALL histbeg("histrac_inst", iim,zx_lon, jjm+1,zx_lat,
    28 cym     .                 1,iim,1,jjm+1,  itau_phy, zjulian, pdtphys,
    29 cym     .                 nhori, nid_tra3)
    30          CALL histbeg_phy("histrac_inst",itau_phy, zjulian, pdtphys,
    31      .                 nhori, nid_tra3)
    32 
    33          call histvert(nid_tra2, "presnivs", "Vertical levels", "mb",
    34      .                 klev, presnivs, nvert)
    35          call histvert(nid_tra3, "presnivs", "presnivs", "mb",
    36      .                 klev, presnivs, nvert)
    37 #endif
    3821
    3922
    … …  
    4831     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    4932     .                "once",  zsto,zout)
    50 #ifdef INCA
    51          CALL histdef(nid_tra, "ps", "Surface pressure", "Pa",
    52      .                iim,jjphy_nb,nhori, 1,1,1,-99, 32,
    53      .                "ave(X)", zsto,zout)
    54 
    55          CALL histdef(nid_tra, "ptrop", "Tropopause pressure", "Pa",
    56      .                iim,jjphy_nb,nhori, 1,1,1,-99, 32,
    57      .                "ave(X)", zsto,zout)
    58 
    59 C   3d FIELDS
    60          CALL histdef(nid_tra, "temp", "Air temperature", "K",
    61      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    62      .                "ave(X)", zsto,zout)
    63 
    64          CALL histdef(nid_tra, "u", "zonal wind component", "m/s",
    65      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    66      .                "ave(X)", zsto,zout)
    67 
    68          CALL histdef(nid_tra, "v", "zonal wind component", "m/s",
    69      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    70      .                "ave(X)", zsto,zout)
    71 
    72          CALL histdef(nid_tra, "h2o", "Specific Humidity", "MMR",
    73      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    74      .                "ave(X)", zsto,zout)
    75 
    76          CALL histdef(nid_tra, "pmid", "Pressure", "Pa",
    77      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    78      .                "ave(X)", zsto,zout)
    79 
    80          CALL histdef(nid_tra, "pdel", "Delta Pressure", "Pa",
    81      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    82      .                "ave(X)", zsto,zout)
    83 
    84 ! MS info
    85 #ifdef INCA_AER
    86          CALL histdef(nid_tra, "airm", "Air mass", "kg",
    87      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    88      .                "ave(X)", zsto,zout)
    89 #endif
    90 #ifdef INCA
    91 #ifdef INCAINFO
    92 
    93 #ifdef INCA_CH4
    94          DO it=1, phtcnt
    95          WRITE(str2,'(i2.2)') it
    96          CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
    97      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    98      .                "ave(X)", zsto,zout)
    99          ENDDO
    100          DO it=1, hetcnt
    101          WRITE(str2,'(i2.2)') it
    102          CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
    103      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    104      .                "ave(X)", zsto,zout)
    105          ENDDO
    106          DO it=1, extcnt
    107          WRITE(str2,'(i2.2)') it
    108          CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
    109      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    110      .                "ave(X)", zsto,zout)
    111          ENDDO
    112 #endif
    113 #ifdef INCA_NMHC
    114          DO it=1, phtcnt
    115          WRITE(str2,'(i2.2)') it
    116          CALL histdef(nid_tra, "j"//str2,"j"//str2, "CM-3 S-1",
    117      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    118      .                "ave(X)", zsto,zout)
    119          ENDDO
    120          DO it=1, hetcnt
    121          WRITE(str2,'(i2.2)') it
    122          CALL histdef(nid_tra, "w"//str2,"w"//str2, "S-1",
    123      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    124      .                "ave(X)", zsto,zout)
    125          ENDDO
    126 
    127          DO it=1, extcnt
    128          WRITE(str2,'(i2.2)') it
    129          CALL histdef(nid_tra, "ext"//str2,"ext"//str2, "CM-3 S-1",
    130      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    131      .                "ave(X)", zsto,zout)
    132          ENDDO
    133 #endif
    134          DO it=1, nfs
    135          WRITE(str2,'(i2.2)') it
    136          CALL histdef(nid_tra, "INV"//str2, "INV"//str2, "CM-3",
    137      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    138      .                "ave(X)", zsto,zout)
    139          ENDDO
    140 
    141 #else
    142 
    143 #ifdef INCA_NMHC
    144          CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
    145      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    146      .                "ave(X)", zsto,zout)
    147          CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
    148      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    149      .                "ave(X)", zsto,zout)
    150          CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
    151      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    152      .                "ave(X)", zsto,zout)
    153          CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
    154      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    155      .                "ave(X)", zsto,zout)
    156          CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
    157      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    158      .                "ave(X)", zsto,zout)
    159          CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
    160      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    161      .                "ave(X)", zsto,zout)
    162 #endif
    163 #ifdef INCA_CH4
    164          CALL histdef(nid_tra, "jO3","jO3", "CM-3 S-1",
    165      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    166      .                "ave(X)", zsto,zout)
    167          CALL histdef(nid_tra, "jNO2","jNO2", "CM-3 S-1",
    168      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    169      .                "ave(X)", zsto,zout)
    170          CALL histdef(nid_tra, "jH2O2","jH2O2", "CM-3 S-1",
    171      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    172      .                "ave(X)", zsto,zout)
    173          CALL histdef(nid_tra, "wHNO3","wHNO3", "S-1",
    174      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    175      .                "ave(X)", zsto,zout)
    176          CALL histdef(nid_tra, "kN2O5", "kN2O5","CM-3 S-1",
    177      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    178      .                "ave(X)", zsto,zout)
    179          CALL histdef(nid_tra, "LghtNO","LghtNO", "CM-3 S-1",
    180      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    181      .                "ave(X)", zsto,zout)
    182 #endif
    183 #endif
    184 
    185 !         DO it=1, grpcnt
    186 !         CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
    187 !     .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    188 !     .                "ave(X)", zsto,zout)
    189 !         ENDDO
    190 #endif
    191 
    192 #ifdef INCA_AER
    193 
    194         CALL histdef(nid_tra2, "scavcoef_st","scavcoef_st", "S-1",
    195      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    196      .                "ave(X)", zsto,zout)
    197         CALL histdef(nid_tra2, "scavcoef_cv","scavcoef_cv", "S-1",
    198      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    199      .                "ave(X)", zsto,zout)
    200         CALL histdef(nid_tra2, "AngstroemComp","AngstroemComp",
    201      .      "angs comp", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    202      .      "ave(X)", zsto,zout)
    203 
    204         CALL histdef(nid_tra2, "TOTAERH2O", "Total aerosol water",
    205      .      "kg/m2", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    206      .      "ave(X)", zsto,zout)
    207 #endif
    208 #endif
    20933         DO it=1,nqmax
    21034C champ 2D
    211 #ifdef INCA
    212          IF ( prt_flag_ts(it) == 0 ) CYCLE
    213          CALL histdef(nid_tra, "Emi_"//solsym(it), "Emi_"//solsym(it),
    214      .           "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    215      .           "ave(X)", zsto,zout)
    216          CALL histdef(nid_tra, "Dep_"//solsym(it), "Dep_"//solsym(it),
    217      .           "cm/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    218      .           "ave(X)", zsto,zout)
    219 #ifdef INCA_AER
    220 
    221          CALL histdef(nid_tra2, "Dry_"//solsym(it), "Dry_"//solsym(it),
    222      .           "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    223      .           "ave(X)", zsto,zout)
    224          
    225         IF  ((it .ge. trmx) .and. (it .le. trnx)) then
    226           CALL histdef(nid_tra2, "Sed_"//solsym(it), "Sed_"//solsym(it),
    227      .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    228      .      "ave(X)", zsto,zout)
    229           CALL histdef(nid_tra2, "Wet_"//solsym(it), "Wet_"//solsym(it),
    230      .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    231      .      "ave(X)", zsto,zout)
    232           CALL histdef(nid_tra2, "WetST_"//solsym(it), "WetST_"//solsym(it),
    233      .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    234      .      "ave(X)", zsto,zout)
    235           CALL histdef(nid_tra2, "WetCV_"//solsym(it), "WetCV_"//solsym(it),
    236      .      "kg/m2/s", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    237      .      "ave(X)", zsto,zout)
    238           CALL histdef(nid_tra2, "Emi_alt_"//solsym(it), "Emi_alt_"//solsym(it),
    239      .      "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    240      .      "ave(X)", zsto,zout)
    241           CALL histdef(nid_tra2, "Load_"//solsym(it), "Load_"//solsym(it),
    242      .      "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    243      .      "ave(X)", zsto,zout)
    244           CALL histdef(nid_tra2, "SConc_"//solsym(it), "SConc_"//solsym(it),
    245      .      "kg/m3", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    246      .      "ave(X)", zsto,zout)
    247           do la=1,las
    248             CALL histdef(nid_tra2, "OD"//cla(la)//"_"//solsym(it), "OD"//cla(la)//"_"//solsym(it),
    249      .        "opt. depth", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    250      .        "ave(X)", zsto,zout)
    251 
    252             CALL histdef(nid_tra2, "OD3D"//cla(la)//"_"//solsym(it), "OD3d_"//cla(la)//"_"//solsym(it),
    253      .        "opt. depth 3D", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    254      .        "ave(X)", zsto,zout)
    255           enddo
    256 
    257           CALL histdef(nid_tra2, "CLOAD_"//solsym(it), "3D LOAD_"//solsym(it),
    258      .      "kg/m2 resp. #/m2", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    259      .      "ave(X)", zsto,zout)
    260 
    261           CALL histdef(nid_tra2, "MD_"//solsym(it), "MD_"//solsym(it),
    262      .      "median diameter [m]", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    263      .      "ave(X)", zsto,zout)
    264 
    265           CALL histdef(nid_tra2, "MDW_"//solsym(it), "MDW_"//solsym(it),
    266      .      "wet median diameter [m]", iim,jjphy_nb,nhori, klev,1,klev, nvert, 32,
    267      .      "ave(X)", zsto,zout)
    268 
    269 
    270 
    271           CALL histdef(nid_tra3, "Inst_Load_"//solsym(it), "Inst_Load_"//solsym(it),
    272      .      "kg/m2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    273      .      "inst(X)", zout,zout)
    274 
    275         endif
    276 
    277 #endif
    278          CALL histdef(nid_tra, solsym(it), solsym(it), "VMR",
    279      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    280      .                "ave(X)", zsto,zout)
    281 #else
    28235         iq=it+2
    28336         iiq=niadv(iq)
    … …  
    30659c---fin Olivia   
    30760
    308 #endif
    30961         ENDDO
    31062
    311 #ifdef INCA
    312          CALL histdef(nid_tra, "O3_column", "O3_column",
    313      .           "DU", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    314      .           "ave(X)", zsto,zout)
    315          CALL histdef(nid_tra, "CO_column", "CO_column",
    316      .           "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    317      .           "ave(X)", zsto,zout)
    318          CALL histdef(nid_tra, "CH4_column", "CH4_column",
    319      .           "10^18 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    320      .           "ave(X)", zsto,zout)
    321          CALL histdef(nid_tra, "NO2_column", "NO2_column",
    322      .           "10^15 CM-2", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    323      .           "ave(X)", zsto,zout)
    324          CALL histdef(nid_tra, "O3_ste", "O3_ste",
    325      .           "CM-2 S-1", iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    326      .           "ave(X)", zsto,zout)
    327          CALL histdef(nid_tra, "O3_prod", "O3_prod", "CM-3 S-1",
    328      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    329      .                "ave(X)", zsto,zout)
    330          CALL histdef(nid_tra, "O3_loss", "O3_loss", "CM-3 S-1",
    331      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    332      .                "ave(X)", zsto,zout)
    333 
    334 !        Special variables for daytime averaging
    335 !        CALL histdef(nid_tra, "day_cnt", "day_cnt", "-",
    336 !    .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    337 !    .                "t_sum(X)", zsto,zout)
    338 !        CALL histdef(nid_tra, "NO_day", "NO_day", "VMR",
    339 !    .                iim,jjm+1,nhori, klev,1,klev,nvert, 32,
    340 !    .                "t_sum(X)", zsto,zout)
    341 
    342 
    343 
    344 !! Ajout Anne
    345 #ifdef INCA_AER
    346 ! for sulfur cycle
    347          CALL histdef(nid_tra, "SO2_p_dmsoh", "SO2_p_dmsoh", "CM-3 S-1",
    348      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    349      .                "ave(X)", zsto,zout)
    350          CALL histdef(nid_tra, "SO2_p_dmsno3", "SO2_p_dmsno3", "CM-3 S-1",
    351      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    352      .                "ave(X)", zsto,zout)
    353          CALL histdef(nid_tra, "SO2_p_h2soh", "SO2_p_h2soh", "CM-3 S-1",
    354      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    355      .                "ave(X)", zsto,zout)
    356          CALL histdef(nid_tra, "SO2_p_dmsooh", "SO2_p_dmsooh", "CM-3 S-1",
    357      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    358      .                "ave(X)", zsto,zout)
    359          CALL histdef(nid_tra, "DMSO_p_dmsoh", "DMSO_p_dmsoh", "CM-3 S-1",
    360      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    361      .                "ave(X)", zsto,zout)
    362          CALL histdef(nid_tra, "ASMSAM_p_dmsooh", "ASMSAM_p_dmsooh", "CM-3 S-1",
    363      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    364      .                "ave(X)", zsto,zout)
    365          CALL histdef(nid_tra, "ASSO4M_p_so2oh", "ASSO4M_p_so2oh", "CM-3 S-1",
    366      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    367      .                "ave(X)", zsto,zout)
    368          CALL histdef(nid_tra, "ASSO4M_p_so2h2o2", "ASSO4M_p_h2o2","CM-3 S-1",
    369      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    370      .                "ave(X)", zsto,zout)
    371          CALL histdef(nid_tra, "ASSO4M_p_so2o3", "ASSO4M_p_so2o3", "CM-3 S-1",
    372      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    373      .                "ave(X)", zsto,zout)
    374 
    375          CALL histdef(nid_tra, "Wet3D_SO2", "Wet3D_SO2",
    376      .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    377      .           "ave(X)", zsto,zout)
    378          CALL histdef(nid_tra, "Wet3D_DMS", "Wet3D_DMS",
    379      .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    380      .           "ave(X)", zsto,zout)
    381          CALL histdef(nid_tra, "Wet3D_HNO3", "Wet3D_HNO3",
    382      .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    383      .           "ave(X)", zsto,zout)
    384          CALL histdef(nid_tra, "Wet3D_H2O2", "Wet3D_H2O2",
    385      .           "kg/m2/s", iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    386      .           "ave(X)", zsto,zout)
    387 
    388          CALL histdef(nid_tra, "PH_HIST", "pH over physics time step", "1",
    389      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    390      .                "ave(X)", zsto,zout)
    391 #endif
    392 
    393 #ifdef INCA_CH4
    394          DO it=1, grpcnt
    395          CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
    396      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    397      .                "ave(X)", zsto,zout)
    398          ENDDO
    399 
    400 #endif
    401 
    402 
    403 #ifdef INCA_NMHC
    404          CALL histdef(nid_tra, "CO2_basprod", "CO2_basprod", "CM-3 S-1",
    405      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    406      .                "ave(X)", zsto,zout)
    407          CALL histdef(nid_tra, "CO2_nmhcprod", "CO2_nmhcprod", "CM-3 S-1",
    408      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    409      .                "ave(X)", zsto,zout)
    410          CALL histdef(nid_tra, "CO2_radicalprod", "CO2_radicalprod", "CM-3 S-1",
    411      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    412      .                "ave(X)", zsto,zout)
    413          CALL histdef(nid_tra, "HNO3_prod", "HNO3_prod", "CM-3 S-1",
    414      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    415      .                "ave(X)", zsto,zout)
    416          CALL histdef(nid_tra, "HNO3_loss", "HNO3_loss", "CM-3 S-1",
    417      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    418      .                "ave(X)", zsto,zout)
    419          CALL histdef(nid_tra, "CO_prod", "CO_prod", "CM-3 S-1",
    420      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    421      .                "ave(X)", zsto,zout)
    422          CALL histdef(nid_tra, "CO_loss", "CO_loss", "CM-3 S-1",
    423      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    424      .                "ave(X)", zsto,zout)
    425 
    426          DO it=1, grpcnt
    427          CALL histdef(nid_tra, grpsym(it), grpsym(it), "VMR",
    428      .                iim,jjphy_nb,nhori, klev,1,klev,nvert, 32,
    429      .                "ave(X)", zsto,zout)
    430          ENDDO
    431 
    432 
    433 #endif
    434 ! Fin ajout Anne
    435 #else
    43663         CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",
    43764     .                iim,jjphy_nb,nhori, 1,1,1, -99, 32,
    … …  
    493120     .                "ave(X)", zsto,zout)
    494121
    495 #endif
    496122c
    497123         CALL histend(nid_tra)
    498 #ifdef INCA_AER
    499          CALL histend(nid_tra2)
    500          CALL histend(nid_tra3)
    501 #endif
    502124         ndex2d = 0
    503125         ndex3d = 0
    504126         ndex = 0
     127#endif

  • LMDZ4/branches/V3_test/libf/phylmd/physiq.F

    r720 r726  
    2525      USE misc_mod, mydebug=>debug
    2626      USE vampir
    27 #ifdef INCA
    28 cym      USE chemshut
    29       USE species_names
    30 #ifdef INCA_CH4
    31 !      USE obs_pos
    32 #endif
    33 #endif
    3427      IMPLICIT none
    3528c======================================================================
    … …  
    19451938           WRITE(lunout,*) 'Appel CHEMINI ...'
    19461939#endif
    1947            CALL chemini( rpi,
     1940           CALL chemini(
    19481941     $                   rg,
    19491942     $                   ra,
    … …  
    30213014
    30223015#ifdef INCA_AER
    3023       call AEROSOL_METEO_CALC(calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs,
    3024      &   prfl,psfl,pctsrf(1,3),airephy,xjour,rlat,rlon)
     3016      call AEROSOL_METEO_CALC(calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs
     3017     &   ,prfl,psfl,pctsrf,airephy,xjour,rlat,rlon,u10m,v10m)
    30253018#endif
    30263019

  • LMDZ4/branches/V3_test/libf/phylmd/phytrac.F

    r704 r726  
    6767      USE iophy
    6868      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
    9870      IMPLICIT none
    9971c======================================================================
    … …  
    165137#ifdef INCA
    166138      REAL flxmass_w(klon,klev)
     139      CHARACTER(len=8) :: solsym(nqmax)
    167140#endif
    168141c      integer iflag_con
    … …  
    230203cAA Pour l'instant seuls les cas du rn et du pb ont ete envisages.
    231204
    232       REAL source(klon)       ! a voir lorsque le flux est prescrit
     205      REAL source(klon,nqmax)       ! a voir lorsque le flux est prescrit
    233206cAA
    234207cAA Pour la source de radon et son reservoir de sol
    … …  
    365338      INTEGER           :: lastgas
    366339      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) 
    409341
    410342      REAL, PARAMETER   :: dry_mass = 28.966
    … …  
    413345      REAL              :: calday
    414346      REAL              :: pdel(klon,klev)
    415       REAL              :: dummy(klon,klev) = 0.
     347      REAL              :: dummy(klon,klev)
    416348#endif
    417349#ifdef INCA_AER
    … …  
    420352c
    421353c======================================================================
     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
    422399         modname='phytrac'
    423400
    … …  
    549526     $                 rneb,         ! for chimiaq
    550527     $                 t_seri,       ! for chimiaq
    551      $                 rh)
     528     $                 rh,
     529     $                 lafin)
    552530! fin changement anne
    553531
    … …  
    555533
    556534        CALL chemmain (tr_seri,    !mmr
    557      $                 nas,        !nas
    558535     $                 nstep,      !nstep
    559536     $                 calday,     !calday
    … …  
    590567     $                 obuf,       !obuf
    591568     $                 iip1,       !nx
    592      $                 jjp1)       !ny
     569     $                 jjp1,       !ny
     570     $                 source,
     571     $                 solsym)
    593572#ifdef INCAINFO
    594573#ifdef INCA_AER
    … …  
    797776C         CALL minmaxqfi(tr_seri(1,1,it),0.,1.e33,'cltracrn it='//itn)
    798777      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
    804779
    805780Cmaf provisoire source / traceur a creer
    806781        DO i=1, klon
    807           source(i) = 0.0 ! pas de source, pour l'instant
     782          source(i,it) = 0.0 ! pas de source, pour l'instant
    808783        ENDDO
    809784C
    810785#endif
    811786          CALL cltrac(pdtphys, coefh,t_seri,
    812      s               tr_seri(1,1,it), source,
     787     s               tr_seri(1,1,it), source(:,it),
    813788     e               paprs, pplay, delp,
    814789     s               d_tr_cl(1,1,it))

  • LMDZ4/branches/V3_test/libf/phylmd/write_histrac.h

    r704 r726  
    22! $Header$
    33!
     4#ifndef INCA
    45      ndex = 0
    56      ndex2d = 0
    … …  
    1314cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,airephy,zx_tmp_2d)     
    1415      CALL histwrite_phy(nid_tra,"aire",itau_w,airephy)
    15 #ifdef INCA
    16 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ps,zx_tmp_2d)
    17       CALL histwrite_phy(nid_tra,"ps",itau_w,ps)
    18 
    19 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ptrop,zx_tmp_2d)
    20       CALL histwrite_phy(nid_tra,"ptrop",itau_w,ptrop)
    21 
    22 C   3d FIELDS
    23 
    24 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,t_seri, zx_tmp_3d)
    25       CALL histwrite_phy(nid_tra,"temp",itau_w,t_seri)
    26 
    27 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,u, zx_tmp_3d)
    28       CALL histwrite_phy(nid_tra,"u",itau_w,u)
    29 
    30 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,v, zx_tmp_3d)
    31       CALL histwrite_phy(nid_tra,"v",itau_w,v)
    32 
    33 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,sh, zx_tmp_3d)
    34       CALL histwrite_phy(nid_tra,"h2o",itau_w,sh)
    35 
    36 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pdel, zx_tmp_3d)
    37       CALL histwrite_phy(nid_tra,"pdel",itau_w,pdel)
    38 
    39 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,pplay, zx_tmp_3d)
    40       CALL histwrite_phy(nid_tra,"pmid",itau_w,pplay)
    41 
    42 ! Ajout Anne
    43 #ifdef INCA_AER
    44 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,airm, zx_tmp_3d)
    45       CALL histwrite_phy(nid_tra,"airm",itau_w,airm)
    46 #endif
    47 
    48 ! Fin ajout Anne
    49 
    50 #ifdef INCA
    51 #ifdef INCAINFO
    52 #ifdef INCA_NMHC
    53       DO it=1, phtcnt
    54       WRITE(str2,'(i2.2)') it
    55 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it),
    56 cym     .     zx_tmp_3d)
    57       CALL histwrite_phy(nid_tra,"j"//str2,itau_w,jrates(:,:,it))
    58       ENDDO
    59 
    60       DO it=1, hetcnt
    61       WRITE(str2,'(i2.2)') it
    62 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it),
    63 cym     .     zx_tmp_3d)
    64       CALL histwrite_phy(nid_tra,"w"//str2,itau_w,hrates(:,:,it))
    65       ENDDO
    66 
    67       DO it=1, extcnt
    68       WRITE(str2,'(i2.2)') it
    69 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it),
    70      .     zx_tmp_3d)
    71       CALL histwrite_phy(nid_tra,"ext"//str2,itau_w,extflx(:,:,it))
    72       ENDDO
    73 #endif
    74 #ifdef INCA_CH4
    75        DO it=1, phtcnt
    76       WRITE(str2,'(i2.2)') it
    77 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,it),
    78 cym     .     zx_tmp_3d)
    79       CALL histwrite_phy(nid_tra,"j"//str2,itau_w,jrates(:,:,it))
    80       ENDDO
    81 
    82       DO it=1, hetcnt
    83       WRITE(str2,'(i2.2)') it
    84 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,it),
    85 cym     .     zx_tmp_3d)
    86       CALL histwrite_phy(nid_tra,"w"//str2,itau_w,hrates(:,:,it))
    87       ENDDO
    88 
    89       DO it=1, extcnt
    90       WRITE(str2,'(i2.2)') it
    91 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,it),
    92 cym   .     zx_tmp_3d)
    93       CALL histwrite_phy(nid_tra,"ext"//str2,itau_w,extflx(:,:,it))
    94       ENDDO
    95 #endif
    96 
    97       DO it=1, nfs
    98       WRITE(str2,'(i2.2)') it
    99 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,invariants(1,1,it),
    100 cym     .     zx_tmp_3d)
    101       CALL histwrite_phy(nid_tra,"INV"//str2,itau_w,invariants(:,:,it))
    102       ENDDO
    103 #else
    104 #ifdef INCA_NMHC
    105 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
    106 cym     .     zx_tmp_3d)
    107       CALL histwrite_phy(nid_tra,"jO3",itau_w,jrates(:,:,2))
    108 
    109 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4),
    110 cym     .     zx_tmp_3d)
    111       CALL histwrite_phy(nid_tra,"jNO2",itau_w,jrates(:,:,4))
    112 
    113 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13),
    114 cym     .     zx_tmp_3d)
    115       CALL histwrite_phy(nid_tra,"jH2O2",itau_w,jrates(:,:,13))
    116 
    117 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1),
    118 cym     .     zx_tmp_3d)
    119       CALL histwrite_phy(nid_tra,"wHNO3",itau_w,hrates(:,:,1))
    120 
    121 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1),
    122 cym     .     zx_tmp_3d)
    123       CALL histwrite_phy(nid_tra,"kN2O5",itau_w,krates(:,:,1))
    124 
    125 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1),
    126 cym     .     zx_tmp_3d)
    127       CALL histwrite_phy(nid_tra,"LghtNO",itau_w,extflx(:,:,1))
    128 #endif
    129 #ifdef INCA_CH4
    130 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,2),
    131 cym     .     zx_tmp_3d)
    132       CALL histwrite_phy(nid_tra,"jO3",itau_w,jrates(:,:,2))
    133 
    134 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,4),
    135 cym     .     zx_tmp_3d)
    136       CALL histwrite_phy(nid_tra,"jNO2",itau_w,jrates(:,:,4))
    137 
    138 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,jrates(1,1,13),
    139 cym     .     zx_tmp_3d)
    140       CALL histwrite_phy(nid_tra,"jH2O2",itau_w,jrates(:,:,13))
    141 
    142 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hrates(1,1,1),
    143 cym     .     zx_tmp_3d)
    144       CALL histwrite_phy(nid_tra,"wHNO3",itau_w,hrates(:,:,1))
    145 
    146 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,krates(1,1,1),
    147 cym     .     zx_tmp_3d)
    148       CALL histwrite_phy(nid_tra,"kN2O5",itau_w,krates(:,:,1))
    149 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,extflx(1,1,1),
    150 cym     .     zx_tmp_3d)
    151       CALL histwrite_phy(nid_tra,"LghtNO",itau_w,extflx(:,:,1))
    152 #endif
    153 #endif
    154 !      DO it=1, grpcnt
    155 !
    156 !cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
    157 !      zx_tmp_3d = zx_tmp_3d * dry_mass / nadv_mass(it)
    158 !      CALL histwrite_phy(nid_tra,grpsym(it),itau_w,zx_tmp_3d,
    159 !     .                                   iim*(jjm+1)*klev,ndex3d)
    160 !      ENDDO
    161 #endif
    162 
    163 #ifdef INCA_AER
    164 
    165       it = id_CIDUSTM
    166 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_st(1,1,it),
    167 cym     .                  zx_tmp_3d)
    168        CALL histwrite_phy(nid_tra2,"scavcoef_st",itau_w,
    169      .                    scavcoef_st(:,:,it))
    170 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,scavcoef_cv(1,1,it),
    171 cym     .                  zx_tmp_3d)
    172        CALL histwrite_phy(nid_tra2,"scavcoef_cv",itau_w,
    173      .                    scavcoef_cv(:,:,it))
    174 
    175 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,angst(1),zx_tmp_2d)
    176        CALL histwrite_phy(nid_tra2,"AngstroemComp",itau_w,angst))
    177 
    178 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,totaerh2o(1,1), zx_tmp_3d)
    179        CALL histwrite_phy(nid_tra2,"TOTAERH2O",itau_w,totaerh2o)
    180 
    181 #endif
    182 #endif
    18316
    18417      DO it=1,nqmax
    18518C champs 2D
    186 #ifdef INCA
    187       IF ( prt_flag_ts(it) == 0 ) CYCLE
    18819
    189 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, eflux(1,it),zx_tmp_2d)
    190       CALL histwrite_phy(nid_tra,"Emi_"//solsym(it),itau_w,eflux(:,it))
    191 
    192 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, dvel(1,it),zx_tmp_2d)
    193       CALL histwrite_phy(nid_tra,"Dep_"//solsym(it),itau_w,dvel(:,it))
    194 #ifdef INCA_AER
    195       call diag(airephy,tr_seri,.false.)
    196 
    197 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,dflux(1,it),zx_tmp_2d)
    198       CALL histwrite_phy(nid_tra2,"Dry_"//solsym(it),itau_w,
    199      .                   dflux(:,it))
    200 
    201       IF  ((it .ge. trmx) .and. (it .le. trnx)) then
    202 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,sflux(1,it),zx_tmp_2d)
    203         CALL histwrite_phy(nid_tra2,"Sed_"//solsym(it),itau_w,
    204      .                     sflux(:,it))
    205        
    206 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,wflux(1,it),zx_tmp_2d)
    207         CALL histwrite_phy(nid_tra2,"Wet_"//solsym(it),itau_w,
    208      .                     wflux(:,it))
    209 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,wsflux(1,it),zx_tmp_2d)
    210         CALL histwrite_phy(nid_tra2,"WetST_"//solsym(it),itau_w,
    211      .                     wsflux(:,it))
    212 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,wcflux(1,it),zx_tmp_2d)
    213         CALL histwrite_phy(nid_tra2,"WetCV_"//solsym(it),itau_w,
    214      .                     wcflux(:,it))
    215        
    216 cym      CALL gr_fi_ecrit(klev, klon,iim,jjm+1,eflux_alt(1,1,it),zx_tmp_3d)
    217         CALL histwrite_phy(nid_tra2,"Emi_alt_"//solsym(it),itau_w,
    218      .                     eflux_alt(:,:,it))
    219        
    220 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,aload(1,it),zx_tmp_2d)
    221         CALL histwrite_phy(nid_tra2,"Load_"//solsym(it),itau_w,
    222      .                     aload(:,it))
    223         CALL histwrite_phy(nid_tra3,"Inst_Load_"//solsym(it),itau_w,
    224      .                     aload(:,it))
    225 
    226 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,sconc(1,it),zx_tmp_2d)
    227         CALL histwrite_phy(nid_tra2,"SConc_"//solsym(it),itau_w,
    228      .                     sconc(:,it))
    229 
    230         do la=1,las
    231 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1,tausum(1,la,it),zx_tmp_2d)
    232           CALL histwrite_phy(nid_tra2,"OD"//cla(la)//"_"//solsym(it),
    233      .                       itau_w,tausum(:,la,it))
    234 cym      CALL gr_fi_ecrit(klev, klon,iim,jjm+1,tau(1,1,la,it),zx_tmp_3d)
    235           CALL histwrite_phy(nid_tra2,"OD3D"//cla(la)//"_"//solsym(it),
    236      .                       itau_w,tau(:,:,la,it))
    237         enddo
    238 
    239 cym      CALL gr_fi_ecrit(klev, klon,iim,jjm+1,md(1,1,it),zx_tmp_3d)
    240         CALL histwrite_phy(nid_tra2,"MD_"//solsym(it),itau_w,
    241      .                     md(:,:,it))
    242 
    243 cym      CALL gr_fi_ecrit(klev, klon,iim,jjm+1,mdw(1,1,it),zx_tmp_3d)
    244         CALL histwrite_phy(nid_tra2,"MDW_"//solsym(it),itau_w,
    245      .                     mdw(:,:,it))
    246 
    247 cym      CALL gr_fi_ecrit(klev, klon,iim,jjm+1,cload(1,1,it),zx_tmp_3d)
    248         CALL histwrite_phy(nid_tra2,"CLOAD_"//solsym(it),
    249      .                     itau_w,cload(:,:,it))
    250 
    251 
    252       endif
    253 
    254 
    255 #endif
    256 C champs 3D
    257 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
    258 
    259        !Prefer vmr to mmr for transported species
    260        if( adv_mass(it) /= 0. ) then
    261 #ifdef INCA_AER
    262          if (it .lt. trmx) then
    263 #endif
    264        zx_tmp_fi3d = tr_seri(:,:,it) * dry_mass / adv_mass(it)
    265 #ifdef INCA_AER
    266          endif
    267 #endif
    268        else
    269 #ifdef INCA_CH4
    270        if ( solsym(it) == 'OX' ) then
    271        zx_tmp_fi3d = tr_seri(:,:,it) * dry_mass / nadv_mass(id_o3)
    272        end if
    273 #endif
    274 #ifdef INCA_NMHC
    275        if ( solsym(it) == 'OX' ) then
    276        zx_tmp_fi3d = tr_seri(:,:,it) * dry_mass / nadv_mass(id_o3)
    277        end if
    278 #endif
    279        end if
    280 
    281        CALL histwrite_phy(nid_tra,solsym(it),itau_w,zx_tmp_fi3d)
    282 #else
    28320
    28421cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,tr_seri(1,1,it),zx_tmp_3d)
    … …  
    30239c---fin Olivia     
    30340     
    304 #endif
    30541      ENDDO
    30642
    307 #ifdef INCA
    308 !#ifdef INCA_CH4
    309 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_col(1), zx_tmp_2d)
    310       CALL histwrite_phy(nid_tra,"O3_column",itau_w,o3_col)
    311 
    312 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, co_col(1), zx_tmp_2d)
    313       CALL histwrite_phy(nid_tra,"CO_column",itau_w,co_col)
    314 
    315 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, ch4_col(1), zx_tmp_2d)
    316       CALL histwrite_phy(nid_tra,"CH4_column",itau_w,ch4_col)
    317 
    318 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, no2_col(1), zx_tmp_2d)
    319       CALL histwrite_phy(nid_tra,"NO2_column",itau_w,no2_col)
    320 
    321 cym      CALL gr_fi_ecrit(1, klon,iim,jjm+1, o3_st_flx(1), zx_tmp_2d)
    322       CALL histwrite_phy(nid_tra,"O3_ste",itau_w,o3_st_flx)
    323 
    324 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,o3_prod(1,1),
    325 cym     .     zx_tmp_3d)
    326       CALL histwrite_phy(nid_tra,"O3_prod",itau_w,o3_prod)
    327 
    328 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,o3_loss(1,1),
    329 cym    .     zx_tmp_3d)
    330       CALL histwrite_phy(nid_tra,"O3_loss",itau_w,o3_loss)
    331 ! Ajout Anne
    332 #ifdef INCA_AER
    333 ! for sulfur cycle
    334 
    335 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsoh(1,1),
    336 cym     .     zx_tmp_3d)
    337       CALL histwrite_phy(nid_tra,"SO2_p_dmsoh",itau_w,SO2_p_dmsoh)
    338 
    339 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsno3(1,1),
    340 cym     .     zx_tmp_3d)
    341       CALL histwrite_phy(nid_tra,"SO2_p_dmsno3",itau_w,SO2_p_dmsno3)
    342 
    343 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_h2soh(1,1),
    344 cym     .     zx_tmp_3d)
    345       CALL histwrite_phy(nid_tra,"SO2_p_h2soh",itau_w,SO2_p_h2soh)
    346 
    347 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,SO2_p_dmsooh(1,1),
    348 cym     .     zx_tmp_3d)
    349       CALL histwrite_phy(nid_tra,"SO2_p_dmsooh",itau_w,SO2_p_dmsooh)
    350 
    351 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,DMSO_p_dmsoh(1,1),
    352 cym     .     zx_tmp_3d)
    353       CALL histwrite_phy(nid_tra,"DMSO_p_dmsoh",itau_w,DMSO_p_dmsoh)
    354 
    355 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASMSAM_p_dmsooh(1,1),
    356 cym     .     zx_tmp_3d)
    357       CALL histwrite_phy(nid_tra,"ASMSAM_p_dmsooh",itau_w,ASMSAM_p_dmsooh)
    358 
    359 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2oh(1,1),
    360 cym     .     zx_tmp_3d)
    361       CALL histwrite_phy(nid_tra,"ASSO4M_p_so2oh",itau_w,ASSO4M_p_so2oh)
    362 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2h2o2(1,1),
    363 cym     .     zx_tmp_3d)
    364       CALL histwrite_phy(nid_tra,"ASSO4M_p_so2h2o2",itau_w,ASSO4M_p_so2h2o2)
    365 
    366 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,ASSO4M_p_so2o3(1,1),
    367 cym     .     zx_tmp_3d)
    368       CALL histwrite_phy(nid_tra,"ASSO4M_p_so2o3",itau_w,ASSO4M_p_so2o3)
    369 
    370 c closing the sulfur budget
    371 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_so2(1,1), zx_tmp_3d)
    372       CALL histwrite_phy(nid_tra,"Wet3D_SO2",itau_w,wet3d_so2)
    373 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_dms(1,1), zx_tmp_3d)
    374       CALL histwrite_phy(nid_tra,"Wet3D_DMS",itau_w,wet3d_dms)
    375 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_hno3(1,1), zx_tmp_3d)
    376       CALL histwrite_phy(nid_tra,"Wet3D_HNO3",itau_w,wet3d_hno3)
    377 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wet3d_h2o2(1,1), zx_tmp_3d)
    378       CALL histwrite_phy(nid_tra,"Wet3D_H2O2",itau_w,wet3d_h2o2)
    379 
    380 
    381 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,PH_HIST(1,1),
    382 cym     .     zx_tmp_3d)
    383       CALL histwrite_phy(nid_tra,"PH_HIST",itau_w,PH_HIST)
    384 #endif
    385 
    386 #ifdef INCA_CH4
    387       DO it=1, grpcnt
    388 
    389 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
    390       zx_tmp_fi3d = nas(:,:,it) * dry_mass / nadv_mass(it)
    391       CALL histwrite_phy(nid_tra,grpsym(it),itau_w,zx_tmp_fi3d)
    392       ENDDO
    393 
    394 #endif
    395 
    396 
    397 #ifdef INCA_NMHC
    398 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_basprod(1,1),
    399 cym     .     zx_tmp_3d)
    400       CALL histwrite_phy(nid_tra,"CO2_basprod",itau_w,CO2_basprod)
    401 
    402 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_nmhcprod(1,1),
    403 cym     .     zx_tmp_3d)
    404       CALL histwrite_phy(nid_tra,"CO2_nmhcprod",itau_w,CO2_nmhcprod)
    405 
    406 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,CO2_radicalprod(1,1),
    407 cym     .     zx_tmp_3d)
    408       CALL histwrite_phy(nid_tra,"CO2_radicalprod",itau_w,
    409      .                   CO2_radicalprod)
    410 
    411 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_prod(1,1),
    412 cym     .     zx_tmp_3d)
    413       CALL histwrite_phy(nid_tra,"HNO3_prod",itau_w,hno3_prod)
    414 
    415 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,hno3_loss(1,1),
    416 cym     .     zx_tmp_3d)
    417       CALL histwrite_phy(nid_tra,"HNO3_loss",itau_w,hno3_loss)
    418 
    419 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_prod(1,1),
    420 cym     .     zx_tmp_3d)
    421       CALL histwrite_phy(nid_tra,"CO_prod",itau_w,co_prod)
    422 
    423 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,co_loss(1,1),
    424 cym     .     zx_tmp_3d)
    425       CALL histwrite_phy(nid_tra,"CO_loss",itau_w,co_loss)
    426 
    427       DO it=1, grpcnt
    428 
    429 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,nas(1,1,it),zx_tmp_3d)
    430       zx_tmp_3d = nas(:,:,it) * dry_mass / nadv_mass(it)
    431       CALL histwrite_phy(nid_tra,grpsym(it),itau_w,zx_tmp_fi3d)
    432       ENDDO
    433 
    434 
    435 #endif
    436 
    437 ! Fin ajout Anne
    438 
    439 !     ... Special section for daytime averaging
    440 ! cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,day_cnt(1,1),
    441 !    .       zx_tmp_3d)
    442 !       CALL histwrite_phy(nid_tra,"day_cnt",itau_w,zx_tmp_3d,
    443 !    .                                  iim*(jjm+1)*klev,ndex3d)
    444 ! cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1,no_daytime(1,1),
    445 !    .       zx_tmp_3d)
    446 !       CALL histwrite_phy(nid_tra,"NO_day",itau_w,zx_tmp_3d,
    447 !    .                                  iim*(jjm+1)*klev,ndex3d)
    448 
    449 !#endif
    450 #else
    45143
    45244C abder
    … …  
    50294
    50395c abder
    504 #endif
    50596
    50697      if (ok_sync) then
    50798         call histsync(nid_tra)
    508 #ifdef INCA_AER
    509          call histsync(nid_tra2)
    510          call histsync(nid_tra3)
    511 #endif
    51299       endif
    513100
     101#endif
    514102
     103
     104
Note: See TracChangeset for help on using the changeset viewer.