Changeset 69 for trunk/NEMO/OPA_SRC/mppini_2.h90

Timestamp:

2004-04-22T14:17:26+02:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX043 : Change the bathymetry file name "bathymetry" to "bathy_level"

File:

• trunk/NEMO/OPA_SRC/mppini_2.h90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/mppini_2.h90

@@ -40 +40 @@
       !!----------------------------------------------------------------------
       !! Local variables
-      CHARACTER (len=25) ::   clexp           ! temporary name
+      CHARACTER (len=25) ::   clexp,       &  ! temporary name
+         &                    clname          !    "        "
+      LOGICAL ::   llbon                      ! check the existence of bathy files
       INTEGER :: ji, jj, jn, jproc, jarea     ! dummy loop indices
       INTEGER ::   inum = 11                  ! temporary logical unit
@@ -98 +100 @@
 
       ! open the file
+      clname = 'bathy_level'
       IF(lwp) WRITE(numout,*) '           bathymetry file in levels'
-      CALL ctlopn( inum, 'bathymetry', 'OLD', 'FORMATTED', 'SEQUENTIAL',   &
-                   1      , numout      , lwp  , 1                          )
-
-      ! read bathymetry file
-      REWIND inum
-      READ(inum,9001) clexp, iim, ijm
-      READ(inum,'(/)')
-      ifreq=40
-      il1=1
-      DO jn=1,jpidta/ifreq+1
+      INQUIRE( FILE=clname, EXIST=llbon )
+      IF( llbon ) THEN
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) '         read level bathymetry in ', clname
+         IF(lwp) WRITE(numout,*)
+         OPEN( UNIT=inum, FILE=clname, FORM='FORMATTED', ACCESS='SEQUENTIAL', STATUS='OLD' )
+
+         ! read bathymetry file
+         REWIND inum
+         READ(inum,9001) clexp, iim, ijm
          READ(inum,'(/)')
-         il2 = MIN( jpidta, il1+ifreq-1 )
-         READ(inum,9002) ( ii, ji = il1, il2, 5 )
-         READ(inum,'(/)')
-         DO jj = jpjdta, 1, -1
-            READ(inum,9003) ij, ( idata(ji,jj), ji = il1, il2 )
-         END DO
-         il1 = il1 + ifreq
-      END DO
-      CLOSE(inum)
-
-9001  FORMAT(1x,a15,2i8)
-9002  FORMAT(3x,13(i3,12x))
-9003  FORMAT(i3,41i3)
+         ifreq=40
+         il1=1
+         DO jn=1,jpidta/ifreq+1
+            READ(inum,'(/)')
+            il2 = MIN( jpidta, il1+ifreq-1 )
+            READ(inum,9002) ( ii, ji = il1, il2, 5 )
+            READ(inum,'(/)')
+            DO jj = jpjdta, 1, -1
+               READ(inum,9003) ij, ( idata(ji,jj), ji = il1, il2 )
+            END DO
+            il1 = il1 + ifreq
+         END DO
+         CLOSE(inum)
+
+9001     FORMAT(1x,a15,2i8)
+9002     FORMAT(3x,13(i3,12x))
+9003     FORMAT(i3,41i3)
+
+      ELSE
+         IF(lwp) WRITE(numout,cform_err)
+         IF(lwp) WRITE(numout,*)'    mppini_2 : unable to read the file', clname
+         nstop = nstop + 1
+      ENDIF
 
       ! land/sea mask over the global/zoom domain