compositing wind and pressure data based on certain OLR values

[Charles Seman]

Jason,

Perhaps the "fndlvl" function may be helpful to define a "mask" which
then could be used with "maskout"?  The maskout web site has an example
of averaging only over land using a "mask"...
http://grads.iges.org/grads/gadoc/gradfuncfndlvl.html
http://grads.iges.org/grads/gadoc/gradfuncmaskout.html

Chuck

Jason Snyder wrote:
> I am trying to make a script that composites pressure values from dates
> when OLR values over a certain area (latitude range: 5S to 5N and
> longitude: 120E and 120W) are below a certain value.  I tried using this
> script and it did not seem to work properly.  What should I do to make
> this script so that I can make this script work properly?  Also how can I
> composite wind data at specific pressure levels such as 850 mb and 500 mb
> based on this script?
> 'sdfopen olr.day.mean.nc'
> 'sdfopen pres.sfc.1993.nc'
> 'set t 6790'
> 'define presa = 0'
> 'define number = 0'
> count = 6790
> while (count < 12267)
>    'set t 'count
>    'define olra = aave(olr,lon=120,lon=240,lat=-5,lat=5)'
>    if (olra < 230)
>       'define presb = pres.2'
>       'define n=1'
>    else
>       'define presb = 0'
>       'define n = 0'
>    endif
> 'define presa = presa + presb'
> 'define number = number + n'
> count = count + 1
> endwhile
> 'define presavg=presa'
>
>
> If someone has any suggestions on what to do with this script to make it
> composite wind/pressure data for days where the average OLR values within
> the area 120 W to 120 E and 15S to 15N, is less than 234, that would be
> greatly appreciated.
>
> Thanks,
>
> Jason
> -Jason
>

--

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
 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 necessarily
reflect any position of the Government or NOAA."