Changeset 3961


Ignore:
Timestamp:
Jul 13, 2021, 12:42:55 PM (3 years ago)
Author:
acozic
Message:

hadv and vadv weren't allocated in inca case
Add modname to help in debug (error message)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/LMDZ-tracers/libf/dyn3d_common/infotrac.F90

    r3960 r3961  
    269269
    270270  IF(readTracersFiles(type_trac, fType, tracers)) CALL abort_gcm(modname,'problem with tracers file(s)',1)
     271  modname = 'infotrac_init'
    271272  CALL msg(fType == 0, 'WARNING: USING DEFAULT VALUES !')
    272273
     
    310311        nqtrue = nbtr + nqo                                          !--- Total number of tracers
    311312        ALLOCATE(ttr(nqtrue)); ttr(1:nqo) = tracers(1:nqo)
     313        ALLOCATE(hadv(nqtrue)); hadv(:) = 0
     314        ALLOCATE(vadv(nqtrue)); vadv(:) = 0
    312315        DO iq = nqo+1, nqtrue
    313316          ttr(iq)%name = solsym(iq)
    314317          ttr(iq)%prnt = tran0
    315318          ttr(iq)%igen = 1
    316           hadv = hadv_inca(iq-nqo)
    317           vadv = vadv_inca(iq-nqo)
     319
     320          hadv(iq) = hadv_inca(iq-nqo)
     321          vadv(iq) = vadv_inca(iq-nqo)
     322
    318323        END DO
    319324        CALL MOVE_ALLOC(FROM=ttr, TO=tracers)
Note: See TracChangeset for help on using the changeset viewer.