Changeset 3961 for LMDZ6/branches/LMDZ-tracers

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:

• LMDZ6/branches/LMDZ-tracers/libf/dyn3d_common/infotrac.F90 (modified) (2 diffs)

LMDZ6/branches/LMDZ-tracers/libf/dyn3d_common/infotrac.F90

@@ -269 +269 @@
 
   IF(readTracersFiles(type_trac, fType, tracers)) CALL abort_gcm(modname,'problem with tracers file(s)',1)
+  modname = 'infotrac_init'
   CALL msg(fType == 0, 'WARNING: USING DEFAULT VALUES !')
 
@@ -310 +311 @@
         nqtrue = nbtr + nqo                                          !--- Total number of tracers
         ALLOCATE(ttr(nqtrue)); ttr(1:nqo) = tracers(1:nqo)
+        ALLOCATE(hadv(nqtrue)); hadv(:) = 0
+        ALLOCATE(vadv(nqtrue)); vadv(:) = 0
         DO iq = nqo+1, nqtrue
           ttr(iq)%name = solsym(iq)
           ttr(iq)%prnt = tran0
           ttr(iq)%igen = 1
-          hadv = hadv_inca(iq-nqo)
-          vadv = vadv_inca(iq-nqo)
+
+          hadv(iq) = hadv_inca(iq-nqo)
+          vadv(iq) = vadv_inca(iq-nqo)
+
         END DO
         CALL MOVE_ALLOC(FROM=ttr, TO=tracers)