Ignore:
Timestamp:
Jun 17, 2022, 4:24:49 PM (2 years ago)
Author:
lguez
Message:

Sync latest trunk changes to branch LMDZ-ECRAD.

Location:
LMDZ6/branches/LMDZ-ECRAD
Files:
2 deleted
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/LMDZ-ECRAD

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/calfis.F

    r2604 r4171  
    2929c    Auteur :  P. Le Van, F. Hourdin
    3030c   .........
    31       USE infotrac, ONLY: nqtot, niadv, tname
     31      USE infotrac, ONLY: nqtot, tracers
    3232      USE control_mod, ONLY: planet_type, nsplit_phys
    3333#ifdef CPP_PHYS
    … …  
    135135c    -----------------
    136136
    137       INTEGER i,j,l,ig0,ig,iq,iiq
     137      INTEGER i,j,l,ig0,ig,iq,itr
    138138      REAL zpsrf(ngridmx)
    139139      REAL zplev(ngridmx,llm+1),zplay(ngridmx,llm)
    … …  
    281281c   ---------------
    282282c
     283      itr=0
    283284      DO iq=1,nqtot
    284           iiq=niadv(iq)
     285         IF(.NOT.tracers(iq)%isAdvected) CYCLE
     286         itr = itr + 1
    285287         DO l=1,llm
    286             zqfi(1,l,iq) = pq(1,1,l,iiq)
    287             ig0          = 2
     288            zqfi(1,l,itr) = pq(1,1,l,iq)
     289            ig0           = 2
    288290            DO j=2,jjm
    289291               DO i = 1, iim
    290                   zqfi(ig0,l,iq)  = pq(i,j,l,iiq)
     292                  zqfi(ig0,l,itr) = pq(i,j,l,iq)
    291293                  ig0             = ig0 + 1
    292294               ENDDO
    293295            ENDDO
    294             zqfi(ig0,l,iq) = pq(1,jjp1,l,iiq)
     296            zqfi(ig0,l,itr) = pq(1,jjp1,l,iq)
    295297         ENDDO
    296298      ENDDO
    … …  
    481483         lafin_split=lafin.and.isplit==nsplit_phys
    482484
    483         CALL call_physiq(ngridmx,llm,nqtot,tname,
     485!      if (planet_type=="earth") then
     486        CALL call_physiq(ngridmx,llm,nqtot,tracers(:)%name,
    484487     &                   debut_split,lafin_split,
    485488     &                   jD_cur,jH_cur_split,zdt_split,
    … …  
    490493     &                   flxwfi,pducov,
    491494     &                   zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)
    492                              
    493 !      if (planet_type=="earth") then
    494 !
    495 !         CALL physiq (ngridmx,
    496 !     .             llm,
    497 !     .             debut_split,
    498 !     .             lafin_split,
    499 !     .             jD_cur,
    500 !     .             jH_cur_split,
    501 !     .             zdt_split,
    502 !     .             zplev,
    503 !     .             zplay,
    504 !     .             zphi,
    505 !     .             zphis,
    506 !     .             presnivs,
    507 !     .             zufi,
    508 !     .             zvfi, zrfi,
    509 !     .             ztfi,
    510 !     .             zqfi,
    511 !     .             flxwfi,
    512 !     .             zdufi,
    513 !     .             zdvfi,
    514 !     .             zdtfi,
    515 !     .             zdqfi,
    516 !     .             zdpsrf,
    517 !     .             pducov)
    518495!
    519496!      else if ( planet_type=="generic" ) then
    … …  
    522499!     .             llm,            !! nlayer
    523500!     .             nqtot,          !! nq
    524 !     .             tname,          !! tracer names from dynamical core (given in infotrac)
     501!     .             tracers(:)%name,!! tracer names from dynamical core (given in infotrac)
    525502!     .             debut_split,    !! firstcall
    526503!     .             lafin_split,    !! lastcall
    … …  
    622599      pdqfi(:,:,:,:)=0.
    623600C
     601      itr = 0
    624602      DO iq=1,nqtot
    625          iiq=niadv(iq)
     603         IF(.NOT.tracers(iq)%isAdvected) CYCLE
     604         itr = itr + 1
    626605         DO l=1,llm
    627606            DO i=1,iip1
    628                pdqfi(i,1,l,iiq)    = zdqfi(1,l,iq)
    629                pdqfi(i,jjp1,l,iiq) = zdqfi(ngridmx,l,iq)
     607               pdqfi(i,1,l,iq)    = zdqfi(1,l,itr)
     608               pdqfi(i,jjp1,l,iq) = zdqfi(ngridmx,l,itr)
    630609            ENDDO
    631610            DO j=2,jjm
    632611               ig0=1+(j-2)*iim
    633612               DO i=1,iim
    634                   pdqfi(i,j,l,iiq) = zdqfi(ig0+i,l,iq)
     613                  pdqfi(i,j,l,iq) = zdqfi(ig0+i,l,itr)
    635614               ENDDO
    636                pdqfi(iip1,j,l,iiq) = pdqfi(1,j,l,iq)
     615               pdqfi(iip1,j,l,iq) = pdqfi(1,j,l,itr)
    637616            ENDDO
    638617         ENDDO

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/calfis_loc.F

    r2604 r4171  
    4545      USE Times
    4646#endif
    47       USE infotrac, ONLY: nqtot, niadv, tname
     47      USE infotrac, ONLY: nqtot, tracers
    4848      USE control_mod, ONLY: planet_type, nsplit_phys
    4949#ifdef CPP_PHYS
    … …  
    154154c    -----------------
    155155
    156       INTEGER i,j,l,ig0,ig,iq,iiq
     156      INTEGER i,j,l,ig0,ig,iq,itr
    157157      REAL,ALLOCATABLE,SAVE :: zpsrf(:)
    158158      REAL,ALLOCATABLE,SAVE :: zplev(:,:),zplay(:,:)
    … …  
    366366c
    367367
     368      itr = 0
    368369      DO iq=1,nqtot
    369          iiq=niadv(iq)
     370         IF(.NOT.tracers(iq)%isAdvected) CYCLE
     371         itr = itr + 1
    370372c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    371373         DO l=1,llm
    … …  
    375377             i=index_i(ig0)
    376378             j=index_j(ig0)
    377              zqfi(ig0,l,iq)  = pq(i,j,l,iiq)
     379             zqfi(ig0,l,itr)  = pq(i,j,l,iq)
    378380           enddo
    379381         ENDDO
    … …  
    731733         lafin_split=lafin.and.isplit==nsplit_phys
    732734
    733         CALL call_physiq(klon,llm,nqtot,tname,
     735        CALL call_physiq(klon,llm,nqtot,tracers(:)%name,
    734736     &                   debut_split,lafin_split,
    735737     &                   jD_cur,jH_cur_split,zdt_split,
    … …  
    10691071C
    10701072!cdir NODEP
     1073      itr = 0
    10711074      DO iq=1,nqtot
    1072          iiq=niadv(iq)
     1075         IF(.NOT.tracers(iq)%isAdvected) CYCLE
     1076         itr = itr + 1
    10731077c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
    10741078         DO l=1,llm
    … …  
    10791083              i=index_i(ig0)
    10801084              j=index_j(ig0)
    1081               pdqfi(i,j,l,iiq) = zdqfi(ig0,l,iq)
    1082               if (i==1) pdqfi(iip1,j,l,iiq) = zdqfi(ig0,l,iq)
     1085              pdqfi(i,j,l,iq) = zdqfi(ig0,l,itr)
     1086              if (i==1) pdqfi(iip1,j,l,iq) = zdqfi(ig0,l,itr)
    10831087            ENDDO
    10841088           
    10851089            IF (is_north_pole_dyn) then
    10861090              DO i=1,iip1
    1087                 pdqfi(i,1,l,iiq)    = zdqfi(1,l,iq)
     1091                pdqfi(i,1,l,iq)    = zdqfi(1,l,itr)
    10881092              ENDDO
    10891093            ENDIF
    … …  
    10911095            IF (is_south_pole_dyn) then
    10921096              DO i=1,iip1
    1093                 pdqfi(i,jjp1,l,iiq) = zdqfi(klon,l,iq)
     1097                pdqfi(i,jjp1,l,iq) = zdqfi(klon,l,itr)
    10941098              ENDDO
    10951099            ENDIF

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/inigeomphy_mod.F90

    r3435 r4171  
    9393  ALLOCATE(boundslat_reg(jjm+1,2))
    9494 
    95   DO i=1,iim
    96    boundslon_reg(i,east)=rlonu(i+1)
    97    boundslon_reg(i,west)=rlonu(i)
     95  ! specific handling of the -180 longitude scalar grid point boundaries
     96  boundslon_reg(1,east)=rlonu(1)
     97  boundslon_reg(1,west)=rlonu(iim)-2*PI
     98  DO i=2,iim
     99   boundslon_reg(i,east)=rlonu(i)
     100   boundslon_reg(i,west)=rlonu(i-1)
    98101  ENDDO
    99102
    … …  
    124127    cufi_glo(1) = cu(1)
    125128    cvfi_glo(1) = cv(1)
    126     boundslonfi_glo(1,north_east)=0
     129    boundslonfi_glo(1,north_east)=PI
    127130    boundslatfi_glo(1,north_east)=PI/2
    128     boundslonfi_glo(1,north_west)=2*PI
     131    boundslonfi_glo(1,north_west)=-PI
    129132    boundslatfi_glo(1,north_west)=PI/2
    130     boundslonfi_glo(1,south_west)=2*PI
     133    boundslonfi_glo(1,south_west)=-PI
    131134    boundslatfi_glo(1,south_west)=rlatv(1)
    132     boundslonfi_glo(1,south_east)=0
     135    boundslonfi_glo(1,south_east)=PI
    133136    boundslatfi_glo(1,south_east)=rlatv(1)
    134137    DO j=2,jjm
    … …  
    141144        boundslonfi_glo(k,north_east)=rlonu(i)
    142145        boundslatfi_glo(k,north_east)=rlatv(j-1)
    143         boundslonfi_glo(k,north_west)=rlonu(i+1)
     146        if (i.eq.1) then
     147          ! special case for the first longitude's west bound
     148          boundslonfi_glo(k,north_west)=rlonu(iim)-2*PI
     149          boundslonfi_glo(k,south_west)=rlonu(iim)-2*PI
     150        else
     151          boundslonfi_glo(k,north_west)=rlonu(i-1)
     152          boundslonfi_glo(k,south_west)=rlonu(i-1)
     153        endif
    144154        boundslatfi_glo(k,north_west)=rlatv(j-1)
    145         boundslonfi_glo(k,south_west)=rlonu(i+1)
    146155        boundslatfi_glo(k,south_west)=rlatv(j)
    147156        boundslonfi_glo(k,south_east)=rlonu(i)
    … …  
    154163    cufi_glo(klon_glo) = cu((iim+1)*jjm+1)
    155164    cvfi_glo(klon_glo) = cv((iim+1)*jjm-iim)
    156     boundslonfi_glo(klon_glo,north_east)= 0
     165    boundslonfi_glo(klon_glo,north_east)= PI
    157166    boundslatfi_glo(klon_glo,north_east)= rlatv(jjm)
    158     boundslonfi_glo(klon_glo,north_west)= 2*PI
     167    boundslonfi_glo(klon_glo,north_west)= -PI
    159168    boundslatfi_glo(klon_glo,north_west)= rlatv(jjm)
    160     boundslonfi_glo(klon_glo,south_west)= 2*PI
     169    boundslonfi_glo(klon_glo,south_west)= -PI
    161170    boundslatfi_glo(klon_glo,south_west)= -PI/2
    162     boundslonfi_glo(klon_glo,south_east)= 0
     171    boundslonfi_glo(klon_glo,south_east)= PI
    163172    boundslatfi_glo(klon_glo,south_east)= -Pi/2
    164173

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/phylmd/callphysiq_mod.F90

    r2604 r4171  
    2020
    2121  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
    22   USE control_mod, ONLY: planet_type
     22  USE control_mod, ONLY: planet_type, ok_dyn_xios
    2323  USE physiq_mod, ONLY: physiq
     24#ifdef CPP_XIOS
     25  USE mod_xios_dyn3dmem, ONLY : dyn3d_ctx_handle
     26  USE xios, ONLY : xios_set_current_context
     27#endif
    2428  IMPLICIT NONE
    2529
    … …  
    9094              zdpsrf_omp)
    9195
     96! switching back to LMDZDYN context
     97#ifdef CPP_XIOS
     98!$OMP MASTER
     99  if (ok_dyn_xios) then
     100     CALL xios_set_current_context(dyn3d_ctx_handle)
     101  endif
     102!$OMP END MASTER
     103#endif 
     104
    92105
    93106END SUBROUTINE call_physiq

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/phylmd/ce0l.F90

    r2665 r4171  
    2929  USE iniphysiq_mod,  ONLY: iniphysiq
    3030  USE mod_const_mpi,  ONLY: comm_lmdz
     31
    3132#ifdef CPP_PARA
    3233  USE mod_const_mpi,  ONLY: init_const_mpi
    33   USE parallel_lmdz,  ONLY: init_parallel, mpi_rank, omp_rank
     34  USE parallel_lmdz,  ONLY: init_parallel, mpi_rank, omp_rank, using_mpi
    3435  USE bands,          ONLY: read_distrib, distrib_phys
    3536  USE mod_hallo,      ONLY: init_mod_hallo
    3637  USE mod_interface_dyn_phys, ONLY: init_interface_dyn_phys
    37 #endif
     38#ifdef CPP_XIOS
     39  USE xios, only: xios_finalize
     40#endif
     41#endif
     42
    3843  USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, kappa, omeg, r, rad, &
    3944                          pi, jmp1
    … …  
    5055  include "comgeom2.h"
    5156  include "iniprint.h"
     57#ifdef CPP_MPI
     58  include 'mpif.h'
     59#endif
     60 
    5261  REAL               :: masque(iip1,jjp1)             !--- CONTINENTAL MASK
    5362  REAL               :: phis  (iip1,jjp1)             !--- GROUND GEOPOTENTIAL
    … …  
    6776  REAL, ALLOCATABLE  :: masktmp(:)
    6877
    69 #ifndef CPP_PARA
     78#ifdef CPP_PARA
     79  integer ierr
     80#else
    7081! for iniphysiq in serial mode
    7182  INTEGER,PARAMETER :: mpi_rank=0
    … …  
    168179    ocemask = ocetmp
    169180    IF(dlat_omask(1)<dlat_omask(jml_omask)) THEN
    170       DO j=1,jjp1; ocemask(:,j) = ocetmp(:,jjp1-j+1); END DO
     181       DO j=1,jjp1
     182          ocemask(:,j) = ocetmp(:,jjp1-j+1)
     183       END DO
    171184    END IF
    172185    DEALLOCATE(ocetmp,lon_omask,lat_omask,dlon_omask,dlat_omask)
    … …  
    239252#ifdef CPP_PARA
    240253  END IF
     254#ifdef CPP_XIOS
     255  CALL xios_finalize
     256#endif
     257#ifdef CPP_MPI
     258  IF (using_mpi) call MPI_FINALIZE(ierr)
     259#endif
    241260#endif
    242261

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

    r3435 r4171  
    3939  USE comvert_mod, ONLY: ap, bp, preff, pressure_exner
    4040  USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn, itau_phy, start_time
     41  USE strings_mod, ONLY: strLower
    4142 
    4243  IMPLICIT NONE
    … …  
    7475  USE exner_hyb_m,    ONLY: exner_hyb
    7576  USE exner_milieu_m, ONLY: exner_milieu
    76   USE infotrac,       ONLY: nqtot, tname
     77  USE infotrac,       ONLY: nqtot, tracers
    7778  USE filtreg_mod
    7879  IMPLICIT NONE
    … …  
    8485! Local variables:
    8586  CHARACTER(LEN=256) :: modname, fmt
    86   INTEGER            :: i, j, l, ji, itau, iday
     87  INTEGER            :: i, j, l, ji, itau, iday, iq
    8788  REAL               :: xpn, xps, time, phystep
    8889  REAL, DIMENSION(iip1,jjp1)       :: psol
    … …  
    145146! Look for ozone tracer:
    146147#ifndef INCA
    147   DO i=1,nqtot; IF(ANY(["O3","o3"]==tname(i))) EXIT; END DO
    148   IF(i/=nqtot+1) THEN
     148  DO iq=1,nqtot; IF(strLower(tracers(iq)%name)=="o3") EXIT; END DO
     149  IF(iq/=nqtot+1) THEN
    149150    CALL regr_lat_time_coefoz
    150151    CALL press_coefoz
    151     CALL regr_pr_o3(p3d, q3d(:,:,:,i))
    152     q3d(:,:,:,i)=q3d(:,:,:,i)*48./ 29.                  !--- Mole->mass fraction         
     152    CALL regr_pr_o3(p3d, q3d(:,:,:,iq))
     153    q3d(:,:,:,iq)=q3d(:,:,:,iq)*48./ 29.                !--- Mole->mass fraction         
    153154  END IF
    154155#endif

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

    r3773 r4171  
    9393  USE conf_phys_m, ONLY: conf_phys
    9494  USE init_ssrf_m, ONLY: start_init_subsurf
     95  USE phys_state_var_mod, ONLY: beta_aridity, delta_tsurf, awake_dens, cv_gen, &
     96       ratqs_inter, rneb_ancien
    9597  !use ioipsl_getincom
    9698  IMPLICIT NONE
    … …  
    119121  INTEGER :: flag_aerosol
    120122  INTEGER :: flag_aerosol_strat
     123  INTEGER :: flag_volc_surfstrat
    121124  LOGICAL :: flag_aer_feedback
    122125  LOGICAL :: flag_bc_internal_mixture
    … …  
    126129  LOGICAL :: filtre_oro=.false.
    127130
     131  INCLUDE "compbl.h"
     132  INCLUDE "alpale.h"
     133 
    128134  deg2rad= pi/180.0
    129135  iml=assert_eq(SIZE(masque,1),SIZE(phis,1),TRIM(modname)//" iml")
    … …  
    138144                   iflag_cldcon,                                        &
    139145                   iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs,            &
    140                    ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan,          &
    141                    aerosol_couple,                                      &
    142                    chemistry_couple, flag_aerosol, flag_aerosol_strat,  &
    143                    flag_aer_feedback,                                   &
    144                    flag_bc_internal_mixture, bl95_b0, bl95_b1, &
     146                   ok_ade, ok_aie, ok_alw, ok_cdnc, ok_volcan, flag_volc_surfstrat,     &
     147                   aerosol_couple, chemistry_couple, flag_aerosol, flag_aerosol_strat,  &
     148                   flag_aer_feedback, flag_bc_internal_mixture, bl95_b0, bl95_b1,       &
    145149                   read_climoz, alp_offset)
    146150  CALL phys_state_var_init(read_climoz)
    … …  
    231235  wake_delta_pbl_TKE(:,:,:)=0
    232236  wake_dens(:)=0
     237  awake_dens = 0.
     238  cv_gen = 0.
    233239  ale_bl = 0.
    234240  ale_bl_trig =0.
    … …  
    240246 
    241247  z0m(:,is_oce) = rugmer(:)
    242   z0m(:,is_ter) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
    243   z0m(:,is_lic) = MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
     248  z0m(:,is_ter) = 0.01 !MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
     249  z0m(:,is_lic) = 0.001 !MAX(1.0e-05,zstd(:)*zsig(:)/2.0)
    244250  z0m(:,is_sic) = 0.001
    245251  z0h(:,:)=z0m(:,:)
    … …  
    273279  CALL fonte_neige_init(run_off_lic_0)
    274280  CALL pbl_surface_init( fder, snsrf, qsurf, tsoil )
     281
     282  IF (iflag_pbl>1 .AND. iflag_wake>=1  .AND. iflag_pbl_split >=1) then
     283     delta_tsurf = 0.
     284     beta_aridity = 0.
     285  end IF
     286
     287  ratqs_inter = 0.002
     288  rneb_ancien = 0.
    275289  CALL phyredem( "startphy.nc" )
    276290

  • LMDZ6/branches/LMDZ-ECRAD/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90

    r3677 r4171  
    1616  USE mod_phys_lmdz_para, ONLY: klon_omp ! number of columns (on local omp grid)
    1717  USE vertical_layers_mod, ONLY : init_vertical_layers
    18   USE infotrac, ONLY: nqtot,nqo,nbtr,tname,ttext,type_trac,&
    19                       niadv,conv_flg,pbl_flg,solsym,&
    20                       nqfils,nqdesc,nqdesc_tot,iqfils,iqpere,&
    21                       ok_isotopes,ok_iso_verif,ok_isotrac,&
    22                       ok_init_iso,niso_possibles,tnat,&
    23                       alpha_ideal,use_iso,iqiso,iso_num,&
    24                       iso_indnum,zone_num,phase_num,&
    25                       indnum_fn_num,index_trac,&
    26                       niso,ntraceurs_zone,ntraciso
     18  USE infotrac, ONLY: nbtr,nqCO2,tracers,isotopes,type_trac,conv_flg,pbl_flg,nqtottr
    2719#ifdef CPP_StratAer
    28   USE infotrac, ONLY: nbtr_bin, nbtr_sulgas, id_OCS_strat, &
     20  USE infotrac_phy, ONLY: nbtr_bin, nbtr_sulgas, id_OCS_strat, &
    2921                      id_SO2_strat, id_H2SO4_strat, id_BIN01_strat
    3022
    … …  
    116108 
    117109!!$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/)
    118 !$OMP PARALLEL DEFAULT(SHARED) &
    119110!       Copy all threadprivate variables in temps_mod
    120 !$OMP COPYIN(annee_ref, day_ini, day_ref, start_time)
     111!$OMP PARALLEL DEFAULT(SHARED) COPYIN(annee_ref,day_ini,day_ref,start_time)
    121112
    122113  ! Initialize physical constants in physics:
    … …  
    146137
    147138  ! Initialize tracer names, numbers, etc. for physics
    148   CALL init_infotrac_phy(nqtot,nqo,nbtr,tname,ttext,type_trac,&
    149                          niadv,conv_flg,pbl_flg,solsym,&
    150                          nqfils,nqdesc,nqdesc_tot,iqfils,iqpere,&
    151                          ok_isotopes,ok_iso_verif,ok_isotrac,&
    152                          ok_init_iso,niso_possibles,tnat,&
    153                          alpha_ideal,use_iso,iqiso,iso_num,&
    154                          iso_indnum,zone_num,phase_num,&
    155                          indnum_fn_num,index_trac,&
    156                          niso,ntraceurs_zone,ntraciso&
    157 #ifdef CPP_StratAer
    158                          ,nbtr_bin,nbtr_sulgas&
    159                          ,id_OCS_strat,id_SO2_strat,id_H2SO4_strat,id_BIN01_strat&
    160 #endif
    161                          )
     139  CALL init_infotrac_phy(type_trac, tracers, isotopes, nqtottr, nqCO2, pbl_flg, conv_flg)
    162140
    163141  ! Initializations for Reprobus
    … …  
    172150!$OMP END PARALLEL
    173151
    174   IF (type_trac == 'inca') THEN
    175 #ifdef INCA
    176      call init_const_lmdz( &
    177           anneeref,dayref, iphysiq,day_step,nday,  &
    178           nbsrf, is_oce,is_sic, is_ter,is_lic, calend, &
    179           config_inca)
    180      call init_inca_para( &
    181           nbp_lon,nbp_lat,nbp_lev,klon_glo,mpi_size, &
    182           distrib_phys,communicator)
    183 #endif
    184   END IF
     152
    185153  IF (type_trac == 'repr') THEN
    186154#ifdef REPROBUS
    … …  
    198166  END IF
    199167
    200   IF (type_trac == 'inca') THEN
     168  IF (type_trac == 'inca' .OR. type_trac == 'inco') THEN
    201169#ifdef INCA
    202      CALL init_inca_dim(klon_omp,nbp_lev,nbp_lon,nbp_lat - 1, &
    203           rlonudyn,rlatudyn,rlonvdyn,rlatvdyn)
     170     CALL init_inca_dim_reg(nbp_lon, nbp_lat - 1, &
     171          rlonudyn, rlatudyn, rlonvdyn, rlatvdyn)
    204172#endif
    205     IF (type_trac == 'repr') THEN
    206 #ifdef REPROBUS
    207        CALL Init_chem_rep_phys(klon_omp,nbp_lev)
    208 #endif
    209     END IF
    210173  END IF
    211174
Note: See TracChangeset for help on using the changeset viewer.