<div dir="ltr"><div class="gmail_default" style="font-size:small">Can you please try wrf_map_overlays instead of wrf_map_overlay? Â wrf_map_overlay is a deprecated function.</div><div class="gmail_default" style="font-size:small">
<br></div><div class="gmail_default" style="font-size:small">Please see the wrf_map_overlays documentation for more information:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style>
<a href="http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml">http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml</a><br></div><div class="gmail_default" style><br></div><div class="gmail_default" style>
Also, there are some cases where the WRF functions don't recognize the gsnXXXX resources, and so it's possible that gsnPaperOrientation is getting ignored.</div><div class="gmail_default" style><br></div><div class="gmail_default" style>
You can try not setting gsnMaximize to True to see if that turns off the rotation of your plot.</div><div class="gmail_default" style><br></div><div class="gmail_default" style>If you continue to have problems, it would help if you could provide the data file needed to run your script.</div>
<div class="gmail_default" style><br></div><div class="gmail_default" style>Thanks,</div><div class="gmail_default" style><br></div><div class="gmail_default" style>--Mary</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On Mon, Aug 4, 2014 at 9:43 PM, æ¨çŽ‰åŽ <span dir="ltr"><<a href="mailto:yangyh@mail.typhoon.gov.cn" target="_blank">yangyh@mail.typhoon.gov.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<u></u>
<div>
<div><div align="left">
<p>Hi,</p>
<p>Â I'd like to specify the orientation of the plot to portrait, but I'm failed even I using the following scripts:</p>
<p>Â Â Â Â Â Â Â Â <a href="mailto:opts_r@gsnMaximize" target="_blank">opts_r@gsnMaximize</a>Â Â Â Â Â Â Â Â Â = True<br>Â Â Â Â Â Â Â Â <a href="mailto:opts_r@gsnPaperOrientation" target="_blank">opts_r@gsnPaperOrientation</a> = "portrait"</p>
<p>or </p>
<p>Â Â Â Â Â Â Â Â <a href="mailto:type@wkOrientation" target="_blank">type@wkOrientation</a> = "portrait"</p>
<p>Attachment is my scripts, please check it.</p>
<p>Thanks in advance.</p>
<p>Yuhua </p>
<p> </p><strong>您于 Tue, 5 Aug 2014 11:12:04 +0800 çš„æ¥ä¿¡ä¸å†™åˆ°: </strong><br><br></div>
<div>
<div dir="ltr">Dear Sir
<div><br></div>
<div>i am trying to read TRMM(tropical rainfall measuring mission data) on NCL version 6.1.2 but as the latitude and longitude are different from geographic values i.e from ncl_filedump i get total  longitude = 1440  latitude = 400, the spatial resolution is 25km .</div>
<div>i need to extract rainfall data at a location lat=1.3911 and long =103.9864 for that i used the formula </div>
<div>trmmlat=</div>
<table style="WIDTH:125pt;BORDER-COLLAPSE:collapse" border="0" cellspacing="0" cellpadding="0" width="167">
<tbody>
<tr style="HEIGHT:15pt" height="20">
<td style="WIDTH:125pt;HEIGHT:15pt" height="20" width="167">((lat<a href="tel:%2B49.875%29%2F0.25" value="+49875025" target="_blank">+49.875)/0.25</a>)+1<br><br>trmmlong=<br>
<table style="WIDTH:137pt;BORDER-COLLAPSE:collapse" border="0" cellspacing="0" cellpadding="0" width="182">
<tbody>
<tr style="HEIGHT:15pt" height="20">
<td style="WIDTH:137pt;HEIGHT:15pt" height="20" width="182">((long+179.875)/0.25)+1<br><br>and then extract the pcp variables using ncl at these trmmlat and trmmlong ,i want to make sure whether my method is correct or not.<br>
<br>regards<br>saurabh</td></tr></tbody></table></td></tr></tbody></table>
<div><br></div></div>
<div class="gmail_extra"><br><br>
<div class="gmail_quote">On Tue, Aug 5, 2014 at 2:00 AM, <span dir="ltr"><<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a>></span> wrote:<br>
<blockquote style="BORDER-LEFT:#ccc 1px solid;MARGIN:0px 0px 0px 0.8ex;PADDING-LEFT:1ex" class="gmail_quote">Send ncl-talk mailing list submissions to<br>Â Â Â Â <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<br>To subscribe or unsubscribe via the World Wide Web, visit<br>Â Â Â Â <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>or, via email, send a message with subject or body 'help' to<br>
    <a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a><br><br>You can reach the person managing the list at<br>    <a href="mailto:ncl-talk-owner@ucar.edu" target="_blank">ncl-talk-owner@ucar.edu</a><br>
<br>When replying, please edit your Subject line so it is more specific<br>than "Re: Contents of ncl-talk digest..."<br><br><br>Today's Topics:<br><br>  1. Re: triple2grid (Dennis Shea)<br>  2. Re: How to create two y axis with a common x axis in bar  plot<br>
   (Mary Haley)<br>  3. Re: (no subject) (Dennis Shea)<br>  4. Re: Running example scripts for beginners (Mary Haley)<br>  5. Re: read separated files (Dennis Shea)<br>  6. Using a loop to find values greater than 0 (calculate   PDD)<br>
   (Lauren Jean Vargo)<br>  7. Re: Using a loop to find values greater than 0   (calculate<br>   PDD) (Dennis Shea)<br>  8. Re: Using a loop to find values greater than 0   (calculate<br>   PDD) (Kyle Griffin)<br>
  9. Re: Strange white lines across plot (Mary Haley)<br> 10. Re: Fwd: Masking a U.S. county from a netcdf file (Mary Haley)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>
Date: Mon, 4 Aug 2014 07:17:29 -0600<br>From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>Subject: Re: [ncl-talk] triple2grid<br>To: Manisha Ganeshan <<a href="mailto:manisha.ganeshan86@gmail.com" target="_blank">manisha.ganeshan86@gmail.com</a>><br>
Cc: NCL-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>Message-ID:<br>Â Â Â Â <CAOF1d_4pLjOOnFbXb-7ju2BQbkXwbx=iCVsuVWjiUM=<a href="mailto:wyR_Efw@mail.gmail.com" target="_blank">wyR_Efw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br><br>Not that I am aware of. Sorry.<br><br><br>On Fri, Aug 1, 2014 at 4:10 PM, Manisha Ganeshan <<br><a href="mailto:manisha.ganeshan86@gmail.com" target="_blank">manisha.ganeshan86@gmail.com</a>> wrote:<br>
<br>> Thanks, Dennis. Is there a function that calculates the mean or saves the<br>> total number of nearest neighbor source points for each grid point in the<br>> destination grid?<br>><br>> Regards,<br>> Nisha<br>
><br>><br>> On Fri, Aug 1, 2014 at 6:03 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>><br>>> It does not calculate the mean. New 'nearest neighbor' values overwrite<br>
>> previous values.<br>>><br>>><br>>> On Fri, Aug 1, 2014 at 3:59 PM, Manisha Ganeshan <<br>>> <a href="mailto:manisha.ganeshan86@gmail.com" target="_blank">manisha.ganeshan86@gmail.com</a>> wrote:<br>
>><br>>>> Hello,<br>>>><br>>>> I have a question about the *triple2grid* function. The description<br>>>> says that the function does not perform interpolation but instead assigns<br>
>>> each unstructured (source) data point to the nearest grid point<br>>>> (destination). How does it treat multiple unstructured data points (source)<br>>>> that have a common nearest grid point (destination)? Does it calculate the<br>
>>> mean value (after ignoring missing values)?<br>>>><br>>>> Thanks,<br>>>> Nisha<br>>>><br>>>> _______________________________________________<br>>>> ncl-talk mailing list<br>
>>> List instructions, subscriber options, unsubscribe:<br>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>>>><br>
>>><br>>><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/81cbd948/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/81cbd948/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 2<br>Date: Mon, 4 Aug 2014 08:29:44 -0600<br>From: Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>Subject: Re: [ncl-talk] How to create two y axis with a common x axis<br>
    in bar  plot<br>To: Ipshita Majhi <<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>><br>Cc: NCL USERS <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <CACNN_CKd2OjK7wbauhLmC01QPRpsCMfMq2aYC1CDzqh=<a href="mailto:UQeVog@mail.gmail.com" target="_blank">UQeVog@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>Ipsita,<br>
<br>In order to overlay a bunch of bar charts, it works best if you can shift<br>each of the bars by a certain bar width, and then you can draw them<br>side-by-side. Â I think you may have used part of bar_11.ncl on our bar<br>
examples page, and this is a good one to work with. However, you want to<br>add some overlay code in there, which will make all the plots part of the<br>same plot.<br><br>The key is to draw each set of bars with a slightly shifted set of X<br>
values, so the bars are drawn next to each other.<br><br>See the attached example and script.<br><br>--Mary<br><br><br><br>On Fri, Aug 1, 2014 at 4:27 PM, Ipshita Majhi <<a href="mailto:ipmajhi@alaska.edu" target="_blank">ipmajhi@alaska.edu</a>> wrote:<br>
<br>> Hi,<br>> I want to plot trend values and monthly means in one plot. Both are an<br>> array of size 12 .<br>> The program gives no error but the plot are overlapping. I am not sure how<br>> to fix it. I realized that if I use the below mentioned method they overlap<br>
> is there a way to have two y axis in bar plots. Since the trends are much<br>> smaller then avgs.<br>><br>> ;*****************************************************<br>> ;Trend and monthly climatology together in one bar plot<br>
> ;*******************************************************<br>><br>> Â Â Â load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>> Â Â Â load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
><br>> Â Â Â months=(/1,2,3,4,5,6,7,8,9,10,11,12/)<br>><br>><br>> monthly_avg=asciiread("~/Documents/PhD_June_2015/NCL_Moonsoon/Data_output/Monthly_avg_air.txt",-1,<br>> "float")<br>><br>
> monthly_Trend=asciiread("~/Documents/PhD_June_2015/NCL_Moonsoon/Data_output/trend_air.txt",<br>> -1, "float")<br>> Â Â Â printVarSummary(monthly_avg)<br>><br>> Â Â Â wks = gsn_open_wks("ps","Trend_Avg")<br>
><br>><br>> Â Â Â sres = True<br>> Â Â Â sres@vpWidthF = 0.7<br>> Â Â Â sres@vpHeightF = 0.5<br>> Â Â Â sres@gsnDraw = True<br>> Â Â Â sres@gsnFrame = False<br>> Â Â Â sres@gsnXYBarChart = True<br>> Â Â Â sres@gsnXYBarChartBarWidth = 0.15 Â Â Â Â Â ; change bar widths<br>
>    sres@tmXBMode      = "Explicit"     ; explicit labels<br>>    sres@tmXBValues     = (/1,2,3,4,5,6,7,8,9,10,11,12/)<br>>    sres@tmXBLabels =<br>> (/"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"/)<br>
> Â Â Â sres@tmXBLabelFontHeightF = 0.0205<br>> Â Â Â sres@tmXTLabelFontHeightF = 0.0205<br>> Â Â Â sres@tmYLLabelFontHeightF = 0.0225<br>> Â Â Â sres@tiMainFontHeightF = 0.025<br>> Â Â Â sres@tiMainFont = "helvetica"<br>
> Â Â Â sres@tiMainString = "All India Rainfall Climatology and Trend<br>> 1880-2012"<br>> Â Â Â sres@tiYAxisString = "(mm)"<br>><br>> Â Â Â plot1 = gsn_csm_xy (wks,months,monthly_avg,sres)<br>
> Â Â Â sres@gsnXYBarChartColors = (/"red"/) Â Â Â Â Â Â Â Â Â Â ;<br>> seperately, not<br>> Â Â Â plot2 = gsn_csm_xy(wks,months,monthly_Trend,sres)<br>> Â Â Â Â Â Â Â Â ; advancing the frame<br>> Â Â Â sres@gsnXYBarChartColors = (/"blue"/) Â Â Â Â Â Â Â Â Â Â ; but<br>
> tweaking where<br>><br>><br>> _______________________________________________<br>> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.html</a><br>
-------------- next part --------------<br>A non-text attachment was scrubbed...<br>Name: bar.png<br>Type: image/png<br>Size: 42598 bytes<br>Desc: not available<br>Url : <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.png" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.png</a><br>
-------------- next part --------------<br>A non-text attachment was scrubbed...<br>Name: bar_11.ncl<br>Type: application/octet-stream<br>Size: 3061 bytes<br>Desc: not available<br>Url : <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.obj" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/b04f1178/attachment-0001.obj</a><br>
<br>------------------------------<br><br>Message: 3<br>Date: Mon, 4 Aug 2014 08:44:26 -0600<br>From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>Subject: Re: [ncl-talk] (no subject)<br>
To: Modise Wiston <<a href="mailto:modise.wiston@postgrad.manchester.ac.uk" target="_blank">modise.wiston@postgrad.manchester.ac.uk</a>><br>Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <<a href="mailto:CAOF1d_4GEj90MZ04-u1vemt3UbnH%2Bvd--oWt1Ax73BAP2poQag@mail.gmail.com" target="_blank">CAOF1d_4GEj90MZ04-u1vemt3UbnH+vd--oWt1Ax73BAP2poQag@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br>
<br>While not WRF, the last example at:<br><a href="http://www.ncl.ucar.edu/Applications/iso.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/iso.shtml</a><br><br>may give you an idea of how to approach the issue.<br>
<br>===<br>Also, please send *all* WRF questions to <a href="mailto:wrfhelp@ucar.edu" target="_blank">wrfhelp@ucar.edu</a><br>The know the WRF model best.<br>You can cc ncl-talk@ucar but, really, wrfhelp should be your primary<br>
contact for assorted WRF related issues.<br><br>Good luck<br><br><br><br>On Sat, Aug 2, 2014 at 12:06 PM, Modise Wiston <<br><a href="mailto:modise.wiston@postgrad.manchester.ac.uk" target="_blank">modise.wiston@postgrad.manchester.ac.uk</a>> wrote:<br>
<br>> Â *To whom it may concern*<br>><br>> Â *Hi all,*<br>><br>> Â *I am trying to determine the cloud droplet number<br>> concentration (QNDROP) from my [WRF-Chem] model simulation. The script [see<br>> the code below] is such that the code loops through the x-y-z (i-j-k)<br>
> planes and finds the max (& min) QNDROP along the vertical (model levels).<br>> The code is such that it exits/stops just after the first cloud band only<br>> -without continuing through all the layers above the first cloud band. This<br>
> also calculates CCN number concentrations below the cloud base.*<br>><br>> Â *How can I make it (or write a function/code) to calculate the height at<br>> which the max QNDROP occurs..? I have been trying to go through several NCL<br>
> examples (functions) but couldn't figure out out how to do this.. (I would<br>> like to plot 'height vs max QNDROP' on the horizontal lat-lon domain) *<br>><br>> Â *Any help is appreciated please..*<br>
> *Thank you,*<br>> *M. Wiston*<br>> *The University of Manchester  *<br>><br>> -----------------------------------------------------------------------------------------------------------<br>>  dt = 6<br>
> Â do ifil = 0, 0; numFILES-48, dt<br>> Â Â a = addfile(FILES(ifil)+".nc","r")<br>> Â Â ;list time step to output data<br>> Â Â times = wrf_user_list_times(a)<br>> Â Â ntimes = dimsizes(times)<br>
><br>> Â Â do it = 0,ntimes-1,2<br>> Â do ispec = 0, nspec-1<br>> Â Â Â Â Â species = speclist(ispec)<br>><br>> Â ;Obtain the variables from the file<br>> Â Â if(isfilevar(a,"QNDROP"))<br>> N = wrf_user_getvar(a,"QNDROP",it)<br>
> Nd = (N*1.2923)/1e6 Â Â Â Â Â Â Â Â Â Â Â ; convert from /kg to /sm^3 then<br>> to /scm^3<br>> Nd@units = "/cm^3" Â Â Â Â Â Â Â Â Â Â Â ; converted units<br>> Â Â end if<br>><br>> Â Â Â Â Â Â ;define an array to hold data<br>
> Â Â dimsNd = dimsizes(Nd)<br>> Â Â print(dimsNd)<br>><br>> Â Â Â Â Â Â ;----------------Determine the domain size (grid points and<br>> vert. levels)----------------;<br>> Â Â ; Â Â Â Â (321 x 251) Â Â Â Â Â Â Â Â Â Â Â ; grid points along "i"<br>
> (hor) by "j" (vert) direction  ;<br>>   ;   (e_we:lat, e_sn:lon)<br>>                     ;<br>>   ;   40 vert levels:(k)           ; model levels<br>>                    ;<br>
> Â Â ; Â Â Print output gives: k = 40 Â Â Â Â ;NB: this is read starting from<br>> zero (i.e. 0-39) Â Â ;<br>> Â Â ; Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â j = 250<br>> Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â ;<br>> Â Â ; Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â i = 320<br>
> Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â ;<br>><br>> Â ;------------------------------------------------------------------------------------------------;<br>> Â Â min_Nd = Â new( (/dimsNd(1),dimsNd(2)/),float )<br>
> Â Â max_Nd = Â new( (/dimsNd(1),dimsNd(2)/),float )<br>> Â Â c_btm = new( (/dimsNd(1),dimsNd(2)/),integer )<br>> Â Â c_top = new( (/dimsNd(1),dimsNd(2)/),integer )<br>><br>> Â Â do j = 0, dimsNd(1)-1 Â Â Â Â Â Â Â Â Â Â Â Â Â ; loop along j<br>
> do i = 0, dimsNd(2)-1 Â Â Â Â Â Â Â Â Â Â Â ; loop along i<br>> Â Â Â ; cldfrc = (i*j)<br>> Â Â Â c_btm(j,i) = -1 Â Â Â Â Â Â Â Â Â Â Â Â Â ;define cloud base and<br>> cloud top limits to confine the execution<br>
> Â Â Â c_top(j,i) = -1<br>> Â Â Â do k = 0, dimsNd(0)-2 Â Â Â Â Â Â Â Â ; loop [vertically] through<br>> levels<br>> Â ; print("k= "+k)<br>> Â if(Nd(k,j,i) .le. cut_pnt .and. Nd((k+1),j,i) .gt. cut_pnt)<br>
> Â Â Â Â c_btm(j,i) = (k+1)<br>> ;[cloud base]<br>> Â Â Â Â ; print("c_btm = "+c_btm(j,i))<br>><br>> Â else if (Nd(k,j,i) .gt. cut_pnt .and. Nd((k+1),j,i) .le. cut_pnt) then<br>> Â Â Â Â c_top(j,i) = k<br>
>  ;[cloud top]<br>> break               ; exit the the process after the first<br>> cloud band<br>>  end if<br>>      end if<br>>    end do<br>>            ; min_Nd(i) = min(Nd(:,i))<br>
> Â Â Â ; max_Nd(i) = max(Nd(:,i))<br>> Â Â Â ; min_Nd(j) = min(Nd(:,j))<br>> Â Â Â ; max_Nd(j) = max(Nd(:,j))<br>> Â Â Â ; min_Nd(j,i) = min(Nd(:,j,i))<br>> Â Â Â ; max_Nd(j,i) = max(Nd(0:3,j,i))<br>><br>> Â Â Â if(c_btm(j,i) .eq. -1 .and. c_top(j,i) .eq. -1)<br>
> Â Â Â Â Â Â max_Nd(j,i) = 0.0<br>><br>> Â Â Â else if (c_btm(j,i) .gt. -1 .and. c_top(j,i) .gt. -1) then<br>> Â Â Â Â Â Â max_Nd(j,i) = max( Nd(c_btm(j,i):c_top(j,i),j,i) )<br>> Â Â Â Â Â Â ;print(min_Nd(j,i))<br>
> Â Â Â Â Â Â ;print("max_Nd = " +max_Nd(j,i))<br>> Â print("| "+sprintf("%5.3f",max_Nd(j,i))+ " | "+sprintf("%2.1f",<br>> c_btm(j,i))+" | "+sprintf("%2.1f", c_btm(j,i))+" |" )<br>
> Â Â Â Â Â Â Â end if<br>> Â Â Â end if<br>> Â Â Â end do<br>> Â Â Â Â end do<br>><br>> Â ;------------------------------------------------------------------------------------------<br>><br>> Â Â Â Â ;;get the CCN number concentration from the file<br>
> Cn1 = wrf_user_getvar(a,"CCN1",it)<br>> Cn2 = wrf_user_getvar(a,"CCN2",it)<br>> Cn3 = wrf_user_getvar(a,"CCN3",it)<br>> Cn4 = wrf_user_getvar(a,"CCN4",it)<br>> Cn5 = wrf_user_getvar(a,"CCN5",it)<br>
> Â Â Â Â Â Cn6 = wrf_user_getvar(a,"CCN6",it)<br>> Â Â Â Â Â z = wrf_user_getvar(a,"z",it)<br>> Â Â Â Â Â w = wrf_user_getvar(a,"wa",it) Â Â Â Â Â Â Â Â Â ; vertical<br>> [updraft..??] velocity<br>
> Â Â Â Â Â w@units = "m/s"<br>><br>> Â p = wrf_user_getvar(a, "p",it) Â Â Â Â Â Â Â Â Â ;extract p and tk to<br>> convert to STP<br>> Â Â Â Â Â Â Â Â Â tk = wrf_user_getvar(a,"tk",it)<br>
> C1 = Cn1*((tk/p)*(P_std/T_std))<br>> C2 = Cn2*((tk/p)*(P_std/T_std))<br>> C3 = Cn3*((tk/p)*(P_std/T_std))<br>> C4 = Cn4*((tk/p)*(P_std/T_std))<br>> C5 = Cn5*((tk/p)*(P_std/T_std))<br>> C6 = Cn6*((tk/p)*(P_std/T_std))<br>
> C1@units = "/cm^3"<br>><br>> C1_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>> C2_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>> C3_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>> C4_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>
> C5_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>> C6_btm = new( (/dimsNd(1),dimsNd(2)/),float )<br>> SS = new( (/dimsNd(1),dimsNd(2)/),float )<br>><br>> do j = 0, dimsNd(1)-1<br>> Â Â Â do i = 0, dimsNd(2)-1<br>
> Â Â Â Â Â if(c_btm(j,i) .eq. -1)<br>> Â Â Â Â Â Â Â C1_btm(j,i) = 0.0<br>> Â Â Â Â Â Â Â C2_btm(j,i) = 0.0<br>> Â Â Â Â Â Â Â C3_btm(j,i) = 0.0<br>> Â Â Â Â Â Â Â C4_btm(j,i) = 0.0<br>> Â Â Â Â Â Â Â C5_btm(j,i) = 0.0<br>
> Â Â Â Â Â Â Â C6_btm(j,i) = 0.0<br>><br>> Â Â Â Â Â else if(c_btm(j,i) .ne. -1) then<br>> Â Â Â Â Â Â Â C1_btm(j,i) = ( C1(c_btm(j,i),j,i) )<br>> Â Â Â Â Â Â Â C2_btm(j,i) = ( C2(c_btm(j,i),j,i) )<br>> Â Â Â Â Â Â Â C3_btm(j,i) = ( C3(c_btm(j,i),j,i) )<br>
> Â Â Â Â Â Â Â C4_btm(j,i) = ( C4(c_btm(j,i),j,i) )<br>> Â Â Â Â Â Â Â C5_btm(j,i) = ( C5(c_btm(j,i),j,i) )<br>> Â Â Â Â Â Â Â C6_btm(j,i) = ( C6(c_btm(j,i),j,i) )<br>><br>> Â Â Â Â Â Â Â ; print("Cn1 = " +C1_btm(j,i))<br>
> Â Â Â Â Â Â Â ; print("Cn2 = " +C2_btm(j,i))<br>> Â Â Â Â Â Â Â ; print("Cn3 = " +C3_btm(j,i))<br>> Â Â Â Â Â Â Â ; print("Cn4 = " +C4_btm(j,i))<br>> Â Â Â Â Â Â Â ; print("Cn5 = " +C5_btm(j,i))<br>
> Â Â Â Â Â Â Â ; print("Cn6 = " +C6_btm(j,i))<br>> Â print("| "+sprintf("%4.3f", C1_btm(j,i))+" | "+sprintf("%4.3f",<br>> C2_btm(j,i))+" | "+sprintf("%4.3f", C3_btm(j,i))+" | "+\<br>
> sprintf("%4.3f", C4_btm(j,i))+" | "+sprintf("%4.3f", C5_btm(j,i))+" |<br>> "+sprintf("%4.3f", C6_btm(j,i))+" |" )<br>> Â Â Â Â end if<br>> Â Â Â Â end if<br>
> Â Â Â Â end do<br>> Â Â Â Â Â Â end do<br>><br>> ;`````````````````````````````````````````````````````````````````````````````````<br>><br>> Â Â Â Â Â Â Â Â Â ; if(max_Nd(j,i) .le. C1_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = 0.01<br>
><br>> Â Â Â Â Â Â Â Â Â ; Â else if(max_Nd(j,i) .gt. C1_btm(j,i) .and.<br>> max_Nd(j,i) .le. C2_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = "0.02%"<br>><br>> Â Â Â Â Â Â Â Â Â ; Â else if(max_Nd(j,i) .gt. C2_btm(j,i) .and.<br>
> max_Nd(j,i) .le. C3_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = "0.1%"<br>><br>> Â Â Â Â Â Â Â Â Â ; Â else if (max_Nd(j,i) .gt. C3_btm(j,i) .and.<br>> max_Nd(j,i) .le. C4_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = "0.2%"<br>
><br>> Â Â Â Â Â Â Â Â Â ; Â else if (max_Nd(j,i) .gt. C4_btm(j,i) .and.<br>> max_Nd(j,i) .le. C5_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = "0.5%"<br>><br>> Â Â Â Â Â Â Â Â Â ; Â else if (max_Nd(j,i) .gt. C5_btm(j,i) .and.<br>
> max_Nd(j,i) .le. C6_btm(j,i))<br>> Â Â Â Â Â Â Â Â Â Â Â Â ; Â SS(j,i) = "1.0%"<br>> Â Â Â Â Â Â Â Â Â ; Â end if<br>> Â Â Â Â Â Â Â Â Â ; Â end if<br>> Â Â Â Â Â Â Â Â Â ; Â end if<br>> Â Â Â Â Â Â Â Â Â ; Â end if<br>
> Â Â Â Â Â Â Â Â Â ; Â end if<br>> Â Â Â Â Â Â Â Â Â ; Â end if<br>><br>> ;-------------------------------------------------------------------------------------------<br>> Â Â Â ;Set some Basic Plot options<br>
>    res                 = True<br>>    res@tiMainString       = "CCN nmbr @: " +times(it)  ; main<br>> field title<br>>    mpres              = True<br>>    pltres               = True<br>
><br>>    res@gsnMaximize          = False<br>>    res@gsnFrame            = False<br>>    res@gsnDraw            = False<br>>    res@cnFillOn             = True<br>
>    res@cnLinesOn           = False    ; turn off contour<br>> lines<br>>    res@cnLineLabelsOn         = False    ; lables the data<br>> on contours/map grids<br>>    res@gsnSpreadColors        = True<br>
>    res@lbOrientation         = "Horizontal"     ;<br>> orientates the linear scale<br>>        res@cnLevelSelectionMode      = "ExplicitLevels"<br>><br>> ;-------------------------------------------------------------------------------------------<br>
>    res@tiYAxisString    = ""<br>>    res@tiXAxisString    = "Nmbr cont'n (/cm^3)"<br>>    res@tiXAxisSide     = "Bottom"<br>><br>>    ;res@tiMainString  = "Droplet nmbr at: " +times(it)<br>
> Â Â Â contour_Nd = wrf_contour(a,wks,max_Nd,res)<br>><br>> Â Â Â Â Â Â Â plot_Nd = wrf_map_overlays(a,wks,(/contour_Nd/),pltres,res)<br>> Â Â ; plots data on a map backgrnd<br>> Â Â Â ; contour_Cn1 = wrf_contour(a,wks,C1_btm,res)<br>
> Â Â Â ; plot_Cn1 = wrf_map_overlays(a,wks,(/contour_Cn1/),pltres,res)<br>> Â Â Â Â ; Â contour_Cn2 = wrf_contour(a,wks,C2_btm,res)<br>> Â Â Â ; Â plot_Cn2 = wrf_map_overlays(a,wks,(/contour_Cn2/),pltres,res)<br>> Â Â Â ; Â contour_Cn3 = wrf_contour(a,wks,C3_btm,res)<br>
> Â Â Â ; Â plot_Cn3 = wrf_map_overlays(a,wks,(/contour_Cn3/),pltres,res)<br>> Â Â Â ; Â contour_Cn4 = wrf_contour(a,wks,C4_btm,res)<br>> Â Â Â ; Â plot_Cn4 = wrf_map_overlays(a,wks,(/contour_Cn4/),pltres,res)<br>> Â Â Â ; Â contour_Cn5 = wrf_contour(a,wks,C5_btm,res)<br>
> Â Â Â ; Â plot_Cn5 = wrf_map_overlays(a,wks,(/contour_Cn5/),pltres,res)<br>> Â Â Â Â ; Â contour_Cn6 = wrf_contour(a,wks,C6_btm,res)<br>> Â Â Â ; Â plot_Cn6 = wrf_map_overlays(a,wks,(/contour_Cn6/),pltres,res)<br>
><br>> Â Â Â Â Â Â Â ; delete(contour_Nd)<br>> Â Â Â ; delete(contour_Cn1)<br>><br>> Â Â Â ; map = wrf_map(wks,a,res)<br>> Â Â Â ; overlay(map,contour_Nd)<br>> Â Â Â ; overlay(map,contour_Cn1)<br>> Â Â Â ; overlay(plot0,plot3)<br>
> Â Â Â ; gsn_panel(wks,plots,(/3,2/),True)<br>> Â Â Â Â ; draw(plot)<br>> Â Â Â ; frame(wks)<br>> Â Â Â Â Â end do<br>> Â Â end do<br>><br>> Â end<br>><br>> ---------------------------------------------------------------------------------------------------------------<br>
><br>> _______________________________________________<br>> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/be554295/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/be554295/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 4<br>Date: Mon, 4 Aug 2014 08:52:43 -0600<br>From: Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>Subject: Re: [ncl-talk] Running example scripts for beginners<br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>><br>Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <<a href="mailto:CACNN_CJySFTTZvaJQbpnCE%2BN3mBi3BncJCF7RxZHxgMmDd41nQ@mail.gmail.com" target="_blank">CACNN_CJySFTTZvaJQbpnCE+N3mBi3BncJCF7RxZHxgMmDd41nQ@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br>
<br>?Zilore,<br><br>WRF data is on a curvilinear grid, which means each data point is<br>represented by a unique lat/lon value. You said that "wrf_map_overlays"<br>didn't work for you, but it would help if you explain how it didn't work.<br>
 Did you get an error message? If so, include the error message. Did you<br>get the wrong plot? If so, indicate this. The wrf_xxxx scripts are tailored<br>for WRF ARW data, so they should work pretty well.<br><br>If you want to use the gsn_csm_xxxx scripts, you have to provide lat/lon<br>
information, because the gsn_csm_xxxx scripts are *not* tailored for WRF<br>data. Otherwise NCL will have no way of knowing what map projection or<br>where your data is supposed to be overlaid on the map.<br><br>If you are trying to use gsn_csm_contour_map and gsn_csm_contour to contour<br>
WRF data, then I suggest visiting this page:<br><br>?<a href="http://www.ncl.ucar.edu/Applications/wrfgsn.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/wrfgsn.shtml</a><br><br>The first example, wrf_gsn_1.ncl, is the simplest one, but it only draws<br>
one field on a map.<br><br>The last example, wrf_gsn_5.ncl, shows how to overlay multiple fields on<br>one plot. Â It shows how to use both gsn_csm_xxx scripts, or wrf_xxxx<br>scripts. Â It's up to you which one you want.<br>
<br>Please look at these examples and try one of them. Â If you continue to have<br>problems, include your script again, but this time include more information<br>about what is wrong.<br><br>Thanks,<br><br>--Mary<br><br><br>
<br>On Fri, Aug 1, 2014 at 4:05 PM, zilore mumba <<a href="mailto:zmumba@yahoo.com" target="_blank">zmumba@yahoo.com</a>> wrote:<br><br>> My sincere apologies for posting very basic issues. I have been trying to<br>
> make headway in learning ncl with very little progress. All the examples I<br>> try do not seem to work. Examples,<br>><br>><br>><br>><br>><br>><br>> wrfout_d01_2014-07-29_00^%00^%<a href="http://00.nc/" target="_blank">00.nc</a><br>
> <<a href="https://www.dropbox.com/s/aeoqmfgrk6pbd99/wrfout_d01_2014-07-29_00%5E%2500%5E%2500.nc" target="_blank">https://www.dropbox.com/s/aeoqmfgrk6pbd99/wrfout_d01_2014-07-29_00%5E%2500%5E%2500.nc</a>><br>> zilore mumba shared from Dropbox<br>
> View on <a href="http://www.dropbox.com/" target="_blank">www.dropbox.com</a><br>> <<a href="https://www.dropbox.com/s/aeoqmfgrk6pbd99/wrfout_d01_2014-07-29_00%5E%2500%5E%2500.nc" target="_blank">https://www.dropbox.com/s/aeoqmfgrk6pbd99/wrfout_d01_2014-07-29_00%5E%2500%5E%2500.nc</a>><br>
> Preview by Yahoo<br>><br>><br>> 1. I have tried the two examples on "wrf_map_overlay" examples found at<br>> wrf_map_overlays<br>> <<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml</a>> both<br>
> don't work for me, only modification being the filename (which works in<br>> grads)<br>> 2. Â I have tried the "Getting_started example 5"found at Getting started<br>> using NCL - example 5 - contour plots over maps<br>
> <<a href="https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml</a>>,<br>> it does not work for me.<br>
> 3. I have tried *overlay_*" examples found at<br>> <a href="https://www.ncl.ucar.edu/Applications/Scripts/overlay_1.ncl" target="_blank">https://www.ncl.ucar.edu/Applications/Scripts/overlay_1.ncl</a>, they do not<br>
> work for me.<br>> I would appreciate if someone can make overlay_1.ncl (below) work for me.<br>> Modifications are a) my filename and longitude and latitude of my domain.<br>> The wrf output file is attached if allowed since it is too large (through<br>
> dropbox). I moved to Windows which does not allow : in file names.<br>> The error I get is:<br>> fatal:Dimension (bottom_top) of file (wrfout_d01_2014-07-29_00:00:00) does<br>> not have an associated coordinate variable<br>
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 16 in<br>> file overlay_1.ncl<br>><br>> ;*************************************************<br>> ; overlay_1.ncl<br>> ;<br>> ; Concepts illustrated:<br>
> ; Â - Overlaying line contours on filled contours<br>> ; Â - Explicitly setting contour levels<br>> ; Â - Selecting a different color map<br>> ;<br>> ;*************************************************<br>
><br>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>><br>> begin<br>> Â a = addfile("wrfout_d01_2014-07-29_00:00:<a href="http://00.nc/" target="_blank">00.nc</a>","r")<br>
> Â temp = a->T(0,{500},:,:)<br>> Â uwnd = a->U(0,{500},:,:)<br>><br>> Â wks = gsn_open_wks("ps","overlay")<br>> Â gsn_define_colormap(wks,"BlueDarkRed18") Â Â ; use the BlueDarkRed18<br>
> colormap<br>><br>> Â res = True<br>> Â res@mpFillOn = False<br>> Â res@mpMaxLatF = Â -4.07 Â Â Â Â Â Â Â Â Â Â Â ; specify the plot domain<br>> Â res@mpMinLatF = -21.9 Â Â Â Â Â Â Â Â Â Â Â ;<br>> Â res@mpMinLonF = Â 17.25 Â Â Â Â Â Â Â Â Â Â Â ;<br>
>  res@mpMaxLonF =  38.44            ;<br>>  sres@gsnFrame = False            ; do not advance the frame<br>>  sres@cnLevelSelectionMode = "ExplicitLevels" ; use explicit levels<br>
> Â sres@cnLevels = ispan(-5,35,5) Â Â Â Â Â Â Â ; set the contour levels<br>><br>> Â plot = gsn_csm_contour_map(wks,temp,res) Â Â ; create the temperature plot<br>> Â plot_ov = gsn_csm_contour(wks,uwnd,sres) Â Â ; create the U-wind plot<br>
> Â overlay(plot,plot_ov) Â Â Â Â Â Â Â Â Â Â Â ; overlay the U-wind plot on<br>> the temperature plot<br>> Â draw(plot) Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â ; draw the temperature plot<br>> (with the U-wind plot overlaid)<br>
>  frame(wks)                  ; advance the frame<br>> end<br>><br>>  res@mpOutlineOn = True               ; turn the map<br>> outline on<br>>  res@gsnDraw    =  False          ; do not draw the plot<br>
><br>><br>><br>><br>> [image: image]<br>> <<a href="https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml</a>><br>
><br>><br>><br>><br>><br>> Getting started using NCL - example 5 - contour plots ov...<br>> <<a href="https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml</a>><br>
> NCL Home > Documentation > Manuals > Getting Started GSUN > Examples > 1 |<br>> 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11<br>> View on <a href="http://www.ncl.ucar.edu/" target="_blank">www.ncl.ucar.edu</a><br>
> <<a href="https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Getting_Started/Examples/gsun05n.shtml</a>><br>> Preview by Yahoo<br>
><br>><br>><br>><br>><br>><br>><br>><br>> wrf_map_overlays<br>> <<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml</a>><br>
> NCL WRFUserARW.ncl functions (graphics routines)<br>> View on <a href="http://www.ncl.ucar.edu/" target="_blank">www.ncl.ucar.edu</a><br>> <<a href="https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml</a>><br>
> Preview by Yahoo<br>><br>><br>><br>> _______________________________________________<br>> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/2718eee8/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/2718eee8/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 5<br>Date: Mon, 4 Aug 2014 08:55:42 -0600<br>From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>Subject: Re: [ncl-talk] read separated files<br>
To: Yuqiang Zhang <<a href="mailto:yuqiangzhang.thu@gmail.com" target="_blank">yuqiangzhang.thu@gmail.com</a>><br>Cc: NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>Message-ID:<br>
    <<a href="mailto:CAOF1d_4_aU1aywjH6CZvfAskDg276GTDvL%2B3c6hhZeDPHYM2RQ@mail.gmail.com" target="_blank">CAOF1d_4_aU1aywjH6CZvfAskDg276GTDvL+3c6hhZeDPHYM2RQ@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br>
<br>"ListSet Type" could be used  but the following may be easier.<br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/ListSetType.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/ListSetType.shtml</a><br>
<br>Also, your dats is type "short"<br><a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml</a><br>
<br>========<br>Please see 'addfiles'<br><br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a><br>
<br>===<br>Â Â scrDir0 = "/d1/zhangg/prec/"<br>Â Â scrFile0= systemfunc("cd "+ scrDir0 +" ; ls<br>TPCAS-CMFD_V0105_B-01_"+year+".nc"<br>Â Â print(srcFile0)<br><br>Â Â na = addfiles(srcDir0+srcFile0, "r")<br>
<br>* Â pp = short2flt( na[:}->prec(:,107,302))*<br><br>* Â printVarSummary(pp)*<br><br><br><br><br><br>On Sat, Aug 2, 2014 at 10:35 AM, Yuqiang Zhang <<a href="mailto:yuqiangzhang.thu@gmail.com" target="_blank">yuqiangzhang.thu@gmail.com</a>><br>
wrote:<br><br>> Hi,<br>><br>><br>><br>> You can use this function ?ListSetType? to read different dimensions data<br>> from multi files.<br>><br>><br>><br>> Regards,<br>><br>> Yuqiang<br>
><br>><br>><br>> *From:* <a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a> [mailto:<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>] *On<br>
> Behalf Of *guozfruit<br>> *Sent:* Saturday, August 02, 2014 11:38 AM<br>> *To:* NCL<br>> *Subject:* [ncl-talk] read separated files<br>><br>><br>><br>> Hi all,<br>><br>><br>><br>> I have 3-hr netcdf data in 12 separated files (ie for 12 months). I'm<br>
> trying to extract one grid preciptation for the whole year. I have the<br>> problem about how to define one variable with different dimensions when I<br>> read data from different files. Could you please help me with it? Thanks.<br>
><br>><br>><br>> As below are the data info, and my script:<br>><br>> dimensions:<br>><br>> lon = 700 ;<br>><br>> lat = 400 ;<br>><br>> time = UNLIMITED ;<br>><br>> short prec(time, lat, lon) ;<br>
><br>> Â Â prec:scale_factor = 0.002f ;<br>><br>> Â Â prec:add_offset = 50.f ;<br>><br>> Â Â prec:_Fill_Value = -32767s ;<br>><br>> Â Â prec:missing_value = -32767s ;<br>><br>> Â Â prec:units = "mm hr-1" ;<br>
><br>> Â Â prec:long_name = "Precipitation rate" ;<br>><br>><br>><br>><br>><br>> Â year = 2008<br>><br>> Â mon = 12<br>><br>> Â nhour = new(mon,integer)<br>><br>><br>><br>
>  if(year .eq. 2008) then<br>><br>>   day_leap   = (/31,29,31,30,31,30,31,31,30,31,30,31/)<br>><br>>   do i = 0,mon-1<br>><br>>    nhour(i) = day_leap(i)*8<br>><br>>   end do<br>><br>
> Â else<br>><br>> Â Â day_nonleap = (/31,28,31,30,31,30,31,31,30,31,30,31/)<br>><br>> Â Â do i = 0,mon-1<br>><br>> Â Â Â nhour(i) = day_nonleap(i)*8<br>><br>> Â Â end do<br>><br>> Â end if<br>
><br>><br>><br>> *pp = Â new((/12,nhour/),float)*<br>><br>> do imon = 1,12<br>><br>> Â Â Â if(imon.lt.10) then<br>><br>> Â Â Â Â mon3 = "0"+imon<br>><br>> Â Â Â else<br>><br>> Â Â Â Â mon3 = ""+imon<br>
><br>>    end if<br>><br>>    scrDir0 = "/d1/zhangg/prec/"<br>><br>>   do ifl = 0,nhour(imon-1)-1<br>><br>>    na  =<br>> addfile(scrDir0+"prec_ITPCAS-CMFD_V0105_B-01_"+year+mon3+".nc","r")<br>
><br>> Â Â * pp(imon-1,ifl) Â = na->prec(:,107,302)*<br>><br>> Â Â end do<br>><br>> Â end do<br>><br>><br>><br>><br>><br>><br>><br>> X<br>><br>> _______________________________________________<br>
> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/5bf2441c/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/5bf2441c/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 6<br>Date: Mon, 4 Aug 2014 15:16:12 +0000<br>From: Lauren Jean Vargo <<a href="mailto:lvargo@unm.edu" target="_blank">lvargo@unm.edu</a>><br>Subject: [ncl-talk] Using a loop to find values greater than 0<br>
    (calculate    PDD)<br>To: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>Message-ID: <<a href="mailto:70D79029-3BC6-4860-8081-E5C280D91F93@unm.edu" target="_blank">70D79029-3BC6-4860-8081-E5C280D91F93@unm.edu</a>><br>
Content-Type: text/plain; charset="Windows-1252"<br><br>Hello,<br><br>I?m having trouble doing a calculation with NCL, ultimately I am trying to calculate the Positive degree-day (PDD) sum. I?ve read in daily temperature data (which is 3D [time,lat,lon]). What I am trying to do is calculate the sum of the temperatures that are greater than 0C for 1 year.<br>
<br>The way I was trying to do this was with a loop, and to first set any temperature less than or equal to zero, just to zero. Next I want to sum all the temperatures over 365 days at each lat & lon point. If there is an easier way to do this using NCL functions, that would be great.<br>
<br>The error message that I am getting is that the ?loop must be scalar". However, I?m not sure how to specify that I want the temperature value to be analyzed when tasC is (time,lat,lon).<br><br>I?ve uploaded the file "tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc? to the ftp account.<br>
<br>I?m running ncl version 6.1.2, and the system is Darwin Kernel Version 13.3.0<br><br><br>Here is the script:<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>
load ?$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br><br>begin<br><br>; Read in near surface air temperature (Ta)<br><br>Â Â Â Â Â Â Â Â a = addfile("tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc","r")<br>
        tasL = a->tas(0:364,:,:)<br>        tasC = tasL - 273.15 ; Convert K to C<br>        print (tasC)<br>        printVarSummary (tasC)<br><br>        do tasC = 0,364<br>            if (tasC .le. 0)<br>
                H = 0<br>            end if<br>        end do<br><br>        PDD_sum = dim_cumsum_n_Wrap(tasC)<br><br>end<br><br><br>Any help would be greatly appreciated.<br><br>
Thanks,<br><br>Lauren Vargo<br>M.S. Candidate<br>The University of New Mexico<br><br>------------------------------<br><br>Message: 7<br>Date: Mon, 4 Aug 2014 10:11:24 -0600<br>From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Using a loop to find values greater than 0<br>Â Â Â Â (calculate PDD)<br>To: Lauren Jean Vargo <<a href="mailto:lvargo@unm.edu" target="_blank">lvargo@unm.edu</a>><br>Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <CAOF1d_5=0d_Fd0HP_1aOnzafkOmiScwouqZhZbD=<a href="mailto:dSzzozwDbw@mail.gmail.com" target="_blank">dSzzozwDbw@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>
The variable, tasC, is an *array *<br><br>Â Â tasC = tasL - 273.15 ; Convert K to C<br>Â Â printVarSummary(tasC)<br><br>You can not use an array as a loop variable ... in any language.<br><br>Â do tasC = 0,364<br><br>could be, say,<br>
<br>Â do day = 0,364<br><br>See:<br><a href="https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml#Loops" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml#Loops</a><br>
<br>Specifically:<br><br>Â do *loop_identifier* = *scalar_start_expr* ,<br><br>*end_expr*<br><br>*T*he "scalar_start_expr , end_expr" are scalar.<br><br><br><br>===<br><br>Â Â Â Â Â Â Â Â Â a = addfile("tas_day_CCSM4_lgm_<br>
<a href="http://r2i1p1_18700101-19001231.nc/" target="_blank">r2i1p1_18700101-19001231.nc</a>","r")<br>Â Â Â Â Â Â Â Â tasL Â = a->tas(0:364,:,:)<br><br>Â Â Â Â Â Â Â Â TIME = cd_calendar(tasL&time,0) Â ; TIME(ntim,6)<br>
        print(TIME)<br><br>        tasC = tasL - 273.15 ; Convert K to C<br>        tasC = where(tasC.gt.0.0, tasc, 0.0)<br>        pdd  = dim_sum_n_Wrap(tasC, 0)    ; (lat,lon)<br>        pdd@long_name = "degree days for year="+toint(TIME(0,0))<br>
        pdd@units = "degC"<br><br><br><br><br>On Mon, Aug 4, 2014 at 9:16 AM, Lauren Jean Vargo <<a href="mailto:lvargo@unm.edu" target="_blank">lvargo@unm.edu</a>> wrote:<br><br>> Hello,<br>><br>
> I?m having trouble doing a calculation with NCL, ultimately I am trying to<br>> calculate the Positive degree-day (PDD) sum. I?ve read in daily temperature<br>> data (which is 3D [time,lat,lon]). What I am trying to do is calculate the<br>
> sum of the temperatures that are greater than 0C for 1 year.<br>><br>> The way I was trying to do this was with a loop, and to first set any<br>> temperature less than or equal to zero, just to zero. Next I want to sum<br>
> all the temperatures over 365 days at each lat & lon point. If there is an<br>> easier way to do this using NCL functions, that would be great.<br>><br>> The error message that I am getting is that the ?loop must be scalar".<br>
> However, I?m not sure how to specify that I want the temperature value to<br>> be analyzed when tasC is (time,lat,lon).<br>><br>> I?ve uploaded the file "tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc? to<br>
> the ftp account.<br>><br>> I?m running ncl version 6.1.2, and the system is Darwin Kernel Version<br>> 13.3.0<br>><br>><br>> Here is the script:<br>><br>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>> load ?$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>
><br>> begin<br>><br>> ; Read in near surface air temperature (Ta)<br>><br>> Â Â Â Â Â Â Â Â a =<br>> addfile("tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc","r")<br>> Â Â Â Â Â Â Â Â tasL = a->tas(0:364,:,:)<br>
> Â Â Â Â Â Â Â Â tasC = tasL - 273.15 ; Convert K to C<br>> Â Â Â Â Â Â Â Â print (tasC)<br>> Â Â Â Â Â Â Â Â printVarSummary (tasC)<br>><br>> Â Â Â Â Â Â Â Â do tasC = 0,364<br>> Â Â Â Â Â Â Â Â Â Â Â Â if (tasC .le. 0)<br>
> Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â H = 0<br>> Â Â Â Â Â Â Â Â Â Â Â Â end if<br>> Â Â Â Â Â Â Â Â end do<br>><br>> Â Â Â Â Â Â Â Â PDD_sum = dim_cumsum_n_Wrap(tasC)<br>><br>> end<br>><br>><br>> Any help would be greatly appreciated.<br>
><br>> Thanks,<br>><br>> Lauren Vargo<br>> M.S. Candidate<br>> The University of New Mexico<br>> _______________________________________________<br>> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/571eb834/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/571eb834/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 8<br>Date: Mon, 4 Aug 2014 11:36:03 -0500<br>From: Kyle Griffin <<a href="mailto:ksgriffin2@wisc.edu" target="_blank">ksgriffin2@wisc.edu</a>><br>Subject: Re: [ncl-talk] Using a loop to find values greater than 0<br>
    (calculate PDD)<br>To: Lauren Jean Vargo <<a href="mailto:lvargo@unm.edu" target="_blank">lvargo@unm.edu</a>><br>Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <CAKt9o=gyy9Wtwyuv3J32VTubo9zrZVmFFh5LAuN5oju5F6=<a href="mailto:8bg@mail.gmail.com" target="_blank">8bg@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>This also sounds like an ideal spot for implementing the 'where' function.<br>
<br><a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/where.shtml" target="_blank">https://www.ncl.ucar.edu/Document/Functions/Built-in/where.shtml</a><br><br>pos_dd = where(tasC.ge.0,tasC,tasC@_FillValue)<br>
total_pos_dd = dim_sum_n_Wrap(pos_dd,0)<br><br>This will assign the missing value of tasC to any location in the array<br>that is greater than or equal to zero and have each point less than zero<br>assigned to the missing value (usually in the _FillValue attribute. If<br>
tasC@_FillValue<br>is not set, you can set it to anything). You can then add up the values<br>over the 0th dimension (time) to get a total at each lat/lon point.<br><br>Hope that helps,<br><br><br>Kyle<br><br>----------------------------------------<br>
Kyle S. Griffin<br>Department of Atmospheric and Oceanic Sciences<br>University of Wisconsin - Madison<br>Room 1421<br>1225 W Dayton St, Madison, WI 53706<br>Email: <a href="mailto:ksgriffin2@wisc.edu" target="_blank">ksgriffin2@wisc.edu</a><br>
<br><br>On Mon, Aug 4, 2014 at 11:11 AM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br><br>> The variable, tasC, is an *array *<br>><br>> Â Â tasC = tasL - 273.15 ; Convert K to C<br>
> Â Â printVarSummary(tasC)<br>><br>> You can not use an array as a loop variable ... in any language.<br>><br>> Â do tasC = 0,364<br>><br>> could be, say,<br>><br>> Â do day = 0,364<br>><br>> See:<br>
><br>> <a href="https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml#Loops" target="_blank">https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclStatements.shtml#Loops</a><br>><br>> Specifically:<br>
><br>> Â do *loop_identifier* = *scalar_start_expr* ,<br>><br>> *end_expr*<br>><br>> *T*he "scalar_start_expr , end_expr" are scalar.<br>><br>><br>><br>> ===<br>><br>> Â Â Â Â Â Â Â Â Â a = addfile("tas_day_CCSM4_lgm_<br>
> <a href="http://r2i1p1_18700101-19001231.nc/" target="_blank">r2i1p1_18700101-19001231.nc</a>","r")<br>> Â Â Â Â Â Â Â Â tasL Â = a->tas(0:364,:,:)<br>><br>> Â Â Â Â Â Â Â Â TIME = cd_calendar(tasL&time,0) Â ; TIME(ntim,6)<br>
>         print(TIME)<br>><br>>         tasC = tasL - 273.15 ; Convert K to C<br>>         tasC = where(tasC.gt.0.0, tasc, 0.0)<br>>         pdd  = dim_sum_n_Wrap(tasC, 0)    ; (lat,lon)<br>
> Â Â Â Â Â Â Â Â pdd@long_name = "degree days for year="+toint(TIME(0,0))<br>> Â Â Â Â Â Â Â Â pdd@units = "degC"<br>><br>><br>><br>><br>> On Mon, Aug 4, 2014 at 9:16 AM, Lauren Jean Vargo <<a href="mailto:lvargo@unm.edu" target="_blank">lvargo@unm.edu</a>> wrote:<br>
><br>>> Hello,<br>>><br>>> I?m having trouble doing a calculation with NCL, ultimately I am trying<br>>> to calculate the Positive degree-day (PDD) sum. I?ve read in daily<br>>> temperature data (which is 3D [time,lat,lon]). What I am trying to do is<br>
>> calculate the sum of the temperatures that are greater than 0C for 1 year.<br>>><br>>> The way I was trying to do this was with a loop, and to first set any<br>>> temperature less than or equal to zero, just to zero. Next I want to sum<br>
>> all the temperatures over 365 days at each lat & lon point. If there is an<br>>> easier way to do this using NCL functions, that would be great.<br>>><br>>> The error message that I am getting is that the ?loop must be scalar".<br>
>> However, I?m not sure how to specify that I want the temperature value to<br>>> be analyzed when tasC is (time,lat,lon).<br>>><br>>> I?ve uploaded the file "tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc? to<br>
>> the ftp account.<br>>><br>>> I?m running ncl version 6.1.2, and the system is Darwin Kernel Version<br>>> 13.3.0<br>>><br>>><br>>> Here is the script:<br>>><br>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>>> load ?$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>
>><br>>> begin<br>>><br>>> ; Read in near surface air temperature (Ta)<br>>><br>>> Â Â Â Â Â Â Â Â a =<br>>> addfile("tas_day_CCSM4_lgm_r2i1p1_18700101-19001231.nc","r")<br>
>> Â Â Â Â Â Â Â Â tasL = a->tas(0:364,:,:)<br>>> Â Â Â Â Â Â Â Â tasC = tasL - 273.15 ; Convert K to C<br>>> Â Â Â Â Â Â Â Â print (tasC)<br>>> Â Â Â Â Â Â Â Â printVarSummary (tasC)<br>>><br>
>> Â Â Â Â Â Â Â Â do tasC = 0,364<br>>> Â Â Â Â Â Â Â Â Â Â Â Â if (tasC .le. 0)<br>>> Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â H = 0<br>>> Â Â Â Â Â Â Â Â Â Â Â Â end if<br>>> Â Â Â Â Â Â Â Â end do<br>
>><br>>> Â Â Â Â Â Â Â Â PDD_sum = dim_cumsum_n_Wrap(tasC)<br>>><br>>> end<br>>><br>>><br>>> Any help would be greatly appreciated.<br>>><br>>> Thanks,<br>>><br>
>> Lauren Vargo<br>>> M.S. Candidate<br>>> The University of New Mexico<br>>> _______________________________________________<br>>> ncl-talk mailing list<br>>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>>><br>><br>><br>> _______________________________________________<br>
> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/17872935/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/17872935/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 9<br>Date: Mon, 4 Aug 2014 10:52:57 -0600<br>From: Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>Subject: Re: [ncl-talk] Strange white lines across plot<br>
To: Michael Notaro <<a href="mailto:mnotaro@wisc.edu" target="_blank">mnotaro@wisc.edu</a>><br>Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>Â Â Â Â <<a href="mailto:CACNN_CJk8ojmqz8L3oqZiTw2ERjJF6S3xmwprTO--fT7jtnTSQ@mail.gmail.com" target="_blank">CACNN_CJk8ojmqz8L3oqZiTw2ERjJF6S3xmwprTO--fT7jtnTSQ@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br>
<br>Hi Michael,<br><br>I really don't have any idea about this one. Â We did fix a bug in V6.1.2<br>that caused thin white lines to show up, but this doesn't look like a case<br>of that.<br><br>Can you provide your VAR3 variable, say in a NetCDF file, so we can try<br>
running your script here?<br><br>--Mary<br><br><br>On Fri, Aug 1, 2014 at 4:10 PM, Michael Notaro <<a href="mailto:mnotaro@wisc.edu" target="_blank">mnotaro@wisc.edu</a>> wrote:<br><br>><br>> I am using NCL 6.1.2. Please see the attached plot<br>
><br>> and note the odd white lines across the figure. Any<br>><br>> idea what caused it and how to remove it?<br>><br>><br>><br>><br>> The plotting lines are listed below.<br>><br>><br>><br>
><br>> If I change the data values in the plot to some<br>><br>> arbitrary value, sometimes the lines go away.<br>><br>><br>><br>><br>> Michael<br>><br>><br>><br>><br>><br>><br>
><br>> wks=gsn_open_wks("ps","bar_revised2")<br>><br>> gsn_define_colormap(wks, "NCV_blue_red")<br>><br>><br>><br>><br>> res = True<br>><br>> res@gsnDraw = False<br>
><br>> res@gsnFrame = False<br>><br>> res@cnRasterModeOn=True<br>><br>> res@cnFillDrawOrder = "PreDraw"<br>><br>> res@tmXTOn = False<br>><br>> res@tmXTBorderOn = False<br>><br>> res@tmYROn = False<br>
><br>> res@tmYRBorderOn = False<br>><br>> res@tmXBOn = False<br>><br>> res@tmXBBorderOn = False<br>><br>> res@tmYLOn = False<br>><br>> res@tmYLBorderOn = False<br>><br>> res@lbLabelFont = 25<br>
><br>> res@lbLabelFontHeightF = 0.012<br>><br>> res@lbBottomMarginF = 0.5<br>><br>> res@lbTopMarginF = 0.25<br>><br>> res@cnFillOn = True<br>><br>> res@cnLinesOn = False<br>><br>> res@cnLevelSelectionMode = "ManualLevels"<br>
><br>><br>><br>><br>> lower=int2flt(floattointeger(min(VAR3)*100.))/100.<br>><br>> upper=int2flt(floattointeger(max(VAR3)*100.))/100.<br>><br>> delta=fabs(lower-upper)<br>><br>><br>><br>
><br>> res@cnLevelSpacingF=0.01<br>><br>><br>><br>><br>> if ( (delta.ge.0.).and.(delta.lt.0.0025) ) then<br>><br>> Â res@cnLevelSpacingF=0.00001<br>><br>> end if<br>><br>> if ( (delta.ge.0.0025).and.(delta.lt.0.025) ) then<br>
><br>> Â res@cnLevelSpacingF=0.0001<br>><br>> end if<br>><br>> if ( (delta.ge.0.025).and.(delta.lt.0.25) ) then<br>><br>> Â res@cnLevelSpacingF=0.001<br>><br>> end if<br>><br>> if ( (delta.ge.0.25).and.(delta.lt.2.5) ) then<br>
><br>> Â res@cnLevelSpacingF=0.01<br>><br>> end if<br>><br>> if ( (delta.ge.2.5).and.(delta.lt.25.0) ) then<br>><br>> Â res@cnLevelSpacingF=0.1<br>><br>> end if<br>><br>> if ( (delta.ge.25.).and.(delta.lt.250) ) then<br>
><br>> Â res@cnLevelSpacingF=1.0<br>><br>> end if<br>><br>> if ( (delta.ge.250.).and.(delta.lt.2500.) ) then<br>><br>> Â res@cnLevelSpacingF=10.0<br>><br>> end if<br>><br>><br>><br>><br>
> res@gsnSpreadColors=True<br>><br>><br>><br>><br>> if ( (lower.ge.0.).and.(upper.ge.0.) ) then<br>><br>> Â res@gsnSpreadColorStart=120<br>><br>> Â res@gsnSpreadColorEnd=2<br>><br>> end if<br>
><br>> if ( (lower.le.0.).and.(upper.le.0.) ) then<br>><br>> Â res@gsnSpreadColorStart=255<br>><br>> Â res@gsnSpreadColorEnd=137<br>><br>> end if<br>><br>><br>><br>><br>> if ( (lower.lt.0.).and.(upper.gt.0.) ) then<br>
><br>> Â ratio=fabs(lower)/fabs(upper)<br>><br>> Â res@gsnSpreadColorStart=255<br>><br>> Â res@gsnSpreadColorEnd=128-floattointeger(128./ratio)<br>><br>> end if<br>><br>><br>><br>><br>> res@cnMinLevelValF=lower<br>
><br>> res@cnMaxLevelValF=upper<br>><br>><br>><br>><br>> res@lbBoxLinesOn=False<br>><br>><br>><br>><br>><br>><br>> plot = gsn_csm_contour(wks, VAR3, res)<br>><br>> _______________________________________________<br>
> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/8f82b7fe/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/8f82b7fe/attachment-0001.html</a><br>
<br>------------------------------<br><br>Message: 10<br>Date: Mon, 4 Aug 2014 10:58:51 -0600<br>From: Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>Subject: Re: [ncl-talk] Fwd: Masking a U.S. county from a netcdf file<br>
To: Mohammadhassan Mohegh <<a href="mailto:mohegh@usc.edu" target="_blank">mohegh@usc.edu</a>><br>Cc: ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>Message-ID:<br>Â Â Â Â <CACNN_CLmsqp6fP8Q50GDM=<a href="mailto:ebG%2B4RZhFgugzK8c04fHkeExH-3A@mail.gmail.com" target="_blank">ebG+4RZhFgugzK8c04fHkeExH-3A@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br><br>Hi,<br><br>You can use a combination of shapefiles and the "gc_inout" function in<br>order to collect values for a given county and get an average.<br><br>
My suggestion is to use the USA_adm2.shp shapefile downloaded from<br><a href="http://gadm.org/country" target="_blank">gadm.org/country</a>, and then look at example "shapefiles_13.ncl" at:<br><br><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex13" target="_blank">http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex13</a><br>
<br>You can also look at the example right after that, which does the same<br>thing for the coastal counties of Germany, "<br>Germany_coastal_counties_DEU_adm.ncl<br><<a href="http://www.ncl.ucar.edu/Applications/Scripts/Germany_coastal_counties_DEU_adm.ncl" target="_blank">http://www.ncl.ucar.edu/Applications/Scripts/Germany_coastal_counties_DEU_adm.ncl</a>><br>
".<br><br>The "shapefiles_13.ncl" example only does the counties of Georgia. Â Let me<br>know if you need help extending it to do all the counties of the US. It<br>would help, however, if you could provide me with your two WRF output files.<br>
<br>--Mary<br><br><br><br>On Thu, Jul 31, 2014 at 5:09 PM, Mohammadhassan Mohegh <<a href="mailto:mohegh@usc.edu" target="_blank">mohegh@usc.edu</a>><br>wrote:<br><br>><br>> Hi<br>><br>> I'm working with WRF model output and I need to average the differences of<br>
> 2 outputs within a county. Is there any way to do this? Are there<br>> shapefiles or netcdf files available for US counties that NCL recognizes?<br>><br>> I have created a netcdf file containing these counties using ARCGIS. The<br>
> problem is that it does not have the same grid and I can't use the mask()<br>> function with it. Does any one know of a<br>><br>> Thanks in advance for your help.<br>><br>><br>><br>> --<br>> Mohammadhassan Mohegh, PhD student<br>
> Astani Dept of Civil and Environmental Engineering<br>> University of Southern California<br>> 3620 South Vermont Ave, KAP227<br>> Los Angeles, CA<br>> <a href="mailto:mohegh@usc.edu" target="_blank">mohegh@usc.edu</a><br>
><br>> _______________________________________________<br>> ncl-talk mailing list<br>> List instructions, subscriber options, unsubscribe:<br>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/8992036d/attachment-0001.html" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140804/8992036d/attachment-0001.html</a><br>
<br>------------------------------<br><br>_______________________________________________<br>ncl-talk mailing list<br><a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br><a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br><br>End of ncl-talk Digest, Vol 129, Issue 4<br>****************************************<br></blockquote></div><br><br clear="all">
<div><br></div>-- <br>
<div dir="ltr"><font style="FONT-FAMILY:arial,helvetica,sans-serif;COLOR:rgb(0,0,0)"><span>With regards</span><span></span></font><font style="FONT-FAMILY:arial,helvetica,sans-serif"><span style="COLOR:rgb(0,0,0)"><span><br>
</span></span></font><font style="FONT-FAMILY:arial,helvetica,sans-serif"><span style="COLOR:rgb(0,0,0)"><span>Saurabh Kumar singh</span><span></span> </span></font><br style="COLOR:rgb(51,204,0)"><br><font style="COLOR:rgb(102,102,102)" size="1">Â <b><i><span style="FONT-FAMILY:Webdings" lang="EN-GB">P</span></i></b><b><i><span style="FONT-FAMILY:Tahoma,sans-serif" lang="EN-GB"> :Â </span></i></b><b><i><span style="FONT-FAMILY:Arial,sans-serif" lang="EN-GB"> </span></i></b><b><i><span style="FONT-FAMILY:Tahoma,sans-serif" lang="EN-GB">Please consider the environment before printing this e-mail</span></i></b> </font><br style="COLOR:rgb(102,102,102)">
<br><br><br></div></div></div></div>
<br>
ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼ï¼
<br>
å称： ä¸å›½æ°”象局上海å°é£Žç ”究所
<br>
地å€ï¼š ä¸å›½ä¸Šæµ·å¸‚å¾æ±‡åŒºè’²è¥¿è·¯166å·
<br>
电è¯ï¼š<a href="tel:%EF%BC%88%EF%BC%8B86%EF%BC%89021-64391966" value="+862164391966" target="_blank">（＋86）021-64391966</a>
<br>
ä¼ çœŸï¼š<a href="tel:%EF%BC%88%EF%BC%8B86%EF%BC%89021-64391966" value="+862164391966" target="_blank">（＋86）021-64391966</a>
<br>
============================
</div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>