Changeset 5251 for LMDZ6/trunk
- Timestamp:
- Oct 22, 2024, 12:31:08 PM (3 months ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 13 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/dynetat0.F90
r5214 r5251 25 25 #endif 26 26 USE iso_params_mod ! tnat_* and alpha_ideal_* 27 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 27 28 28 29 IMPLICIT NONE … … 146 147 ix = strIdx(newHNO3, var); IF(ix /= 0) oldVar = oldHNO3(ix) !--- REPROBUS HNO3 exceptions 147 148 #endif 148 #ifdef INCA 149 IF (CPPKEY_INCA) THEN 149 150 IF(var == 'O3') oldVar = 'OX' !--- DEAL WITH INCA OZONE EXCEPTION 150 #endif 151 END IF 151 152 !-------------------------------------------------------------------------------------------------------------------------- 152 153 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 121 121 USE CHEM_REP, ONLY: Init_chem_rep_trac 122 122 #endif 123 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 123 124 IMPLICIT NONE 124 125 !============================================================================================================================== … … 147 148 ! Local variables 148 149 INTEGER, ALLOCATABLE :: hadv(:), vadv(:) !--- Horizontal/vertical transport scheme number 149 #ifdef INCA150 150 INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA 151 151 vad (:), vadv_inca(:), pbl_flg_inca(:) 152 152 CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:) !--- Tracers names for INCA 153 153 INTEGER :: nqINCA 154 #endif155 154 CHARACTER(LEN=2) :: suff(9) !--- Suffixes for schemes of order 3 or 4 (Prather) 156 155 CHARACTER(LEN=3) :: descrq(30) !--- Advection scheme description tags … … 199 198 SELECT CASE(type_trac) 200 199 CASE('inca', 'inco') 201 #ifndef INCA 200 IF (.NOT. CPPKEY_INCA) THEN 202 201 CALL abort_gcm(modname, 'You must add cpp key INCA and compile with INCA code', 1) 203 #endif 202 END IF 204 203 CASE('repr') 205 204 #ifndef REPROBUS … … 232 231 IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac)) THEN !=== FOUND OLD STYLE INCA "traceur.def" 233 232 !--------------------------------------------------------------------------------------------------------------------------- 234 #ifdef INCA 233 IF (CPPKEY_INCA) THEN 235 234 nqo = SIZE(tracers) - nqCO2 236 235 CALL Init_chem_inca_trac(nqINCA) !--- Get nqINCA from INCA … … 263 262 IF(setGeneration(tracers)) CALL abort_gcm(modname,'See above',1) !- SET FIELDS %iGeneration, %gen0Name 264 263 DEALLOCATE(had, hadv_inca, vad, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca) 265 #endif 264 END IF 266 265 !--------------------------------------------------------------------------------------------------------------------------- 267 266 ELSE !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE) … … 271 270 nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O') !--- Passed to phytrac 272 271 nqCO2 = COUNT( [type_trac == 'inco', type_trac == 'co2i'] ) 273 #ifdef INCA 272 IF (CPPKEY_INCA) THEN 274 273 nqINCA = COUNT(tracers(:)%component == 'inca') 275 #endif 274 END IF 276 275 IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "hadv"', 1) 277 276 IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "vadv"', 1) … … 401 400 CALL msg('niso = '//TRIM(int2str(niso)), modname) 402 401 CALL msg('ntiso = '//TRIM(int2str(ntiso)), modname) 403 #ifdef INCA 402 IF (CPPKEY_INCA) THEN 404 403 CALL msg('nqCO2 = '//TRIM(int2str(nqCO2)), modname) 405 404 CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname) 406 #endif 405 END IF 407 406 t => tracers 408 407 CALL msg('Information stored in '//TRIM(modname)//': ', modname) -
LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90
r5214 r5251 26 26 #endif 27 27 USE iso_params_mod ! tnat_* and alpha_ideal_* 28 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 28 29 29 30 IMPLICIT NONE … … 171 172 ix = strIdx(newHNO3, var); IF(ix /= 0) oldVar = oldHNO3(ix) !--- REPROBUS HNO3 exceptions 172 173 #endif 173 #ifdef INCA 174 IF (CPPKEY_INCA) THEN 174 175 IF(var == 'O3') oldVar = 'OX' !--- DEAL WITH INCA OZONE EXCEPTION 175 #endif 176 END IF 176 177 !-------------------------------------------------------------------------------------------------------------------------- 177 178 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 46 46 xios_set_current_context, & 47 47 using_xios 48 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 48 49 49 50 IMPLICIT NONE … … 1519 1520 endif 1520 1521 1521 #ifdef INCA 1522 IF (CPPKEY_INCA) THEN 1522 1523 IF (ANY(type_trac == ['inca','inco'])) THEN 1523 1524 CALL finalize_inca … … 1529 1530 !$OMP END MASTER 1530 1531 ENDIF 1531 #endif 1532 END IF 1532 1533 #ifdef REPROBUS 1533 1534 if (type_trac == 'repr') CALL finalize_reprobus … … 1576 1577 !$OMP END MASTER 1577 1578 1578 #ifdef INCA 1579 IF (CPPKEY_INCA) THEN 1579 1580 IF (ANY(type_trac == ['inca','inco'])) THEN 1580 1581 CALL finalize_inca … … 1586 1587 !$OMP END MASTER 1587 1588 ENDIF 1588 #endif 1589 END IF 1589 1590 #ifdef REPROBUS 1590 1591 if (type_trac == 'repr') CALL finalize_reprobus … … 1749 1750 !$OMP END MASTER 1750 1751 1751 #ifdef INCA 1752 IF (CPPKEY_INCA) THEN 1752 1753 IF (ANY(type_trac == ['inca','inco'])) THEN 1753 1754 CALL finalize_inca … … 1760 1761 ENDIF 1761 1762 1762 #endif 1763 END IF 1763 1764 #ifdef REPROBUS 1764 1765 if (type_trac == 'repr') CALL finalize_reprobus … … 1863 1864 !$OMP END MASTER 1864 1865 1865 #ifdef INCA 1866 IF (CPPKEY_INCA) THEN 1866 1867 IF (ANY(type_trac == ['inca','inco'])) THEN 1867 1868 CALL finalize_inca … … 1874 1875 ENDIF 1875 1876 1876 #endif 1877 END IF 1877 1878 #ifdef REPROBUS 1878 1879 if (type_trac == 'repr') CALL finalize_reprobus -
LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90
r4050 r5251 77 77 USE infotrac, ONLY: nqtot, tracers 78 78 USE filtreg_mod 79 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 79 80 IMPLICIT NONE 80 81 !------------------------------------------------------------------------------- … … 145 146 !******************************************************************************* 146 147 ! Look for ozone tracer: 147 #ifndef INCA 148 IF (.NOT. CPPKEY_INCA) THEN 148 149 DO iq=1,nqtot; IF(strLower(tracers(iq)%name)=="o3") EXIT; END DO 149 150 IF(iq/=nqtot+1) THEN … … 151 152 CALL press_coefoz 152 153 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 154 155 END IF 155 #endif 156 END IF 156 157 q3d(iip1,:,:,:)=q3d(1,:,:,:) 157 158 -
LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90
r4389 r5251 38 38 USE phyaqua_mod, ONLY: iniaqua 39 39 USE comconst_mod, ONLY: omeg, rad 40 #ifdef INCA41 40 USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic 42 41 #ifdef CPP_PARA … … 45 44 #endif 46 45 USE mod_phys_lmdz_omp_data, ONLY: klon_omp 47 #endif48 46 USE ioipsl_getin_p_mod, ONLY: getin_p 49 47 USE slab_heat_transp_mod, ONLY: ini_slab_transp_geom 48 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 50 49 IMPLICIT NONE 51 50 … … 167 166 168 167 IF (ANY(type_trac == ['inca','inco'])) THEN 169 #ifdef INCA 168 IF (CPPKEY_INCA) THEN 170 169 CALL init_inca_dim_reg(nbp_lon, nbp_lat - 1, & 171 170 rlonudyn, rlatudyn, rlonvdyn, rlatvdyn) 172 #endif 171 END IF 173 172 END IF 174 173 -
LMDZ6/trunk/libf/phy_common/mod_phys_lmdz_mpi_data.F90
r4602 r5251 3 3 ! 4 4 MODULE mod_phys_lmdz_mpi_data 5 USE lmdz_mpi 5 USE lmdz_mpi 6 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 6 7 7 8 INTEGER,SAVE :: ii_begin … … 174 175 CALL Print_module_data 175 176 176 #ifdef INCA 177 IF (CPPKEY_INCA) THEN 177 178 CALL init_inca_para( & 178 179 nbp_lon,nbp_lat,nbp_lev,klon_glo,mpi_size, & 179 180 distrib,nvertex, grid_type,communicator) 180 #endif 181 END IF 181 182 182 183 END SUBROUTINE Init_phys_lmdz_mpi_data -
LMDZ6/trunk/libf/phy_common/physics_distribution_mod.F90
r4389 r5251 17 17 USE CHEM_REP, ONLY : Init_chem_rep_phys 18 18 #endif 19 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 19 20 20 21 IMPLICIT NONE … … 33 34 CALL init_dimphy(klon_omp,nbp_lev) 34 35 35 #ifdef INCA 36 IF (CPPKEY_INCA) THEN 36 37 CALL Init_inca_dim(klon_omp,nbp_lev) 37 #endif 38 END IF 38 39 39 40 #ifdef REPROBUS -
LMDZ6/trunk/libf/phylmd/infotrac_phy.F90
r5237 r5251 138 138 USE CHEM_REP, ONLY: Init_chem_rep_trac 139 139 #endif 140 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 140 141 IMPLICIT NONE 141 142 !============================================================================================================================== … … 164 165 ! Local variables 165 166 INTEGER, ALLOCATABLE :: hadv(:), vadv(:) !--- Horizontal/vertical transport scheme number 166 #ifdef INCA167 167 INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA 168 168 vad (:), vadv_inca(:), pbl_flg_inca(:) 169 169 CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:) !--- Tracers names for INCA 170 170 INTEGER :: nqINCA 171 #endif172 171 #ifdef CPP_StratAer 173 172 CHARACTER(LEN=maxlen), ALLOCATABLE :: tnames(:) … … 224 223 SELECT CASE(type_trac) 225 224 CASE('inca', 'inco') 226 #ifndef INCA 225 IF (.NOT. CPPKEY_INCA) THEN 227 226 CALL abort_physic(modname, 'You must add cpp key INCA and compile with INCA code', 1) 228 #endif 227 END IF 229 228 CASE('repr') 230 229 #ifndef REPROBUS … … 267 266 IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac) .AND. lInit) THEN !=== FOUND OLD STYLE INCA "traceur.def" 268 267 !--------------------------------------------------------------------------------------------------------------------------- 269 #ifdef INCA 268 IF (CPPKEY_INCA) THEN 270 269 nqo = SIZE(tracers) - nqCO2 271 270 CALL Init_chem_inca_trac(nqINCA) !--- Get nqINCA from INCA … … 298 297 IF(setGeneration(tracers)) CALL abort_physic(modname,'See below',1) !- SET FIELDS %iGeneration, %gen0Name 299 298 DEALLOCATE(had, hadv_inca, vad, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca) 300 #endif 299 END IF 301 300 !--------------------------------------------------------------------------------------------------------------------------- 302 301 ELSE !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE) … … 306 305 nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O') !--- Passed to phytrac 307 306 nqCO2 = COUNT( [type_trac == 'inco', type_trac == 'co2i'] ) 308 #ifdef INCA 307 IF (CPPKEY_INCA) THEN 309 308 nqINCA = COUNT(tracers(:)%component == 'inca') 310 #endif 309 END IF 311 310 IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "hadv"', 1) 312 311 IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "vadv"', 1) … … 418 417 CALL msg('niso = '//TRIM(int2str(niso)), modname) 419 418 CALL msg('ntiso = '//TRIM(int2str(ntiso)), modname) 420 #ifdef INCA 419 IF (CPPKEY_INCA) THEN 421 420 CALL msg('nqCO2 = '//TRIM(int2str(nqCO2)), modname) 422 421 CALL msg('nqINCA = '//TRIM(int2str(nqINCA)), modname) 423 #endif 422 END IF 424 423 t => tracers 425 424 CALL msg('Information stored in '//TRIM(modname)//': ', modname) -
LMDZ6/trunk/libf/phylmd/oasis.F90
r4754 r5251 122 122 USE mod_phys_lmdz_mpi_data, ONLY: klon_mpi_para_nb 123 123 use config_ocean_skin_m, only: activate_ocean_skin 124 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 124 125 125 126 ! Local variables … … 356 357 357 358 358 #ifdef INCA 359 IF (CPPKEY_INCA) THEN 359 360 IF (dms_cycle_cpl .OR. n2o_cycle_cpl) THEN 360 361 CALL init_inca_oasis(inforecv(idr_ocedms:idr_ocen2o)) 361 362 ENDIF 362 #endif 363 END IF 363 364 364 365 !************************************************************************************ -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r5208 r5251 109 109 110 110 111 #ifdef INCA112 111 USE geometry_mod, ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo 113 112 USE time_phylmdz_mod, ONLY: ndays 114 113 USE infotrac_phy, ONLY: nqCO2 115 #endif116 114 #ifdef REPROBUS 117 115 USE chem_rep, ONLY: Init_chem_rep_xjour, d_q_rep, d_ql_rep, d_qi_rep, & … … 120 118 USE strataer_emiss_mod, ONLY: strataer_emiss_init 121 119 #endif 122 #if defined INCA || defined REPROBUS123 120 USE time_phylmdz_mod, ONLY: annee_ref, day_ini, day_ref, start_time 124 121 USE vertical_layers_mod, ONLY: aps, bps, ap, bp 125 #endif 126 122 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 127 123 128 124 #ifdef CPP_RRTM … … 1243 1239 REAL, dimension(klon, klev) :: cldfrarad ! fraction nuageuse 1244 1240 1245 #ifdef INCA1246 1241 REAL :: calday, zxsnow_dummy(klon) 1247 1242 ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA … … 1250 1245 REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_cginca 1251 1246 REAL, DIMENSION(klon,klev,nbands) :: init_ccminca 1252 #endif1253 1247 REAL, DIMENSION(klon,nbtr) :: init_source 1254 1248 … … 2178 2172 ! 2179 2173 IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 2180 #ifdef INCA 2174 IF (CPPKEY_INCA) THEN 2181 2175 CALL VTe(VTphysiq) 2182 2176 CALL VTb(VTinca) … … 2188 2182 config_inca) 2189 2183 2190 CALL init_inca_geometry( & 2184 CALL init_inca_geometry( & 2191 2185 longitude, latitude, & 2192 2186 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 2196 2190 CALL chemini( pplay, & 2197 2191 nbp_lon, nbp_lat, & … … 2251 2245 piz_aero(:,:,:,:) = init_pizinca 2252 2246 cg_aero(:,:,:,:) = init_cginca 2253 ! 2247 ! 2254 2248 2255 2249 2256 2250 CALL VTe(VTinca) 2257 2251 CALL VTb(VTphysiq) 2258 #endif 2252 END IF 2259 2253 ENDIF 2260 2254 ! … … 4181 4175 4182 4176 IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 4183 #ifdef INCA 4177 IF (CPPKEY_INCA) THEN 4184 4178 CALL VTe(VTphysiq) 4185 4179 CALL VTb(VTinca) … … 4234 4228 CALL VTe(VTinca) 4235 4229 CALL VTb(VTphysiq) 4236 #endif 4230 END IF 4237 4231 ENDIF !type_trac = inca or inco 4238 4232 IF (type_trac == 'repr') THEN … … 4525 4519 4526 4520 IF (aerosol_couple.AND.config_inca=='aero') THEN 4527 #ifdef INCA 4521 IF (CPPKEY_INCA) THEN 4528 4522 CALL radlwsw_inca & 4529 4523 (chemistry_couple, kdlon,kflev,dist, rmu0, fract, solaire, & … … 4545 4539 cldtaupirad, & 4546 4540 topswai_aero, solswai_aero) 4547 #endif 4541 END IF 4548 4542 ELSE 4549 4543 ! … … 5530 5524 ! 5531 5525 IF (ANY(type_trac == ['inca','inco'])) THEN 5532 #ifdef INCA 5526 IF (CPPKEY_INCA) THEN 5533 5527 CALL VTe(VTphysiq) 5534 5528 CALL VTb(VTinca) … … 5550 5544 CALL VTe(VTinca) 5551 5545 CALL VTb(VTphysiq) 5552 #endif 5546 END IF 5553 5547 ENDIF 5554 5548 … … 5808 5802 IF (using_xios) THEN 5809 5803 5810 #ifdef INCA 5811 IF (type_trac == 'inca') THEN 5812 IF (is_omp_master .AND. grid_type==unstructured) THEN 5804 IF (CPPKEY_INCA) THEN 5805 IF (type_trac == 'inca') THEN 5806 IF (is_omp_master .AND. grid_type==unstructured) THEN 5813 5807 CALL finalize_inca 5814 5808 ENDIF 5815 5809 ENDIF 5816 #endif 5810 END IF 5817 5811 5818 5812 ! close xios physiq context (call LMDZ) -
LMDZ6/trunk/libf/phylmd/tracinca_mod.F90
r4500 r5251 62 62 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 63 63 USE aero_mod, ONLY : naero_grp 64 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 64 65 IMPLICIT NONE 65 66 … … 145 146 END DO 146 147 147 #ifdef INCA 148 IF (config_inca == 'aero' .OR. config_inca == 'chem') THEN 149 zpmfu(:,:)=pmfu(:,:) 148 IF (CPPKEY_INCA) THEN 149 IF (config_inca == 'aero' .OR. config_inca == 'chem') THEN 150 zpmfu(:,:)=pmfu(:,:) 150 151 ELSE IF (config_inca == 'aeNP') THEN 151 152 zpmfu(:,:)=upwd(:,:) … … 156 157 pplay,pdel,prfl,pmflxr,psfl, & 157 158 pmflxs,zpmfu,itop_con,ibas_con, & 158 pphi,cell_area,nstep,rneb,t_seri, & 159 pphi,cell_area,nstep,rneb,t_seri, & 159 160 rh,tau_aero,piz_aero,cg_aero, & 160 161 rfname,ccm,lafin) 161 #endif 162 163 164 #ifdef INCA 162 END IF 163 164 165 IF (CPPKEY_INCA) THEN 165 166 CALL chemmain (tr_seri, & !mmr 166 167 nstep, & !nstep … … 198 199 nbp_lat, & !ny 199 200 source ) 200 #endif 201 END IF 201 202 202 203 CALL VTe(VTinca) -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r5244 r5251 109 109 110 110 111 #ifdef INCA112 111 USE geometry_mod, ONLY: longitude, latitude, boundslon, boundslat, ind_cell_glo 113 112 USE time_phylmdz_mod, ONLY: ndays 114 113 USE infotrac_phy, ONLY: nqCO2 115 #endif116 114 #ifdef REPROBUS 117 115 USE chem_rep, ONLY: Init_chem_rep_xjour, d_q_rep, d_ql_rep, d_qi_rep, & … … 120 118 USE strataer_emiss_mod, ONLY: strataer_emiss_init 121 119 #endif 122 #if defined INCA || defined REPROBUS123 120 USE time_phylmdz_mod, ONLY: annee_ref, day_ini, day_ref, start_time 124 121 USE vertical_layers_mod, ONLY: aps, bps, ap, bp 125 #endif126 122 127 123 … … 441 437 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra 442 438 USE output_physiqex_mod, ONLY: output_physiqex 439 USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_INCA 443 440 444 441 … … 1349 1346 REAL, dimension(klon, klev) :: cldfrarad ! fraction nuageuse 1350 1347 1351 #ifdef INCA1352 1348 REAL :: calday, zxsnow_dummy(klon) 1353 1349 ! set de variables utilisees pour l'initialisation des valeurs provenant de INCA … … 1356 1352 REAL, DIMENSION(klon,klev,naero_grp,nbands) :: init_cginca 1357 1353 REAL, DIMENSION(klon,klev,nbands) :: init_ccminca 1358 #endif1359 1354 REAL, DIMENSION(klon,nbtr) :: init_source 1360 1355 … … 2340 2335 ! 2341 2336 IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 2342 #ifdef INCA 2337 IF (CPPKEY_INCA) THEN 2343 2338 CALL VTe(VTphysiq) 2344 2339 CALL VTb(VTinca) … … 2350 2345 config_inca) 2351 2346 2352 CALL init_inca_geometry( & 2347 CALL init_inca_geometry( & 2353 2348 longitude, latitude, & 2354 2349 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 2358 2353 CALL chemini( pplay, & 2359 2354 nbp_lon, nbp_lat, & … … 2413 2408 piz_aero(:,:,:,:) = init_pizinca 2414 2409 cg_aero(:,:,:,:) = init_cginca 2415 ! 2410 ! 2416 2411 2417 2412 2418 2413 CALL VTe(VTinca) 2419 2414 CALL VTb(VTphysiq) 2420 #endif 2415 END IF 2421 2416 ENDIF 2422 2417 ! … … 5635 5630 5636 5631 IF (ANY(type_trac == ['inca','inco'])) THEN ! ModThL 5637 #ifdef INCA 5632 IF (CPPKEY_INCA) THEN 5638 5633 CALL VTe(VTphysiq) 5639 5634 CALL VTb(VTinca) … … 5688 5683 CALL VTe(VTinca) 5689 5684 CALL VTb(VTphysiq) 5690 #endif 5685 END IF 5691 5686 ENDIF !type_trac = inca or inco 5692 5687 IF (type_trac == 'repr') THEN … … 5978 5973 5979 5974 IF (aerosol_couple.AND.config_inca=='aero') THEN 5980 #ifdef INCA 5975 IF (CPPKEY_INCA) THEN 5981 5976 CALL radlwsw_inca & 5982 5977 (chemistry_couple, kdlon,kflev,dist, rmu0, fract, solaire, & … … 5998 5993 cldtaupirad, & 5999 5994 topswai_aero, solswai_aero) 6000 #endif 5995 END IF 6001 5996 ELSE 6002 5997 ! … … 7109 7104 ! 7110 7105 IF (ANY(type_trac == ['inca','inco'])) THEN 7111 #ifdef INCA 7106 IF (CPPKEY_INCA) THEN 7112 7107 CALL VTe(VTphysiq) 7113 7108 CALL VTb(VTinca) … … 7129 7124 CALL VTe(VTinca) 7130 7125 CALL VTb(VTphysiq) 7131 #endif 7126 END IF 7132 7127 ENDIF 7133 7128 … … 7440 7435 IF (using_xios) THEN 7441 7436 7442 #ifdef INCA 7443 IF (type_trac == 'inca') THEN 7444 IF (is_omp_master .AND. grid_type==unstructured) THEN 7437 IF (CPPKEY_INCA) THEN 7438 IF (type_trac == 'inca') THEN 7439 IF (is_omp_master .AND. grid_type==unstructured) THEN 7445 7440 CALL finalize_inca 7446 7441 ENDIF 7447 7442 ENDIF 7448 #endif 7443 END IF 7449 7444 7450 7445 ! close xios physiq context (call LMDZ)
Note: See TracChangeset
for help on using the changeset viewer.