Changeset 22 for readTracFiles_mod.f90
- Timestamp:
- Sep 9, 2022, 5:20:56 PM (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
readTracFiles_mod.f90
r21 r22 175 175 END SELECT 176 176 !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 177 177 178 178 179 IF(ALL([2,3] /= fType)) RETURN … … 613 614 !------------------------------------------------------------------------------------------------------------------------------ 614 615 TYPE(trac_type), ALLOCATABLE, INTENT(INOUT) :: tr(:) !--- Tracer derived type vector 615 INTEGER, ALLOCATABLE :: iy(:), iz(:) 616 INTEGER :: ig, ng, iq, jq, ip, nq, n, ix(SIZE(tr)), k 616 !------------------------------------------------------------------------------------------------------------------------------ 617 TYPE(trac_type), ALLOCATABLE :: tr2(:) 618 INTEGER, ALLOCATABLE :: iy(:), iz(:) 619 INTEGER :: ig, ng, iq, jq, ip, nq, n, ix(SIZE(tr)), k 620 INTEGER :: it 621 ! tr2 is introduced in order to cope with a bug in gfortran 4.8.5 compiler 617 622 !------------------------------------------------------------------------------------------------------------------------------ 618 623 nq = SIZE(tr) … … 620 625 iq = strIdx(tr(:)%name, addPhase('H2O', ip)) 621 626 IF(iq == 0) CYCLE 622 tr = [tr(iq), tr(1:iq-1), tr(iq+1:nq)] 627 tr2 = tr(:) 628 tr = [tr2(iq), tr2(1:iq-1), tr2(iq+1:nq)] 623 629 END DO 624 630 IF(lSortByGen) THEN
Note: See TracChangeset
for help on using the changeset viewer.