Discussion:
CONTOUR: Data coordinate system not established error
(too old to reply)
v***@gmail.com
2014-07-15 03:09:38 UTC
Permalink
Hello,

I am trying to plot a temperature contour vs log pressure and latitude. I can successfully produce a plot with a log y axis or with a linear reversed y axis (high pressure at the bottom). However when I try to use both keywords I get the above error. Below is a copy of the code in question:

contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],yrange=[1000,0],/ylog,/fill, $
Levels=level, C_Color=col1

The following two versions do work -

contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],yrange=[1000,0],/fill, $
Levels=level, C_Color=col1

contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],/ylog,/fill, $
Levels=level, C_Color=col1


Any ideas what might be going wrong?
David Fanning
2014-07-15 03:40:09 UTC
Permalink
Post by v***@gmail.com
contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],yrange=[1000,0],/ylog,/fill, $
Levels=level, C_Color=col1
The following two versions do work -
contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],yrange=[1000,0],/fill, $
Levels=level, C_Color=col1
contour, zt(*,*,itime), lat,zp(*,*,itime), color=0, background=255, $
Xtitle='Latitude', Ytitle='Pressure [Pa]', $
Title = 'Zonally Averaged Surface Temperature at Ls='+string(ls_year(itime)), $
thick = 6, Xrange = [-90,90],/ylog,/fill, $
Levels=level, C_Color=col1
Any ideas what might be going wrong?
I would say what is wrong is you are trying to take the Log(0), which is
undefined. Try setting your YRange=[1,1000]. That will work better. :-)

Cheers,

David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.idlcoyote.com/
Sepore ma de ni thue. ("Perhaps thou speakest truth.")
Loading...