Changeset 886

Timestamp:

02/08/10 10:47:55 (15 years ago)

Author:

bellier

Message:

added optional argument fillvalue to fliodefv
for external Netcdf _FillValue attribute

File:

• IOIPSL/trunk/src/fliocom.f90 (modified) (9 diffs)

IOIPSL/trunk/src/fliocom.f90

@@ -205 +205 @@
 !! SUBROUTINE fliodefv &
 !! & (f_i,v_n,[v_d],v_t, &
-!! &  axis,standard_name,long_name,units,valid_min,valid_max)
+!! &  axis,standard_name,long_name,units, &
+!! &  valid_min,valid_max,fillvalue)
 !!
 !! INPUT
@@ -227 +228 @@
 !! (C) axis,standard_name,long_name,units : Attributes
 !!     (axis should be used only for coordinates)
-!! (R) valid_min,valid_max : Attributes
+!! (R) valid_min,valid_max,fillvalue : Attributes
 !!--------------------------------------------------------------------
   MODULE PROCEDURE &
@@ -1287 +1288 @@
 SUBROUTINE fliodv_r0d &
  & (f_i,v_n,v_t, &
- &  axis,standard_name,long_name,units,valid_min,valid_max)
+ &  axis,standard_name,long_name,units,valid_min,valid_max,fillvalue)
 !---------------------------------------------------------------------
   IMPLICIT NONE
@@ -1296 +1297 @@
   CHARACTER(LEN=*),OPTIONAL,INTENT(IN) :: &
  & axis,standard_name,long_name,units
-  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max
+  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max,fillvalue
 !---------------------------------------------------------------------
   CALL flio_udv &
  &  (f_i,0,v_n,(/0/),v_t, &
- &   axis,standard_name,long_name,units,valid_min,valid_max)
+ &   axis,standard_name,long_name,units,valid_min,valid_max,fillvalue)
 !------------------------
 END SUBROUTINE fliodv_r0d
@@ -1306 +1307 @@
 SUBROUTINE fliodv_rnd &
  & (f_i,v_n,v_d,v_t, &
- &  axis,standard_name,long_name,units,valid_min,valid_max)
+ &  axis,standard_name,long_name,units,valid_min,valid_max,fillvalue)
 !---------------------------------------------------------------------
   IMPLICIT NONE
@@ -1316 +1317 @@
   CHARACTER(LEN=*),OPTIONAL,INTENT(IN) :: &
  & axis,standard_name,long_name,units
-  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max
+  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max,fillvalue
 !---------------------------------------------------------------------
   CALL flio_udv &
  &  (f_i,SIZE(v_d),v_n,v_d,v_t, &
- &   axis,standard_name,long_name,units,valid_min,valid_max)
+ &   axis,standard_name,long_name,units,valid_min,valid_max,fillvalue)
 !------------------------
 END SUBROUTINE fliodv_rnd
@@ -1326 +1327 @@
 SUBROUTINE flio_udv &
  & (f_i,n_d,v_n,v_d,v_t, &
- &  axis,standard_name,long_name,units,valid_min,valid_max)
+ &  axis,standard_name,long_name,units,valid_min,valid_max,fillvalue)
 !---------------------------------------------------------------------
   IMPLICIT NONE
@@ -1336 +1337 @@
   CHARACTER(LEN=*),OPTIONAL,INTENT(IN) :: &
  & axis,standard_name,long_name,units
-  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max
+  REAL,OPTIONAL,INTENT(IN) :: valid_min,valid_max,fillvalue
 !-
   INTEGER :: f_e,m_k,i_v,i_rc,ii,idd
@@ -1460 +1461 @@
         CALL ipslerr (2,'fliodefv', &
    &      'Variable '//TRIM(v_n),'attribute valid_max', &
+   &      'not supported for this external type')
+      END SELECT
+    ENDIF
+    IF (PRESENT(fillvalue)) THEN
+      SELECT CASE (m_k)
+      CASE(flio_i1,flio_i2)
+        i_rc = NF90_PUT_ATT(f_e,i_v,'_FillValue',NINT(fillvalue,KIND=i_2))
+      CASE(flio_i4)
+        i_rc = NF90_PUT_ATT(f_e,i_v,'_FillValue',NINT(fillvalue,KIND=i_4))
+      CASE(flio_r4)
+        i_rc = NF90_PUT_ATT(f_e,i_v,'_FillValue',REAL(fillvalue,KIND=r_4))
+      CASE(flio_r8)
+        i_rc = NF90_PUT_ATT(f_e,i_v,'_FillValue',REAL(fillvalue,KIND=r_8))
+      CASE DEFAULT
+        CALL ipslerr (2,'fliodefv', &
+   &      'Variable '//TRIM(v_n),'attribute fillvalue', &
    &      'not supported for this external type')
       END SELECT