Dear Simon,
I think, the first thing you have to decide is how you want to go about the different resolutions of your data set.
You can regrid the topography to the ¼ degree, or leave as is, which would mean that you would get a lot of identical values for all the 1 minute locations
that are in the same ¼ degree grid cell.
You could also downscale your 1/4 degree precip to the 1 minute, but I guess, that might be more difficult.
To achieve this your data sets need to be on proper grids; I suggest you have proper lat and lon axis. Then ferret knows how to deal with averaging and such.
For the technical side look at defining grids and regridding in ferret in the manual and mailing list. Or send some more information on the data sets – then
somebody might be able to give you a bit more detail.
You should be able to do the masking as soon you have the data on grids ferret can understand – so when they are on a lat/lon grids – whether you average, downscale
or not.
‘hope that was of some help,
Best wishes,
Joerg
From: owner-ferret_users@xxxxxxxx [mailto:owner-ferret_users@xxxxxxxx]
On Behalf Of simon EK
Sent: 16 October 2015 10:36
To: Andre Paim; Ansley Manke; Ferret users mailing list
Subject: Re: [ferret_users] ferret - repeat command - range = 70:80:.25 -error
I have a daily rainfall dataset on .25x.25 grid and
I have a topography dataset with 1 minute resolution.
I would like to study the impact of altitude on rainfall. Hence I would like to get spacial averages of rainfall over different altitude regions, say
0 to 100 meters
100 to 500 meters
500 to 1000 meters
and above 1000 meters.
I have seen certain ferret scripts using masking, etc.
Can you pl. help me?
Regards
Simon
On Fri, Sep 18, 2015 at 8:40 AM, simon EK <simonambat@xxxxxxxxx> wrote:
Thank you so much for the timely help.
The solution works well without any error.
1. repeat.jnl
===================================================
use rfclim_1901_2013.nc
set region/x=66.4E:100.1E/y=6.4N:38.6N/t=01-JAN 12:00 to 31-DEC 17:48
! Repeat on 10.5deg Latitude
repeat/x=76.25:79.75:.25 (go test.jnl)
===================================================
2. test.jnl
===============================
let rf = rf_clim[y=10.5N, l=1:365 ]
plot rf
frame/format=gif/file="rfclim-10-5x`x`.gif"
===============================
But there is a problem in creation of the graph files.
It seems that x assumes only integers in creating file name. The files created are:
=================
rfclim-10-5x76.gif
rfclim-10-5x77.gif
rfclim-10-5x78.gif
rfclim-10-5x79.gif
rfclim-10-5x80.gif
================
Value of x ranges from 76.25 to 79.75 in steps of .25. The gif images are created, but files names are not created for decimal values.
How to resolve this, any idea?
On Fri, Sep 18, 2015 at 1:50 AM, Andre Paim <paim.oceano@xxxxxxxxx> wrote:
Dear Simon,
Why don't you just use X in your repeat? Something like:
repeat/x=8:30:.25 (go test.jnl)
================================================
let rf = rf_clim[ y=10.5N, l=1:365 ]
plot rf
frame/format=gif/file="rfclim-10-5x`x`.gif"
================================================
That might solve your problem.
On 17 September 2015 at 13:50, simon EK <simonambat@xxxxxxxxx> wrote:
hi all,
I have grided rainfall dataset with .25x.25 deg resolution.
I want to plot a graph of daily rainfall for each grid point along the lonitude with .25 deg interval and I tried using repeat command. But I get error.
1. repeat.jnl
===============================================
use rfclim_1901_2013.nc
set region/x=66.4E:100.1E/y=6.4N:38.6N/t=01-JAN 12:00 to 31-DEC 17:48
! Repeat on 77.5deg Longitude
repeat/name=a/range=8:30:.25 (go test.jnl)
===============================================
2. test.jnl
================================================
let rf = rf_clim[x=`a`, y=10.5N, l=1:365 ]
plot rf
frame/format=gif/file="rfclim-10-5x`a`.gif"
====================================================
While running, I get the following error:
yes? go repeat-plot-long.jnl
use rfclim_1901_2013.nc
set region/x=66.4E:100.1E/y=6.4N:38.6N/t=01-JAN 12:00 to 31-DEC 17:48
! Repeat on 77.5deg Longitude
repeat/name=a/range=8:30:.25 (go test.jnl)
**ERROR: command syntax: .25
repeat/name=a/range=8:30:.25 (go test.jnl)
Command file, command group, or REPEAT execution aborted
**ERROR: invalid command: Delta value on REPEAT must be nonzero
Can somebody help me in resolving the issue?