<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">http://www.ncl.ucar.edu/Applications/Scripts/spec_1.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>;************************************************<br>begin<br>;***************************************************** <br>;read in data <br>;***************************************************** <br>f0 = addfile ("/../1WAH_daily_precip_b387clim_Bangladesh_100_1988_2007_land_only.nc", "r");<br>f1 = addfile("/../APHRO_1988_2007_correct_daily_precip.nc", "r")<br>f2 = addfile("/../GPCC_1988_2007_correct_daily_precip.nc","r")<br><br>;***************************************************** <br>; parameters<br>;***************************************************** <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,1.0,1)<br>fldmean1 = wgt_areaave_Wrap(PPT_obs1,1.0,1.0,1)<br>fldmean2= wgt_areaave_Wrap(PPT_obs2,1.0,1.0,1)<br><br>;printVarSummary(time)<br><br>printVarSummary(PPT_mod)<br>printVarSummary(PPT_obs1)<br>printVarSummary(PPT_obs2)<br><br>;************************************************<br>; set function arguments<br>;************************************************<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>;************************************************<br>; calculate spectrum<br>;************************************************<br> spec = specx_anal(fldmean2,d,sm,pct)<br>;************************************************<br>; plotting<br>;************************************************<br> wks = gsn_open_wks("eps","spec") ; 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,spec@spcx,res) ; create plot <br>;***********************************************<br>end<br></span></p><p class="MsoNormal"><span><br></span></p>
</div>