<div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif;font-size:small;color:rgb(0,0,0)">Dear NCL users</div><div class="gmail_default" style="font-family:georgia,serif;font-size:small;color:rgb(0,0,0)"><br></div><div class="gmail_default" style="font-family:georgia,serif;font-size:small;color:rgb(0,0,0)">I am using NCL 6.4.0 on ubuntu. I need to write a variable to netcdf file. I did it successfully using Method 1 given on ncl homepage. I have little problem that when I open nc file in grads, the time is not realistic. I mean its like;</div><div class="gmail_default" style="font-family:georgia,serif;font-size:small;color:rgb(0,0,0)"><br></div><div class="gmail_default"><div class="gmail_default" style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small">Time values set: 2320:6:25:0 2320:6:25:0</div><div style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small"> </div><div style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small">while the variable has the time from 190101 to 201312</div><div style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small"><br></div><div style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small">Can someone please help me how can add time dim as that attached with the define variable? My code is given below</div><div style="color:rgb(0,0,0);font-family:georgia,serif;font-size:small"><br></div><div><div><font color="#ff0000">;====================================================</font></div><div><font color="#ff0000">begin</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000"> diri  = "./"</font></div><div><font color="#ff0000"> fili  = "<a href="http://precip.mon.total.v7.nc">precip.mon.total.v7.nc</a>"</font></div><div><font color="#ff0000"> f   = addfile(diri+fili, "r")</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000"> prc  = flt2dble(f->precip(:,{12:33},{30:60}))</font></div><div><font color="#ff0000"> pmsg  = prc@_FillValue   ; convenience</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">;*********************************</font></div><div><font color="#ff0000">; plot parameters</font></div><div><font color="#ff0000">;*********************************</font></div><div><font color="#ff0000"> dimprc = dimsizes(prc)</font></div><div><font color="#ff0000"> ntim  = dimprc(0)</font></div><div><font color="#ff0000"> nlat  = dimprc(1)</font></div><div><font color="#ff0000"> mlon  = dimprc(2)</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000"> tim = f->time</font></div><div><font color="#ff0000"> utc_date = cd_calendar(tim, 0)  </font></div><div><font color="#ff0000"> year  = tointeger(utc_date(:,0))</font></div><div><font color="#ff0000"> yrStrt = year(0)</font></div><div><font color="#ff0000"> yrLast = year(ntim-1)</font></div><div><font color="#ff0000"> nyear  = yrLast-yrStrt+1</font></div><div><font color="#ff0000">;print(year + " " + yrStrt + " "+ yrLast +" "+ nyear)</font></div><div><font color="#ff0000"> yyyymm = yyyymm_time(yrStrt, yrLast, "integer") </font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000"> yrfrac = (/ yyyymm_to_yyyyfrac(yyyymm, 0.0) /)</font></div><div><span style="white-space:pre"><font color="#ff0000"> </font></span></div><div><font color="#ff0000"> prc&time = (/ yyyymm /)</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000">;write to nc file</font></div><div><font color="#ff0000">;================================================================</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">  system("/bin/rm -f <a href="http://prc_3mon.nc">prc_3mon.nc</a>") ;remove any pre-existing file</font></div><div><font color="#ff0000">  ncdf = addfile("<a href="http://prc_3mon.nc">prc_3mon.nc</a>" ,"c") ; open output netCDF file</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000">; create global attributes of the file (optional)</font></div><div><font color="#ff0000">;================================================================</font></div><div><font color="#ff0000">    fAtt        = True    ; assign file attributes</font></div><div><font color="#ff0000">    fAtt@title     = "GPCC Rainfall"</font></div><div><font color="#ff0000">    fAtt@source_file  = "<a href="http://precip.mon.total.v7.nc">precip.mon.total.v7.nc</a>"</font></div><div><font color="#ff0000">    fAtt@Conventions  = "None"</font></div><div><font color="#ff0000">    fAtt@creation_date = systemfunc ("date")</font></div><div><font color="#ff0000">    fileattdef( ncdf, fAtt )      ; copy file attributes</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000">; make time an UNLIMITED dimension;recommended for most applications</font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">    filedimdef(ncdf,"time",-1,True) </font></div><div><font color="#ff0000"><br></font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000">; output variables directly; NCL will call appropriate functions</font></div><div><font color="#ff0000">; to write the meta data associated with each variable</font></div><div><font color="#ff0000">;===============================================================</font></div><div><font color="#ff0000">    ncdf->spi = spi             ; 3D        </font></div><div><font color="#ff0000">end</font></div></div><div><br></div><div><br></div><div>Regards</div></div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><i>Dr. MUHAMMAD AFZAAL KARORI</i><br>Pakistan Meteorological Department<br>Sector H-8/2, Pitras Bukhari Road,<br>Islamabad - 44000<br>Pakistan <br><br><br></div></div></div></div></div>
</div>