Re: Plot Press vs Latitude

Adam,

I want plot the difference between model and observation, but when I did
it, appear the message:

(0) gsn_csm_pres_hgt: Fatal: The first dimension of the input data must
(0) have a coordinate variable called 'lev.'
(0) Cannot create plot.
fatal:Illegal right-hand side type for assignment
fatal:Execute: Error occurred at or near line 50 in file teste1.ncl

I can't understand. What did I wrong? My files are attach.

Thanks,

Guilherme.

2013/1/14 Adam Phillips <asphilli@ucar.edu>

> Hi Guilherme,
> Your level array runs from 1000->100, so I am guessing that your vertical
> coordinate is in units of hPa. Simply assign the vertical coordinate
> dimension unit attribute like this:
> u&lev@units = "hPa"
>
> and your script should be fine. Note that as you are setting your contour
> levels from -70:70 you might want to convert your u array to Celsius from
> Kelvin.. Adam
>
> On 01/14/2013 11:08 AM, Guilherme Martins wrote:
>
> Hi users,
>
> When I run my script (in attach and the file .nc) I get the message
> below:
>
> (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'
> (0) Cannot create plot.
> fatal:Illegal right-hand side type for assignment
> fatal:Execute: Error occurred at or near line 41 in file teste.ncl
>
> My file is output from MCGA model
>
> And I saw that vertical level is:
>
> double lev ( lev )
> long_name : generic
> units : level
> axis : Z
>
> How can I fix this?
>
> Thanks,
>
> Guilherme Martins.
>
> --
> *Guilherme Martins*
> http://guilhermemartins.hd1.com.br/
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> --
> ______________________________________________________________
> Adam Phillips asphilli@ucar.edu
> NCAR/Climate and Global Dynamics Division (303) 497-1726
> P.O. Box 3000
> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-- 
*Guilherme Martins*
http://guilhermemartins.hd1.com.br/

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