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