time step when using filwgts_lanczos

Hi everyone,

I'm trying to filter some satellite data, which have different time
units and time steps.
My original NCL script was to band-pass filter daily data, that did have
a data for each day.
(it works well).
Now, I have to deal with weekly data, but the time unit remains being
the day.
Another data set I have to filter has a 3-day time step, the unit still
being the day.

So I don't know how to "tell" the filtering function all these different
time information.
I need to filter between 10 and 50 days.
When I put liminf = 10 and limp = 50, the filter doesn't seem to
understand I want to band-pass filter between 10 and 50 days.

Do you know how I can do that ?

Hereafter you can find the beginning of my script.

Thank you very much for any help !

Cecile

;=============================================================
  type_filt = 2
  liminf = 10
  limsup = 50
  date1 = 20000101
  date2 = 20051231
;=============================================================
  if ((data.eq."SSH").and.(expid.eq."TP_ERS"))then
    new_units = "days since 1992-10-14" ; because ut_calendar is
sensitive to the unit case !!!
  end if

  if ((data.eq."SSH").and.(expid.eq."TPJ_JPL"))then
    new_units = "days since 1992-01-01 12:00:00"
    date1 = 20000102
    date2 = 20051231
  end if

  print("Reading field")
  r=addfile(indir+infile,"r")
  time = r->time

  time_at_units = new_units
  printVarSummary(time)
  date = ut_calendar(time,2)
  print(date)
  i1 = ind(date.eq.date1)
  i2 = ind(date.eq.date2)

  fld1= r->$fldname1$(i1:i2,{south:north},{west:east})

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 28 2008 - 16:42:24 MDT