Changeset 5251 for LMDZ6/trunk


Ignore:
Timestamp:
Oct 22, 2024, 12:31:08 PM (3 months ago)
Author:
abarral
Message:

Wrap uses of cpp key INCA
Add INCA wrapper

Location:
LMDZ6/trunk/libf
Files:
13 edited
1 copied

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3d/dynetat0.F90

    r5214 r5251  
    2525#endif
    2626  USE iso_params_mod   ! tnat_* and alpha_ideal_*
     27  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    2728
    2829  IMPLICIT NONE
     
    146147    ix = strIdx(newHNO3, var); IF(ix /= 0) oldVar = oldHNO3(ix)                          !--- REPROBUS HNO3 exceptions
    147148#endif
    148 #ifdef INCA
     149IF (CPPKEY_INCA) THEN
    149150    IF(var == 'O3') oldVar = 'OX'                                                        !--- DEAL WITH INCA OZONE EXCEPTION
    150 #endif
     151END IF
    151152    !--------------------------------------------------------------------------------------------------------------------------
    152153    IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr .AND. .NOT.lSkip) THEN                !=== REGULAR CASE: AVAILABLE VARIABLE
  • LMDZ6/trunk/libf/dyn3d_common/infotrac.F90

    r5237 r5251  
    121121   USE CHEM_REP, ONLY: Init_chem_rep_trac
    122122#endif
     123   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    123124   IMPLICIT NONE
    124125!==============================================================================================================================
     
    147148! Local variables
    148149   INTEGER, ALLOCATABLE :: hadv(:), vadv(:)                          !--- Horizontal/vertical transport scheme number
    149 #ifdef INCA
    150150   INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA
    151151                           vad (:), vadv_inca(:),  pbl_flg_inca(:)
    152152   CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:)                   !--- Tracers names for INCA
    153153   INTEGER :: nqINCA
    154 #endif
    155154   CHARACTER(LEN=2)      ::   suff(9)                                !--- Suffixes for schemes of order 3 or 4 (Prather)
    156155   CHARACTER(LEN=3)      :: descrq(30)                               !--- Advection scheme description tags
     
    199198   SELECT CASE(type_trac)
    200199      CASE('inca', 'inco')
    201 #ifndef INCA
     200IF (.NOT. CPPKEY_INCA) THEN
    202201         CALL abort_gcm(modname, 'You must add cpp key INCA and compile with INCA code', 1)
    203 #endif
     202END IF
    204203      CASE('repr')
    205204#ifndef REPROBUS
     
    232231   IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac)) THEN      !=== FOUND OLD STYLE INCA "traceur.def"
    233232   !---------------------------------------------------------------------------------------------------------------------------
    234 #ifdef INCA
     233IF (CPPKEY_INCA) THEN
    235234      nqo = SIZE(tracers) - nqCO2
    236235      CALL Init_chem_inca_trac(nqINCA)                               !--- Get nqINCA from INCA
     
    263262      IF(setGeneration(tracers)) CALL abort_gcm(modname,'See above',1) !- SET FIELDS %iGeneration, %gen0Name
    264263      DEALLOCATE(had, hadv_inca, vad, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca)
    265 #endif
     264END IF
    266265   !---------------------------------------------------------------------------------------------------------------------------
    267266   ELSE                                                              !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE)
     
    271270   nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O')     !--- Passed to phytrac
    272271   nqCO2  =      COUNT( [type_trac == 'inco', type_trac == 'co2i'] )
    273 #ifdef INCA
     272IF (CPPKEY_INCA) THEN
    274273   nqINCA =      COUNT(tracers(:)%component == 'inca')
    275 #endif
     274END IF
    276275   IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "hadv"', 1)
    277276   IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "vadv"', 1)
     
    401400   CALL msg('niso   = '//TRIM(int2str(niso)),   modname)
    402401   CALL msg('ntiso  = '//TRIM(int2str(ntiso)),  modname)
    403 #ifdef INCA
     402IF (CPPKEY_INCA) THEN
    404403   CALL msg('nqCO2  = '//TRIM(int2str(nqCO2)),  modname)
    405404   CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname)
    406 #endif
     405END IF
    407406   t => tracers
    408407   CALL msg('Information stored in '//TRIM(modname)//': ', modname)
  • LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90

    r5214 r5251  
    2626#endif
    2727  USE iso_params_mod   ! tnat_* and alpha_ideal_*
     28  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    2829
    2930  IMPLICIT NONE
     
    171172    ix = strIdx(newHNO3, var); IF(ix /= 0) oldVar = oldHNO3(ix)                          !--- REPROBUS HNO3 exceptions
    172173#endif
    173 #ifdef INCA
     174IF (CPPKEY_INCA) THEN
    174175    IF(var == 'O3') oldVar = 'OX'                                                        !--- DEAL WITH INCA OZONE EXCEPTION
    175 #endif
     176END IF
    176177    !--------------------------------------------------------------------------------------------------------------------------
    177178    IF(NF90_INQ_VARID(fID, var, vID) == NF90_NoErr .AND. .NOT.lSkip) THEN                !=== REGULAR CASE: AVAILABLE VARIABLE
  • LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.F90

    r5246 r5251  
    4646         xios_set_current_context, &
    4747         using_xios
     48   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    4849
    4950  IMPLICIT NONE
     
    15191520     endif
    15201521
    1521 #ifdef INCA
     1522IF (CPPKEY_INCA) THEN
    15221523     IF (ANY(type_trac == ['inca','inco'])) THEN
    15231524        CALL finalize_inca
     
    15291530!$OMP END MASTER
    15301531     ENDIF
    1531 #endif
     1532END IF
    15321533#ifdef REPROBUS
    15331534     if (type_trac == 'repr') CALL finalize_reprobus
     
    15761577!$OMP END MASTER
    15771578
    1578 #ifdef INCA
     1579IF (CPPKEY_INCA) THEN
    15791580          IF (ANY(type_trac == ['inca','inco'])) THEN
    15801581             CALL finalize_inca
     
    15861587!$OMP END MASTER
    15871588          ENDIF
    1588 #endif
     1589END IF
    15891590#ifdef REPROBUS
    15901591          if (type_trac == 'repr') CALL finalize_reprobus
     
    17491750!$OMP END MASTER
    17501751
    1751 #ifdef INCA
     1752IF (CPPKEY_INCA) THEN
    17521753             IF (ANY(type_trac == ['inca','inco'])) THEN
    17531754                CALL finalize_inca
     
    17601761             ENDIF
    17611762
    1762 #endif
     1763END IF
    17631764#ifdef REPROBUS
    17641765             if (type_trac == 'repr') CALL finalize_reprobus
     
    18631864!$OMP END MASTER
    18641865
    1865 #ifdef INCA
     1866IF (CPPKEY_INCA) THEN
    18661867  IF (ANY(type_trac == ['inca','inco'])) THEN
    18671868     CALL finalize_inca
     
    18741875  ENDIF
    18751876
    1876 #endif
     1877END IF
    18771878#ifdef REPROBUS
    18781879  if (type_trac == 'repr') CALL finalize_reprobus
  • LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

    r4050 r5251  
    7777  USE infotrac,       ONLY: nqtot, tracers
    7878  USE filtreg_mod
     79  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    7980  IMPLICIT NONE
    8081!-------------------------------------------------------------------------------
     
    145146!*******************************************************************************
    146147! Look for ozone tracer:
    147 #ifndef INCA
     148IF (.NOT. CPPKEY_INCA) THEN
    148149  DO iq=1,nqtot; IF(strLower(tracers(iq)%name)=="o3") EXIT; END DO
    149150  IF(iq/=nqtot+1) THEN
     
    151152    CALL press_coefoz
    152153    CALL regr_pr_o3(p3d, q3d(:,:,:,iq))
    153     q3d(:,:,:,iq)=q3d(:,:,:,iq)*48./ 29.                !--- Mole->mass fraction         
     154    q3d(:,:,:,iq)=q3d(:,:,:,iq)*48./ 29.                !--- Mole->mass fraction
    154155  END IF
    155 #endif
     156END IF
    156157  q3d(iip1,:,:,:)=q3d(1,:,:,:)
    157158
  • LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90

    r4389 r5251  
    3838  USE phyaqua_mod, ONLY: iniaqua
    3939  USE comconst_mod, ONLY: omeg, rad
    40 #ifdef INCA
    4140  USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic
    4241#ifdef CPP_PARA
     
    4544#endif
    4645  USE mod_phys_lmdz_omp_data, ONLY: klon_omp
    47 #endif
    4846  USE ioipsl_getin_p_mod, ONLY: getin_p
    4947  USE slab_heat_transp_mod, ONLY: ini_slab_transp_geom
     48  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    5049  IMPLICIT NONE
    5150
     
    167166
    168167  IF (ANY(type_trac == ['inca','inco'])) THEN
    169 #ifdef INCA
     168IF (CPPKEY_INCA) THEN
    170169     CALL init_inca_dim_reg(nbp_lon, nbp_lat - 1, &
    171170          rlonudyn, rlatudyn, rlonvdyn, rlatvdyn)
    172 #endif
     171END IF
    173172  END IF
    174173
  • LMDZ6/trunk/libf/phy_common/mod_phys_lmdz_mpi_data.F90

    r4602 r5251  
    33!
    44MODULE mod_phys_lmdz_mpi_data
    5   USE lmdz_mpi
     5  USE lmdz_mpi
     6  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    67
    78  INTEGER,SAVE :: ii_begin
     
    174175    CALL Print_module_data
    175176
    176 #ifdef INCA
     177IF (CPPKEY_INCA) THEN
    177178       CALL init_inca_para( &
    178179            nbp_lon,nbp_lat,nbp_lev,klon_glo,mpi_size, &
    179180            distrib,nvertex, grid_type,communicator)
    180 #endif
     181END IF
    181182   
    182183  END SUBROUTINE Init_phys_lmdz_mpi_data
  • LMDZ6/trunk/libf/phy_common/physics_distribution_mod.F90

    r4389 r5251  
    1717  USE CHEM_REP, ONLY : Init_chem_rep_phys
    1818#endif
     19  USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    1920
    2021  IMPLICIT NONE
     
    3334    CALL init_dimphy(klon_omp,nbp_lev)
    3435
    35 #ifdef INCA
     36IF (CPPKEY_INCA) THEN
    3637       CALL Init_inca_dim(klon_omp,nbp_lev)
    37 #endif
     38END IF
    3839
    3940#ifdef REPROBUS
  • LMDZ6/trunk/libf/phylmd/infotrac_phy.F90

    r5237 r5251  
    138138   USE CHEM_REP, ONLY: Init_chem_rep_trac
    139139#endif
     140   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    140141   IMPLICIT NONE
    141142!==============================================================================================================================
     
    164165! Local variables
    165166   INTEGER, ALLOCATABLE :: hadv(:), vadv(:)                          !--- Horizontal/vertical transport scheme number
    166 #ifdef INCA
    167167   INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA
    168168                           vad (:), vadv_inca(:),  pbl_flg_inca(:)
    169169   CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:)                   !--- Tracers names for INCA
    170170   INTEGER :: nqINCA
    171 #endif
    172171#ifdef CPP_StratAer
    173172   CHARACTER(LEN=maxlen), ALLOCATABLE :: tnames(:)
     
    224223   SELECT CASE(type_trac)
    225224      CASE('inca', 'inco')
    226 #ifndef INCA
     225IF (.NOT. CPPKEY_INCA) THEN
    227226         CALL abort_physic(modname, 'You must add cpp key INCA and compile with INCA code', 1)
    228 #endif
     227END IF
    229228      CASE('repr')
    230229#ifndef REPROBUS
     
    267266   IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac) .AND. lInit) THEN  !=== FOUND OLD STYLE INCA "traceur.def"
    268267   !---------------------------------------------------------------------------------------------------------------------------
    269 #ifdef INCA
     268IF (CPPKEY_INCA) THEN
    270269      nqo = SIZE(tracers) - nqCO2
    271270      CALL Init_chem_inca_trac(nqINCA)                               !--- Get nqINCA from INCA
     
    298297      IF(setGeneration(tracers)) CALL abort_physic(modname,'See below',1) !- SET FIELDS %iGeneration, %gen0Name
    299298      DEALLOCATE(had, hadv_inca, vad, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca)
    300 #endif
     299END IF
    301300   !---------------------------------------------------------------------------------------------------------------------------
    302301   ELSE                                                              !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE)
     
    306305   nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O')     !--- Passed to phytrac
    307306   nqCO2  =      COUNT( [type_trac == 'inco', type_trac == 'co2i'] )
    308 #ifdef INCA
     307IF (CPPKEY_INCA) THEN
    309308   nqINCA =      COUNT(tracers(:)%component == 'inca')
    310 #endif
     309END IF
    311310   IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "hadv"', 1)
    312311   IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "vadv"', 1)
     
    418417   CALL msg('niso   = '//TRIM(int2str(niso)),   modname)
    419418   CALL msg('ntiso  = '//TRIM(int2str(ntiso)),  modname)
    420 #ifdef INCA
     419IF (CPPKEY_INCA) THEN
    421420   CALL msg('nqCO2  = '//TRIM(int2str(nqCO2)),  modname)
    422421   CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname)
    423 #endif
     422END IF
    424423   t => tracers
    425424   CALL msg('Information stored in '//TRIM(modname)//': ', modname)
  • LMDZ6/trunk/libf/phylmd/oasis.F90

    r4754 r5251  
    122122    USE mod_phys_lmdz_mpi_data, ONLY: klon_mpi_para_nb
    123123    use config_ocean_skin_m, only: activate_ocean_skin
     124    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    124125
    125126! Local variables
     
    356357
    357358
    358 #ifdef INCA
     359IF (CPPKEY_INCA) THEN
    359360    IF (dms_cycle_cpl .OR. n2o_cycle_cpl) THEN
    360361       CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o))
    361362    ENDIF
    362 #endif
     363END IF
    363364 
    364365!************************************************************************************
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r5208 r5251  
    109109
    110110
    111 #ifdef INCA
    112111    USE geometry_mod,      ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo
    113112    USE time_phylmdz_mod,  ONLY: ndays
    114113    USE infotrac_phy,      ONLY: nqCO2
    115 #endif
    116114#ifdef REPROBUS
    117115    USE chem_rep, ONLY: Init_chem_rep_xjour, d_q_rep, d_ql_rep, d_qi_rep, &
     
    120118    USE strataer_emiss_mod, ONLY: strataer_emiss_init
    121119#endif
    122 #if defined INCA || defined REPROBUS
    123120    USE time_phylmdz_mod,    ONLY: annee_ref, day_ini, day_ref, start_time
    124121    USE vertical_layers_mod, ONLY: aps, bps, ap, bp
    125 #endif
    126 
     122    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    127123
    128124#ifdef CPP_RRTM
     
    12431239    REAL, dimension(klon, klev) :: cldfrarad   ! fraction nuageuse
    12441240
    1245 #ifdef INCA
    12461241    REAL :: calday, zxsnow_dummy(klon)
    12471242    ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA
     
    12501245    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_cginca
    12511246    REAL, DIMENSION(klon,klev,nbands) :: init_ccminca
    1252 #endif
    12531247    REAL, DIMENSION(klon,nbtr) :: init_source
    12541248
     
    21782172       !
    21792173       IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL
    2180 #ifdef INCA
     2174IF (CPPKEY_INCA) THEN
    21812175          CALL VTe(VTphysiq)
    21822176          CALL VTb(VTinca)
     
    21882182          config_inca)
    21892183
    2190           CALL init_inca_geometry( & 
     2184          CALL init_inca_geometry( &
    21912185               longitude, latitude, &
    21922186               boundslon, boundslat, &
    2193                cell_area, ind_cell_glo) 
    2194 
    2195           if (grid_type==unstructured) THEN 
     2187               cell_area, ind_cell_glo)
     2188
     2189          if (grid_type==unstructured) THEN
    21962190             CALL chemini(  pplay, &
    21972191                  nbp_lon, nbp_lat, &
     
    22512245          piz_aero(:,:,:,:) = init_pizinca
    22522246          cg_aero(:,:,:,:) = init_cginca
    2253 !         
     2247!
    22542248
    22552249
    22562250          CALL VTe(VTinca)
    22572251          CALL VTb(VTphysiq)
    2258 #endif
     2252END IF
    22592253       ENDIF
    22602254       !
     
    41814175
    41824176    IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL
    4183 #ifdef INCA
     4177IF (CPPKEY_INCA) THEN
    41844178       CALL VTe(VTphysiq)
    41854179       CALL VTb(VTinca)
     
    42344228       CALL VTe(VTinca)
    42354229       CALL VTb(VTphysiq)
    4236 #endif
     4230END IF
    42374231    ENDIF !type_trac = inca or inco
    42384232    IF (type_trac == 'repr') THEN
     
    45254519
    45264520       IF (aerosol_couple.AND.config_inca=='aero') THEN
    4527 #ifdef INCA
     4521IF (CPPKEY_INCA) THEN
    45284522          CALL radlwsw_inca  &
    45294523               (chemistry_couple, kdlon,kflev,dist, rmu0, fract, solaire, &
     
    45454539               cldtaupirad, &
    45464540               topswai_aero, solswai_aero)
    4547 #endif
     4541END IF
    45484542       ELSE
    45494543          !
     
    55305524    !
    55315525    IF (ANY(type_trac == ['inca','inco'])) THEN
    5532 #ifdef INCA
     5526IF (CPPKEY_INCA) THEN
    55335527       CALL VTe(VTphysiq)
    55345528       CALL VTb(VTinca)
     
    55505544       CALL VTe(VTinca)
    55515545       CALL VTb(VTphysiq)
    5552 #endif
     5546END IF
    55535547    ENDIF
    55545548
     
    58085802       IF (using_xios) THEN
    58095803
    5810 #ifdef INCA
    5811           IF (type_trac == 'inca') THEN 
    5812              IF (is_omp_master .AND. grid_type==unstructured) THEN 
     5804IF (CPPKEY_INCA) THEN
     5805          IF (type_trac == 'inca') THEN
     5806             IF (is_omp_master .AND. grid_type==unstructured) THEN
    58135807                CALL finalize_inca
    58145808             ENDIF
    58155809          ENDIF
    5816 #endif
     5810END IF
    58175811
    58185812! close xios physiq context (call LMDZ)
  • LMDZ6/trunk/libf/phylmd/tracinca_mod.F90

    r4500 r5251  
    6262    USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
    6363    USE aero_mod, ONLY : naero_grp
     64    USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    6465    IMPLICIT NONE
    6566   
     
    145146    END DO
    146147 
    147 #ifdef INCA
    148     IF (config_inca == 'aero' .OR. config_inca == 'chem') THEN 
    149        zpmfu(:,:)=pmfu(:,:)       
     148IF (CPPKEY_INCA) THEN
     149    IF (config_inca == 'aero' .OR. config_inca == 'chem') THEN
     150       zpmfu(:,:)=pmfu(:,:)
    150151    ELSE IF (config_inca == 'aeNP') THEN
    151152       zpmfu(:,:)=upwd(:,:)
     
    156157         pplay,pdel,prfl,pmflxr,psfl,    &
    157158         pmflxs,zpmfu,itop_con,ibas_con,  &
    158          pphi,cell_area,nstep,rneb,t_seri, &     
     159         pphi,cell_area,nstep,rneb,t_seri, &
    159160         rh,tau_aero,piz_aero,cg_aero,   &
    160161         rfname,ccm,lafin)
    161 #endif
    162 
    163 
    164 #ifdef INCA
     162END IF
     163
     164
     165IF (CPPKEY_INCA) THEN
    165166    CALL chemmain (tr_seri, &   !mmr
    166167         nstep,      & !nstep
     
    198199         nbp_lat,    & !ny
    199200         source )
    200 #endif
     201END IF
    201202   
    202203    CALL VTe(VTinca)
  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r5244 r5251  
    109109
    110110
    111 #ifdef INCA
    112111    USE geometry_mod,      ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo
    113112    USE time_phylmdz_mod,  ONLY: ndays
    114113    USE infotrac_phy,      ONLY: nqCO2
    115 #endif
    116114#ifdef REPROBUS
    117115    USE chem_rep, ONLY: Init_chem_rep_xjour, d_q_rep, d_ql_rep, d_qi_rep, &
     
    120118    USE strataer_emiss_mod, ONLY: strataer_emiss_init
    121119#endif
    122 #if defined INCA || defined REPROBUS
    123120    USE time_phylmdz_mod,    ONLY: annee_ref, day_ini, day_ref, start_time
    124121    USE vertical_layers_mod, ONLY: aps, bps, ap, bp
    125 #endif
    126122
    127123
     
    441437       reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra
    442438       USE output_physiqex_mod, ONLY: output_physiqex
     439       USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA
    443440
    444441
     
    13491346    REAL, dimension(klon, klev) :: cldfrarad   ! fraction nuageuse
    13501347
    1351 #ifdef INCA
    13521348    REAL :: calday, zxsnow_dummy(klon)
    13531349    ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA
     
    13561352    REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_cginca
    13571353    REAL, DIMENSION(klon,klev,nbands) :: init_ccminca
    1358 #endif
    13591354    REAL, DIMENSION(klon,nbtr) :: init_source
    13601355
     
    23402335       !
    23412336       IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL
    2342 #ifdef INCA
     2337IF (CPPKEY_INCA) THEN
    23432338          CALL VTe(VTphysiq)
    23442339          CALL VTb(VTinca)
     
    23502345          config_inca)
    23512346
    2352           CALL init_inca_geometry( & 
     2347          CALL init_inca_geometry( &
    23532348               longitude, latitude, &
    23542349               boundslon, boundslat, &
    2355                cell_area, ind_cell_glo) 
    2356 
    2357           if (grid_type==unstructured) THEN 
     2350               cell_area, ind_cell_glo)
     2351
     2352          if (grid_type==unstructured) THEN
    23582353             CALL chemini(  pplay, &
    23592354                  nbp_lon, nbp_lat, &
     
    24132408          piz_aero(:,:,:,:) = init_pizinca
    24142409          cg_aero(:,:,:,:) = init_cginca
    2415 !         
     2410!
    24162411
    24172412
    24182413          CALL VTe(VTinca)
    24192414          CALL VTb(VTphysiq)
    2420 #endif
     2415END IF
    24212416       ENDIF
    24222417       !
     
    56355630
    56365631    IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL
    5637 #ifdef INCA
     5632IF (CPPKEY_INCA) THEN
    56385633       CALL VTe(VTphysiq)
    56395634       CALL VTb(VTinca)
     
    56885683       CALL VTe(VTinca)
    56895684       CALL VTb(VTphysiq)
    5690 #endif
     5685END IF
    56915686    ENDIF !type_trac = inca or inco
    56925687    IF (type_trac == 'repr') THEN
     
    59785973
    59795974       IF (aerosol_couple.AND.config_inca=='aero') THEN
    5980 #ifdef INCA
     5975IF (CPPKEY_INCA) THEN
    59815976          CALL radlwsw_inca  &
    59825977               (chemistry_couple, kdlon,kflev,dist, rmu0, fract, solaire, &
     
    59985993               cldtaupirad, &
    59995994               topswai_aero, solswai_aero)
    6000 #endif
     5995END IF
    60015996       ELSE
    60025997          !
     
    71097104    !
    71107105    IF (ANY(type_trac == ['inca','inco'])) THEN
    7111 #ifdef INCA
     7106IF (CPPKEY_INCA) THEN
    71127107       CALL VTe(VTphysiq)
    71137108       CALL VTb(VTinca)
     
    71297124       CALL VTe(VTinca)
    71307125       CALL VTb(VTphysiq)
    7131 #endif
     7126END IF
    71327127    ENDIF
    71337128
     
    74407435       IF (using_xios) THEN
    74417436
    7442 #ifdef INCA
    7443           IF (type_trac == 'inca') THEN 
    7444              IF (is_omp_master .AND. grid_type==unstructured) THEN 
     7437IF (CPPKEY_INCA) THEN
     7438          IF (type_trac == 'inca') THEN
     7439             IF (is_omp_master .AND. grid_type==unstructured) THEN
    74457440                CALL finalize_inca
    74467441             ENDIF
    74477442          ENDIF
    7448 #endif
     7443END IF
    74497444
    74507445! close xios physiq context (call LMDZ)
Note: See TracChangeset for help on using the changeset viewer.