<div dir="ltr"><div>I am not familiar with ferret's graphics.</div><div>===</div><div><br></div><div>I speculate that <br></div><div> printVarSummary(clice_clim)</div><div> printVarSummary(clice_clim_sm)</div><div>would be <br></div><div> clice_clim(366)</div><div> clice_clim_sm(366)</div><div><br></div><div>and ...</div><div> printVarSummary(<b> </b>clice_anom )</div><div>would be</div><div> clice_anom(10957)</div><div><br></div><div>All one dimensional. <br></div><div><b> </b></div><div><b> </b></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jun 25, 2021 at 5:17 AM Ushnanshu Dutta via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.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">Hi,<div><br></div><div>I want to save smooth climatology using ncl. Calculation part is okay. But unfortunately it is getting saved in X-axis (as revealed by ferret). I want it on L-axis.</div><div> My code is as follows,</div><div><br></div><div><b>------------------------------------------------------</b></div><div><div><b> 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/diagnostics_cam.ncl" <br><br><br>begin<br><br>a= addfile("<a href="http://qici_era5.nc" target="_blank">qici_era5.nc</a>","r")<br>b= addfile("/home/SSPMRES/hazra/ushnanshu/paper2/data/<a href="http://olr_noaa_1981_2010.nc" target="_blank">olr_noaa_1981_2010.nc</a>","r")<br><br>clice = doubletofloat(a->qici(:,{1},{1}))<br><br>olr = b->OLR(:,{1},{1})<br>dummy = olr(:)<br><br></b></div><div><b> ymdStrt = 19810101 ; start yyyymmdd<br> ymdLast = 20101231 ; last <br> yrStrt = ymdStrt/10000<br> yrLast = ymdLast/10000<br> <br> nhar = 4 ; number of fourier comp <br><br> ;***********************************************************<br>; Read user specified time and create required yyyyddd <br>;***********************************************************<br> <br> time = b->TIME ; all times on file<br> ymd = cd_calendar(time, -2) ; yyyymmdd<br> iStrt = ind(ymd.eq.ymdStrt) ; index start<br> iLast = ind(ymd.eq.ymdLast) ; index last <br> delete(time)<br> delete(ymd)<br><br>;***********************************************************<br>; Read user specified time and create required yyyyddd <br>;***********************************************************<br> time = b->TIME(iStrt:iLast) ; time:units = "hours since"<br> TIME = cd_calendar(time, 0) ; type float <br> year = floattointeger( TIME(:,0) )<br> month = floattointeger( TIME(:,1) )<br> day = floattointeger( TIME(:,2) ) <br> ddd = day_of_year(year, month, day) <br> yyyyddd = year*1000 + ddd ; needed for input<br><br>;***********************************************************<br>; Compute daily climatology: cloud ICE raw and then 'smoothed' <br>;***********************************************************<br> clice_clim = clmDayT(clice, yyyyddd) ; daily climatology at CI box<br> printVarSummary(clice_clim) <br><br>;***********************************************************<br>; Compute smoothed daily climatology using 'nhar' harmonics<br>;***********************************************************<br> clice_clim_sm = smthClmDayT(clice_clim, nhar)<br> printVarSummary(clice_clim_sm)<br><br>;***********************************************************<br>; Compute daily anomalies using smoothed climatologies<br>;***********************************************************<br> <br> clice_anom = calcDayAnomT(clice, yyyyddd, clice_clim_sm) <br> clice_anom@long_name = "cldice Anomalies of CI box "<br><br>clice_clim_sm!0=time<br><br>copy_VarCoords(dummy,clice_anom)<br> printVarSummary(clice_clim_sm)<br><br>;------------------------------------------------------------------------------------<br>;------------------------------------------------------------------------------------<br><br>setfileoption("nc","Format","NetCDF4Classic")<br>;---------Writing the output---------------------<br> diro = "./"<br> filo = "anomalies_iwp_CI_avg_era5.nc"<br> system("/bin/rm -f "+diro+filo) ; rm any pre-existing file<br> ;filedimdef(filo,"tnew",-1,True) <br> f_select = addfile(diro+filo, "c")<br><br> f_select->clice_anom =clice_anom<br> f_select->clice_clim =clice_clim_sm<br>end</b><br></div><div><br></div><div><br></div><div><br></div><div>Output file is aving like this. I want it in the L dimension (i.e time) instead of I dimesnion</div><div> 1> ./anomalies_iwp_CI_avg_era5.nc (default)<br> name title I J K L M N<br> CLICE_ANOM<br> cldice Anomalies of CI box ... ... ... 1:10957 ... ...<br> CLICE_CLIM<br> Daily Climatology 1:366 ... ... ... ... ...<br></div><div><br></div><div><br></div><div><br></div><div>please help</div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><div><b><span style="font-family:garamond,serif">Ushnanshu Dutta</span><br></b></div><b><span style="font-family:garamond,serif">Research Fellow<br></span></b></div><b><span style="font-family:garamond,serif">Indian Institute of Tropical Meteorology<br></span></b></div><b><span style="font-family:garamond,serif">Pune,411008<br></span></b></div><b><span style="font-family:garamond,serif">Mob: 9405022799</span></b><br></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>