<div dir="ltr"><div>I have taken a 'quick' look at your script. I do not 'know' the model (ie, data) you are using. Your script was a bit confusing. :-)<br><br>--<br></div><div>Attached:<br><br></div><div>[1] <br>A very basic script that illustrates plotting the SWVP variable *and* using the linint2 and area_hi3lores interpolation functions on this field. This is fairly simple. <br><br>[2]<br></div><div>I took a look at your original script. I was a bit confused. I made some changes. I know the numbers are *not correct*. Again I do not know your data. I computed some sample weighted vertical sums but not the weighted vertical averages [ sum(X*DP)/sum(DP) ] . I did not look at or make any changes to your plot codes. You will have to change these. This should give you a start. The key is for you to to correctly compute the layer thickness(es) [dp]. I think your surface pressures must be (PP+p ???).<br><br>[3]<br></div><div>When you send questions to ncl-talk, it is best to indicate the version of NCL you are using. Further, clean scripts are appreciated. <br></div><div><br></div><div>Good luck<br></div><div>D<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 29, 2016 at 10:10 AM, Paolina Bongioannini Cerlini <span dir="ltr"><<a href="mailto:paolina.cerlini@unipg.it" target="_blank">paolina.cerlini@unipg.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Sorry sorry, I forgot to say that in BM-256x256x64_3km_SWVP-write.<wbr>ncl script you need to read before<br>
file with 2D output and then 3D output, since there are different variables stored.<br>
I put both types of files in the incoming ftp dir.<br>
:)<br>
<span class="im HOEnZb"><br>
Paolina<br>
----- Original Message -----<br>
From: "Dennis Shea" <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>><br>
To: "paolina cerlini" <<a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>><br>
Cc: "ncl-talk" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
</span><span class="im HOEnZb">Sent: Thursday, 29 September, 2016 04:19:58<br>
</span><div class="HOEnZb"><div class="h5">Subject: Re: [ncl-talk] Fwd: regrid and plot<br>
<br>
What is needed is one or two sample netCDF files:<br>
<br>
EG:<br>
<br>
RCE_IDEAL_256x256x64_3000m_<wbr>10s_90days_uv0_256_0004237920.<wbr>nc<br>
RCE_...<br>
<br>
THX<br>
<br>
On Wed, Sep 28, 2016 at 10:47 AM, Paolina Bongioannini Cerlini <<br>
<a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>> wrote:<br>
<br>
> Hello,<br>
> thank you for suggesting me to upload the .ncl .ps<br>
> output.txt files in you server.<br>
> Paolina<br>
><br>
><br>
> ----- Original Message -----<br>
> From: "Dennis Shea" <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>><br>
> To: "paolina cerlini" <<a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>><br>
> Cc: "ncl-talk" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
> Sent: Wednesday, 28 September, 2016 18:30:02<br>
> Subject: Re: [ncl-talk] Fwd: regrid and plot<br>
><br>
> Could you please put the file onto our ftp area?<br>
><br>
> ftp <a href="http://ftp.cgd.ucar.edu" rel="noreferrer" target="_blank">ftp.cgd.ucar.edu</a><br>
> anonymous<br>
> your_email<br>
> cd incoming<br>
> prompt<br>
> put ...<br>
> quit<br>
><br>
> Please let ncl-talk know ehen the file transfer has been successfully<br>
> completed.<br>
><br>
> THX<br>
><br>
><br>
> On Wed, Sep 28, 2016 at 10:23 AM, Paolina Bongioannini Cerlini <<br>
> <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>> wrote:<br>
><br>
> ><br>
> ><br>
> > ----- Forwarded Message -----<br>
> > From: "paolina cerlini" <<a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a>><br>
> > To: "ncl-talk" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
> > Sent: Wednesday, 28 September, 2016 17:24:51<br>
> > Subject: regrid and plot<br>
> ><br>
> > Hello NCL people,<br>
> ><br>
> > I have on output from a model and I want to read it, to elaborate it,<br>
> > then to re-grid it on a lower resolution grid, then re-plot it.<br>
> ><br>
> > I succeeded in doing it but<br>
> > the plot is not doing the right axes and titles.<br>
> > Maybe it is some trivial error about Levels.<br>
> > I attach the script and output both on the screen and .ps<br>
> > Thank for your help.<br>
> ><br>
> > Paolina<br>
> ><br>
> > --<br>
> > Paolina Bongioannini Cerlini Ph.D.<br>
> > University of Perugia<br>
> > CIRIAF -room 8<br>
> > 06125 - Perugia<br>
> ><br>
> > Tel.<a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
> > e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
> > skype: pcerlini<br>
> > --<br>
> > --<br>
> > Paolina Bongioannini Cerlini Ph.D.<br>
> > University of Perugia<br>
> > CIRIAF -room 8<br>
> > 06125 - Perugia<br>
> ><br>
> > Tel.<a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
> > e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
> > skype: pcerlini<br>
> > --<br>
> ><br>
> > ______________________________<wbr>_________________<br>
> > ncl-talk mailing list<br>
> > <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
> > List instructions, subscriber options, unsubscribe:<br>
> > <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
> ><br>
> ><br>
> --<br>
> Paolina Bongioannini Cerlini Ph.D.<br>
> University of Perugia<br>
> CIRIAF -room 8<br>
> 06125 - Perugia<br>
><br>
> Tel.<a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
> e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
> skype: pcerlini<br>
> --<br>
> ______________________________<wbr>_________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
><br>
--<br>
Paolina Bongioannini Cerlini Ph.D.<br>
University of Perugia<br>
CIRIAF -room 8<br>
06125 - Perugia<br>
<br>
Tel.<a href="tel:%2B39%20075%205853576" value="+390755853576">+39 075 5853576</a><br>
e-mail: <a href="mailto:paolina.cerlini@unipg.it">paolina.cerlini@unipg.it</a><br>
skype: pcerlini<br>
--<br>
</div></div></blockquote></div><br></div>