Changeset 128 for trunk/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2004-07-07T14:38:57+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE080 : Now the previous bathymetry level file (bathy_level) which was in ASCII format must be converted into NetCDF format (bathy_level.nc) to be read

File:

• trunk/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -250 +250 @@
       !!      ntopo=-1 :   rectangular channel or bassin with a bump 
       !!      ntopo= 0 :   flat rectangular channel or basin 
-      !!      ntopo= 1 :   mbathy is read in 'bathy_level' file
+      !!      ntopo= 1 :   mbathy is read in 'bathy_level.nc' NetCDF file
       !!                   bathy  is read in 'bathy_meter.nc' NetCDF file
       !!      C A U T I O N : mbathy will be modified during the initializa-
@@ -372 +372 @@
          !                                            ! =============== !
 
-         clname = 'bathy_level'                          ! level bathymetry
+         clname = 'bathy_level.nc'                       ! Level bathymetry
          INQUIRE( FILE=clname, EXIST=llbon )
          IF( llbon ) THEN
@@ -378 +378 @@
             IF(lwp) WRITE(numout,*) '         read level bathymetry in ', clname
             IF(lwp) WRITE(numout,*)
-            OPEN( UNIT=inum, FILE=clname, FORM='FORMATTED', ACCESS='SEQUENTIAL', STATUS='OLD' )
-
-         REWIND(inum)
-         READ  (inum,9101) clexp, iim, ijm
-         READ  (inum,'(/)')
-         ifreq = 40
-         il1   = 1
-         DO jn = 1, jpidta/ifreq+1
-            READ(inum,'(/)')
-            il2 = MIN( jpidta, il1+ifreq-1 )
-            READ(inum,9201) ( ii, ji = il1, il2, 5 )
-            READ(inum,'(/)')
-            DO jj = jpjdta, 1, -1
-               READ(inum,9202) ij, ( idta(ji,jj), ji = il1, il2 )
-            END DO
-            il1 = il1 + ifreq
-         END DO
-         CLOSE(inum)
-
-9101     FORMAT(1x,a15,2i8)
-9201     FORMAT(3x,13(i3,12x))
-9202     FORMAT(i3,41i3)
+            itime = 1
+            ipi = jpidta
+            ipj = jpjdta
+            ipk = 1
+            zdt = rdt
+            CALL flinopen( clname, 1, jpidta, 1, jpjdta, .FALSE.,   &
+                           ipi, ipj, ipk, zlamt, zphit, zdept, itime, istep, zdate0, zdt, inum )
+            CALL flinget( inum, 'Bathy_level', jpidta, jpjdta, 1,   &
+                          itime, 1, 1, 1, jpidta, 1, jpjdta, zdta(:,:) )
+            idta(:,:) = zdta(:,:)
+            CALL flinclo( inum )
 
          ELSE