Re: nc_add_recs.m

[Dave Allured <dave.allured@xxxxxxxx>]

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