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
-- email signature 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>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Nov 2 16:48:53 2013
This archive was generated by hypermail 2.1.8 : Mon Nov 11 2013 - 09:45:34 MST