Re: terrain vertical cross sectional

The Y coordinate array of the data array you pass to gsn_csm_pres_hgt must be a pressure array, and not a height array. That's why the error message is complaining about the units of "the first dimension", which refers to the leftmost dimension.

The gsn_csm_pres_hgt script is basically the same as gsn_csm_contour, except for the following:

  - It expects a Y coordinate pressure array in pascals, millibars, or hecto-pascals
  - It reverses the Y axis and uses a log scale
  - It puts height values on the right Y axis

I recommend using gsn_csm_contour if you are not interested in the above.

--Mary

On Nov 10, 2013, at 9:58 PM, Ahmed Lasheen <ahmed4kernel@gmail.com> wrote:

> Hello All
> I am trying to overlay terrain over vertical cross sectional of zonal wind. I am trying to use the ETOPO5 data. I first tried to plot it alone using the command plot = gsn_csm_pres_hgt(wks,elv,res), But i got the following error
> (0) gsn_csm_pres_hgt: Fatal: The coordinate array for the first dimension of the input data must be in Pascals, Hecto-pascals, or millibars
> (0) and it must contain the attribute 'units' set to one of the following strings (depending on your units):
> (0) 'mb' 'Mb' 'MB' 'millibar' 'millibars' 'MILLIBARS' 'hybrid_sigma_pressure' 'Pa' 'pa' 'PA' 'Pascals' 'pascals' 'PASCALS' 'hpa' 'hPa' 'HPA' 'hecto-pascals' 'HECTO-PASCALS'
>
> It looks like a problem related to the data itself.
> is there is any solution or other way to do that?
> I have attached the data and script for the vertical cross section of zonal wind that I am trying to plot the terrain over.
> thanks in advance
> ---------------
> Ahmed L. Shaaban | Graduate Student
> Department of Atmospheric and Environmental science(DEAS).
> State University of New York, University (SUNY) at Albany
> 1400 Washington Avenue, Albany, NY 12222
> ----------------
>
> <terrain_vertical_cross_section.ncl>_______________________________________________
> 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 Mon Nov 11 16:00:10 2013