Changeset 5251 for LMDZ6/trunk/libf/phylmd/infotrac_phy.F90
- Timestamp:
- Oct 22, 2024, 12:31:08 PM (5 days ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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)
Note: See TracChangeset
for help on using the changeset viewer.