Oops, here is the attachment.
--Dave
Dave Allured wrote:
This is easy to do in NCL. See attached NCL script. This demo makes a
new Netcdf file with a 3D variable as you specified. By changing to
"append mode" in the addfile statement, you could add the same kind of
array to a pre-existing Netcdf file. Will an NCL solution work for you?
--Dave A.
CU/CIRES Climate Diagnostics Center (CDC)
NOAA/ESRL/PSD, Climate Analysis Branch (CAB)
http://www.cdc.noaa.gov/
Boulder, Colorado
303-497-6825
David Wang wrote:
Hello,
i've been scratching my head how to add a three-dimensional (time,
lat, lon) variable into a netcdf file where time is the unlimited
dimension and both lat and lon are dimensions of size 1. i work with
matlab (nc_add_recs.m in snctools) but matlab doesn't allow trailing
singleton dimensions. so while nc_add_recs.m (and nc_varput.m as well
presumably) requires an input of rank 3, i'm not able to provide it an
array of size (ntime, 1, 1). it's a pain in the neck. if there is any
workaround out there, please give me a hint.
thanks,
david
===============================================================================
To unsubscribe netcdfgroup, visit:
http://www.unidata.ucar.edu/mailing-list-delete-form.html
===============================================================================
;-----------------------------------------------------------------------------
;
; singleton.ncl - Demo script to create a netcdf variable dimensioned (n,1,1).
; 2007-aug-03 Original version. By Dave Allured.
;
; This version names the three dimensions, but does not define the
; coordinate arrays. Creation of coordinate arrays is documented on
; the NCL website.
;
;-----------------------------------------------------------------------------
begin
n = 3
dims = (/ n, 1, 1 /) ; create a simple 3D array
data = new (dims, "float")
data(:,0,0) = (/ 11, 22, 33 /)
data!0 = "time" ; name the three dimensions
data!1 = "lat"
data!2 = "lon"
print (data)
outfile = addfile ("out.nc", "c") ; create a new Netcdf file
time_name = "time" ; special to make unlimited dimension
time_dummy_size = 0
time_unlimited = True
filedimdef (outfile, time_name, time_dummy_size, time_unlimited)
outfile->var = data ; write data array and companion
; dimensions to the output file
end
