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

Timestamp:

2007-01-04T10:45:24+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_bugfix_001 : SM + CT : - add an initialization of iom_file(kiomid)%luld(kiv) field

• change "offset" to "add_offset" characters
• use a loop instead of the MINLOC function to find the minimum value of an array
• allow to write/read headers of DIMG files over 2 records (iom_rstdimg.F90)

File:

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

trunk/NEMO/OPA_SRC/IOM/iom_nf90.F90

@@ -60 +60 @@
       INTEGER            ::   if90id   ! nf90 identifier of the opened file
       INTEGER            ::   idmy     ! dummy variable
+      INTEGER            ::   jl       ! loop variable
       !---------------------------------------------------------------------
 
@@ -109 +110 @@
       ! =============
       IF( istop == nstop ) THEN   ! no error within this routine
-         kiomid = MINLOC(iom_file(:)%nfid, dim = 1)
+!does not work with some compilers         kiomid = MINLOC(iom_file(:)%nfid, dim = 1)
+         kiomid = 0
+         DO jl = jpmax_files, 1, -1
+            IF( iom_file(jl)%nfid == 0 )   kiomid = jl
+         ENDDO
          iom_file(kiomid)%name   = TRIM(cdname)
          iom_file(kiomid)%nfid   = if90id
@@ -174 +179 @@
          iom_file(kiomid)%ndims(kiv)  = i_nvd
          CALL iom_nf90_check(NF90_Inquire_Variable(if90id, ivarid, dimids = idimid(1:i_nvd)), clinfo)   ! dimensions ids
-         DO ji = 1, i_nvd   ! dimensions size
+         iom_file(kiomid)%luld(kiv) = .FALSE.   ! default value
+         DO ji = 1, i_nvd                       ! dimensions size
             CALL iom_nf90_check(NF90_Inquire_Dimension(if90id, idimid(ji), len = iom_file(kiomid)%dimsz(ji,kiv)), clinfo)   
             IF( idimid(ji) == iom_file(kiomid)%iduld ) iom_file(kiomid)%luld(kiv) = .TRUE.   ! unlimited dimension? 
          END DO
-         !---------- Deal with scale_factor and offset
+         !---------- Deal with scale_factor and add_offset
          llok = NF90_Inquire_attribute(if90id, ivarid, 'scale_factor') == nf90_noerr
          IF( llok) THEN
@@ -185 +191 @@
             iom_file(kiomid)%scf(kiv) = 1.
          END IF
-         llok = NF90_Inquire_attribute(if90id, ivarid, 'offset') == nf90_noerr
+         llok = NF90_Inquire_attribute(if90id, ivarid, 'add_offset') == nf90_noerr
          IF( llok ) THEN
-            CALL iom_nf90_check(NF90_GET_ATT(if90id, ivarid, 'offset', iom_file(kiomid)%ofs(kiv)), clinfo)
+            CALL iom_nf90_check(NF90_GET_ATT(if90id, ivarid, 'add_offset', iom_file(kiomid)%ofs(kiv)), clinfo)
          ELSE
             iom_file(kiomid)%ofs(kiv) = 0.