Changeset 4393
- Timestamp:
- Jan 23, 2023, 5:39:48 PM (22 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/traclmdz_mod.F90
r4389 r4393 96 96 USE indice_sol_mod 97 97 USE print_control_mod, ONLY: lunout 98 USE strings_mod, ONLY: strLower99 98 100 99 ! Input variables … … 177 176 IF(.NOT.(tracers(iq)%isInPhysics)) CYCLE 178 177 it = it+1 179 SELECT CASE( strLower(tracers(iq)%name))180 CASE(" rn");id_rn = it ! radon181 CASE(" pb");id_pb = it ! plomb182 CASE(" aga");id_aga = it ! Age of stratospheric air183 CASE(" be","be7");id_be = it ! Recherche du Beryllium 7184 CASE("o3" );id_o3 = it ! Recherche de l'ozone185 CASE("pcsat" );id_pcsat = it186 CASE("pcocsat" );id_pcocsat= it187 CASE("pcq" );id_pcq = it188 CASE("pcs0" );id_pcs0 = it189 CASE("pcos0" );id_pcos0 = it190 CASE("pcq0" );id_pcq0 = it178 SELECT CASE(tracers(iq)%name) 179 CASE("RN"); id_rn = it ! radon 180 CASE("PB"); id_pb = it ! plomb 181 CASE("Aga","AGA"); id_aga = it ! Age of stratospheric air 182 CASE("Be","BE","Be7","BE7"); id_be = it ! Recherche du Beryllium 7 183 CASE("o3","O3"); id_o3 = it ! Recherche de l'ozone 184 CASE("pcsat", "Pcsat"); id_pcsat = it 185 CASE("pcocsat","Pcocsat"); id_pcocsat= it 186 CASE("pcq", "Pcq"); id_pcq = it 187 CASE("pcs0", "Pcs0"); id_pcs0 = it 188 CASE("pcos0", "Pcos0"); id_pcos0 = it 189 CASE("pcq0", "Pcq0"); id_pcq0 = it 191 190 CASE DEFAULT 192 191 WRITE(lunout,*) 'This is an unknown tracer in LMDZ : ', trim(tracers(iq)%name) 193 192 END SELECT 194 193 195 SELECT CASE( strLower(tracers(iq)%name))196 CASE(" pb") !--- RomP >>> profil initial de PB210194 SELECT CASE(tracers(iq)%name) 195 CASE("PB") !--- RomP >>> profil initial de PB210 197 196 OPEN(ilesfil2,file='prof.pb210',status='old',iostat=irr2) 198 197 IF(irr2 == 0) THEN … … 205 204 WRITE(lunout,*)'prof.pb210 does not exist: use restart values' 206 205 END IF 207 CASE(" aga")206 CASE("Aga","AGA") 208 207 radio(id_aga) = .FALSE. 209 208 aerosol(id_aga) = .FALSE. … … 217 216 lev_1p5km=klev/2 218 217 END IF 219 CASE(" be","be7")218 CASE("Be","BE","Be7","BE7") 220 219 ALLOCATE( srcbe(klon,klev) ) 221 220 radio(id_be) = .TRUE. … … 234 233 WRITE(lunout,*)'Prof. Be7 does not exist: use restart values' 235 234 END IF 236 CASE("o3" )!--- Parametrisation par la chimie de Cariolle235 CASE("o3","O3") !--- Parametrisation par la chimie de Cariolle 237 236 CALL alloc_coefoz !--- Allocate ozone coefficients 238 237 CALL press_coefoz !--- Read input pressure levels 239 CASE("pcs0"," pcos0","pcq0")238 CASE("pcs0","Pcs0", "pcos0","Pcos0", "pcq0","Pcq0") 240 239 conv_flg(it)=0 !--- No transport by convection for this tracer 241 240 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.