Re: ocean depth versus time (or lat/lon), emphasising upper ocean

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Thu, 25 Jun 2009 06:25:21 -0600 (MDT)

Mark,

[FYI: I saw your attachments in the original message.]

The reason you're seeing the linear axis is because the depth
values that you used for "fo" can't be used by NCL (which I
don't understand), so hence it reverts back to using
linear coordinates.

NCL doesn't like the jump from using coordinates that are
spaced by 100 to coordinates that are spaced by 1000.

In fact, there are several warning mesages that are spit out as a
result of this, and it states it's reverting back to a linear axis:

   warning:_NhlCreateSplineCoordApprox: Attempt to create spline
       approximation for Y axis failed: consider adjusting trYTensionF
       value
   warning:IrTransSetValues: error creating spline approximation for
       trYCoordPoints; defaulting to linear

Just for grins, I added one extra value to your "yo" array, 1500,
and this got rid of the errors and gave me an irregular axis.

yo = (/10.,100.,200.,300.,400.,500.,600.,700.,800.,900.,1000.,1500., \
        2000.,3000.,4000.,5000.,6000./)

See the attached "test.png".

--Mary

On Thu, 25 Jun 2009 Mark.Collier_at_csiro.au wrote:

> Hi,
> it appears that my attachments could not attach. I have downloaded my ncl script and file to the top level of ftp.dar.csiro.au,macftp,macftp if someone can help.
> Regards,
> Mark.
>
> ________________________________________
> From: Collier, Mark (CMAR, Aspendale)
> Sent: Thursday, 25 June 2009 1:11 PM
> To: ncl-talk_at_ucar.edu
> Subject: ocean depth versus time (or lat/lon), emphasising upper ocean
>
> Hi,
> I am trying to plot ocean depth vs time. for a number of ocean models (data has been supplied to me from these different modelling groups) and I want to emphasis the upper ocean within one plot (not 2 plots on top of each other as is commonly done) so that I can take advantage of the panelling tools.
>
> When I generate the plot without
>
> res_at_gsnYAxisIrregular2Linear = False
>
> it gives me a good representation of what I want, however, a bit too much emphasis on the very upper data because of the high resolution near the ocean surface.
>
> And so I thought of regridding the data to a depth coordinate that would give me something closer to what I want, and the idea is to do this for all models so that I have a consistent depth axis.
>
> However, when I do this it gives me a perfectly linear depth axis without any emphasis on the first 1000m. I suspect some attribute is missing from the coordinate.
>
> I have attached my script (just for one ocean model) and left it the code to plot the reasonable scaling for the depth axis which can be turned on by swapping the 3 lines where it says "here". When the existing code
> is used depth goes the wrong way and doesn't respond to the resource trYReverse but that is of secondary concern at the moment.
>
> Thanks,
> Mark.
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Thu Jun 25 2009 - 06:25:21 MDT

This archive was generated by hypermail 2.2.0 : Thu Jun 25 2009 - 06:30:15 MDT