Changeset 4359 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Nov 30, 2022, 6:23:36 PM (22 months ago)

Author:

dcugnet

Message:

Fix for a gfortran specific byg. In the following expression:

nam = [(TRIM(name(i))TRIM(int2str(nsrf,2)), i=1, SIZE(name))]

where nam(:) is an ALLOCATABLE string, the gfortran compiled executable crashes, but should not:

"Fortran runtime error: Index '0' of dimension 1 of array 'name' below lower bound of 1"

File:

• LMDZ6/trunk/libf/phylmdiso/phyetat0_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmdiso/phyetat0_mod.F90

@@ -766 +766 @@
 !-----------------------------------------------------------------------------
   INTEGER :: nsrf, i
-  CHARACTER(LEN=maxlen), ALLOCATABLE :: nam(:)
-  CHARACTER(LEN=maxlen) :: tname, des
+  CHARACTER(LEN=maxlen) :: nam(SIZE(name)), tname, des
   IF(SIZE(field,3)>99) CALL abort_physic("phyetat0", "Too much sub-cells", 1)
   DO nsrf = 1, SIZE(field,3)
-    nam = [(TRIM(name(i))//TRIM(int2str(nsrf,2)), i=1, SIZE(name))]
+    DO i = 1, SIZE(name); nam(i) = TRIM(name(i))//TRIM(int2str(nsrf,2)); END DO
     des = TRIM(descr)//" srf:"//int2str(nsrf,2)
     lFound = phyetat0_get21(field(:,:,nsrf), nam, TRIM(des), default, tname)