Re: Encountering problem with filedimdef

From: Hobbs, Will R (3244-CalTech) <William.R.Hobbs_at_nyahnyahspammersnyahnyah>
Date: Tue Dec 28 2010 - 18:10:12 MST

Andy

If you are defining a record dimension (i.e. an unlimited dimension, time in your case) it must always be the first (leftmost) coordinate.

If you want time as an unlimited dimension, then set

  dim_names = (/ "time","lon", "lat", "lev"/)
  dim_sizes = (/-1, nlon , nlat , nlev /) ;dimension size of an unlimited coord should be -1, regardless of actual length
  dimUnlim = (/ True, False , False, False /)
  filedimdef( ncf, dim_names, dim_sizes, dimUnlim )

If it's more important to you that time is the rightmost dimension, then just set it to size 'ntim'

 dim_names = (/ "lon", "lat", "lev","time" /)
  dim_sizes = (/-1, nlon , nlat , nlev,ntim /) ;dimension size of an unlimited coord should be -1, regardless of actual length
  dimUnlim = (/ False , False, False, False /)
  filedimdef( ncf, dim_names, dim_sizes, dimUnlim )

Hope that helps,

Will

**********************************************************
Will Hobbs, Ph.D. William.R.Hobbs@jpl.nasa.gov
Jet Propulsion Laboratory
4800 Oak Grove Dr. office: 300-324a
M/S 300-323 phone: (818) 354-0466
Pasadena, CA 91109 fax: (818) 354-0966
**********************************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 28 18:10:54 2010

This archive was generated by hypermail 2.1.8 : Tue Jan 04 2011 - 09:16:52 MST