Changeset 5234 for LMDZ6/trunk
- Timestamp:
- Sep 30, 2024, 5:55:57 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/infotrac_phy.F90
r5216 r5234 371 371 CALL abort_physic(modname, 'problem with the computation of nqtottr', 1) 372 372 373 #ifdef CPP_StratAer 374 IF (type_trac == 'coag') THEN 375 nbtr_bin = COUNT([(tracers(iq)%name(1:3)=='BIN', iq=1, nqtot)]) 376 nbtr_sulgas = COUNT([(tracers(iq)%name(1:3)=='GAS', iq=1, nqtot)]) 377 tnames = PACK(tracers(:)%name, MASK=tracers(:)%isInPhysics) 378 id_BIN01_strat = strIdx(tnames, 'BIN01' ) 379 id_OCS_strat = strIdx(tnames, 'GASOCS' ) 380 id_SO2_strat = strIdx(tnames, 'GASSO2' ) 381 id_H2SO4_strat = strIdx(tnames, 'GASH2SO4') 382 id_TEST_strat = strIdx(tnames, 'GASTEST' ) 383 END IF 384 #endif 385 373 386 !=== DISPLAY THE RESULTS 374 387 IF(.NOT.is_master) RETURN … … 405 418 #ifdef CPP_StratAer 406 419 IF (type_trac == 'coag') THEN 407 nbtr_bin = COUNT([(tracers(iq)%name(1:3)=='BIN', iq=1, nqtot)])408 nbtr_sulgas = COUNT([(tracers(iq)%name(1:3)=='GAS', iq=1, nqtot)])409 tnames = PACK(tracers(:)%name, MASK=tracers(:)%isInPhysics)410 id_BIN01_strat = strIdx(tnames, 'BIN01' )411 id_OCS_strat = strIdx(tnames, 'GASOCS' )412 id_SO2_strat = strIdx(tnames, 'GASSO2' )413 id_H2SO4_strat = strIdx(tnames, 'GASH2SO4')414 id_TEST_strat = strIdx(tnames, 'GASTEST' )415 420 CALL msg('nbtr_bin ='//TRIM(int2str(nbtr_bin )), modname) 416 421 CALL msg('nbtr_sulgas ='//TRIM(int2str(nbtr_sulgas )), modname)
Note: See TracChangeset
for help on using the changeset viewer.