currently SET data sets:1> ./rainfall.topo.nc (default)
name title I J K L
MEAN_RAINFALL mean rainfall 1:360 1:180 ... ...
TOPO etopo 1:360 1:180 ... ...
let mask = IF (TOPO GE 100 AND topo LT 1000) THEN 1let rf = MEAN_RAINFALL*mask
Here, I got only one level which is from 100 to 1000m. I could not able to define that as a function of K dimension.
Or Is it any possible way to deal with this situation. Please suggest.
Regards
Pratik
Hi,If I am understanding correctly, you should do this manually by listing each coordinates' rain and elevation value to a text file and then reading it later. Following are few steps:(1) let m_rain=if mean_rainfall ge 0 then mean_rainfall else (-1)! To make bad values or missing values as -1. Later use -1 as bad(2) list/clobber/file=Elevation_rain.dat/nohead/form=(f9.3,x,f9.3,x,f9.3,x,f9.3,) topo,m_rain,x[gx=topo],y[gy=topo]! ^ Last x and y are just to keep track of lat and long in case required in future.(3) Come out of ferret session and restart new session and read the file by : file/var=ele,m_r,lon,lat Elevation_rain.dat ; sh da ;By this you can do rain vs elevation analysis.On Mon, May 27, 2019 at 6:49 AM Pratik Kad <pratikkad17@xxxxxxxxx> wrote:Dear All,I need to save rainfall with respect to elevation. I tried with masking by elevation (Etopo) data; also I faced rang issue during visualisation.How can I take z-axis as elevation topography for rainfall parameter from attached data?Thanks in advance.Regards
PratikRegards,G Satyesh, India.