<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>