<div dir="ltr">On Mon, Jun 30, 2014 at 9:37 AM, Anil Kumar <span dir="ltr"><<a href="mailto:anilku@oceanmanager.com" target="_blank">anilku@oceanmanager.com</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><div><p class="MsoNormal">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Hi Grads users,<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">I am awaiting an useful comments or direction on the below case. Please lead me, to get all the maximum values from the small – small areas of every 2 degree in Lat and lon say for example. Below is my complete query.<u></u><u></u></span></p>
<div><div class="h5"><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">My Requirement</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> : I have a larger dataset (.NC file) 1440*720 points. I want a new destination grid with linear interpolation having all the all area maximum values using amax() of small-small regions. And ultimately, need to query the destination grid for output.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> </span></p></div></div></div></div></blockquote><div><br></div><div>I am afraid the regridding method you are looking for is not implemented by lterp() or re(). This is a good example of when GrADS is not the best tool for the job at hand. But if you insist, a very inefficient approach would be</div>
<div><br></div><div>1) define a 2x2 global grid (say, use lterp with linear interpolation),</div><div>2) loop over every 2x2 region in the 0.25x0.5 deg grid, compute the 2 x2 area max with amax(), and use "set defval" to replace the corresponding value in the previous defined variable in 1). As I said, this is very inefficient.</div>
<div><br></div><div>If you are up for some C/Fortran programming, you can use the opengrads IPC extension (which is available with your Windows build):</div><div><br></div><div> <a href="http://opengrads.org/doc/udxt/libipc/">http://opengrads.org/doc/udxt/libipc/</a><br>
</div><div><br></div><div>to write your variables to disk,</div><div><br></div><div>ga-> ipc_save u_velocity u.bin</div><div>ga-> ipc_save v_velocity v.bin</div><div><br></div><div>Then write a C/Fortran program that read these files (u.bin,v.bin) and compute the amax at the 2x2 grid you require, writing output files (u2x2.bin.v2x2.bin). </div>
<div><br></div><div>ga-> ! my_amaxregrid </div><div><br></div><div>Consult the IPC documentation for the description of the transfer file format (very similar to the old GrADS 1.x UDF.) You can finally read these back with ipc_load() extension:</div>
<div><br></div><div>ga-> d ipc_load("u2x2.bin");ipc_load("v2x2.bin")</div><div><br></div><div>Notice that this is a poor-man version of the file-based UDF that existed in GrADS 1.x. Of course, you can accomplish the same more efficiently by writing a genuine OpenGrADS extension with no transfer files involved.</div>
<div><br></div><div> Good Luck,</div><div><br></div><div> Arlindo</div><div><br></div><div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple"><div><div><div class="h5"><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u></span></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">What Am doing so far : <u></u><u></u></span></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">I have a NetCDF file having below CTL Info.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dset D:\WeatherGribFiles\2014\06\22\<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">undef 1.26765e+30<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dtype netcdf<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">xdef 1440 linear 74.16 0.24999<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">ydef 720 linear -90 0.24999<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">zdef 1 linear 1 1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">tdef 1 linear 00Z23JUN2014 1mn<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">vars 2<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">endvars<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">I want a new destination grid ; which would contain the area maximum values as points. Destination grid should contain maximum value from every 8 points in x-axis and y axis. For this, I have created below, CTL file named as <b>rtofs_glo_2ds_f024_daily_prog.ctl</b>:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dset D:\WeatherGribFiles\2014\06\22\<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">index D:\WeatherGribFiles\2014\06\22\<span style="color:red">rtofs_glo_2ds_f024_daily_prog_DUMMY.nc </span><u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:10pt;font-family:Calibri,sans-serif">title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">undef 1.26765e+30 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dtype netcdf <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">xdef 180 linear 74.16 2.0 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">ydef 90 linear -90 2.0</span><span style="font-size:10pt;font-family:Calibri,sans-serif"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">zdef 1 linear 1 1 <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">tdef 1 linear 00Z23JUN2014 1mn<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">options template <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">vars 2 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s) <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s) <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">endvars</span><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Below is the command I use to display the data from destination grid <b>rtofs_glo_2ds_f024_daily_prog_DUMMY.nc </b>defined in above CTL file. Below is the command set using.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">reinit<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">sdfopen D:/WeatherGribFiles/2014/06/22/<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">open D:/WeatherGribFiles/2014/06/22/rtofs_glo_2ds_f024_daily_prog2.ctl<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set parea 0 11 0 8.5<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set vpage 0 11 0 8.5<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set lon -74.16 434<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set lat -90 90<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set cthick 6<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set arrscl 0.10 <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">define val = mag(lterp(u_velocity.1*1.95,lat.2,amax),lterp(v_velocity.1*1.95,lat.2,amax))<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">d lterp(u_velocity.1*1.95,lat.2,amax);lterp(v_velocity.1*1.95,lat.2,amax);val<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">When I change the destination grid dimensions slightly in <b>rtofs_glo_2ds_f024_daily_prog.ctl </b>file<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">xdef 135 linear 74.16 2.5 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">ydef 68 linear -90 2.5</span><span style="font-size:10pt;font-family:Calibri,sans-serif"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">and checks for the output, It varies at the same point i.e. if I look at certain point before changing and after changing grid dimensions, the result comes of different magnitude.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">What I want is, When I increase or decrease the destination grid points in “<b>rtofs_glo_2ds_f024_daily_prog.ctl”, </b>and the maximum value should be at its place. <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">How can I achieve this ????? <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
</div></div><div><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)">Thanks & Regards<u></u><u></u></span></b></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)">Anil Kumar – </span></b><b><span style="font-size:11pt;color:rgb(33,88,104)">Sr. </span></b><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)">Software Developer<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(0,161,218)">OceanManager</span></b><b><i><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)"> </span></i></b><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(0,161,218)">Inc.</span></b><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(0,155,210)"> </span></b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:black">| </span><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><a href="http://www.oceanmanager.com/" target="_blank">www.oceanmanager.com</a></span><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(49,132,155)"><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;color:rgb(237,125,49)">Tel: </span><span style="font-size:11pt;color:rgb(237,125,49)"><a href="tel:%2B%2091-172-5026090" value="+911725026090" target="_blank">+ 91-172-5026090</a>, 5026091, 4646070<u></u><u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(0,161,218)">OceanManager- Marine Software<u></u><u></u></span></b></p></div><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<div><div style="border-style:solid none none;border-top-color:rgb(225,225,225);border-top-width:1pt;padding:3pt 0in 0in"><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif"> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>] <b>On Behalf Of </b>Anil Kumar<br>
<b>Sent:</b> Friday, June 27, 2014 12:55 PM<br><b>To:</b> <a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>; 'GrADS Users Forum'<br><b>Subject:</b> Re: [gradsusr] Help on calculating Area Maximum from NetCDF file using Grads<u></u><u></u></span></p>
</div></div><div><div class="h5"><p class="MsoNormal"><u></u> <u></u></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Hi Jennifer and All,<u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Am seeking help on calculating area maximum with linear interpolation using Grads. Earlier I was using below methodology to calculate area average over small areas defined. It was working good but now I need to calculate amax using lterp in grads <a href="http://www.iges.org/grads/gadoc/gradfunclterp.html" target="_blank">but in documentation I cannot find that ltrep uses amax as argument</a> ; neither it is giving correct result nor giving error.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">My Requirement</span></b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"> : I have a larger dataset (.NC file) 1440*720 points. I want a new destination grid with linear interpolation having all the all area maximum values using amax() of small-small regions. And ultimately, need to query the destination grid for output.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">What Am doing so far : <u></u><u></u></span></b></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">I have a NetCDF file having below CTL Info.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dset D:\WeatherGribFiles\2014\06\22\<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">undef 1.26765e+30<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dtype netcdf<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">xdef 1440 linear 74.16 0.24999<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">ydef 720 linear -90 0.24999<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">zdef 1 linear 1 1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">tdef 1 linear 00Z23JUN2014 1mn<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">vars 2<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">endvars<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">I want a new destination grid ; which would contain the area maximum values as points. Destination grid should contain maximum value from every 8 points in x-axis and y axis. For this, I have created below, CTL file named as <b>rtofs_glo_2ds_f024_daily_prog.ctl</b>:<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dset D:\WeatherGribFiles\2014\06\22\<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">index D:\WeatherGribFiles\2014\06\22\<span style="color:red">rtofs_glo_2ds_f024_daily_prog_DUMMY.nc </span><u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:10pt;font-family:Calibri,sans-serif">title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">undef 1.26765e+30 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">dtype netcdf <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">xdef 180 linear 74.16 2.0 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">ydef 90 linear -90 2.0</span><span style="font-size:10pt;font-family:Calibri,sans-serif"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">zdef 1 linear 1 1 <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">tdef 1 linear 00Z23JUN2014 1mn<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">options template <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">vars 2 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s) <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s) <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">endvars</span><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u><u></u></span></p><p class="MsoNormal">
<span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Below is the command I use to display the data from destination grid <b>rtofs_glo_2ds_f024_daily_prog_DUMMY.nc </b>defined in above CTL file. Below is the command set using.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">reinit<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">sdfopen D:/WeatherGribFiles/2014/06/22/<a href="http://rtofs_glo_2ds_f024_daily_prog.nc" target="_blank">rtofs_glo_2ds_f024_daily_prog.nc</a><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">open D:/WeatherGribFiles/2014/06/22/rtofs_glo_2ds_f024_daily_prog2.ctl<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set parea 0 11 0 8.5<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set vpage 0 11 0 8.5<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set lon -74.16 434<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set lat -90 90<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set cthick 6<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">set arrscl 0.10 <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">define val = mag(lterp(u_velocity.1*1.95,lat.2,amax),lterp(v_velocity.1*1.95,lat.2,amax))<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif">d lterp(u_velocity.1*1.95,lat.2,amax);lterp(v_velocity.1*1.95,lat.2,amax);val<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">*****************************************************************<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">When I change the destination grid dimensions slightly in <b>rtofs_glo_2ds_f024_daily_prog.ctl </b>file<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">xdef 135 linear 74.16 2.5 <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;font-family:Calibri,sans-serif;background-color:yellow">ydef 68 linear -90 2.5</span><span style="font-size:10pt;font-family:Calibri,sans-serif"> <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">and checks for the output, It varies at the same point i.e. if I look at certain point before changing and after changing grid dimensions, the result comes of different magnitude.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">What I want is, When I increase or decrease the destination grid points in “<b>rtofs_glo_2ds_f024_daily_prog.ctl”, </b>and the maximum value should be at its place. <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">How can I achieve this ????? <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)">Thanks for any guidance.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(31,73,125)"><u></u> <u></u></span></p><div><p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)">Thanks & Regards<u></u><u></u></span></b></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(33,88,104)">Anil Kumar</span></b><b><span style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(0,161,218)"><u></u><u></u></span></b></p>
</div></div></div></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><font face="'comic sans ms', sans-serif" size="6" color="#006600"><span style="font-family:Helvetica;font-size:12px"></span><span></span><span></span>Arlindo da Silva</font><br>
<font size="4"><i><a href="mailto:dasilva@alum.mit.edu" target="_blank">dasilva@alum.mit.edu</a></i></font>
</div></div>