Changeset 4233


Ignore:
Timestamp:
Aug 31, 2022, 6:40:27 PM (21 months ago)
Author:
fhourdin
Message:

Commission pour Jean-Yves

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/misc/readTracFiles_mod.f90

    r4194 r4233  
    175175  END SELECT
    176176  !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
     177
    177178
    178179  IF(ALL([2,3] /= fType)) RETURN
     
    613614!------------------------------------------------------------------------------------------------------------------------------
    614615  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
    617622!------------------------------------------------------------------------------------------------------------------------------
    618623  nq = SIZE(tr)
     
    620625    iq = strIdx(tr(:)%name, addPhase('H2O', ip))
    621626    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)]
    623629  END DO
    624630  IF(lSortByGen) THEN
Note: See TracChangeset for help on using the changeset viewer.