problem in regridig of data

["Yogesh K. Tiwari" <ktiwari@xxxxxxxxxxxxxxx>]

Hello Ferret Users,

I caaled a script in ferret which counts the 
number of observations and then regrids the data

The problem is if I shade the data or try to save 
it ,it gets hanged for hours and hours.

Some times few data works and sometime another 
file contains same of data dosen't works.

Currently the file which I am using and got hanged 
has 40 mb data size.

If any person can help me, if the ferret hangs 
without giving any error then what may be the 
problem and how we can decide the check points.

I called script like that :

*********************************

yes? set data ecmwf_co2_jan2003.nc  !

 File name is ./ecmwf_co2_jan2003.nc

yes? set data tm3airs_new_January_2003.nc

 File name is ./tm3airs_new_January_2003.nc

yes?sh dat
     currently SET data sets:
    1> ./ecmwf_co2_jan2003.nc
 name       I         J         K         L
 LON       1:360     ...       ...       ...
 LAT       1:180     ...       ...       ...
 CO2_T     1:360     1:180     ...       ...


    2> ./tm3airs_new_January_2003.nc  (default)
 name       I         J         K         L

 LON       ...       ...       ...       1:119641
 LAT       ...       ...       ...       1:119641
 TM3CO2    ...       ...       ...       1:119641



yes? ! define destination axes from lat and lon in 
ecmwf_co2_jan2003


yes?      define axis/x/units=longitudes/from_data 
xlon=lon[d=1]

yes?      define axis/y/units=latitudes/from_data 
 ylat=lat[d=1]

yes?      define grid/x=xlon/y=ylat gcom

yes?      define symbol radius = 0.25

yes?      define symbol cutoff = 2.0

yes?      let lon_1D   = XSEQUENCE(lon[d=2])
yes?      let lat_1D   = XSEQUENCE(lat[d=2])
yes?      let co2_1D   = XSEQUENCE(tm3co2[d=2])

yes?      let obs_count = 
scatgrid_nobs_xy(lon_1D,lat_1D,x[gx=gcom],y[gy=gcom])


yes? !griding CO2 using the mask if there is 
atleast 10 observations in one grid

yes?      let co2_point = IF obs_count GE 10 THEN\
...? 
scat2gridgauss_xy(lon_1D,lat_1D,co2_1D,\
...? 
x[gx=gcom],y[gy=gcom],($radius),($radius),($cutoff),($cutoff))\
...?                      ELSE 1e35
 !-> DEFINE VARIABLE co2_point = IF obs_count GE 
10 THEN 
scat2gridgauss_xy(lon_1D,lat_1D,co2_1D, 
x[gx=gcom],y[gy=gcom],0.25,0.25,2.0,2.0) 
           ELSE 1e35


yes? shade co2_point
 File name is ./default.spk
 File name is /pool/ia64/ferret/v5.70/ppl/default.spk

*************************************************************

The ferret dosen't comes out from here after 
waiting of three-four hours.

Many thanks in advance for help,

Regards,
Yogesh

--

===========================================
Yogesh K. Tiwari,
Max-Planck Institute for Biogeochemistry,
Beutenberg Campus, Hans-Knoell-Strasse 10,
D-07745 Jena,
Germany

Office   : 0049 3641 576 376
Home     : 0049 3641 223 163
Fax      : 0049 3641 577 300
Cell     : 0049 1520 4591 008
e-mail   : yogesh.tiwari@bgc-jena.mpg.de
===========================================