<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">qici_era5.nc</a>","r")<br>b= addfile("/home/SSPMRES/hazra/ushnanshu/paper2/data/<a href="http://olr_noaa_1981_2010.nc">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" class="gmail_signature" data-smartmail="gmail_signature"><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>