time-slice

Hello,

I'm trying to write a simple script to output a new .nc file at each
timestep. The data is 3 dimensional (time-365, lat-354, lon-720). I am
able to take a timeslice one at a time with the following script. The
output has no metadata. How can I change this? Also, what would a "do
loop" look like that would output a new file(365 files) with different
names for every time (0:364)?

Thanks,
Nate

  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
  a =
addfile("/home/mikl6340/MsTMIP_Model/DLEM/BG1_DLEM_v2_daily_GPP_2000.nc","r")
  b = a->GPP(0,:,:)

  c = ndtooned(b)
  system("/bin/rm -f day1.nc") ; remove any pre-exist file
  cout = addfile("day1.nc", "c") ; new netCDF file
  cout@title = "day1"
  cout->c = c

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Mar 14 10:29:42 2013