Most of the changes are intended to help to eventually remove the constraints about the tracers assumptions, in particular water tracers.
IF(tracers(iq)%isH2Ofamily) CYCLE
and DO it=1,nqtottr; iq=itr_indice(it) loops are replaced with:
it = 0 DO iq = 1, nqtot IF(.NOT.tracers(iq)%isAdvected .OR. tracers(iq)%isH2Ofamily) CYCLE it = it+1
it = 0 DO iq = 1, nqtot
IF(.NOT.tracers(iq)%isAdvected .OR. tracers(iq)%isH2Ofamily) CYCLE it = it+1
DO iq=1,nqpere loops are replaced with: DO iq=1,nqtot IF(tracers(iq)%parent/='air') CYCLE
DO iq=1,nqpere loops are replaced with: DO iq=1,nqtot
IF(tracers(iq)%parent/='air') CYCLE