<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>