Re: sea level pressure contours

Hi Alexander,
Are you getting any error messages that are being returned by the
plotting functions? I would guess that what is going on is that you are
losing your latitude/longitude coordinate variable information in this
coding:

pres = in->MSL_GDS0_SFC(108,:,:)
var2= pres/100

Since you are dividing pres by 100 NCL is not sure that you want to keep
the coordinate variable info, and is not passing the coordinate
variables from the right hand sign of the equal sign to the left hand
side. Add a couple printVarSummary statements to see what I mean:

pres = in->MSL_GDS0_SFC(108,:,:)
printVarSummary(pres)
var2= pres/100
printVarSummary(var2)

Assuming I'm right, try changing your coding to this:

pres = in->MSL_GDS0_SFC(108,:,:)
var2 = pres
var2 = (/ pres/100 /)

or

pres = in->MSL_GDS0_SFC(108,:,:)
pres = pres/100
.....
plot2 = gsn_csm_contour(wks,pres,res2)

Hope that helps. If not, let ncl-talk know.
Adam

On 03/19/2012 12:22 PM, Alexander Semenov wrote:
> Hi!
>
> I'm trying to overlay 2 variables: sea level pressure and temperature
> from ERA reanalysis.
>
> For some reason, sea level pressure field (contours) is completely
> messed up. I attached the plot and the code.
>
> Appreciate your help
>
> --
> regards
>
> *******************************************************
> *Alexander Semenov*
> *
> *
> PhD Student - Research Assistant
> International Arctic Research Center
> Department of atmospheric sciences
> University of Alaska Fairbanks
>
> 930 Koyukuk dr, 408c3
> Fairbanks, AK, USA, 99775
> work phone +19074742672
>
>
>
> _______________________________________________
> 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
Received on Mon Mar 19 12:47:31 2012