Changeset 2051 for trunk/LMDZ.GENERIC/libf

Timestamp:

Dec 13, 2018, 2:13:43 PM (6 years ago)

Author:

jleconte

Message:

12/12/2018 == JL

• Correct a bug from commit 2032 in callcorrk: automatic arrays must be allocated

only if they haven't been before (callcorrk is called twice in physiq when cloud cover is not uniform).

File:

• trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -248 +248 @@
          ! now that L_NGAUSS has been initialized (by sugas_corrk)
          ! allocate related arrays
-         ALLOCATE(dtaui(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
-         ALLOCATE(dtauv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
-         ALLOCATE(cosbv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
-         ALLOCATE(cosbi(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
-         ALLOCATE(wbari(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
-         ALLOCATE(wbarv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
-         ALLOCATE(tauv(L_NLEVRAD,L_NSPECTV,L_NGAUSS))
-         ALLOCATE(taucumv(L_LEVELS,L_NSPECTV,L_NGAUSS))
-         ALLOCATE(taucumi(L_LEVELS,L_NSPECTI,L_NGAUSS))
-         ALLOCATE(taugsurf(L_NSPECTV,L_NGAUSS-1))
-         ALLOCATE(taugsurfi(L_NSPECTI,L_NGAUSS-1))
+         if(.not.allocated(dtaui)) ALLOCATE(dtaui(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
+         if(.not.allocated(dtauv)) ALLOCATE(dtauv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
+         if(.not.allocated(cosbv)) ALLOCATE(cosbv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
+         if(.not.allocated(cosbi)) ALLOCATE(cosbi(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
+         if(.not.allocated(wbari)) ALLOCATE(wbari(L_NLAYRAD,L_NSPECTI,L_NGAUSS))
+         if(.not.allocated(wbarv)) ALLOCATE(wbarv(L_NLAYRAD,L_NSPECTV,L_NGAUSS))
+         if(.not.allocated(tauv)) ALLOCATE(tauv(L_NLEVRAD,L_NSPECTV,L_NGAUSS))
+         if(.not.allocated(taucumv)) ALLOCATE(taucumv(L_LEVELS,L_NSPECTV,L_NGAUSS))
+         if(.not.allocated(taucumi)) ALLOCATE(taucumi(L_LEVELS,L_NSPECTI,L_NGAUSS))
+         if(.not.allocated(taugsurf)) ALLOCATE(taugsurf(L_NSPECTV,L_NGAUSS-1))
+         if(.not.allocated(taugsurfi)) ALLOCATE(taugsurfi(L_NSPECTI,L_NGAUSS-1))
 
          if((igcm_h2o_vap.eq.0) .and. varactive)then