Changeset 960 for LMDZ4/trunk/libf/dyn3d


Ignore:
Timestamp:
May 27, 2008, 6:46:50 PM (17 years ago)
Author:
lsce
Message:
  • Ajoute du parametre config_inca dans conf_gcm.F config_inca='none'(sans INCA, par defaut) config_inca='chem'(avec INCA config chemie) config_inca='aero'(avec INCA config aerosol)
  • Menage parmis les cles CPP INCA
  • Enleve le calcul d'omega dans calfis.F et active le calcul correspondant dans advtrac.F(avant uniquement pour INCA).

JG

Location:
LMDZ4/trunk/libf/dyn3d
Files:
10 edited

Legend:

Unmodified
Added
Removed

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

    r762 r960  
    44c
    55c
    6 #ifdef INCA
    76      SUBROUTINE advtrac(pbaru,pbarv ,
    87     *                   p,  masse,q,iapptrac,teta,
    98     *                  flxw,
    109     *                  pk)
    11 #else
    12       SUBROUTINE advtrac(pbaru,pbarv ,
    13      *                   p,  masse,q,iapptrac,teta,
    14      *                  pk)
    15 #endif
    1610c     Auteur :  F. Hourdin
    1711c
    … …  
    4842      REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
    4943      REAL pk(ip1jmp1,llm)
    50 #ifdef INCA
    51       REAL               :: flxw(ip1jmp1,llm)
    52 #endif
     44      REAL flxw(ip1jmp1,llm)
    5345
    5446c-------------------------------------------------------------
    … …  
    124116        CALL groupe( massem, pbaruc,pbarvc, pbarug,pbarvg,wg )
    125117
    126 #ifdef INCA
    127118      ! ... Flux de masse diaganostiques traceurs
    128119      flxw = wg / FLOAT(iapp_tracvl)
    129 #endif
    130120
    131121c  test sur l'eventuelle creation de valeurs negatives de la masse
    … …  
    209199        else if (iadv(iq).eq.20) then
    210200            call pentes_ini (q(1,1,iq),wg,massem,pbarug,pbarvg,0)
    211 #ifdef INCA
    212        do iiq = iq+1, iq+3
    213 c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
    214           q(:,:,iiq)=q(:,:,iiq)*1
    215        enddo
    216 #endif
    217201
    218202c   ----------------------------------------------------------------
    … …  
    228212           call  prather(q(1,1,iq),wg,massem,pbarug,pbarvg,
    229213     s                     n,dtbon)
    230 #ifdef INCA
    231        do iiq = iq+1, iq+9
    232 c         q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
    233          q(:,:,iiq)=q(:,:,iiq)*1
    234        enddo
    235 #endif
     214
    236215c   ----------------------------------------------------------------
    237216c   Schemas PPM Lin et Rood

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

    r762 r960  
    44c
    55c
    6 #ifdef INCA
    76            SUBROUTINE caladvtrac(q,pbaru,pbarv ,
    87     *                   p ,masse, dq ,  teta,
    9      *                   flxw,
    10      *                   pk)
    11 #else
    12             SUBROUTINE caladvtrac(q,pbaru,pbarv ,
    13      *                   p ,masse, dq ,  teta,
    14      *                   pk)
    15 #endif
    16 
     8     *                   flxw, pk)
    179c
    1810      IMPLICIT NONE
    … …  
    3931      REAL p( ip1jmp1,llmp1),q( ip1jmp1,llm,nqmx),dq( ip1jmp1,llm,2 )
    4032      REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
    41 #ifdef INCA
    4233      REAL               :: flxw(ip1jmp1,llm)
    43 #endif
    4434
    4535c  ..................................................................
    … …  
    6959c   advection
    7060
    71 #ifdef INCA
    72       CALL advtrac( pbaru,pbarv,
    73      *             p,  masse,q,iapptrac, teta,
    74      .             flxw,
    75      .             pk)
    76 #else
    77       CALL advtrac( pbaru,pbarv,
    78      *             p,  masse,q,iapptrac, teta,
    79      .             pk)
    80 #endif
     61        CALL advtrac( pbaru,pbarv,
     62     *       p,  masse,q,iapptrac, teta,
     63     .       flxw, pk)
    8164c
    8265

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

    r956 r960  
    2222     $                  pdteta,
    2323     $                  pdq,
    24      $                  pw,
    25 #ifdef INCA
    2624     $                  flxw,
    27 #endif
    2825     $                  clesphy0,
    2926     $                  pdufi,
    … …  
    121118      REAL pdq(iip1,jjp1,llm,nqmx)
    122119c
    123       REAL pw(iip1,jjp1,llm)
    124 
    125120      REAL pps(iip1,jjp1)
    126121      REAL pp(iip1,jjp1,llmp1)
    … …  
    151146      REAL pcvgu(ngridmx,llm), pcvgv(ngridmx,llm)
    152147      REAL pcvgt(ngridmx,llm), pcvgq(ngridmx,llm,2)
    153 c
    154       REAL pvervel(ngridmx,llm)
    155148c
    156149      REAL zdufi(ngridmx,llm),zdvfi(ngridmx,llm)
    … …  
    169162      REAL PVteta(ngridmx,ntetaSTD)
    170163c
    171 #ifdef INCA
    172       REAL flxw(iip1,jjp1,llm)
    173       REAL flxwfi(ngridmx,llm)
    174 #endif
     164      REAL flxw(iip1,jjp1,llm)  ! Flux de masse verticale sur la grille dynamique
     165      REAL flxwfi(ngridmx,llm)  ! Flux de masse verticale sur la grille physiq
    175166c
    176167     
    … …  
    328319
    329320c   ....  Calcul de la vitesse  verticale  ( en Pa*m*s  ou Kg/s )  ....
    330 c
    331       DO l=1,llm
    332         pvervel(1,l)=pw(1,1,l) * g /apoln
    333         ig0=2
    334        DO j=2,jjm
    335            DO i = 1, iim
    336               pvervel(ig0,l) = pw(i,j,l) * g * unsaire(i,j)
    337               ig0 = ig0 + 1
    338            ENDDO
    339        ENDDO
    340         pvervel(ig0,l)=pw(1,jjp1,l) * g /apols
    341       ENDDO
     321c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux
     322c    de masse est calclue dans advtrac.F 
     323c      DO l=1,llm
     324c        pvervel(1,l)=pw(1,1,l) * g /apoln
     325c        ig0=2
     326c       DO j=2,jjm
     327c           DO i = 1, iim
     328c              pvervel(ig0,l) = pw(i,j,l) * g * unsaire(i,j)
     329c              ig0 = ig0 + 1
     330c           ENDDO
     331c       ENDDO
     332c        pvervel(ig0,l)=pw(1,jjp1,l) * g /apols
     333c      ENDDO
    342334
    343335c
    … …  
    442434     $           ntetaSTD,rtetaSTD,PVteta)
    443435c
    444 #ifdef INCA
     436c On change de grille, dynamique vers physiq, pour le flux de masse verticale
    445437      CALL gr_dyn_fi(llm,iip1,jjp1,ngridmx,flxw,flxwfi)
    446 #endif
    447 
    448438
    449439c-----------------------------------------------------------------------
    … …  
    470460     .             ztfi,
    471461     .             zqfi,
    472      .             pvervel,
    473 #ifdef INCA
    474462     .             flxwfi,
    475 #endif
    476463     .             zdufi,
    477464     .             zdvfi,

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

    r956 r960  
    507507       offline = .FALSE.
    508508       CALL getin('offline',offline)
     509
     510!Config  Key  = config_inca
     511!Config  Desc = Choix de configuration de INCA
     512!Config  Def  = none
     513!Config  Help = Choix de configuration de INCA :
     514!Config         'none' = sans INCA
     515!Config         'chem' = INCA avec calcul de chemie
     516!Config         'aero' = INCA avec calcul des aerosols
     517      config_inca = 'none'
     518      CALL getin('config_inca',config_inca)
    509519
    510520
    … …  
    541551      write(lunout,*)' tauyy = ', tauyy
    542552      write(lunout,*)' offline = ', offline
     553      write(lunout,*)' config_inca = ', config_inca
    543554
    544555      RETURN
    … …  
    653664       offline = .FALSE.
    654665       CALL getin('offline',offline)
    655       write(lunout,*)' offline = ', offline
    656 
     666
     667!Config  Key  = config_inca
     668!Config  Desc = Choix de configuration de INCA
     669!Config  Def  = none
     670!Config  Help = Choix de configuration de INCA :
     671!Config         'none' = sans INCA
     672!Config         'chem' = INCA avec calcul de chemie
     673!Config         'aero' = INCA avec calcul des aerosols
     674      config_inca = 'none'
     675      CALL getin('config_inca',config_inca)
    657676
    658677      write(lunout,*)' #########################################'
    … …  
    688707      write(lunout,*)' tauy = ', tauy
    689708      write(lunout,*)' offline = ', offline
     709      write(lunout,*)' config_inca = ', config_inca
    690710c
    691711      RETURN

  • LMDZ4/trunk/libf/dyn3d/control.h

    r956 r960  
    1414     &              iperiod,iapp_tracvl,iconser,iecri,idissip,iphysiq , &
    1515     &              periodav,iecrimoy,dayref,anneeref,                  &
    16      &              raz_date,offline,ip_ebil_dyn
     16     &              raz_date,offline,ip_ebil_dyn,config_inca
    1717
    1818      INTEGER   nday,day_step,iperiod,iapp_tracvl,iconser,iecri,        &
    … …  
    2121      REAL periodav
    2222      logical offline
    23 
     23      CHARACTER*4 config_inca
    2424!-----------------------------------------------------------------------

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

    r822 r960  
    2929#include "indicesol.h"
    3030#include "advtrac.h"
    31 cym#include "dimphy.h"
     31#include  "control.h"
    3232      REAL :: masque(iip1,jjp1)
    3333      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
    … …  
    4343      mmt_adj(:,:,:,:) = 1
    4444
     45      IF (config_inca /= 'none') THEN
    4546#ifdef INCA
    46       call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
    47       print *, 'nbtrac =' , nbtrac
     47         call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
    4848#endif
     49         print *, 'nbtrac =' , nbtrac
     50      END IF
    4951
    5052      CALL Init_Phys_lmdz(iim,jjp1,llm,nqmx-2,1,(jjm-1)*iim+2)

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

    r822 r960  
    521521#include "advtrac.h"
    522522#include "temps.h"
     523#include "control.h"
    523524
    524525      INTEGER nq, l
    … …  
    623624#endif
    624625
    625 #ifdef INCA
     626      IF (config_inca /= 'none') THEN
    626627! Ajout Anne pour lecture valeurs traceurs dans un fichier start_trac.nc
    627       ierr_file = NF_OPEN ("start_trac.nc", NF_NOWRITE,nid_trac)
    628       IF (ierr_file .NE.NF_NOERR) THEN
    629         write(6,*)' Pb d''ouverture du fichier start_trac.nc'
    630         write(6,*)' ierr = ', ierr_file
    631       ENDIF
    632 #endif
     628         ierr_file = NF_OPEN ("start_trac.nc", NF_NOWRITE,nid_trac)
     629         IF (ierr_file .NE.NF_NOERR) THEN
     630            write(6,*)' Pb d''ouverture du fichier start_trac.nc'
     631            write(6,*)' ierr = ', ierr_file
     632         ENDIF
     633      END IF
    633634
    634635      IF(nq.GE.1) THEN
    635        do iq=1,nq
    636 #ifdef INCA
     636      do iq=1,nq
     637
     638         IF (config_inca == 'none') THEN
     639            ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
     640            IF (ierr .NE. NF_NOERR) THEN
     641               PRINT*, "Variable  tname(iq) n est pas definie"
     642               CALL abort
     643            ENDIF
     644#ifdef NC_DOUBLE
     645            ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
     646#else
     647            ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
     648#endif
     649        ELSE ! config_inca = 'chem' ou 'aero'
    637650! lecture de la valeur du traceur dans start_trac.nc
    638          if (ierr_file .ne. 2) then
    639             ierr = NF_INQ_VARID (nid_trac, tname(iq), nvarid_trac)
    640             IF (ierr .NE. NF_NOERR) THEN
    641                PRINT*, tname(iq),"est absent de start_trac.nc"
    642                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
    643                IF (ierr .NE. NF_NOERR) THEN
    644                   PRINT*, "Variable ", tname(iq)," n est pas definie"
    645                   CALL abort
    646                ENDIF
    647 #ifdef NC_DOUBLE
    648                ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
    649 #else
    650                ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
    651 #endif
    652 
    653             ELSE
    654                PRINT*, tname(iq), "est present dans start_trac.nc"
     651           IF (ierr_file .ne. 2) THEN
     652             ierr = NF_INQ_VARID (nid_trac, tname(iq), nvarid_trac)
     653             IF (ierr .NE. NF_NOERR) THEN
     654                PRINT*, tname(iq),"est absent de start_trac.nc"
     655                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
     656                IF (ierr .NE. NF_NOERR) THEN
     657                   PRINT*, "Variable ", tname(iq)," n est pas definie"
     658                   CALL abort
     659                ENDIF
     660#ifdef NC_DOUBLE
     661                ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
     662#else
     663                ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
     664#endif
     665               
     666             ELSE
     667                PRINT*, tname(iq), "est present dans start_trac.nc"
    655668#ifdef NC_DOUBLE
    656669               ierr = NF_GET_VAR_DOUBLE(nid_trac, nvarid_trac, trac_tmp)
    … …  
    658671               ierr = NF_GET_VAR_REAL(nid_trac, nvarid_trac, trac_tmp)
    659672#endif
    660                IF (ierr .NE. NF_NOERR) THEN
    661                   PRINT*, "Lecture echouee pour", tname(iq)
    662                   CALL abort
    663                ENDIF
    664                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
    665                IF (ierr .NE. NF_NOERR) THEN
    666                   PRINT*, "Variable ", tname(iq)," n est pas definie"
    667                   CALL abort
    668                ENDIF
    669 #ifdef NC_DOUBLE
    670                ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,trac_tmp)
    671 #else
    672                ierr = NF_PUT_VAR_REAL (nid,nvarid,trac_tmp)
    673 #endif
    674 
    675           ENDIF
     673                IF (ierr .NE. NF_NOERR) THEN
     674                   PRINT*, "Lecture echouee pour", tname(iq)
     675                   CALL abort
     676                ENDIF
     677                ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
     678                IF (ierr .NE. NF_NOERR) THEN
     679                   PRINT*, "Variable ", tname(iq)," n est pas definie"
     680                   CALL abort
     681                ENDIF
     682#ifdef NC_DOUBLE
     683                ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,trac_tmp)
     684#else
     685                ierr = NF_PUT_VAR_REAL (nid,nvarid,trac_tmp)
     686#endif
     687               
     688             ENDIF ! IF (ierr .NE. NF_NOERR)
    676689! fin lecture du traceur
    677           ELSE  ! si il n'y a pas de fichier start_trac.nc
     690          ELSE                  ! si il n'y a pas de fichier start_trac.nc
    678691!             print *, 'il n y a pas de fichier start_trac'
    679692             ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
    … …  
    687700             ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
    688701#endif
    689           ENDIF
    690 ! endif ifdef INCA
    691 #endif
    692          
    693 #ifndef INCA
    694         ierr = NF_INQ_VARID(nid, tname(iq), nvarid)
    695         IF (ierr .NE. NF_NOERR) THEN
    696            PRINT*, "Variable  tname(iq) n est pas definie"
    697            CALL abort
    698         ENDIF
    699 #ifdef NC_DOUBLE
    700           ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,q(1,1,iq))
    701 #else
    702           ierr = NF_PUT_VAR_REAL (nid,nvarid,q(1,1,iq))
    703 #endif
    704 ! endif ifndef INCA
    705 #endif
     702          ENDIF ! (ierr_file .ne. 2)
     703       END IF   ! config_inca
     704     
    706705      ENDDO
    707706      ENDIF

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

    r956 r960  
    220220!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    221221
     222      IF (config_inca /= 'none') THEN
    222223#ifdef INCA
    223224      call init_const_lmdz(nbtrac,anneeref,dayref,iphysiq,day_step,nday)
    224225      call init_inca_para(iim,jjm+1,klon,1,klon_mpi_para_nb,0)
    225226#endif
    226 
     227      END IF
    227228c
    228229c
    … …  
    259260      endif
    260261
     262      IF (config_inca /= 'none') THEN
    261263#ifdef INCA
    262       call init_inca_dim(klon,llm,iim,jjm,
    263      $     rlonu,rlatu,rlonv,rlatv)
    264 #endif
     264         call init_inca_dim(klon,llm,iim,jjm,
     265     $        rlonu,rlatu,rlonv,rlatv)
     266#endif
     267      END IF
    265268
    266269

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

    r762 r960  
    2727#include "dimensions.h"
    2828#include "advtrac.h"
    29 
     29#include "control.h"
    3030c   local
    3131      character*3 descrq(30)
    … …  
    5555      descrq(30)='PRA'
    5656
     57      IF (config_inca /= 'none') THEN
    5758#ifdef INCA
    58 
    59       CALL init_transport(
    60      $     hadv_flg,
    61      $     vadv_flg,
    62      $     conv_flg,
    63      $     pbl_flg,
    64      $     tracnam)
     59         CALL init_transport(
     60     $        hadv_flg,
     61     $        vadv_flg,
     62     $        conv_flg,
     63     $        pbl_flg,
     64     $        tracnam)
    6565#endif
     66      END IF
    6667
    6768c-----------------------------------------------------------------------
    … …  
    109110c    dans fichier traceur.def
    110111c------------------------------------------------------------------
    111 #ifdef INCA
     112      IF (config_inca == 'none') THEN
     113         print*,'ouverture de traceur.def'
     114         open(90,file='traceur.def',form='formatted',status='old',
     115     s        iostat=ierr)
     116         if(ierr.eq.0) then
     117            print*,'ouverture de traceur.def ok'
     118            read(90,*) nq
     119            print*,'nombre de traceurs ',nq
     120            if (nq.gt.nqmx) then
     121               print*,'nombre de traceurs trop important'
     122               print*,'verifier traceur.def'
     123               stop
     124            endif
     125C     
     126            do iq=1,nq
     127               read(90,999) hadv(iq),vadv(iq),tnom(iq)
     128            end do
     129            close(90) 
     130            PRINT*,'lecture de traceur.def :'   
     131            do iq=1,nq
     132               write(*,*) hadv(iq),vadv(iq),tnom(iq)
     133            end do       
     134         else
     135            print*,'pb ouverture traceur.def'
     136            print*,'ATTENTION on prend des valeurs par defaut'
     137            nq = 4
     138            hadv(1) = 14
     139            vadv(1) = 14
     140            tnom(1) = 'H2Ov'
     141            hadv(2) = 10
     142            vadv(2) = 10
     143            tnom(2) = 'H2Ol'
     144            hadv(3) = 10
     145            vadv(3) = 10
     146            tnom(3) = 'RN'
     147            hadv(4) = 10
     148            vadv(4) = 10
     149            tnom(4) = 'PB'
     150         ENDIF
     151         PRINT*,'Valeur de traceur.def :'
     152         do iq=1,nq
     153            write(*,*) hadv(iq),vadv(iq),tnom(iq)
     154         end do       
     155      ELSE  ! config_inca='aero' ou 'chem'
    112156C le module de chimie fournit les noms des traceurs
    113157C et les schemas d'advection associes.
    114       tnom(1)='H2Ov'
    115       tnom(2)='H2Ol'
    116       nq=nbtrac+2
    117        if (nq.gt.nqmx) then
    118        print*,'nombre de traceurs incompatible INCA/LMDZT'
    119        stop
    120        endif
    121       do iq =3,nq
    122       tnom(iq)=tracnam(iq-2)
    123       end do
    124       do iq =1,nq
    125       hadv(iq)= hadv_flg(iq)
    126       vadv(iq)= vadv_flg(iq)
    127       end do
    128 #else
    129       print*,'ouverture de traceur.def'
    130       open(90,file='traceur.def',form='formatted',status='old',
    131      s     iostat=ierr)
    132       if(ierr.eq.0) then
    133         print*,'ouverture de traceur.def ok'
    134         read(90,*) nq
    135         print*,'nombre de traceurs ',nq
    136         if (nq.gt.nqmx) then
    137           print*,'nombre de traceurs trop important'
    138           print*,'verifier traceur.def'
    139           stop
    140         endif
    141 C
    142         do iq=1,nq
    143           read(90,999) hadv(iq),vadv(iq),tnom(iq)
    144         end do
    145         close(90) 
    146         PRINT*,'lecture de traceur.def :'   
    147         do iq=1,nq
    148           write(*,*) hadv(iq),vadv(iq),tnom(iq)
    149         end do       
    150       else
    151         print*,'pb ouverture traceur.def'
    152         print*,'ATTENTION on prend des valeurs par defaut'
    153         nq = 4
    154         hadv(1) = 14
    155         vadv(1) = 14
    156         tnom(1) = 'H2Ov'
    157         hadv(2) = 10
    158         vadv(2) = 10
    159         tnom(2) = 'H2Ol'
    160         hadv(3) = 10
    161         vadv(3) = 10
    162         tnom(3) = 'RN'
    163         hadv(4) = 10
    164         vadv(4) = 10
    165         tnom(4) = 'PB'
    166       ENDIF
    167       PRINT*,'Valeur de traceur.def :'
    168       do iq=1,nq
    169         write(*,*) hadv(iq),vadv(iq),tnom(iq)
    170       end do       
    171 
    172 #endif
     158         tnom(1)='H2Ov'
     159         tnom(2)='H2Ol'
     160         nq=nbtrac+2
     161         if (nq.gt.nqmx) then
     162            print*,'nombre de traceurs incompatible INCA/LMDZT'
     163            stop
     164         endif
     165         do iq =3,nq
     166            tnom(iq)=tracnam(iq-2)
     167         end do
     168         do iq =1,nq
     169            hadv(iq)= hadv_flg(iq)
     170            vadv(iq)= vadv_flg(iq)
     171         end do
     172      END IF ! config_inca
     173
    173174c a partir du nom court du traceur et du schema d'advection au detemine le nom long.
    174175        iiq=0

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

    r956 r960  
    149149
    150150      REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm)
    151 #ifdef INCA
    152       REAL :: flxw(ip1jmp1,llm)
    153 #endif
     151      REAL :: flxw(ip1jmp1,llm)  ! flux de masse verticale
    154152
    155153c+jld variables test conservation energie
    … …  
    293291      IF( forward. OR . leapf )  THEN
    294292
    295 c
    296 #ifdef INCA
    297              CALL caladvtrac(q,pbaru,pbarv,
    298      *                      p, masse, dq,  teta,
    299      .             flxw,
    300      .             pk)
    301 #else
    302              CALL caladvtrac(q,pbaru,pbarv,
    303      *                      p, masse, dq,  teta,
    304      .             pk)
    305 #endif
    306 
     293         CALL caladvtrac(q,pbaru,pbarv,
     294     *        p, masse, dq,  teta,
     295     .        flxw, pk)
     296         
    307297         IF (offline) THEN
    308298Cmaf stokage du flux de masse pour traceurs OFF-LINE
    … …  
    383373        CALL calfis( nq, lafin ,rdayvrai,time  ,
    384374     $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
    385      $               du,dv,dteta,dq,w,
    386 #ifdef INCA
     375     $               du,dv,dteta,dq,
    387376     $               flxw,
    388 #endif
    389377     $               clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi  )
    390378
Note: See TracChangeset for help on using the changeset viewer.