Changeset 5609 for trunk/NEMOGCM/TOOLS/SIREN/src/iom_rstdimg.f90

Timestamp:

2015-07-17T17:42:15+02:00 (9 years ago)

Author:

jpaul

Message:

commit changes/bugfix/... for SIREN; see ticket #1580

File:

• trunk/NEMOGCM/TOOLS/SIREN/src/iom_rstdimg.f90 (modified) (3 diffs)

trunk/NEMOGCM/TOOLS/SIREN/src/iom_rstdimg.f90

@@ -1058 +1058 @@
       INTEGER(i4), DIMENSION(ip_maxdim)            :: il_start
       INTEGER(i4), DIMENSION(ip_maxdim)            :: il_count
+
       REAL(dp),    DIMENSION(:,:,:,:), ALLOCATABLE :: dl_value
 
@@ -1204 +1205 @@
       ENDIF
 
+      ! force to change _FillValue to avoid mistake 
+      ! with dummy zero _FillValue
+      IF( td_var%d_fill == 0._dp )THEN
+         CALL var_chg_FillValue(td_var)
+      ENDIF
+
+      ! use scale factor and offset
+      WHERE( td_var%d_value(:,:,:,:) /= td_var%d_fill )
+         td_var%d_value(:,:,:,:) = &
+         &  td_var%d_value(:,:,:,:)*td_var%d_scf + td_var%d_ofs
+      END WHERE
+
    END SUBROUTINE iom_rstdimg__read_var_value
    !-------------------------------------------------------------------
@@ -1660 +1673 @@
          ! change FillValue to 0.
          CALL var_chg_FillValue(td_file%t_var(ji),0._dp)
+
+         ! use scale factor and offset
+         WHERE( td_file%t_var(ji)%d_value(:,:,:,:) /= &
+         &      td_file%t_var(ji)%d_fill )
+            td_file%t_var(ji)%d_value(:,:,:,:) = &
+            &  (td_file%t_var(ji)%d_value(:,:,:,:)-td_file%t_var(ji)%d_ofs) /&
+            &    td_file%t_var(ji)%d_scf
+         END WHERE
 
          cl_name(ji)  = TRIM(td_file%t_var(ji)%c_name)