Changeset 13009 for NEMO/trunk/src/OCE/LBC

Timestamp:

2020-06-03T08:35:19+02:00 (4 years ago)

Author:

smasson

Message:

trunk: bugfix, introduced in [12933], when jpnij = 1, see 2477

File:

• NEMO/trunk/src/OCE/LBC/lib_mpp.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/LBC/lib_mpp.F90

@@ -1126 +1126 @@
             CALL ctl_opn( inum,'ocean.output', 'APPEND', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
             WRITE(inum,*)
-            idg = MAX( INT(LOG10(REAL(jpnij-1,wp))) + 1, 4 )        ! how many digits to we need to write? min=4, max=9
-            WRITE(clfmt, "('(a,i', i1, '.', i1, ')')") idg, idg     ! '(a,ix.x)'
+            idg = MAX( INT(LOG10(REAL(MAX(1,jpnij-1),wp))) + 1, 4 )        ! how many digits to we need to write? min=4, max=9
+            WRITE(clfmt, "('(a,i', i1, '.', i1, ')')") idg, idg            ! '(a,ix.x)'
             WRITE(inum,clfmt) ' ===>>> : see E R R O R in ocean.output_', narea - 1
          ENDIF
@@ -1232 +1232 @@
       IF( PRESENT( karea ) ) THEN
          IF( karea > 1 ) THEN
-            idg = MAX( INT(LOG10(REAL(jpnij-1,wp))) + 1, 4 )        ! how many digits to we need to write? min=4, max=9
-            WRITE(clfmt, "('(a,a,i', i1, '.', i1, ')')") idg, idg   ! '(a,a,ix.x)'
+            ! Warning: jpnij is maybe not already defined when calling ctl_opn -> use mppsize instead of jpnij
+            idg = MAX( INT(LOG10(REAL(MAX(1,mppsize-1),wp))) + 1, 4 )      ! how many digits to we need to write? min=4, max=9
+            WRITE(clfmt, "('(a,a,i', i1, '.', i1, ')')") idg, idg          ! '(a,a,ix.x)'
             WRITE(clfile, clfmt) TRIM(clfile), '_', karea-1
          ENDIF