<div dir="ltr"><div>Hi,</div><div>I have computed climatology using smooth climatology function. <br></div><div>While saving the file the time is coming in "lev" dimension instead of time dimension like this</div><div>   1> ./<a href="http://cfsr_clim.nc">cfsr_clim.nc</a>  (default)<br> name     title                             I         J         K         L         M         N<br> <b>QNET     Daily Climatology                1:80      1:60      1:366     ...       ...       ...</b></div><div>I want 1:366 to be in "L"</div><div>any possible wayout?</div><div><b>below is the script</b></div><div><br></div><div><br></div><div>load "/usr/local/ncl/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "/usr/local/ncl/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "/usr/local/ncl/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "/usr/local/ncl/lib/ncarg/nclscripts/csm/shea_util.ncl"<br><br>begin<br><br><br><br>b    = addfile("<a href="http://cfsr_qnet_new.nc">cfsr_qnet_new.nc</a>","r")<br><br>cqnet1= (b->QNET)<br>cqnet1!0="time"<br>cqnet1!1="lat"<br>cqnet1!2="lon"<br><br><br>f1=addfile("/run/media/iitm/VEERU/manuscritp1/final_2019/<a href="http://landsea.nc">landsea.nc</a>","r")<br>lsdata1 = f1->LSMASK<br>;printVarSummary(lsdata)<br>lsm1  = landsea_mask(lsdata1,cqnet1&lat,cqnet1&lon)<br>   cqnet = mask(cqnet1,lsm1.eq.1,False)<br>copy_VarCoords(cqnet1,cqnet)<br><br><br>   ymdStrt = 20000101                         ; start yyyymmdd<br>   ymdLast = 20091231                         ; 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>   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: <br>;***********************************************************<br>   cqnet_clim = clmDayTLL(cqnet, yyyyddd)     ; daily climatology at CI box<br><br>   printVarSummary(cqnet_clim)            <br><br>;***********************************************************<br>; Compute smoothed daily climatology using 'nhar' harmonics<br>;***********************************************************<br>   cqnet_clim_sm = smthClmDayTLL(cqnet_clim, nhar)<br>   printVarSummary(cqnet_clim_sm)<br><br>tim=cqnet_clim_sm&year_day<br><br>printVarSummary(tim)<br>cqnet_clim_sm!0 = "tim"<br><br>setfileoption("nc","Format","NetCDF4Classic")<br><br>fnc          = "<a href="http://cfsr_clim.nc">cfsr_clim.nc</a>"<br>  system("/bin/rm -f "+fnc)               ; rm any pre-existing file<br>  fout         = addfile(fnc, "c")   ; new netCDF file<br>  fout->qnet    = cqnet_clim_sm<br>end</div><div><br>-- <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"></span></b><b><span style="font-family:garamond,serif"></span></b><br></div></div></div></div></div></div></div></div></div></div></div>