<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:14pt"><div><span>Regarding calculating the vertical integrated moisture flux</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>Hi grad users</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>I am running the script to calculate the moisture flux for the levels from surface to 850,500,700, and 300</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>I have used the vint and the -hdvg(qu,qv) but I am getting the error messges</span></div><div style="color:
rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>Here are the file </span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>sdfopen slp.2006.nc</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>sdfopen shum.2006.nc</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>sdfopen uwnd.2006.nc</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>sdfopen
vwnd.2006.nc</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>Now I set the dimensions of the data as follows\</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>set x 1 144</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>set y 1 73</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>set t 1 365</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color:
transparent;"><span>set z 1</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>define a=shum.2*uwnd.3</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>define b=shum.2*vwnd.4</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>difine qu=vint(press(z=1),a,850)</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><span>define qv=vint(pres(z=1),b,850)</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size:
18.66px; font-style: normal; background-color: transparent;"><span>Then I define the vertical integrated moisture as follow</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><br></div><span><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;">VIMF=<font size="3">
<font face="verdana, helvetica, sans-serif">divergence of qu and qv</font></font></div><div style="color: rgb(0, 0, 0); font-family: verdana, helvetica, sans-serif; font-size: 16px; font-style: normal; background-color: transparent;"><font face="Verdana" size="3"> =<font size="5">hdivg(qu,qv)</font></font></div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;"><font size="5">Here I get the following error </font></div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;">two few levels and othe more</div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;"><br></div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color:
transparent;"><br></div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;">then I tried using the slp to be the surface pressure I end with almost the same errors</div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;">I there any one to help to fix this</div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;">regards</div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;">kaikombo</div><div style="color: rgb(0, 0, 0); font-family: Verdana; font-size: 24px; font-style: normal; background-color: transparent;"><br></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 16px; font-style: normal; background-color:
transparent;"><font size="3">
</font><span style='font-family: "Times New Roman","serif"; font-size: 14pt; mso-fareast-font-family: "Times New Roman";'> -</span></div><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><font size="3">
<br></font></div></span><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><br></div><div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"><div style="color: rgb(0, 0, 0); font-family: times new roman, new york, times, serif; font-size: 18.66px; font-style: normal; background-color: transparent;"><font face="Arial" size="2">On Wednesday, June 25, 2014 3:18 AM, "gradsusr-request@gradsusr.org" <gradsusr-request@gradsusr.org> wrote:<br> </font></div><div style="font-family: times new roman, new york, times, serif; font-size: 14pt;"><div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 12pt;"> <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. Re: gradsusr Digest, Vol 52, Issue 77 (kombo
kai)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>Date: Wed, 25 Jun 2014 03:14:49 -0700<br>From: kombo kai <<a href="mailto:kaikombo@yahoo.com" ymailto="mailto:kaikombo@yahoo.com">kaikombo@yahoo.com</a>><br>Subject: Re: [gradsusr] gradsusr Digest, Vol 52, Issue 77<br>To: "<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>" <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Message-ID:<br> <<a href="mailto:1403691289.4643.YahooMailNeo@web121005.mail.ne1.yahoo.com" ymailto="mailto:1403691289.4643.YahooMailNeo@web121005.mail.ne1.yahoo.com">1403691289.4643.YahooMailNeo@web121005.mail.ne1.yahoo.com</a>><br>Content-Type: text/plain; charset="iso-8859-1"<br><br>Hi grads users<br>I have got four files and I want to calculate the <br>(1) vertical integrated
moisture flux <br>(2) the divergence of the? vertical integrated moisture flux<br>the .nc file are? pres.sfc.2006.nc, shum.2006.nc (1000-300mb),uwnd.2006.nc (1000-50mb) and vwnd.2006.nc (1000-50mb)<br><br>I used the following algorithm to calculate the moisture flux<br>'set x 1 144'<br>'set y 1 73'<br>'set t 1 365'<br>'set lev 1'<br>'define qu=vint(pres,shum.2*uwnd.3,300)'<br>'define qv=vint((pres,shum.2*vwnd.4,300)'<br>'define vi=-hdivg(qu,qv)'<br>when I ran these command in the grads HDF I got the following error <br>error invalid expression 'define qu=vint(pres,shum.2*uwnd.3,300)<br>Error from UNIT? too few levels in default file <br>operation error? =>error? from vint function error occurred at column 1<br><br>Then I changed the file for surface pressure and used the slp.2006.nc but the same case happen again <br>please ay one can help<br>cheers<br>kaikombo<br><br> <br><br><br>On Monday, June 23, 2014 7:26 AM, "<a
href="mailto:gradsusr-request@gradsusr.org" ymailto="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>" <<a href="mailto:gradsusr-request@gradsusr.org" ymailto="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a>> wrote:<br> <br><br><br>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. Re: Use of ^ Character in Descriptor File (Frank Colby)<br>? 2. symbol lookup error (Claudia Wanzeler)<br>? 3. Re: symbol lookup error (Jennifer Adams)<br>? 4. Re: Help on Areal average using Grads (Anil Kumar)<br><br><br>----------------------------------------------------------------------<br><br>Message: 1<br>Date: Mon, 23 Jun 2014 08:22:57 -0400<br>From: Frank Colby <<a href="mailto:Frank_Colby@uml.edu" ymailto="mailto:Frank_Colby@uml.edu">Frank_Colby@uml.edu</a>><br>Subject: Re: [gradsusr] Use of ^ Character in Descriptor File<br>To: <a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>Message-ID: <<a
href="mailto:53A81C21.7010604@uml.edu" ymailto="mailto:53A81C21.7010604@uml.edu">53A81C21.7010604@uml.edu</a>><br>Content-Type: text/plain; charset="iso-8859-1"<br><br>I just looked at my ctl files (windows 7) and they all use the ^ to <br>default to the current directory.? I have to run grads from that <br>directory to make this work.? I use a command window, cd to the <br>directory and then run grads from the command line.<br><br>Frank Colby<br>UMass Lowell<br><br>On 6/22/2014 10:20 AM, Arlindo da Silva wrote:<br>> Sephen,<br>><br>>? ? In the past "^" expansion worked on windows, but since I am not a <br>> windows user I have to rely on users to report such problems.? In <br>> fact, the model.ctl that I use for tests has "^" in it and all windows <br>> builds passes 100% of the tests.<br>><br>>? ? Can you send me the full name of the directory where the data <br>> resides? I want to make sure there is no glitch with the
handling of <br>> specific directory names when expanding "^". Please send me the 2 <br>> versions of the ctl file: one with the "^" and the other with full <br>> pathname and I will take a look.<br>><br>>? ? Arlindo<br>><br>><br>><br>> On Fri, Jun 20, 2014 at 9:23 AM, Stephen McMillan <br>> <<a href="mailto:smcmillan@planalytics.com" ymailto="mailto:smcmillan@planalytics.com">smcmillan@planalytics.com</a> <mailto:<a href="mailto:smcmillan@planalytics.com" ymailto="mailto:smcmillan@planalytics.com">smcmillan@planalytics.com</a>>> wrote:<br>><br>>? ? This question is for those familiar with MS Windows...<br>><br>>? ? The GrADS documentation states that "If a ^ character is placed in<br>>? ? front of /data_filename/, then /data_filename/ is assumed to be<br>>? ? relative to the path of the descriptor file. If you are using the<br>>? ? ^ character in the DSET
entry, then the descriptor file and the<br>>? ? data file may be moved to a new directory without changing any<br>>? ? entries in the data descriptor file, provided their relative paths<br>>? ? remain the same."<br>><br>>? ? However, I have not been able to get this to work.? My descriptor<br>>? ? (.ctl) and data file (.bin) have exactly the same names, apart<br>>? ? from the extensions, and reside in the same folder.? Here's the<br>>? ? applicable "dset" line in my descriptor file:<br>>? ? dset ^tfit.air-nino+soi+amo.lt12.jan51-02.bin<br>><br>>? ? If I include the full pathname in dset, it works fine.? However, I<br>>? ? want to be able to move the files to a new folder without having<br>>? ? to edit the dset entry.? Is there an equivalent for Windows?<br>><br>>? ? I'm using GrADS version 2.0.1.oga.1 on a Windows 7 pc.<br>><br>>?
? Stephen McMillan<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>><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>> Arlindo da Silva<br>> /<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>><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>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a
href="http://gradsusr.org/pipermail/gradsusr/attachments/20140623/a491e448/attachment-0001.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140623/a491e448/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 2<br>Date: Mon, 23 Jun 2014 10:37:42 -0300<br>From: Claudia Wanzeler <<a href="mailto:claudia.wanzeler@gmail.com" ymailto="mailto:claudia.wanzeler@gmail.com">claudia.wanzeler@gmail.com</a>><br>Subject: [gradsusr] symbol lookup error<br>To: <a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>Message-ID:<br>??? <CAN=Q7vnov468bzbMjpqXLqZcOVcUSP3e7UCeQ=<a href="mailto:_ZdKQ_L_OCPw@mail.gmail.com" ymailto="mailto:_ZdKQ_L_OCPw@mail.gmail.com">_ZdKQ_L_OCPw@mail.gmail.com</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>Dear Grads Users,<br><br>I installed the Grads version 2.0.1 on my computer grads. However, when I<br>ran my
script, I am using the NCEP/NCAR Reanalysis data, appeared on the<br>screen the following error message (in red):<br><br>Config: v2.0.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5<br>opendap-grids,stn geotiff shapefile<br>Issue 'q config' command for more detailed configuration information<br>Landscape mode? ('n' for portrait):<br>GX Package Initialization: Size = 11 8.5<br>ga-> run teste_script.gs<br>No hardcopy metafile open<br>All files closed; all defined objects released;<br>All GrADS attributes have been reinitialized<br>1950 365 1<br>grads: symbol lookup error: grads: undefined symbol: nc_get_chunk_cache<br><br><br>Could anyone help me on this?<br><br>Best wishes<br><br><br>-- <br>Claudia? da Costa<br>INPA/UEA<br>Manaus, Brazil<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140623/eeb5bdff/attachment-0001.html"
target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140623/eeb5bdff/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 3<br>Date: Mon, 23 Jun 2014 09:47:11 -0400<br>From: Jennifer Adams <<a href="mailto:jma@cola.iges.org" ymailto="mailto:jma@cola.iges.org">jma@cola.iges.org</a>><br>Subject: Re: [gradsusr] symbol lookup error<br>To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Message-ID: <<a href="mailto:D37713B1-93D6-4F8A-A146-D286ECD6E99F@cola.iges.org" ymailto="mailto:D37713B1-93D6-4F8A-A146-D286ECD6E99F@cola.iges.org">D37713B1-93D6-4F8A-A146-D286ECD6E99F@cola.iges.org</a>><br>Content-Type: text/plain; charset="us-ascii"<br><br>What is the output from 'q config'? What operating system are you running on? Can you also provide the output from 'ncdump -hs' on your data file? <br>--Jennifer<br><br>On Jun 23,
2014, at 9:37 AM, Claudia Wanzeler wrote:<br><br>> Dear Grads Users,<br>> <br>> I installed the Grads version 2.0.1 on my computer grads. However, when I ran my script, I am using the NCEP/NCAR Reanalysis data, appeared on the screen the following error message (in red):<br>> <br>> Config: v2.0.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile<br>> Issue 'q config' command for more detailed configuration information<br>> Landscape mode? ('n' for portrait):? <br>> GX Package Initialization: Size = 11 8.5 <br>> ga-> run teste_script.gs<br>> No hardcopy metafile open<br>> All files closed; all defined objects released;<br>> All GrADS attributes have been reinitialized<br>> 1950 365 1<br>> grads: symbol lookup error: grads: undefined symbol: nc_get_chunk_cache<br>> <br>> <br>> Could anyone help me on this?<br>> <br>> Best wishes<br>> <br>>
<br>> -- <br>> Claudia? da Costa<br>> INPA/UEA<br>> Manaus, Brazil<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>Jennifer M. Adams<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>111 Research Hall, Mail Stop 2B3<br>George Mason University<br>4400 University Drive<br>Fairfax, VA 22030 <br><br><br><br><br><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140623/dc29a2c9/attachment-0001.html"
target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140623/dc29a2c9/attachment-0001.html</a><br><br>------------------------------<br><br>Message: 4<br>Date: Mon, 23 Jun 2014 20:02:43 +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 Areal average using Grads<br>To: "'GrADS Users Forum'" <<a href="mailto:gradsusr@gradsusr.org" ymailto="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>Message-ID: <013f01cf8ef0$014c8140$03e583c0$@oceanmanager.com><br>Content-Type: text/plain; charset="utf-8"<br><br>Dear Jennifer,<br><br><br><br>Referencing your below reply, Am trying to use amax() with linear interpolation to the same problem (need area maximum now) But am not getting the maximum values in destination grid. I have gone throu? the documentation <a
href="http://www.iges.org/grads/gadoc/gradfunclterp.htmlbut" target="_blank">http://www.iges.org/grads/gadoc/gradfunclterp.htmlbut</a>am not getting any clue that we can use amax() function with the ltrep. But when I use it not getting the maximum value?s plot in grads. Below is the command used and CTL file.<br><br><br><br>reinit<br><br>sdfopen D:/WeatherGribFiles/2014/06/11/rtofs_glo_2ds_f024_daily_prog.nc<br><br>open D:/WeatherGribFiles/2014/06/11/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>d lterp(u_velocity.1*1.95,lat.2,amax);lterp(v_velocity.1*1.95,lat.2,amax)<br><br><br><br>CTL file : <br><br><br><br>dset D:\WeatherGribFiles\2014\06\22\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 67 linear -90 3.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>As I change the X def, ydef entries for the dummy destination grid size, value changes but I expect it to be max every time. Please lead me.<br><br><br><br>Thanks & Regards<br><br>Anil Kumar <br><br>OceanManager Inc. <br><br><br><br>From: <a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" ymailto="mailto:gradsusr-bounces@gradsusr.org">gradsusr-bounces@gradsusr.org</a>] On
Behalf Of Jennifer Adams<br>Sent: Thursday, May 01, 2014 9:09 PM<br>To: GrADS Users Forum<br>Subject: Re: [gradsusr] Help on Areal average using Grads<br><br><br><br>Hi, Anil -- <br><br>You need to create a new descriptor file that describes a grid you would like to be your 'destination grid'; presumably it would have grid points that correspond to every 8th grid point in your original grid. It can be a dummy descriptor -- use 'dset ^foo.dat' and 'options template' so GrADS won't care if foo.dat doesn't exist. Open both descriptors, then use lterp to regrid using the aave method:<br><br><br><br>? 'd lterp(u.1,lat.2,aave); lterp(v.1,lat.2,aave)'<br><br><br><br>--Jennifer<br><br><br><br><br><br><br><br><br><br>On May 1, 2014, at 9:58 AM, Anil Kumar wrote:<br><br><br><br><br><br>Hi Grad Users,<br><br><br><br>I have NetCDF file having Ocean Current Data (attached for reference). Am using below command to display u_velocity and v_velocity and
generating vector plot using SKIP() function from Grads.<br><br><br><br>reinit<br><br>sdfopen D:/rtofs_glo_2ds_f024_daily_prog.nc <<a href="http://rtofs_glo_2ds_f024_daily_prog.nc/" target="_blank">http://rtofs_glo_2ds_f024_daily_prog.nc/</a>> <br><br>set vpage 0 11 0 8.5<br><br>set parea 0 11 0 8.5<br><br>set lon 74 434<br><br>set lat -70 70<br><br>set mproj latlon<br><br>set display color white<br><br>set datawarn off<br><br>set grads off<br><br>set timelab off<br><br>set xlab off<br><br>set ylab off<br><br>set frame off<br><br>set t 1<br><br>set cint 40<br><br>set display color white<br><br>set background 0 <br><br>set gxout vector<br><br>set cthick 6<br><br>set ccolor 1<br><br>d skip(u_velocity*1.95,8,8);skip(v_velocity*1.95,8,8)<br><br>printim D:\output.png x1061 y587 -t 0<br><br><br><br>output is attached and unfortunately this is not what i need exactly. It is just plotting every 8th point. I want to plot areal average of every 8 points and
want exactly the same plot.<br><br><br><br><image001.png><br>?<br><br><br><br>NetCDF file has 4320 LON points? and 2160? LAT points. The ?display? command is powerful enough to iterate along the longitude and latitude. But aave() <<a href="http://www.iges.org/grads/gadoc/gradfuncaave.html" target="_blank">http://www.iges.org/grads/gadoc/gradfuncaave.html</a>>? function needs fixed X and Y dimensions, no STRIDE option for X and Y.<br><br><br><br>How can I achieve this using Grads??? Ultimately I need global vector output.<br><br><br><br>Please help me on this. Thanks in advance.<br><br><br><br>-Anil Kumar<br><br><br><br><br><br><rtofs_glo_2ds_f024_daily_prog.nc><output.png>_______________________________________________<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>Jennifer M. Adams<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>111 Research Hall, Mail Stop 2B3<br>George Mason University<br>4400 University Drive<br>Fairfax, VA 22030 <br><br><br><br><br><br><br><br><br><br><br><br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140623/6135e9c6/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140623/6135e9c6/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 52, Issue 77<br>****************************************<br>-------------- next part --------------<br>An HTML attachment was scrubbed...<br>URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20140625/343d380f/attachment.html" target="_blank">http://gradsusr.org/pipermail/gradsusr/attachments/20140625/343d380f/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 52, Issue 83<br>****************************************<br><br><br></div>
</div> </div> </div> </div></body></html>