Hi Helen,
Looking at your resource list, the white area of the plot should
correspond to areas that have missing data, is that the case?
(a print(u(0,{.01},:,0)) should confirm that.
Also, I think you want to change the units of lev_p in the top script
from "hPa" to "Pa". That might alter the way vinth2p works, as you are
specifying the units of your input pressure array incorrectly.
The way I see it, I don't think your alterations to gsn_csm.ncl should
have done anything to the plot, because all you are doing is altering
the YR axis labels.
Adam
Helen Parish wrote:
> I cannot get the top section of my figures to plot correctly. There
> seems to be a mysterious gap in the middle of the top of the plot (the
> white part) which I cannot seem to eliminate. My routine does an
> interpolation to pressure levels, followed by a plot of the variable
> versus pressure and height. Since I have changed the pressure level
> structure in CAM to 50 levels, going from the ground to around 95 km, I
> had to estimate the correspondence between the new pressure levels and
> height.
>
> In order to do this I have adjusted the routine "gsn_csm.ncl". I enclose
> the section of gsn_csm.ncl which I have modified (in the file
> "gsnpart.ncl"), showing the adjusted pressure and height levels.
>
> I have also enclosed a typical anomalous plot, and the ncl script
> ("test.ncl")I used to generate the plot (containing the levels I am
> interpolating onto).
>
> Does anyone know what is causing this problem?.
>
> Thanks,
> Helen.
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-- -------------------------------------------------------------- Adam Phillips asphilli_at_ucar.edu National Center for Atmospheric Research tel: (303) 497-1726 ESSL/CGD/CAS fax: (303) 497-1333 P.O. Box 3000 Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli _______________________________________________ ncl-talk mailing list ncl-talk_at_ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talkReceived on Fri Jun 13 2008 - 13:05:07 MDT
This archive was generated by hypermail 2.2.0 : Fri Jun 13 2008 - 16:48:10 MDT