Dear Jonathan and Praaven,
OK I think it works now. Nevertheless, I had to
replace: let your_var_glev = your_var + 0 * x[gx=temp[d=1]] + 0 *
y[gy=temp[d=1]] by: let your_var_glev = temp[d=1] - temp[d=1] +
your_var + 0 * x[gx=temp[d=1]] + 0 * y[gy=temp[d=1]] because with the first
instruction the 3D structure of levitus isn’t replicated in
your_var_glev.
Thanks a lot. Have a nice
weekend. Mario
----- Original Message -----
Sent: Friday, March 11, 2005 6:43
PM
Subject: Re: representation of a single
point
Mario,
What you need to do is open levitus_climatology
and use it's X and Y axes. Here are more explicit instructions:
use levitus_climatology
use your_data ! A Z only file with your_var in it
let your_var_glev = your_var + 0 * x[gx=temp[d=1]] + 0 * y[gy=temp[d=1]]
At this point you have replicated your data along the levitus X and Y
axes and now you can proceed to create a mask.
--
Jon
Mario Germano wrote:
Thanks Jonathan,
Sorry but TEMP(ZAXLEVITR, YAXLEVITR, XAXLEVITR)
in levitus_climatology.cdf is a 3D variable.
If I do the same with my TEMP I end up with the
following message:
"ERROR: dimensions improperly specified: must
be a 2D region"
How can I handle this ?
Thanks again.
Mario
-----
Original Message -----
Sent:
Friday, March 11, 2005 5:09 PM
Subject:
Re: representation of a single point
Both,
Praveen has supplied half the solution.
I think Mario also needs help converting his 1D data into a 3D dataset
first. Here's the complete solution:
use levitus_climatology
! 1) Create a profile
! 2) Replicate that profile to be a 3D dataset based on the X and Y axes from the temp variable
! 3) Mask the 3D dataset so that only the original profile has valid data
let prof = temp[x=220,y=20]
let prof_glev = prof + 0 * x[gx=temp] + 0 * y[gy=temp]
let masked_prof_glev = IF x[gx=temp] GT 219 AND x[gx=temp] LT 221 AND y[gy=temp] GT 19 AND y[gy=temp] LT 21 THEN prof_glev
! Now you can use the shade command with your profile data
shade masked_prof_glev[k=1]
go land
-- Jon
Praveen V K wrote:
hi Mario
I tried your case with netcdf file
it is very easy. I think You may already know this one
==================================================================
yes? use levitus_climatology
yes? let sst = if X[gx=temp] gt 60 AND X[gx=temp] lt 75 AND Y[gy=temp] gt 10 AND Y[gy=temp] lt 15 then temp
yes? fill/k=1 sst
yes? go land
==================================================================
The above one will give a square plot at 60E-75E and 10N-15N and a white
area elsewhere
good luck
Praveen
On Fri, 11 Mar 2005, Mario Germano wrote:
Hi Ferreters,
I have a file of temperatures at a given location (X_0, Y_0) and at several
levels. I'd like to know how to display the temperature value of a given
level (say K=K_0) over a large domain in a SHADE-like fashion. The result
would be a coloured 1-km square at location (X_0,Y_0) and a white area
elsewhere (see attached figure).
Many thanks for your replies.
Mario.
|