<div dir="ltr">Dear fellow NCL users,<div><br></div><div><br></div><div>I would like to extract the "date" when any grid point in a domain exceeds a threshold value using NCL (daily time steps). I tried to create a netcdf file of daily vorticity using the script below.</div><div><br></div><div>I would like to add a line in this script that will extract the date (day month year) whenever a grid point from this domain is exceeded:</div><div> LAT : 7.5N to 27.5N</div><div> LON : 110E to 135E</div><div><br></div><div>the threshold value that I am using is 3e-5 per second.</div><div><br></div><div><br></div><div>I am attaching the netcdf file and the image in this email.<div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;color:rgb(34,34,34);font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/file/d/0B9faET7Bc2o8RWRXcTdKZEx0YlE/view?usp=drive_web" target="_blank" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0;border:none"><img style="vertical-align: bottom; border: none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:rgb(17,85,204);text-decoration:none;vertical-align:bottom">test_2003_30.png</span></a><img src="//ssl.gstatic.com/ui/v1/icons/common/x_8px.png" style="opacity: 0.55; cursor: pointer; float: right; position: relative; top: -1px; display: none;"></div></div><div><br></div><div>Is this possible in NCL?</div><div>Is there any function in NCL that can extract values from every gridpoint in a specified domain?</div><div><br></div><div><br></div><div>I'll appreciate if anyone can point me to the right direction/functions to implement this in NCL.</div><div><br></div><div>Many thanks in advance!</div><div><br></div><div><br></div><div><div>begin</div><div>;************************************************</div><div>; variable and file handling</div><div>;************************************************</div><div> ufile = addfile("uwind_may_2003_700hPa.nc","r") ; open netcdf file</div><div> u = ufile->uwnd(:,0,:,:) ; pull u off file</div><div> lat = ufile->lat</div><div> lon = ufile->lon</div><div> vfile = addfile("vwind_may_2003_700hPa.nc","r")</div><div> v = vfile->vwnd(:,0,:,:) ; pull v off file</div><div>;************************************************</div><div>; calculate vorticity on a Gaussian Grid</div><div>;************************************************ </div><div> ;vrt = u ; retain coordinates</div><div> ;vrt = uv2vrG_Wrap(u,v)</div><div> vr = uv2vr_cfd (u,v,lat,lon, 3)</div><div> copy_VarMeta(u,vr)</div><div> vr@long_name = "vorticity"</div><div> vr@units = "per second"</div><div>;***********************************************</div><div>; save to a netcdf file</div><div>;***********************************************</div><div> ncdf = addfile("vort_700hPa.nc","c")</div><div> fAtt = True</div><div> fAtt@title = "Vorticity at 700hPa"</div><div> fAtt@source_file = "daily winds from NCEP-NCAR Reanaysis"</div><div> fAtt@Conventions = "None"</div><div> fAtt@creation_date = systemfunc("date")</div><div> fileattdef(ncdf,fAtt) ; copy file attributes</div><div> filedimdef(ncdf,"time",-1,True)</div><div> ncdf->vr = vr</div><div>end</div></div><div><br></div><div><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><font face="comic sans ms, sans-serif"><b><u>Lyndz</u></b></font></div></div></div></div></div></div></div></div>
</div></div>