Re: non conventional missing value type problem

From: Will Hobbs <Will.Hobbs_at_nyahnyahspammersnyahnyah>
Date: Sat Nov 02 2013 - 19:32:56 MDT

Sam

  1. NCL does not recognise the 'mssing_value' attribute; it only recognises the '_FillValue' attribute.

2. In your case, since there is no existing '_FillValue' attribute, you can just set U@_FillValue to whatever value/datatype is appropriate. If this were not the case, you would just need to delete the '_FillValue' attribute first.

Will

From: "Sam McClatchie (NOAA Federal)" <sam.mcclatchie@noaa.gov<mailto:sam.mcclatchie@noaa.gov>>
Reply-To: "sam.mcclatchie@noaa.gov<mailto:sam.mcclatchie@noaa.gov>" <sam.mcclatchie@noaa.gov<mailto:sam.mcclatchie@noaa.gov>>
Date: Sunday, 3 November 2013 7:03 AM
To: NCL maillist <ncl-talk@ucar.edu<mailto:ncl-talk@ucar.edu>>
Subject: non conventional missing value type problem

Hello

I am working with a netcdf of global drifter data <
http://www.aoml.noaa.gov/phod/dac/dac_meanvel.php>
that appears to have a non-conventional missing value.

The variables are mainly type double, but for some reason they appear to have coded the missing value as integer type, unless I read this dump incorrectly:

Variable: f
Type: file
filename: drifter_monthlymeans
path: drifter_monthlymeans.nc
   file global attributes:
      title : drifter_monthlymeans
      description : Climatological Monthly Means of Drifter Data
      version : Version 2.04
      date : Last updated on 24-June-2013
      date2 : File created 24-Jun-2013 13:49:57
      datadate : Climatology updated using data through 31-December-2012
      institution : NOAA/AOML & NOAA/PMEL
      author : Lumpkin and Johnson (2013) J. Geophys. Res., doi:10.1002/jgrc.20210
      CreationDate : 2013/06/24 13:53:04
   dimensions:
      Lat = 317
      Lon = 720
      Time = 12
   variables:
      double Lat ( Lat )
         long_name : Latitude (negative = south)
         units : degrees
         valid_min : -73
         valid_max : 85

      double Lon ( Lon )
         long_name : Longitude (negative = west)
         units : degrees
         valid_min : -180
         valid_max : 179.5

      double Time ( Time )
         long_name : Climatological Month (0=mid-Jan, 1=mid-Feb, etc.)
         calendar : noleap
         valid_min : 0
         valid_max : 11

      double U ( Lat, Lon, Time )
         long_name : Zonal Velocity
         units : meters per second
         valid_min : -2.8581
         valid_max : 4.3286
         missing_value : ( 9, 9, 9, ., 9, 9, 9 )

I don't think rewriting the missing value is an option here, as described in < http://www.ncl.ucar.edu/Support/talk_archives/2011/2628.html>

I tried Dennis' option for turning off default behaviour for missing values posted here < http://www.ncl.ucar.edu/Support/talk_archives/2009/1870.html>

NCL's "setfileoption" allows the user to turn off NCL's default behavior
 http://www.ncl.ucar.edu/Document/Functions/Built-in/setfileoption.shtml
      setfileoption ("nc", "/MissingToFillValue"/, False) ; turn off
default behavior

but I'm obviously still not solving the problembecause my code still gives an error on loading
:
;****************
;; load data
;***************
 setfileoption ("nc", "/MissingToFillValue/", False) ; turn off

  f = addfile("../data/global_drifter_current_climatology/drifter_monthlymeans.nc", "r")
  names = getfilevarnames(f)
  U = f->U(:,:,:)
  V = f->V(:,:,:)
  SST = f->SST(:,:,:)
  eU = f->eU(:,:,:)
  eV = f->eV(:,:,:)
  eSST = f->eSST(:,:,:)
  N = f->N(:,:)
  Latitude = f->Lat
  Longitude = f->Lon
  Month = f->Time

exit
;******************

Error is:
smc@caiman:~/Dropbox/dynamic_data/projects2013/visualization/ncl$ ncl map_global_drifter_currents_monthly_climatology.ncl
 Copyright (C) 1995-2012 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 6.1.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
warning:["NclFile.c":1657]:FileSetFileOption: /MissingToFillValue/ is not a recognized file option for format nc
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (U) types differ: not adding virtual _FillValue attribute
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (V) types differ: not adding virtual _FillValue attribute
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (SST) types differ: not adding virtual _FillValue attribute
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (eU) types differ: not adding virtual _FillValue attribute
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (eV) types differ: not adding virtual _FillValue attribute
warning:NetOpenFile: MissingToFillValue option set True, but missing_value attribute and data variable (eSST) types differ: not adding virtual _FillValue attribute

OK, so I'm obviously not doing something correctly. Can someone please point out my error?

Thanks and Best Fishes

Sam

--
Sam McClatchie,
Supervisory oceanographer, Fisheries oceanography
Southwest Fisheries Science Center, NOAA,
8901 La Jolla Shores Dr.
La Jolla, CA 92037-1509, U.S.A.
Office: 858 546 7083, Cellular:  858 752 8495
Research home page<http://www.fishocean.info>
[cid:part6.08050204.00000803@noaa.gov]


_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

CALCOFI 63rdLogo150b.png
Received on Sat Nov 2 19:33:12 2013

This archive was generated by hypermail 2.1.8 : Mon Nov 11 2013 - 09:45:33 MST