[Thread Prev][Thread Next][Index]

RE: [ferret_users] vland and vfland - global meridional sections

Thanks Ansley for the information.


I rethought what I was trying to do, and although it might be nice to have an “averaged” bathymetry over my zonally averaged plots, it really doesn’t lend much to the finished product. So I’ll leave it alone for now.


The use of vfland and vland do nicely with my section plots though. They are a great tool.


Best from down under,



From: owner-ferret_users@xxxxxxxx [mailto:owner-ferret_users@xxxxxxxx] On Behalf Of Ansley Manke
Sent: Tuesday, 29 September 2015 3:05 AM
To: ferret_users@xxxxxxxx
Subject: Re: [ferret_users] vland and vfland - global meridional sections


Hi Pearse,
Yes, the arguments to those scripts are simple numbers or ranges, not expressions.

To see what the scripts are doing, use

   yes? go/help vland
   yes? go/help vfland

It would be quite simple to do the same thing as vfland, using the x-average of a topography variable.  In vland.jnl, the script opens a topo dataset, and does this plot command.

   PLOT/NOLAB/overlay/VLIM=($ppl$ymax):($ppl$ymin)/SET  -1*fv_rosedepth[($fv_xyregion)]

but in defining the XY region, the script requires a range in one direction and a scalar in the other.  Instead, it sounds as if you want something along these lines:

   USE etopo20
   PLOT/NOLAB/overlay/VLIM=($ppl$ymax):($ppl$ymin)/SET  -1*rose[x=@AVE,y=-90:90]

vfland.jnl sets up some variables and symbols, and calls the "fill_between" script.  Here, I'd suggest you watch it in action for a single value of X, to see how you could set things up to call it with an x-averaged topography variable.  Make your plot calls:

   yes? fill my_variable  

then turn on "mode verify" so you can see the commands as they're run.

   yes? set mode verify:always
   yes? go vfland 20 gray x=150 y=-90:90

What you're looking for in that output, is how the "GO fill_between" script is called. 

On 9/28/2015 12:17 AM, Pearse Buchanan wrote:

Hi fellow ferreters,


I’m trying to make some meridional section plots over various regions of the ocean. This includes a global zonal average, as well as zonal averages for the Indian, pacific and Atlantic oceans with various variables (please see the attached pretty plots).


I came across the vland and vfland scripts during my searching, and was attempting to use them. However, it seems that you must define an exact vertical/horizontal section, i.e. no averages!


This is how I tried to get it to work originally, with no success:

go vfland 20 grey x=@ave y=-90:90

Command file, command group, or REPEAT execution aborted

**ERROR: error: data is not a line


It seems the x=@ave argument is not recognised. Is there an alternative?


Any help would be appreciated.





Pearse J. Buchanan

PhD Candidate / CSIRO-UTAS Quantitative Marine Science

Institute for Marine and Antarctic Studies (IMAS), University of Tasmania




~~~ “Ocean deoxygenation: A Palaeo-Modelling Perspective” ~~~

Parts of the Ocean are predicted to lose oxygen in the coming century. But why?

Looking into the past, how does oxygen in the ocean change across climate transitions?


University of Tasmania Electronic Communications Policy (December, 2014).
This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement