<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 border="0" cellpadding="0" cellspacing="0" width="167" style="border-collapse:collapse;width:125pt">
<tbody><tr height="20" style="height:15pt">
<td height="20" width="167" style="height:15pt;width:125pt">((lat+49.875)/0.25)+1<br><br>trmmlong=<br><table border="0" cellpadding="0" cellspacing="0" width="182" style="border-collapse:collapse;width:137pt"><tbody><tr height="20" style="height:15pt">
<td height="20" width="182" style="height:15pt;width:137pt">((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 class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send ncl-talk mailing list submissions to<br>
<a href="mailto:ncl-talk@ucar.edu">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">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">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">shea@ucar.edu</a>><br>
Subject: Re: [ncl-talk] triple2grid<br>
To: Manisha Ganeshan <<a href="mailto:manisha.ganeshan86@gmail.com">manisha.ganeshan86@gmail.com</a>><br>
Cc: NCL-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CAOF1d_4pLjOOnFbXb-7ju2BQbkXwbx=iCVsuVWjiUM=<a href="mailto:wyR_Efw@mail.gmail.com">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">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">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">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">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">ipmajhi@alaska.edu</a>><br>
Cc: NCL USERS <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CACNN_CKd2OjK7wbauhLmC01QPRpsCMfMq2aYC1CDzqh=<a href="mailto:UQeVog@mail.gmail.com">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">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">shea@ucar.edu</a>><br>
Subject: Re: [ncl-talk] (no subject)<br>
To: Modise Wiston <<a href="mailto:modise.wiston@postgrad.manchester.ac.uk">modise.wiston@postgrad.manchester.ac.uk</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CAOF1d_4GEj90MZ04-u1vemt3UbnH%2Bvd--oWt1Ax73BAP2poQag@mail.gmail.com">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">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">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">haley@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Running example scripts for beginners<br>
To: zilore mumba <<a href="mailto:zmumba@yahoo.com">zmumba@yahoo.com</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CACNN_CJySFTTZvaJQbpnCE%2BN3mBi3BncJCF7RxZHxgMmDd41nQ@mail.gmail.com">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">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">shea@ucar.edu</a>><br>
Subject: Re: [ncl-talk] read separated files<br>
To: Yuqiang Zhang <<a href="mailto:yuqiangzhang.thu@gmail.com">yuqiangzhang.thu@gmail.com</a>><br>
Cc: NCL <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CAOF1d_4_aU1aywjH6CZvfAskDg276GTDvL%2B3c6hhZeDPHYM2RQ@mail.gmail.com">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">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">ncl-talk-bounces@ucar.edu</a> [mailto:<a href="mailto:ncl-talk-bounces@ucar.edu">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">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">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID: <<a href="mailto:70D79029-3BC6-4860-8081-E5C280D91F93@unm.edu">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">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">lvargo@unm.edu</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CAOF1d_5=0d_Fd0HP_1aOnzafkOmiScwouqZhZbD=<a href="mailto:dSzzozwDbw@mail.gmail.com">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">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">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">lvargo@unm.edu</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CAKt9o=gyy9Wtwyuv3J32VTubo9zrZVmFFh5LAuN5oju5F6=<a href="mailto:8bg@mail.gmail.com">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">ksgriffin2@wisc.edu</a><br>
<br>
<br>
On Mon, Aug 4, 2014 at 11:11 AM, Dennis Shea <<a href="mailto:shea@ucar.edu">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">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">haley@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Strange white lines across plot<br>
To: Michael Notaro <<a href="mailto:mnotaro@wisc.edu">mnotaro@wisc.edu</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CACNN_CJk8ojmqz8L3oqZiTw2ERjJF6S3xmwprTO--fT7jtnTSQ@mail.gmail.com">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">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">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">mohegh@usc.edu</a>><br>
Cc: ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CACNN_CLmsqp6fP8Q50GDM=<a href="mailto:ebG%2B4RZhFgugzK8c04fHkeExH-3A@mail.gmail.com">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">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">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">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="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif"><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>