1 |
edhill |
1.1 |
|
2 |
|
|
Examples: NetCDF & MatLAB |
3 |
|
|
========================== |
4 |
|
|
|
5 |
|
|
These directories contain examples of MatLAB commands for reading, |
6 |
|
|
creating, and modifying NetCDF files. All the scripts use Chuck |
7 |
|
|
Denham's MatLAB package which can be obtained at: |
8 |
|
|
|
9 |
|
|
"MexCDF and NetCDF Toolbox For Matlab-5&6" |
10 |
|
|
http://woodshole.er.usgs.gov/staffpages/cdenham/public_html/MexCDF/nc4ml5.html |
11 |
|
|
|
12 |
|
|
and must first be installed. |
13 |
|
|
|
14 |
|
|
|
15 |
|
|
The NetCDF files generated here are intended for use with Ingrid: |
16 |
|
|
|
17 |
|
|
http://ingrid.mit.edu/ |
18 |
|
|
|
19 |
|
|
As a result, they approximately follow the simple COARDS convention |
20 |
|
|
which is described at: |
21 |
|
|
|
22 |
|
|
http://www.unidata.ucar.edu/packages/netcdf/conventions.html |
23 |
|
|
http://ferret.wrc.noaa.gov/noaa_coop/coop_cdf_profile.html |
24 |
|
|
|
25 |
|
|
|
26 |
|
|
The key items for the use of NetCDF files within Ingrid are: |
27 |
|
|
|
28 |
|
|
1) In order to work with the Ingrid "plot land" option, dimensions |
29 |
|
|
should have units of "degree_east" and "degree_north" as, for |
30 |
|
|
example, specified by: |
31 |
|
|
|
32 |
|
|
nc{'lon'}.units = 'degree_east'; |
33 |
|
|
nc{'lat'}.units = 'degree_north'; |
34 |
|
|
|
35 |
|
|
The actual dimension names (eg. "lat" or "Y") are unimportant. |
36 |
|
|
|
37 |
|
|
2) It is generally a good idea to set the "gridtype" to 1 or 0 to |
38 |
|
|
indicate that the coordinate "wraps" or not. The gridtype |
39 |
|
|
attribute allows Ingrid to, for instance, plot longitudinal data |
40 |
|
|
over "wrapped" ranges such as -100 to 500deg. The attribute can |
41 |
|
|
be set using commands such as: |
42 |
|
|
|
43 |
|
|
nc{'X'}.gridtype = ncint(1); |
44 |
|
|
nc{'Y'}.gridtype = ncint(0); |
45 |
|
|
|
46 |
|
|
3) Ingrid does seem to correctly handle the "_FillValue" attribute |
47 |
|
|
and it can be set using a MatLAB command such as: |
48 |
|
|
|
49 |
|
|
nc{ idu }.FillValue_ = ncdouble(-9999.0); |
50 |
|
|
|
51 |
|
|
[ Note how the "_" character is moved to the end of the attribute |
52 |
|
|
name since MatLAB cannot correctly handle variable names with a |
53 |
|
|
leading underscore. ] |