[Thread Prev][Thread Next][Index]

regridding problem



Regridding over years
        Fellow ferreters, I'm having a problem regridding from our orignial time
coordinate to a new one. The original data represents a series of six month aug-jan
runs from 1991-2001. The new time coordinate extends from aug 15 1991 to jul 14 2001
When we do the transformation the first august point is dropped. Here's is the code
we are using:

! NOAA/PMEL TMAP
 ! FERRET v5.51  
 ! Linux 2.4.3-12smp - 03/05/03
 ! 05-Sep-03 13:52     

can var/all
can dat/all
define sym metafile = AM2p12b_clim 
define sym psfile = AM2p12b_clim.ps
define sym varb = precip
define sym epath = /archive/bgs/FMS/ATMOS/jakarta/N45L24
define sym expert = am2p12b_smip
define sym mon_beg = aug
define sym mon_end = jan
define sym year1_monthIC = 1991($mon_beg)
define sym year = 1991($mon_beg)2001($mon_end)
define sym factor = 86400
! define sym factor = 1
use "($epath)/($expert)_($year1_monthIC)_01/history/ts/10yr/($year).($varb).nc"
let num_ens = 1
set mode metafile ($metafile) 
set region/z=200
let ($varb)_ens_mn = ($varb)[d=1]/num_ens*($factor)
def ax/t=14-($mon_beg)-1991:15-jul-2001:1/units=months/calendar="julian" t2 ! Full year axis
let ($varb)_ens_mn2 = ($varb)_ens_mn[gt=t2@nrst]
plot/x=120E/y=5N/symbol=30 ($varb)_ens_mn[l=1:12]
plot/over/x=120E/y=5N/symbol=88 ($varb)_ens_mn2[l=1:25]
quit

        Any thoughts on why the first point is being dropped. Any help would be 
most appreciative. 

-- 


--------------------------------------------------------------------------------

     _/_/_/_/ _/_/_/_/ _/_/_/   _/        Robert G. Smith
    _/       _/       _/    _/ _/         Computer Assistant
   _/  _/_/ _/_/_/   _/    _/ _/          Geophysical Fluid Dynamics Laboratory
  _/    _/ _/       _/    _/ _/           P.O. Box 308
 _/_/_/_/ _/       _/_/_/   _/_/_/_/      Princeton, NJ  08542

--------------------------------------------------------------------------------
             Electronic Mail                             Telephone
       Internet Address: Bob.Smith@noaa.gov             Commercial: (609) 452-6593
   PUGL - http://www.gfdl.gov/~bgs/pugl/pugl.html                                                   FAX: (609) 987-5063
================================================================================



[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement