<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I could figure the required process out. This time followed the example 4 at <a href="http://www.ncl.ucar.edu/Applications/Scripts/spec_4.ncl">http://www.ncl.ucar.edu/Applications/Scripts/spec_4.ncl</a>. However, I would still like to know how can I plot three spectrums of three colours in one single plot. Your help is always useful.<br><br></div>Thank you for putting up wonderful resources.<br><br></div>Best regards,<br></div>Ruksana <br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 5, 2016 at 5:30 PM, Ruksana Abedin <span dir="ltr"><<a href="mailto:ruksana.abedin@gmail.com" target="_blank">ruksana.abedin@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">
<p class="MsoNormal">Hi,</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I have followed the
<a href="http://www.ncl.ucar.edu/Applications/Scripts/spec_1.ncl" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/Scripts/spec_1.<wbr>ncl</a> script to plot a
spectrum for my data but the specs are very difficult to visualize clearly in
the plot. Data is daily precipitation from model and two observations for 1988-2007. <br></p><p class="MsoNormal"><br></p><p class="MsoNormal">1. How to change the x and y axis scales to make it
more visible. <br></p><p class="MsoNormal">2. Can I have the cycles per day instead of months? How to
change this aspect? <br></p><p class="MsoNormal">3. Can I plot three spectrums for three different files in a panel together?<br></p>
<p class="MsoNormal"><br></p><p class="MsoNormal">The variable summary, script and the resultant plot are attached.
Your time and efforts to help me out is highly appreciated.<span> <br></span></p><p class="MsoNormal"><span><br></span></p><p class="MsoNormal"><span><br></span></p><p class="MsoNormal"><span>Variable: PPT_mod<br>Type: float<br>Total Size: 25056000 bytes<br> 6264000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 7200] x [lat | 30] x [lon | 29]<br>Coordinates:<br> time: [52950.5..60149.5]<br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 10<br> z1 : 0<br> standard_name : precipitation_flux<br> long_name : TOTAL PRECIPITATION RATE KG/M2/S<br> units : kg m-2 s-1<br> _FillValue : 1e+20<br> missing_value : 1e+20<br> stash_item : 216<br> stash_section : 5<br> field_code : 90<br> cell_method : time: mean<br><br>Variable: PPT_obs1<br>Type: float<br>Total Size: 25056000 bytes<br> 6264000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 7200] x [lat | 30] x [lon | 29]<br>Coordinates:<br> time: [19880101.5..20071230.5]<br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br> long_name : daily precipitation analysis interpolated onto 0.5deg grids [mm/day]<br> units : mm/day<br> _FillValue : -99.9<br> missing_value : -99.9<br> level_description : Earth surface<br><br></span></p><p class="MsoNormal"><span>Variable: PPT_obs2<br>Type: float<br>Total Size: 25056000 bytes<br> 6264000 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [time | 7200] x [lat | 30] x [lon | 29]<br>Coordinates:<br> time: [19880101.5..20071230.5]<br> lat: [15.25..29.75]<br> lon: [ 83.. 97]<br>Number Of Attributes: 5<br> long_name : full data daily product version 1 precipitation per grid<br> units : mm/day<br> code : 20<br> _FillValue : -99999.99<br> missing_value : -99999.99</span></p><p class="MsoNormal"><span><br></span></p><p class="MsoNormal"><span><br>;*****************************<wbr>*******************<br>begin<br>;*****************************<wbr>************************ <br>;read in data <br>;*****************************<wbr>************************ <br>f0 = addfile ("/../1WAH_daily_precip_<wbr>b387clim_Bangladesh_100_1988_<wbr><a href="http://2007_land_only.nc">2007_land_only.nc</a>", "r");<br>f1 = addfile("/../APHRO_1988_2007_<wbr><a href="http://correct_daily_precip.nc">correct_daily_precip.nc</a>", "r")<br>f2 = addfile("/../GPCC_1988_2007_<wbr><a href="http://correct_daily_precip.nc">correct_daily_precip.nc</a>","r")<br><br>;*****************************<wbr>************************ <br>; parameters<br>;*****************************<wbr>************************ <br><br>PPT_mod = f0->precip(:,0,:,:)<br>PPT_mod = PPT_mod*60*60*24<br>PPT_obs1 = f1->precip<br>PPT_obs2 = f2->precip<br><br>fldmean0 = wgt_areaave_Wrap(PPT_mod,1.0,<wbr>1.0,1)<br>fldmean1 = wgt_areaave_Wrap(PPT_obs1,1.0,<wbr>1.0,1)<br>fldmean2= wgt_areaave_Wrap(PPT_obs2,1.0,<wbr>1.0,1)<br><br>;printVarSummary(time)<br><br>printVarSummary(PPT_mod)<br>printVarSummary(PPT_obs1)<br>printVarSummary(PPT_obs2)<br><br>;*****************************<wbr>*******************<br>; set function arguments<br>;*****************************<wbr>*******************<br>; detrending opt: 0=>remove mean 1=>remove mean and detrend<br> d = 0<br>; smoothing periodogram: (0 <= sm <= ??.) should be at least 3 and odd<br> sm = 7<br>; percent tapered: (0.0 <= pct <= 1.0) 0.10 common. <br> pct = 0.10<br>;*****************************<wbr>*******************<br>; calculate spectrum<br>;*****************************<wbr>*******************<br> spec = specx_anal(fldmean2,d,sm,pct)<br>;*****************************<wbr>*******************<br>; plotting<br>;*****************************<wbr>*******************<br> wks = gsn_open_wks("eps","spec") <wbr> ; Opens a eps file <br><br> res = True ; no plot mods desired<br> res@tiMainString = "Precipitation" ; title<br> res@tiXAxisString = "Frequency (cycles/month)" ; xaxis<br> res@tiYAxisString = "Variance" ; yaxis<br><br> plot=gsn_csm_xy(wks,spec@frq,<wbr>spec@spcx,res) ; create plot <br>;*****************************<wbr>******************<br>end<br></span></p><p class="MsoNormal"><span><br></span></p>
</div>
</blockquote></div><br></div>