<html><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:14pt"><div><span>How to draw Time series in grads using .nc file<br> (dipti kalyanshetti)</span></div><div style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>I think easiest way is calculate the area average of the .nc file for the time to wand to draw the time series.</span></div><div style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>second you have to write these area averaged data points as a text file. then you can use the R , excel or any software to plot your time series</span></div><div style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif; font-size: 18.66px; font-style: normal;
background-color: transparent;"><span>rgd</span></div><div style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>kaikombo</span></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: verdana, helvetica, sans-serif; font-size: 14pt;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> On Wednesday, July 16, 2014 8:37 AM, "gradsusr-request@gradsusr.org" <gradsusr-request@gradsusr.org> wrote:<br> </font> </div> <br><br> <div>Send gradsusr mailing list submissions to<br> <a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><br>To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>or, via email, send a message with subject or body 'help' to<br> <a href="mailto:gradsusr-request@gradsusr.org" ymailto="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br><br>You can reach the person managing the list at<br> <a href="mailto:gradsusr-owner@gradsusr.org" ymailto="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br><br>When replying, please edit your Subject line so it is more specific<br>than "Re: Contents of gradsusr digest..."<br><br><br>Today's Topics:<br><br> 1. how to draw plot tangential and radial wind speed in grads<br> (mannem venkat)<br> 2. How to draw Time series in grads using .nc file<br> (dipti kalyanshetti)<br> 3. Re: How to draw Time series in
grads using .nc file<br> (mannem venkat)<br> 4. Fit Function in the GrADS extension lilbrary (Michael Greve)<br> 5. Re: Help on calculating Area Maximum from NetCDF file using<br> Grads (Anil Kumar)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>Date: Wed, 16 Jul 2014 11:38:26 +0530<br>From: mannem venkat <<a href="mailto:reddyyvu@gmail.com" ymailto="mailto:reddyyvu@gmail.com">reddyyvu@gmail.com</a>><br>Subject: [gradsusr] how to draw plot tangential and radial wind speed<br> in grads<br>To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Message-ID:<br> <<a href="mailto:CAOmoL22QhaRZgCG4ukCH3hR6wVEeP6QzmjUTkuaNjPRGxCap4g@mail.gmail.com"
ymailto="mailto:CAOmoL22QhaRZgCG4ukCH3hR6wVEeP6QzmjUTkuaNjPRGxCap4g@mail.gmail.com">CAOmoL22QhaRZgCG4ukCH3hR6wVEeP6QzmjUTkuaNjPRGxCap4g@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>Hi<br>I am working on modelling activity.<br>I want to plot the vertical cross section of tangential and radial wind<br>speed from WRF model output data.<br>Anybody can send me the related program to plot these two winds in grads.<br><br>thanking you<br><br><br>-- <br>Yours<br><br>M. Venkatramireddy<br>Senior Research Fellow,<br>SARC-ISRO<br>Dept. of Physcis,<br>Yogi Vemana University,<br>Kadapa-516003<br>Andhra pradesh ,India<br>cell:9492515473<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/b08e4264/attachment-0001.html"
target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/b08e4264/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 2<br>Date: Wed, 16 Jul 2014 15:26:27 +0530<br>From: dipti kalyanshetti <<a href="mailto:diptikalyanshetti@gmail.com" ymailto="mailto:diptikalyanshetti@gmail.com">diptikalyanshetti@gmail.com</a>><br>Subject: [gradsusr] How to draw Time series in grads using .nc file<br>To: <a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>Message-ID:<br> <CAGRW+JoaBHhVLr3UD5kcK_j2URpSqkkuwk_ESg0tHFm+s=<a href="mailto:S2ew@mail.gmail.com" ymailto="mailto:S2ew@mail.gmail.com">S2ew@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>Hi,<br><br>I am new to Grads. I want to plot the time series of rainfall using<br>imd_NEM_1901_2000Series.nc file.<br><br>This file contains 1901 to 2000 rainfall data.<br><br>I
used sdfopen file to open above file. Then I used 'display rf' or 'd rf'<br>, but getting only one value as result.<br><br>I want to get graph as rainfall on one axis and years on other.<br><br>Please someone help me in this.<br><br>Thanks & Regards,<br>Dipti<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/5a531396/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/5a531396/attachment-0001.html</a><br>-------------- next part --------------<br>A non-text attachment was scrubbed...<br>Name: imd_NEM_1901_2000Series.nc<br>Type: application/x-netcdf<br>Size: 4188 bytes<br>Desc: not available<br>Url : <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/5a531396/attachment-0001.nc"
target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/5a531396/attachment-0001.nc</a><br><br>------------------------------<br><br>Message: 3<br>Date: Wed, 16 Jul 2014 16:54:41 +0530<br>From: mannem venkat <<a href="mailto:reddyyvu@gmail.com" ymailto="mailto:reddyyvu@gmail.com">reddyyvu@gmail.com</a>><br>Subject: Re: [gradsusr] How to draw Time series in grads using .nc<br> file<br>To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>>,<br> <a href="mailto:diptikalyanshetti@gmail.com" ymailto="mailto:diptikalyanshetti@gmail.com">diptikalyanshetti@gmail.com</a><br>Message-ID:<br> <CAOmoL22CRXt-G5ZZapZGWTyjWYmuQ+<a href="mailto:corP_VYx9kOvXF9xn-gQ@mail.gmail.com" ymailto="mailto:corP_VYx9kOvXF9xn-gQ@mail.gmail.com">corP_VYx9kOvXF9xn-gQ@mail.gmail.com</a>><br>Content-Type: text/plain;
charset="utf-8"<br><br>hello dipti<br>i see your file<br>you using<br>sdfopen *.nc<br>set t 1 100 (mean 1-1901 2-1902 3-1903--------------------)<br>d rf<br><br><br><br>On Wed, Jul 16, 2014 at 3:26 PM, dipti kalyanshetti <<br><a href="mailto:diptikalyanshetti@gmail.com" ymailto="mailto:diptikalyanshetti@gmail.com">diptikalyanshetti@gmail.com</a>> wrote:<br><br>> Hi,<br>><br>> I am new to Grads. I want to plot the time series of rainfall using<br>> imd_NEM_1901_2000Series.nc file.<br>><br>> This file contains 1901 to 2000 rainfall data.<br>><br>> I used sdfopen file to open above file. Then I used 'display rf' or 'd rf'<br>> , but getting only one value as result.<br>><br>> I want to get graph as rainfall on one axis and years on other.<br>><br>> Please someone help me in this.<br>><br>> Thanks & Regards,<br>> Dipti<br>><br>> _______________________________________________<br>> gradsusr
mailing list<br>> <a href="mailto:gradsusr@gradsusr.org" ymailto="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>><br><br><br>-- <br>Yours<br><br>M. Venkatramireddy<br>Senior Research Fellow,<br>SARC-ISRO<br>Dept. of Physcis,<br>Yogi Vemana University,<br>Kadapa-516003<br>Andhra pradesh ,India<br>cell:9492515473<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/07d945f8/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/07d945f8/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 4<br>Date: Wed, 16 Jul 2014 10:35:18 -0400<br>From: Michael Greve <<a href="mailto:mgreve@planalytics.com"
ymailto="mailto:mgreve@planalytics.com">mgreve@planalytics.com</a>><br>Subject: [gradsusr] Fit Function in the GrADS extension lilbrary<br>To: <a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>Message-ID:<br> <<a href="mailto:CAOimPqw6Qjpv_iBVdmwZWDgpBZnmTGg50uwTiEab55iUnTDYVw@mail.gmail.com" ymailto="mailto:CAOimPqw6Qjpv_iBVdmwZWDgpBZnmTGg50uwTiEab55iUnTDYVw@mail.gmail.com">CAOimPqw6Qjpv_iBVdmwZWDgpBZnmTGg50uwTiEab55iUnTDYVw@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>All,<br> Does anyone know what the bracketed numbers are in the raw score<br>regression equation in the output from the "fit" function from Ben-Jei<br>Tsuang's collection? Also does anyone know if the number in the<br>parentheses is actually the residual, i.e (X1 + 0.05...)? If anyone could<br>help, that would be great.<br><br>-- <br>Michael E. Greve
Jr.<br>Planalytics, Inc.<br><br>***************************************************<br>The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. <br><br>If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.<br>***************************************************<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/a6378d92/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/a6378d92/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 5<br>Date: Wed, 16 Jul 2014 21:15:55
+0530<br>From: "Anil Kumar" <<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a>><br>Subject: Re: [gradsusr] Help on calculating Area Maximum from NetCDF<br> file using Grads<br>To: "'GrADS Users Forum'" <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Cc: <a href="mailto:rrajanv@oceanmanager.com" ymailto="mailto:rrajanv@oceanmanager.com">rrajanv@oceanmanager.com</a>, 'Mahesh Gulati'<br> <<a href="mailto:maheshg@oceanmanager.com" ymailto="mailto:maheshg@oceanmanager.com">maheshg@oceanmanager.com</a>><br>Message-ID: <01bd01cfa10d$0aee2620$20ca7260$@oceanmanager.com><br>Content-Type: text/plain; charset="utf-8"<br><br>Dear Arlindo and All Grads Users,<br><br> <br><br>As directed, I have developed below Script, defined a global grid u4x4 and calculated Amax
and used ?set defval? command. When I exported the defined grid to netcdf file. The output was vertical straight lines. Can anyone please point out what is wrong while defining ?set defval? in script? <br><br> <br><br>Any help please.<br><br> <br><br>'reinit'<br><br>'open D:\WeatherGribFiles\2014\06\11\rtofs_glo_2ds_f048_daily_prog.ctl'<br><br> <br><br>* original grid dimensions<br><br>xmin = 74.16<br><br>xmax = 434.0<br><br>ymin = -90<br><br>ymax = 90<br><br> <br><br>* defining "u4x4" global grid for u<br><br>'define u4x4 = re(u_velocity,2)'<br><br> <br><br>* defining var for further use<br><br>'define var = u_velocity'<br><br> <br><br>xinc = 2<br><br>yinc = 1<br><br>xnow = xmin<br><br>ynow = ymin<br><br> <br><br>******* jj and jj_max are the latitude of Global grid u4x4 ****************<br><br>jj = 1<br><br>jj_max = 90<br><br>count = 1<br><br> <br><br>while (xnow < xmax)<br><br>xtempMax = xnow+xinc;ytempMax = ynow+yinc<br><br> <br><br>'d
amax(var,lon='xnow',lon='xtempMax',lat='ynow',lat='ytempMax')'<br><br>maxval=subwrd(result,4)<br><br> <br><br> jj = 1<br><br> while (jj <= jj_max)<br><br> 'set defval u4x4 'count' 'jj' 'maxval<br><br> jj = jj + 1<br><br> endwhile<br><br>xnow = xnow+xinc;ynow = ynow+yinc<br><br>count=count+1<br><br>endwhile<br><br> <br><br> <br><br>*************** Write a netcdf file ****************<br><br>'set sdfwrite a_grid.nc'<br><br>'sdfwrite u4x4'<br><br> <br><br> <br><br>Thanks & Regards<br><br>Anil Kumar <br><br>----------------------------------------<br><br>---------- Forwarded message ----------<br>From: Arlindo da Silva <<a href="mailto:dasilva@alum.mit.edu" ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a> <mailto:<a href="mailto:dasilva@alum.mit.edu" ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>> ><br>Date: Thu, Jul 10, 2014 at 7:48 PM<br>Subject: Re: [gradsusr] Help on
calculating Area Maximum from NetCDF file using Grads<br>To: Anil Kumar <<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a> <mailto:<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a>> >, GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>> ><br><br><br><br>On Thu, Jul 10, 2014 at 5:40 AM, Anil Kumar <<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a> <mailto:<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a>> > wrote:<br><br>Dear Arlindo,<br><br> <br><br>Thanks for your detailed reply. I have a few queries
regarding the task I need to accomplish.<br><br> <br><br>1) If Grads is not the best tool to accomplish my task, can you please suggest me the another tool, if you know, to accomplish the job ?<br><br> <br><br>Probably C or Fortran. You could also accomplish some of that with python with pygrads. All depends on how much programming experience you have.<br><br> <br><br>2) Secondly, the approach which you have already shared with me as below<br><br> <br><br>1) define a 2x2 global grid (say, use lterp with linear interpolation),<br><br>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.<br><br> <br><br>Should I need to create 2x2 global grid first physically first and then need to replace its values. If so, I don?t find any hint on
google, <<a href="http://gradsuser.org/" target="_blank">http://gradsuser.org/</a>> gradsuser.org to create a new grid. What is the procedure to create grid ?<br><br> <br><br> <br><br>What I said was to define a 2x2 grid, for example<br><br> <br><br>define u2x2 = re(u,2)<br><br> <br><br> I really don't have the time to walk you thru the details as it would require me to actually write the code for you. If you can program in C or fortran I'd go that route. If you follow the UDF approach that I outlined you can still do grads for your plotting. And please do not write directly to me, I cannot afford to provide individualized grads support.<br><br> <br><br> Good Luck,<br><br> <br><br> Arlindo<br><br> <br><br> <br><br> <br><br> <br><br>Thanks again for your help.<br><br> <br><br> <br><br>Thanks & Regards<br><br>Anil Kumar <br><br>OceanManager Inc. | <<a
href="http://www.oceanmanager.com/" target="_blank">http://www.oceanmanager.com/</a>> www.oceanmanager.com<br><br>OceanManager - Marine Software<br><br> <br><br>From: <mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>> <a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto: <mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>> <a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Arlindo da Silva<br>Sent: Tuesday, July 01, 2014 7:12 AM<br>To: GrADS Users Forum<br><br><br>Subject: Re: [gradsusr] Help on calculating Area Maximum from NetCDF file using Grads<br><br> <br><br>On Mon, Jun 30, 2014 at 9:37
AM, Anil Kumar <<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a> <mailto:<a href="mailto:anilku@oceanmanager.com" ymailto="mailto:anilku@oceanmanager.com">anilku@oceanmanager.com</a>> > wrote:<br><br>Hi Grads users,<br><br> <br><br>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.<br><br> <br><br>My Requirement : 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.<br><br> <br><br> <br><br>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<br><br> <br><br>1) define a 2x2 global grid (say, use lterp with linear interpolation),<br><br>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.<br><br> <br><br>If you are up for some C/Fortran programming, you can use the opengrads IPC extension (which is available with your Windows build):<br><br> <br><br> <a href="http://opengrads.org/doc/udxt/libipc/" target="_blank">http://opengrads.org/doc/udxt/libipc/</a><br><br> <br><br>to write your variables to disk,<br><br> <br><br>ga-> ipc_save u_velocity u.bin<br><br>ga-> ipc_save v_velocity v.bin<br><br> <br><br>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). <br><br> <br><br>ga-> ! my_amaxregrid <br><br> <br><br>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:<br><br> <br><br>ga-> d ipc_load("u2x2.bin");ipc_load("v2x2.bin")<br><br> <br><br>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.<br><br> <br><br> Good Luck,<br><br> <br><br> Arlindo<br><br> <br><br> <br><br> <br><br> <br><br>What Am doing so far : <br><br>I have a NetCDF file having below CTL Info.<br><br> <br><br>*****************************************************************<br><br>dset D:\WeatherGribFiles\2014\06\22\ <<a
href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC<br><br>undef 1.26765e+30<br><br>dtype netcdf<br><br>xdef 1440 linear 74.16 0.24999<br><br>ydef 720 linear -90 0.24999<br><br>zdef 1 linear 1 1<br><br>tdef 1 linear 00Z23JUN2014 1mn<br><br>vars 2<br><br>u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s)<br><br>v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s)<br><br>endvars<br><br>*****************************************************************<br><br> <br><br>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
rtofs_glo_2ds_f024_daily_prog.ctl:<br><br> <br><br>*****************************************************************<br><br>dset D:\WeatherGribFiles\2014\06\22\ <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>index D:\WeatherGribFiles\2014\06\22\rtofs_glo_2ds_f024_daily_prog_DUMMY.nc <br><br>title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC <br><br>undef 1.26765e+30 <br><br>dtype netcdf <br><br>xdef 180 linear 74.16 2.0 <br><br>ydef 90 linear -90 2.0 <br><br>zdef 1 linear 1 1 <br><br>tdef 1 linear 00Z23JUN2014 1mn<br><br>options template <br><br>vars 2 <br><br>u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s) <br><br>v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s)
<br><br>endvars<br><br>*****************************************************************<br><br>Below is the command I use to display the data from destination grid rtofs_glo_2ds_f024_daily_prog_DUMMY.nc defined in above CTL file. Below is the command set using.<br><br> <br><br>*****************************************************************<br><br>reinit<br><br>sdfopen D:/WeatherGribFiles/2014/06/22/ <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>open D:/WeatherGribFiles/2014/06/22/rtofs_glo_2ds_f024_daily_prog2.ctl<br><br>set parea 0 11 0 8.5<br><br>set vpage 0 11 0 8.5<br><br>set lon -74.16 434<br><br>set lat -90 90<br><br>set cthick 6<br><br>set arrscl 0.10 <br><br>define val = mag(lterp(u_velocity.1*1.95,lat.2,amax),lterp(v_velocity.1*1.95,lat.2,amax))<br><br>d
lterp(u_velocity.1*1.95,lat.2,amax);lterp(v_velocity.1*1.95,lat.2,amax);val<br><br> <br><br>*****************************************************************<br><br> <br><br>When I change the destination grid dimensions slightly in rtofs_glo_2ds_f024_daily_prog.ctl file<br><br> <br><br>xdef 135 linear 74.16 2.5 <br><br>ydef 68 linear -90 2.5 <br><br> <br><br>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.<br><br>What I want is, When I increase or decrease the destination grid points in ?rtofs_glo_2ds_f024_daily_prog.ctl?, and the maximum value should be at its place. <br><br> <br><br>How can I achieve this ????? <br><br> <br><br> <br><br>Thanks & Regards<br><br>Anil Kumar ? Sr. Software Developer<br><br>OceanManager Inc. | <<a href="http://www.oceanmanager.com/"
target="_blank">http://www.oceanmanager.com/</a>> www.oceanmanager.com<br><br>Tel: <tel:%2B%2091-172-5026090> + 91-172-5026090, 5026091, 4646070<br><br>OceanManager- Marine Software<br><br> <br><br>From: <mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>> <a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto: <mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>> <a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On Behalf Of Anil Kumar<br>Sent: Friday, June 27, 2014 12:55 PM<br>To: <mailto:<a href="mailto:jma@cola.iges.org" ymailto="mailto:jma@cola.iges.org">jma@cola.iges.org</a>> <a
href="mailto:jma@cola.iges.org" ymailto="mailto:jma@cola.iges.org">jma@cola.iges.org</a>; 'GrADS Users Forum'<br>Subject: Re: [gradsusr] Help on calculating Area Maximum from NetCDF file using Grads<br><br> <br><br>Hi Jennifer and All,<br><br> <br><br>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">http://www.iges.org/grads/gadoc/gradfunclterp.html</a>> but in documentation I cannot find that ltrep uses amax as argument ; neither it is giving correct result nor giving error.<br><br> <br><br>My Requirement : 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.<br><br> <br><br>What Am doing so far : <br><br>I have a NetCDF file having below CTL Info.<br><br> <br><br>*****************************************************************<br><br>dset D:\WeatherGribFiles\2014\06\22\ <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC<br><br>undef 1.26765e+30<br><br>dtype netcdf<br><br>xdef 1440 linear 74.16 0.24999<br><br>ydef 720 linear -90 0.24999<br><br>zdef 1 linear 1 1<br><br>tdef 1 linear 00Z23JUN2014 1mn<br><br>vars 2<br><br>u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s)<br><br>v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity
(m/s)<br><br>endvars<br><br>*****************************************************************<br><br> <br><br>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 rtofs_glo_2ds_f024_daily_prog.ctl:<br><br> <br><br>*****************************************************************<br><br>dset D:\WeatherGribFiles\2014\06\22\ <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>index D:\WeatherGribFiles\2014\06\22\rtofs_glo_2ds_f024_daily_prog_DUMMY.nc <br><br>title 2D Sfc Daily Prognostic 00Z22jun2014: Forecast, downloaded Jun 22 15:20 UTC <br><br>undef 1.26765e+30 <br><br>dtype netcdf <br><br>xdef 180 linear 74.16 2.0 <br><br>ydef 90 linear -90 2.0
<br><br>zdef 1 linear 1 1 <br><br>tdef 1 linear 00Z23JUN2014 1mn<br><br>options template <br><br>vars 2 <br><br>u_velocity=>u_velocity 1 t,z,y,x eastward_sea_water_velocity (m/s) <br><br>v_velocity=>v_velocity 1 t,z,y,x northward_sea_water_velocity (m/s) <br><br>endvars<br><br>*****************************************************************<br><br>Below is the command I use to display the data from destination grid rtofs_glo_2ds_f024_daily_prog_DUMMY.nc defined in above CTL file. Below is the command set using.<br><br> <br><br>*****************************************************************<br><br>reinit<br><br>sdfopen D:/WeatherGribFiles/2014/06/22/ <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> rtofs_glo_2ds_f024_daily_prog.nc<br><br>open D:/WeatherGribFiles/2014/06/22/rtofs_glo_2ds_f024_daily_prog2.ctl<br><br>set parea 0 11
0 8.5<br><br>set vpage 0 11 0 8.5<br><br>set lon -74.16 434<br><br>set lat -90 90<br><br>set cthick 6<br><br>set arrscl 0.10 <br><br>define val = mag(lterp(u_velocity.1*1.95,lat.2,amax),lterp(v_velocity.1*1.95,lat.2,amax))<br><br>d lterp(u_velocity.1*1.95,lat.2,amax);lterp(v_velocity.1*1.95,lat.2,amax);val<br><br> <br><br>*****************************************************************<br><br> <br><br>When I change the destination grid dimensions slightly in rtofs_glo_2ds_f024_daily_prog.ctl file<br><br> <br><br>xdef 135 linear 74.16 2.5 <br><br>ydef 68 linear -90 2.5 <br><br> <br><br>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.<br><br>What I want is, When I increase or decrease the destination grid points in ?rtofs_glo_2ds_f024_daily_prog.ctl?, and the maximum value should be at its place. <br><br>
<br><br>How can I achieve this ????? <br><br> <br><br>Thanks for any guidance.<br><br> <br><br>Thanks & Regards<br><br>Anil Kumar<br><br><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org" ymailto="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><br><br><br><br> <br><br>-- <br>Arlindo da Silva<br> <mailto:<a href="mailto:dasilva@alum.mit.edu" ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>> <a href="mailto:dasilva@alum.mit.edu" ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a> <br><br><br><br><br><br> <br><br>-- <br>Arlindo da Silva<br> <mailto:<a href="mailto:dasilva@alum.mit.edu"
ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>> <a href="mailto:dasilva@alum.mit.edu" ymailto="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a> <br><br><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org" ymailto="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> <br><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140716/05668621/attachment.html"
target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140716/05668621/attachment.html</a><br><br>------------------------------<br><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" ymailto="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><br>End of gradsusr Digest, Vol 53, Issue 36<br>****************************************<br><br><br></div> </div> </div> </div> </div></body></html>