[gradsusr] Meteorological field in circular area

[davide.sacchetti at arpal.gov.it]

you could create a mask with a circular shape:

'define mask=maskout(lon,-1)'
loop over i, j in region of interest
'set defval mask i j 1'
check if it's a circle (d mask)
'set gxout stat'
'd field*mask'
it's all
hoping it helps
bye bye
Davide
  _____  

From: Rodrigo de Souza Barreto Mathias [mailto:rodrigowind at yahoo.com.br]
To: gradsusr at gradsusr.org [mailto:gradsusr at gradsusr.org]
Sent: Sat, 25 Feb 2012 06:05:32 +0100
Subject: [gradsusr] Meteorological field in circular area




Good morning!
I wanna calculate a meteorological field in a circular area with radius R starting  at a point (lat, lon). For example: average temperature in a circular area with 200 km of radius from a point 20N/020W. I don't know a specific function in Grads to do what I want. Could you help me with some function or a script using a mask?


Thanks a lot!
Rodrigo Mathias  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120227/a45bd9d2/attachment-0003.html