<div dir="ltr"><div>The following is untested. However, it gives you the idea.</div><div>It is your responsibility to look up documentation and look at the printed output.</div><div><br></div><div><br></div><div><div><br></div><div> fils = systemfunc ("ls precip.*.nc") ; file paths<br> a = addfiles (fils, "r")<br> printVarSummary (a)<br><br>;************************************************<br>; Read the file<br>;************************************************<br><br> ymdStrt = 19790101 ; climatology start year/month/day<br> ymdLast = 20151231 ; " last year/month/day <br> time = a[:]->time<br> yyyymmdd= cd_calendar(time, -2) ; note -2<br>print(yyyymmdd)<br> ntStrt = ind(yyyymmdd.eq.ymdStrt) ; start time index<br> ntLast = ind(yyyymmdd.eq.ymdLast) ; last time index<br>print("ntStrt="+ntStrt+" ntLast="+ntLast)<br> prc = a[:]->precip(ntStrt:ntLast,:,:) ; read precip from all files<br><br> printVarSummary (prc)<br> printMinMax (prc,0)<br></div><div> print("=====")</div><div> ; you decide what you need<br></div><div> prcMonth = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml"><b>calculate_monthly_values(</b></a>prc, "avg", 0, False) ; monthly mean<br></div><div> ;;prcMonth = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml"><b>calculate_monthly_values(</b></a>prc, "sum", 0, False) ; monthly total<br></div><div><br></div><div> printVarSummary (prcMonth)<br> printMinMax (prcMonth,0)<br> print("=====")</div><div> opt = True ; examine data distribute [information only]<br> opt@PrintStat = True<br> stat_prc = stat_dispersion(prcMonth, opt )<br><br> prcClm = <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml"><b>clmMonTLL</b></a>(prcMonth) ; Compute monthly climatology<br> printVarSummary (prcClm)<br> printMinMax (prcClm,0)</div><div><br></div><div>====</div><div>Add the plot code.</div><div>Look carefully at: <a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml</a></div><div><br></div><div>EG: The following may be inappropriate for your dataset but this is the idea.<br></div><div><pre> res@cnLevelSelectionMode = "ExplicitLevels"
res@cnLevels = (/ 0.01, 0.02, 0.04, 0.08, 0.16, \
0.32, 0.64, 0.96/)
res@cnFillColors = (/"white","cyan", "green","yellow",\
"darkorange","red","magenta","purple",\
"black"/)
</pre></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Oct 27, 2019 at 2:48 PM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>[0]<br></div><div>As noted on ncl-talk many time the first rule of data processing is *Look at your data!"</div><div>The temperature data were <b>monthly means</b> [K]. <br></div><div>The precipitation data are<b> daily</b> totals [mm/day]. How do I know? I looked at the README file at the URL. You should do that.<br></div><div><br></div><div>[1]</div><div>a =<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" target="_blank"><b> addfiles</b></a> (fils, "r")</div><div><br></div><div>You have the following which is appropriate for an <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" target="_blank"><b>addfile</b></a> reference</div><div> time = a->time </div><div><br></div><div>However, you are using <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" target="_blank"><b>addfiles</b></a> It should be:<br></div><div> time = a<b>[:]</b>->time ; <b>[:]</b> syntax is appropriate for an <b>addfiles</b> reference</div><div><br></div><div>[2] <br></div><div>Did you read the documentation for <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml" target="_blank"><b>cd_calendar </b></a>as I suggested? The whole point of using <b>cd_calendar</b> is to return 'human-readable' start and end times from which the actual start/stop time index values could be determined.. Your script uses the actual time values [ 692496 , 1016808 ] mixed with human-readable yyyymm. Further, you used <br></div><div><br></div><div> yyyymm = cd_calendar(time, <b>-1</b>) ; -1 returns YYYYMM<br></div><div><br></div><div>The precipitation files contain [as noted] daily values. You must use a human readable time like YYYYMMDD [eg: 19790101]. The us</div><div><br></div><div> yyyymmdd = cd_calendar(time, <b>-2</b>) ; -2 returns YYYYMMDD</div><div><br></div><div>This converts the actual time values [eg: 692496] to 'human-readable' yyyymmdd [year=>yyyy, month=>mm, day=dd><b> 692496=>197901</b>01] values. </div><div><br></div><div><br></div><div>===</div><div>Also, your script will need to be changed to accommodate the precipitation variable plot values <br></div><div>There are many examples of how to do this.</div><div><br></div><div>eg: <a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Applications/HiResPrc.shtml</b></a></div><div>===<br></div><div>In a previous email, you mentioned that you converted degrees Kelvin to C.</div><div>It is best to self-document the variable with its new units.<br></div><div><br></div><div> tmp = tmp-273.15</div><div> tmp@units = "degC" ; this will be automatically plotted <br></div><div>===</div><div>Good Luck<br></div><div><br></div><div><b><br></b></div><div><b><br></b></div><div><b><br></b></div><div><b></b><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 26, 2019 at 8:19 AM zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr">Dear all,<div>I would like to have climatology precipitation like temperature climatology discussed through previous emails. However, there are precipitation files for each year from 1979- 2015 (<a href="ftp://ftp.cdc.noaa.gov/Datasets/cpc_global_precip/" target="_blank">ftp://ftp.cdc.noaa.gov/Datasets/cpc_global_precip/</a>). In the first step, I tried to merge all those files as one, and then use Dennis script to plot a map. I faced some errors and need your advice please. the script I am running has been attached.<span style="color:rgb(80,0,80)">Many thanks in advance,</span></div></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr"><span style="color:rgb(80,0,80)">Best regards</span>n Sat, Oct 26, 2019 at 3:20 PM zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I just noticed that I could use <span style="color:rgb(0,0,0);white-space:pre-wrap">res@gsnRightString = " deg C" .</span><div><font color="#000000"><span style="white-space:pre-wrap">Once again, many thanks for your guide.</span></font></div><div><font color="#000000"><span style="white-space:pre-wrap">Best wishes,<br></span></font><div><span style="color:rgb(0,0,0);white-space:pre-wrap"><br></span></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 25, 2019 at 11:07 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>I forgot to mention that plots can be done in 'portrait' [default] or 'landscape' mode.</div><div><br></div><div>You should experiment to see which is appropriate for your needs.<br></div><div><br></div><div>resP@gsnPaperOrientation = "landscape" ; "portrait" is default</div><div><br></div><div> do nmo=0,11 ; loop over the months<br> res@gsnLeftString = months(nmo)<br> plot(nmo) = gsn_csm_contour_map(wks,tmpClm(nmo,:,:), res) ; create plot<br> end do<br><br> if (resP@gsnPaperOrientation .eq."landscape") then<br> gsn_panel(wks,plot,(/3,4/),resP) ; 3 rows x 4 columns<br> else<br> gsn_panel(wks,plot,(/4,3/),resP) ; portrait ; 4 rows x 3 columns<br> end if<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 25, 2019 at 12:15 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hello,</div><div><br></div><div>We do not usually do this but this is offline from ncl-talk.<br></div><div><br></div><div>[1] <br></div><div>Conventional subscripts start at 1 for Fortran and Matlab. Hence, for 12 months [nmos=12]</div><div> Fortran: do nmo=1,nmos <br></div><div> Matlab: for 1:nmos or 1:1:nmos [I am not a Matlab user so this is a guess.]<br></div><div>NCL 'conventional' subscripts start at 0 [like C/C++/IDL/Python]</div><div> NCL: do nmo=0,ntim-1</div><div><br></div><div>[2] Dave illustrated 'coordinate subscripting' which uses the <span style="color:rgb(0,0,255)"><b>{</b></span>...<span style="color:rgb(0,0,255)"><b>}</b></span> syntax.</div><div><br></div><div>The script I am attaching uses conventional subscripting. Just like Matlab/Fortran except it uses a range fro 0 to 11 [12 elements].</div><div><br></div><div>[3] Please carefully read the documentation for <a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml" target="_blank"><b>clmMonTLL</b></a>, <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml" target="_blank"><b>cd_calendar</b></a>,<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" target="_blank"><b> ind</b></a></div><div><br></div><div>[4] Look at the output from '<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml" target="_blank"><b>printVarSummary</b></a>'. USE <b>printVarSummary </b>frequently.<br></div><div> Note the dimensions and added attributes: <br></div><div> time_op_ncl : Climatology: 29 years<br> info : function clmMonTLL: contributed.ncl</div><div><br></div><div>[5] NCL offer a large number of color tables [palettes]:</div><div> <a href="http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml</b></a></div><div> Using a color map that emphasizes features can be very useful.</div><div><br></div><div>[6] You must invest the time to learn any new language. <br></div><div>Karin Meier-Fliesher and Michael Bottinger [DKRZ] wrote a wonderful tutorial. <br></div><div><br></div><div><a href="http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/" target="_blank"><b>http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/</b></a></div><div><br></div><div>I suggest you read it.</div><div><br></div><div>Good luck</div><div>==============<br></div><div><b>%></b> ncl zoe_jacobs.hgcn_cams.ncl</div><div><br></div><div>will produce a png file.<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Oct 25, 2019 at 11:16 AM Dave Allured - NOAA Affiliate via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>1. In your example, please notice that the month loop has a stride of 3. Therefore it is making plots for only four months: January, April, July, and October.</div><div><br></div><div>If you want all 12 months, remove the stride. Then make sure the array "plot" is dimensioned 12. I think you can interchangeably use either a 1-D or 2-D plot array (i.e. 3*4), your choice.</div><div><br></div><div>The secondary index "i" is used to write the four plots into positions 0, 1, 2, 3 in the "plot" graphics array. If you are making 12 plots, you do not really need to use a secondary index.</div><div><br></div><div>Panel plotting uses a graphics array containing multiple plots. Please study basic examples and documentation for making panel plots.</div><div><br></div><div>2. That data file has full coordinates. Therefore you can use either conventional or coordinate subscripting, or mixed, your choice.</div><div><br></div>I suggest using one of the date functions with coordinate subscripting, to index the time subset that you want. Something like this:<br><br> time_units = f->air&time@units<br> time1 = cd_inv_calendar (year1, 1, 1, 0, 0, 0, time_units, 0)<br> time2 = cd_inv_calendar (year2, 12, 31, 23, 0, 0, time_units, 0)</div><div dir="ltr"> air_subset = f->air({time1:time2},:,:)<br><div><br></div><div>Use printVarSummary to ensure that the subset has the dimension sizes that you expect. Then you can proceed to compute the climatology.</div><div><br></div><div><br></div><div>On Fri, Oct 25, 2019 at 10:14 AM zoe jacobs via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><font face="arial, sans-serif">Hi all,</font><div><font face="arial, sans-serif">regarding climo_3.ncl (<a href="https://www.ncl.ucar.edu/Applications/Scripts/climo_3.ncl" target="_blank">https://www.ncl.ucar.edu/Applications/Scripts/climo_3.ncl</a>) ,</font></div><div><font face="arial, sans-serif">I have 2 questions:</font></div><div><font face="arial, sans-serif">1. I need to plot all months on one panel (say 3*4), and cannot understand the logic behind the below loop, which used in the climo_3.ncl script :</font></div><div><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font style="background-color:rgb(241,194,50)" face="arial, sans-serif">i = -1 ; Climatologies
do nmo=0,11,3 ; loop over the months
i = i+1
res@gsnCenterString = months(nmo)+":"+time(0)/100 +"-"+ time(ntim-1)/100
plot(i) = gsn_csm_contour_map(wks,prcClm(nmo,:,:), res) ; create plot
end do</font></pre><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font style="background-color:rgb(241,194,50)" face="arial, sans-serif">How does it work??/</font></pre><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font face="arial, sans-serif">2. I would like to show climotology temperature from 1987- 2015. Data which I am using is<span style="background-color:rgb(180,167,214)"> <a href="http://air.mon.mean.nc" target="_blank">air.mon.mean.nc</a> </span><span style="background-color:rgb(255,255,255)"> ( </span><a href="https://www.esrl.noaa.gov/psd/data/gridded/data.ghcncams.html" target="_blank">https://www.esrl.noaa.gov/psd/data/gridded/data.ghcncams.html</a>) <span style="background-color:rgb(255,255,255)"> and it is </span>using conventional subscripts. I am not familiar with conventional subscripts. So how can I convert coordinate subscripting to conventional subscripts?!!</font></pre><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font face="arial, sans-serif">Please kindly advice me .</font></pre><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font face="arial, sans-serif">Many thanks in advance,</font></pre><pre style="color:rgb(0,0,0);white-space:pre-wrap"><font face="arial, sans-serif">Best regards,</font></pre></div></div></blockquote></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>
</blockquote></div>