<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:14pt"><div><span>Hi grads 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 have got four files and I want to calculate the </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>(1) 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>(2) the divergence of 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>the .nc file are pres.sfc.2006.nc, shum.2006.nc (1000-300mb),uwnd.2006.nc (1000-50mb) and vwnd.2006.nc (1000-50mb)</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><br></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 used the following algorithm to calculate the 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>'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)'</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>when I ran these command in the grads HDF I got the following error </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>error invalid expression 'define qu=vint(pres,shum.2*uwnd.3,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>Error from UNIT too few levels in default 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>operation error =>error from vint
function error occurred at column 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><br></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 changed the file for surface pressure and used the slp.2006.nc but the same case happen again </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>please ay one can help</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>cheers</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>kaikombo</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><br></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><br></span></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <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;"> <div dir="ltr"> <font face="Arial" size="2"> On Monday, June 23, 2014 7:26 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. 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.html" target="_blank">http://www.iges.org/grads/gadoc/gradfunclterp.html</a>but 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><br><br></div> </div> </div> </div> </div></body></html>