q defval/set defval

[Charles Seman]

Elena,

Please find attached a script "gcss5_diff.gs" which has been used to
difference two netCDF files.

The script checks each vertical level for each time level...  currently
it does (x,y) slab differences, but could be adapted to do
point-by-point differences if loops over y and x were added (following
the technique for the z loop)...  If you have a question on how to do
this, please send me an email, OK?

Thanks,
Chuck

Elena A. Yegorova wrote:

>Hello GrADS users,
>
>I am trying to use "q defval" and "set defval" to compare two grids
>pixel-by-pixel and create a new grid based on the results. This computation
>takes a very long time. And my computation is more involved than using
>arithmetic formulas of the two grids like grid1-grid2.
>Is there a more efficient way for creating a new grid with a unique value
>for each grid point (other than set defval)?
>Is there a more efficient way to loop through grid points (other than q
>defval)?
>
>Thanks for you help.
>
>Elena
>
>

--

********************************************************************
 Charles Seman                                charles.seman at noaa.gov
 U.S. Department of Commerce / NOAA / OAR
 Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
 201 Forrestal Road                              fax: (609) 987-5063
 Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not reflect
any position of the Government or NOAA."

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: gcss5_diff.gs
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20051024/f48dac2c/attachment.pl