Changeset 1488 for trunk/NEMO/OPA_SRC/IOM/iom_nf90.F90

Timestamp:

2009-07-16T11:19:46+02:00 (15 years ago)

Author:

smasson

Message:

Addition of clobber and chunksize when opening NetCDF files, see ticket:374

File:

• trunk/NEMO/OPA_SRC/IOM/iom_nf90.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/IOM/iom_nf90.F90

@@ -49 +49 @@
       !! ** Purpose : open an input file with NF90
       !!---------------------------------------------------------------------
-      CHARACTER(len=*)       , INTENT(inout)           ::   cdname    ! File name
-      INTEGER                , INTENT(  out)           ::   kiomid    ! nf90 identifier of the opened file
-      LOGICAL                , INTENT(in   )           ::   ldwrt     ! read or write the file?
-      LOGICAL                , INTENT(in   )           ::   ldok      ! check the existence 
-      INTEGER, DIMENSION(2,5), INTENT(in   ), OPTIONAL ::   kdompar   ! domain parameters: 
-
-      CHARACTER(LEN=100) ::   clinfo   ! info character
-      CHARACTER(LEN=100) ::   cltmp    ! temporary character
-      INTEGER            ::   iln      ! lengths of character
-      INTEGER            ::   istop    ! temporary storage of nstop
-      INTEGER            ::   if90id   ! nf90 identifier of the opened file
-      INTEGER            ::   idmy     ! dummy variable
-      INTEGER            ::   jl       ! loop variable
+      CHARACTER(len=*)       , INTENT(inout)           ::   cdname      ! File name
+      INTEGER                , INTENT(  out)           ::   kiomid      ! nf90 identifier of the opened file
+      LOGICAL                , INTENT(in   )           ::   ldwrt       ! read or write the file?
+      LOGICAL                , INTENT(in   )           ::   ldok        ! check the existence 
+      INTEGER, DIMENSION(2,5), INTENT(in   ), OPTIONAL ::   kdompar     ! domain parameters: 
+
+      CHARACTER(LEN=100) ::   clinfo           ! info character
+      CHARACTER(LEN=100) ::   cltmp            ! temporary character
+      INTEGER            ::   iln              ! lengths of character
+      INTEGER            ::   istop            ! temporary storage of nstop
+      INTEGER            ::   if90id           ! nf90 identifier of the opened file
+      INTEGER            ::   idmy             ! dummy variable
+      INTEGER            ::   jl               ! loop variable
+      INTEGER            ::   ichunk           ! temporary storage of nn_chunksz
+      INTEGER            ::   imode            ! creation mode flag: NF90_CLOBBER or NF90_NOCLOBBER
+      LOGICAL            ::   llclobber        ! local definition of ln_clobber
       !---------------------------------------------------------------------
 
       clinfo = '                    iom_nf90_open ~~~  '
       istop = nstop   ! store the actual value of nstop
-      !
-      IF( ldok ) THEN      ! Open existing file...
+      IF( nn_chunksz > 0 ) THEN   ;   ichunk = nn_chunksz
+      ELSE                        ;   ichunk = NF90_SIZEHINT_DEFAULT
+      ENDIF
+      !
+      llclobber = ldwrt .AND. ln_clobber
+      IF( ldok .AND. .NOT. llclobber ) THEN      ! Open existing file...
          !                 ! =============
          IF( ldwrt ) THEN  ! ... in write mode
             IF(lwp) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in WRITE mode'
-            CALL iom_nf90_check(NF90_OPEN(TRIM(cdname), NF90_WRITE  , if90id), clinfo)
-            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy    ), clinfo)
+            CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_WRITE  , if90id, chunksize = ichunk ), clinfo)
+            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy                          ), clinfo)
          ELSE              ! ... in read mode
             IF(lwp) WRITE(numout,*) TRIM(clinfo)//' open existing file: '//TRIM(cdname)//' in READ mode'
-            CALL iom_nf90_check(NF90_OPEN(TRIM(cdname), NF90_NOWRITE, if90id), clinfo)
-         ENDIF
-      ELSE                 ! the file does not exist
+            CALL iom_nf90_check(NF90_OPEN( TRIM(cdname), NF90_NOWRITE, if90id, chunksize = ichunk ), clinfo)
+         ENDIF
+      ELSE                                       ! the file does not exist (or we overwrite it)
          !                 ! =============
          iln = INDEX( cdname, '.nc' )
@@ -86 +93 @@
             ENDIF
             IF(lwp) WRITE(numout,*) TRIM(clinfo)//' create new file: '//TRIM(cdname)//' in WRITE mode'
-            CALL iom_nf90_check(NF90_CREATE( TRIM(cdname), NF90_NOCLOBBER, if90id ), clinfo)
-            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy          ), clinfo)
+
+            IF( llclobber ) THEN   ;   imode = NF90_CLOBBER 
+            ELSE                   ;   imode = NF90_NOCLOBBER 
+            ENDIF
+            CALL iom_nf90_check(NF90_CREATE( TRIM(cdname), imode, if90id, chunksize = ichunk ), clinfo)
+            CALL iom_nf90_check(NF90_SET_FILL( if90id, NF90_NOFILL, idmy                     ), clinfo)
             ! define dimensions
             CALL iom_nf90_check(NF90_DEF_DIM( if90id, 'x', kdompar(1,1)  , idmy ), clinfo)