From bernd.becker at METOFFICE.GOV.UK Mon Feb 2 11:09:16 2009 From: bernd.becker at METOFFICE.GOV.UK (Bernd Becker) Date: Mon, 2 Feb 2009 16:09:16 +0000 Subject: file check on IO Message-ID: Hello, if I wanted to catch trying to open a non existing file, what do I need? 'sdfopen 'filename * if sdfopen fails, try something else. I tried rc='sdfopen ' filename and other things but to avail. Thanks! Bernd -- Bernd Becker The Monthly Outlook Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.com - http://www.metoffice.com From axl419 at YAHOO.COM Mon Feb 2 20:10:03 2009 From: axl419 at YAHOO.COM (See Hai Ooi) Date: Mon, 2 Feb 2009 17:10:03 -0800 Subject: Error in creating station data using FORTRAN In-Reply-To: <49871C11.2050009@noaa.gov> Message-ID: Dear Messrs Zilore Mumba, Charles Seman and Arlindo, By trial and error and with a hint from Mr Zilore Mumba, I succeeded to resolve the problem myself as shown below:- ............................................................. D:\tfplot>stnmap -i a03.ctl Name of binary data set: d:\tfplot\a03.dat Number of times in the data set: 1 Number of surface variables: 4 Number of level dependent variables: 0 Starting scan of station data binary file. Binary data file open: d:\tfplot\a03.dat Processing time = 1 Time = 1 has stn count = 8 Max reports per time: 8 reports at t = 1 Max data elements in largest report: 4 Version 2 Station map file created: d:\tfplot\a03.map stnmap: WARNING!! This stnmap file can only be accessed by GrADS Version 1.9.0-rc1 stnmap: WARNING!! However, GrADS Version 1.9.0-rc1 can read both versions stnmap: COMMENT -- use the -1 command line option to create a map for older GrADS versions ............................................................... I am pleased to share my adapted and modified work, as attached, with you and the user group. Mr Arlindo, as I did a little bit of modification and addition to the FORTRAN program, you may find this useful to be included in your cookbook. For your information, I use Digital Visual FORTRAN and GrADS v1.9-rc1 (window) in my work. Best regards. OOI See-hai Malaysia Get your new Email address! Grab the Email name you've always wanted before someone else does! http://mail.promotions.yahoo.com/newdomains/aa/ -------------- next part -------------- A non-text attachment was scrubbed... Name: a03.ctl Type: application/octet-stream Size: 320 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090202/ff3b2616/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: a03.dat Type: application/octet-stream Size: 380 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090202/ff3b2616/attachment-0001.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: a03.map Type: application/octet-stream Size: 32 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090202/ff3b2616/attachment-0002.obj -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: a03.txt Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090202/ff3b2616/attachment.txt -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: stndata.for Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090202/ff3b2616/attachment.pl From am_kareem2001 at YAHOO.COM Tue Feb 3 05:38:42 2009 From: am_kareem2001 at YAHOO.COM (ahmed kareem) Date: Tue, 3 Feb 2009 02:38:42 -0800 Subject: correlation errors Message-ID: Dear Grads user?????????? ? I am using correlations scrip, correlating rainfall time series and SSTs. I am always getting this error. ? ? No hardcopy metafile open All files closed; all defined objects released All GrADS attributes have been reinitialized Open Error:? Can't open description file Data Request Error:? Invalid variable name ? Variable 'data3.1' not found in file 1 ? Error ocurred at column 1 DEFINE error:? Invalid expression. SET DFILE error:? file 2 not open Data Request Error:? File number out of range ? Variable = sst.2 ? Error ocurred at column 1 DEFINE error:? Invalid expression. Syntax Error:? Invalid Operand ? 'cru' not a variable or function name Error from TCORR:? Error getting grids Operation Error:? Error from tcorr function ? Error ocurred at column 1 DISPLAY error:? Invalid expression ? Expression = tcorr(cru,sst,t=1,t=30) Cannot plot color bar: No shading information Syntax Error:? Invalid Operand ? 'cru' not a variable or function name Error from TCORR:? Error getting grids Operation Error: ?Error from tcorr function ? Error ocurred at column 1 DISPLAY error:? Invalid expression ? Expression = tcorr(cru,sst,t=1,t=30) successfully done !! ? I am using Grads version 1.8SL9 under windows. ? Attached are the files 1 ? Rainfall files (stn.ctl and mam9.bin) 2- SSTs file (sst6190.ctl and sst6190.oct) 3- The correlation scrip (comb2.gs) ? ? Thanks in advance ? ? ? ? ??????????? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.oct Type: application/octet-stream Size: 1922400 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: COMB2.GS Type: application/octet-stream Size: 1280 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment-0001.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.ctl Type: application/octet-stream Size: 254 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment-0002.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: stn.ctl Type: application/octet-stream Size: 362 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment-0003.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: MAM9.BIN Type: application/octet-stream Size: 720 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/0a9b9043/attachment-0004.obj From j.burn at GEOS.COM Tue Feb 3 06:10:42 2009 From: j.burn at GEOS.COM (James Burn) Date: Tue, 3 Feb 2009 12:10:42 +0100 Subject: Precipitation 3/6 h Message-ID: Hi Fil When looping (with _count) to plot a set of GFS 3hrly precip charts I use: if (math_fmod(_count,2)) 'define APCPsfc6hr=APCPsfc' 'set t '_count-1 'define APCPsfcnew=APCPsfc6hr-APCPsfc' 'd APCPsfcnew' 'set t '_count else 'd APCPsfc' endif Which works nicely. I've not sorted the 3 hourly precip in a meteogram yet, but I imagine something similar would work. Good luck James On Fri, 23 Jan 2009 15:00:06 +0100, Fil wrote: >Hi everyone, > >I'm trying to plot a precipitation meteogram from GFS 3h data but, for >example, +33h timestep has precipitation from +30h to +33h but +36h timestep >also has precipitation from +30h to +33h and +33h to +36h. So what I wanted >to know is how can I subtract precipitation like 30h = 30h - 27h, 36h = 36h >- 33h, 42h = 42h - 39h, etc, so that +30h/+36h/+42h/etc timestep doesn't >include precipitation from +27h/+33h/+39h/etc timestep. > >Also how can I sum precipitation from +3h to 180h? I'm trying to draw the >total precipitation above the bars graph and I tried something like: > >totalprec=sum(APCPsfc,t=2,t=180) >'draw string x y 'totalprec > >But it doesn't work at all. > >Hope someone can help me. > >Cheers. From am_kareem2001 at YAHOO.COM Tue Feb 3 07:26:26 2009 From: am_kareem2001 at YAHOO.COM (ahmed kareem) Date: Tue, 3 Feb 2009 04:26:26 -0800 Subject: dimension error Message-ID: ? ? Dear Grads user?????????? I am using correlations scrip, correlating rainfall time series and SSTs. I am always getting this error. ? ? No hardcopy metafile open All files closed; all defined objects released All GrADS attributes have been reinitialized Data Request Error: Invalid grid coordinates ? World coordinates convert to non-integer grid coordinates ??? Variable = cru? Dimension = 3 Error from TCORR:? Error getting grids Operation Error:? Error from tcorr function ? Error ocurred at column 1 DISPLAY error:? Invalid expression ? Expression = tcorr(cru,sst,t=1,t=30) Cannot plot color bar: No shading information Data Request Error: Invalid grid coordinates ? World coordinates convert to non-integer grid coordinates ??? Variable = cru? Dimension = 3 Error from TCORR:? Error getting grids Operation Error:? Error from tcorr function ? Error ocurred at column 1 DISPLAY error:? Invalid expression ? Expression = tcorr(cru,sst,t=1,t=30) successfully done !! I am using Grads version 1.8SL9 under windows. ? Attached are the files 1 ? Rainfall files (stn.ctl and mam9.bin) 2- SSTs file (sst6190.ctl and sst6190.oct) 3- The correlation scrip (comb2.gs) ? ? Thanks in advance ? ? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.oct Type: application/octet-stream Size: 1922400 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: MAM9.BIN Type: application/octet-stream Size: 720 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment-0001.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.ctl Type: application/octet-stream Size: 254 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment-0002.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: stn.ctl Type: application/octet-stream Size: 362 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment-0003.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: COMB2.GS Type: application/octet-stream Size: 1274 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/342e8adc/attachment-0004.obj From am_kareem2001 at YAHOO.COM Tue Feb 3 14:18:21 2009 From: am_kareem2001 at YAHOO.COM (ahmed kareem) Date: Tue, 3 Feb 2009 11:18:21 -0800 Subject: dimension error Message-ID: Hello I am using Grads version 1.8SL9 under windows xp.?A correlations scrip, correlating rainfall time series and SSTs grids. I am always getting this error. ? ? Data Request Error: Invalid grid coordinates World coordinates convert to non-integer grid coordinates Variable = cru? Dimension = 3 Error from TCORR:? Error getting grids Operation Error:? Error from tcorr function Error ocurred at column 1 DISPLAY error:? Invalid expression ? Expression = tcorr(cru,sst,t=1,t=30) ? ? Attached files are : 1 ? Rainfall files (stn.ctl and mam9.bin) 2- SSTs file (sst6190.ctl and sst6190.oct) 3- The correlation scrip (comb2.gs) ? ? Thanks in advance ? Ahmed M. Abdel Kareem Phd student -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.oct Type: application/octet-stream Size: 1922400 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: MAM9.BIN Type: application/octet-stream Size: 720 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment-0001.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: COMB2.GS Type: application/octet-stream Size: 1274 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment-0002.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.ctl Type: application/octet-stream Size: 254 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment-0003.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: stn.ctl Type: application/octet-stream Size: 362 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/07710a43/attachment-0004.obj From ousmane at IRI.COLUMBIA.EDU Tue Feb 3 14:19:14 2009 From: ousmane at IRI.COLUMBIA.EDU (Ousmane Ndiaye) Date: Tue, 3 Feb 2009 14:19:14 -0500 Subject: using fish_chi extension for daily velocity potential Message-ID: Hi, I am using the function fish_chi to find the velocity potential from the horizontal wind components. The function works great but unfortunately the velocity potential that is defined is on the grid in longitude from 0 to 360 and is not periodic. This is an issue because I would like to display the data over the African continent only, i.e. from longitude 30W to 60E and latitude from 30S to 30N. I have tried to solve this problem by making the input data to the fish_chi function be periodic, i.e. longitude from 0 to 720. When I do this I get the error message: Error from fish: rc = 1000 Operation Error: Error from fish function Operation Error: Error from fish_chi function Error ocurred at column 1 DEFINE error: Invalid expression. Any suggestions on how I can overcome this issue would be greatly appreciated. Or anyone who has the Re-analysis daily velocity potential available that he want to share with. Ousmane From arlindo.dasilva at GMAIL.COM Tue Feb 3 16:08:59 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Tue, 3 Feb 2009 16:08:59 -0500 Subject: using fish_chi extension for daily velocity potential In-Reply-To: <24103394.252641233688754725.JavaMail.root@iripogo1.ldeo.columbia.edu> Message-ID: On Tue, Feb 3, 2009 at 2:19 PM, Ousmane Ndiaye wrote: > Hi, > I am using the function fish_chi to find the velocity potential from the > horizontal wind components. The function works great but unfortunately the > velocity potential that is defined is on the grid in longitude from 0 to 360 > and is not periodic. This is an issue because I would like to display the > data over the African continent only, i.e. from longitude 30W to 60E and > latitude from 30S to 30N. I have tried to solve this problem by making the > input data to the fish_chi function be periodic, > i.e. longitude from 0 to 720. When I do this I get the error message: > > Error from fish: rc = 1000 > Operation Error: Error from fish function > Operation Error: Error from fish_chi function > Error ocurred at column 1 > DEFINE error: Invalid expression. > > Any suggestions on how I can overcome this issue would be greatly > appreciated. Or anyone who has the Re-analysis daily velocity potential > available that he want to share with. Try this: open model set lon -180 180 set lev 200 chi = fish_chi(ua,va) set lon -30 60 set lat -50 60 d chi You should get the attached plot. The trick is to set the domain to [180W,180E] before calling the function. Arlindo > > Ousmane > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090203/3c49c2bd/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: chi.png Type: image/png Size: 77482 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090203/3c49c2bd/attachment.png From anushridutta01 at GMAIL.COM Wed Feb 4 02:20:15 2009 From: anushridutta01 at GMAIL.COM (anushri dutta) Date: Wed, 4 Feb 2009 12:50:15 +0530 Subject: error in ctl file Message-ID: Dear Grads User, Hi, here I am facing a problem, that i want to write a binary data for the 12 month by data taux97-98.cdf, through fwrite function, but when I pot the data through ctl file, it display for a single month, not for all the months, and the error which it showing is ga-> open taux.ctl Scanning description file: taux.ctl Data file tauxobs.grd is open as file 1 LON set to 64.5 160 LAT set to -15 10 LEV set to 1 1 Time values set: 1997:1:1:0 1997:1:1:0 ga-> set t 1 Time values set: 1997:1:1:0 1997:1:1:0 ga-> d taux Contouring: -0.08 to 0.08 interval 0.02 ga-> set t 2 Time values set: 1997:2:1:0 1997:2:1:0 ga-> d taux Low Level I/O Error: Read error on data file Data file name = tauxobs.grd Error reading 192 bytes at location 313344 Data Request Error: Error for variable 'taux' Error ocurred at column 1 DISPLAY error: Invalid expression Expression = taux same as for the 3 to 12 month, I am not understanding why it is so, Please tell me, how do I achieve this, with regards anushri -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/5766b744/attachment.html From sreeds82 at YAHOO.COM Wed Feb 4 03:55:35 2009 From: sreeds82 at YAHOO.COM (sreenivas desamsetti) Date: Wed, 4 Feb 2009 00:55:35 -0800 Subject: Findout SD Message-ID: Dear Sir/Madam, I have two different files with the same dimension. I wish to compute the Standard Deviation (SD) for each grid point basing on the time series of each grid point. Please let me know if any script available to compute the SD. Thanking you D Srinivas Desamsetti Srinivas Research Scholar Department of Meteorology and Oceanography Andhra University Visakhapatnam - 530 003 91 9948530618 (M) sreeds82 at yahoo.com sreeds82 at rediffmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/57d040cd/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/jpeg Size: 11090 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090204/57d040cd/attachment.jpe From bernd.becker at METOFFICE.GOV.UK Wed Feb 4 04:03:35 2009 From: bernd.becker at METOFFICE.GOV.UK (Bernd Becker) Date: Wed, 4 Feb 2009 09:03:35 +0000 Subject: Findout SD In-Reply-To: <905777.14810.qm@web36402.mail.mud.yahoo.com> Message-ID: 'set gxout stat' play with the statistics output device. Good luck, Bernd. On Wed, 2009-02-04 at 00:55 -0800, sreenivas desamsetti wrote: > > > Dear Sir/Madam, > I have two different > files with the same > dimension. > I wish to compute the > Standard Deviation (SD) > for each grid point > basing on the time > series of each grid > point. > Please let me know if > any script available to > compute the SD. > Thanking you > D Srinivas > > > Desamsetti Srinivas > Research Scholar > Department of > Meteorology and > Oceanography > Andhra University > Visakhapatnam - 530 003 > 91 9948530618 (M) > sreeds82 at yahoo.com > sreeds82 at rediffmail.com > > > > > -- Bernd Becker The Monthly Outlook Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.com - http://www.metoffice.com From dbcl at CEH.AC.UK Wed Feb 4 04:03:17 2009 From: dbcl at CEH.AC.UK (Douglas Clark) Date: Wed, 4 Feb 2009 09:03:17 +0000 Subject: error in ctl file Message-ID: Anushri, It's difficult to know what the problem is unless you give more information about how you created the file! There's potentially quite a lot of information that you need to provide, but the following at least would be useful: * the script or commands that you used to perform the fwrite, and what GrADS output to the screen when you did this * the ctl file you created for the output * the results of the command "ncdump -h taux97-98.cdf" (if you have the ncdump utility) In the meantime, is it possible that you only wrote out 1 time of the netCDF file? How big is the file tauxobs.grd (e.g. ls -l in Unix)? If the grid is of size nx*ny points, one year of monthly data should give a file of size nx*ny*12*4 bytes. If your file is smaller, you haven't written all times, so repeat the fwrite stage after 'set t 1 12' (or whatever is appropriate). Try this first! Doug >>> anushridutta01 at GMAIL.COM 04/02/2009 07:20 >>> Dear Grads User, Hi, here I am facing a problem, that i want to write a binary data for the 12 month by data taux97-98.cdf, through fwrite function, but when I pot the data through ctl file, it display for a single month, not for all the months, and the error which it showing is ga-> open taux.ctl Scanning description file: taux.ctl Data file tauxobs.grd is open as file 1 LON set to 64.5 160 LAT set to -15 10 LEV set to 1 1 Time values set: 1997:1:1:0 1997:1:1:0 ga-> set t 1 Time values set: 1997:1:1:0 1997:1:1:0 ga-> d taux Contouring: -0.08 to 0.08 interval 0.02 ga-> set t 2 Time values set: 1997:2:1:0 1997:2:1:0 ga-> d taux Low Level I/O Error: Read error on data file Data file name = tauxobs.grd Error reading 192 bytes at location 313344 Data Request Error: Error for variable 'taux' Error ocurred at column 1 DISPLAY error: Invalid expression Expression = taux same as for the 3 to 12 month, I am not understanding why it is so, Please tell me, how do I achieve this, with regards anushri -- This message (and any attachments) is for the recipient only. NERC is subject to the Freedom of Information Act 2000 and the contents of this email and any reply you make may be disclosed by NERC unless it is exempt from release under the Act. Any material supplied to NERC may be stored in an electronic records management system. From anushridutta01 at GMAIL.COM Wed Feb 4 04:47:10 2009 From: anushridutta01 at GMAIL.COM (anushri dutta) Date: Wed, 4 Feb 2009 15:17:10 +0530 Subject: error in ctl file In-Reply-To: Message-ID: Respected Sir, the command which I am using is the 'reinit' 'sdfopen taux.cdf' 'set gxout fwrite' 'set fwrite tauxobs.grd' 'set x 129.5 320.5' 'set y 121.5 171.5' 'set z 1' 'set t 1 24' 'd taux' 'disable fwrite' and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, and t=24, so, where is the error in file i am not getting it, that is it display the variable for only time t=1, not for the other time, please help me. with redgards anushri On 2/4/09, Douglas Clark wrote: > > Anushri, > > It's difficult to know what the problem is unless you give more information > about how you created the file! > > There's potentially quite a lot of information that you need to provide, > but the following at least would be useful: > * the script or commands that you used to perform the fwrite, and what > GrADS output to the screen when you did this > * the ctl file you created for the output > * the results of the command "ncdump -h taux97-98.cdf" (if you have the > ncdump utility) > > In the meantime, is it possible that you only wrote out 1 time of the > netCDF file? How big is the file tauxobs.grd (e.g. ls -l in Unix)? If the > grid is of size nx*ny points, one year of monthly data should give a file of > size nx*ny*12*4 bytes. If your file is smaller, you haven't written all > times, so repeat the fwrite stage after 'set t 1 12' (or whatever is > appropriate). Try this first! > > Doug > > > >>> anushridutta01 at GMAIL.COM 04/02/2009 07:20 >>> > Dear Grads User, > Hi, > here I am facing a problem, that i want to write a binary data for > the 12 month by data taux97-98.cdf, through fwrite function, > but when I pot the data through ctl file, > it display for a single month, not for all the months, > and the error which it showing is > > ga-> open taux.ctl > Scanning description file: taux.ctl > Data file tauxobs.grd is open as file 1 > LON set to 64.5 160 > LAT set to -15 10 > LEV set to 1 1 > Time values set: 1997:1:1:0 1997:1:1:0 > ga-> set t 1 > Time values set: 1997:1:1:0 1997:1:1:0 > ga-> d taux > Contouring: -0.08 to 0.08 interval 0.02 > ga-> set t 2 > Time values set: 1997:2:1:0 1997:2:1:0 > ga-> d taux > Low Level I/O Error: Read error on data file > Data file name = tauxobs.grd > Error reading 192 bytes at location 313344 > Data Request Error: Error for variable 'taux' > Error ocurred at column 1 > DISPLAY error: Invalid expression > Expression = taux > > > same as for the 3 to 12 month, > > I am not understanding why it is so, > Please tell me, how do I achieve this, > > with regards > > anushri > > > -- > This message (and any attachments) is for the recipient only. NERC > is subject to the Freedom of Information Act 2000 and the contents > of this email and any reply you make may be disclosed by NERC unless > it is exempt from release under the Act. Any material supplied to > NERC may be stored in an electronic records management system. > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/aea88054/attachment.html From madhusoodanan at GMAIL.COM Wed Feb 4 04:54:30 2009 From: madhusoodanan at GMAIL.COM (Madhusoodanan M.S.) Date: Wed, 4 Feb 2009 15:24:30 +0530 Subject: error in ctl file In-Reply-To: <8511ebae0902040147u67f4ceb4o6a96544852e025aa@mail.gmail.com> Message-ID: Anushri, your output (tauxobs.grd) should have 192*51*24*4 bytes (your grid size is 1degx1deg). there could be errors in the .ctl file of tauxobs.grd. cheers, Madhu 2009/2/4 anushri dutta > Respected Sir, > the command which I am using is the > > 'reinit' > 'sdfopen taux.cdf' > 'set gxout fwrite' > 'set fwrite tauxobs.grd' > 'set x 129.5 320.5' > 'set y 121.5 171.5' > 'set z 1' > 'set t 1 24' > 'd taux' > 'disable fwrite' > > and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, and > t=24, > so, where is the error in file i am not getting it, that is it display the > variable for only time t=1, not for the other time, > please help me. > > with redgards > > anushri > > > On 2/4/09, Douglas Clark wrote: >> >> Anushri, >> >> It's difficult to know what the problem is unless you give more >> information about how you created the file! >> >> There's potentially quite a lot of information that you need to provide, >> but the following at least would be useful: >> * the script or commands that you used to perform the fwrite, and what >> GrADS output to the screen when you did this >> * the ctl file you created for the output >> * the results of the command "ncdump -h taux97-98.cdf" (if you have the >> ncdump utility) >> >> In the meantime, is it possible that you only wrote out 1 time of the >> netCDF file? How big is the file tauxobs.grd (e.g. ls -l in Unix)? If the >> grid is of size nx*ny points, one year of monthly data should give a file of >> size nx*ny*12*4 bytes. If your file is smaller, you haven't written all >> times, so repeat the fwrite stage after 'set t 1 12' (or whatever is >> appropriate). Try this first! >> >> Doug >> >> >> >>> anushridutta01 at GMAIL.COM 04/02/2009 07:20 >>> >> Dear Grads User, >> Hi, >> here I am facing a problem, that i want to write a binary data for >> the 12 month by data taux97-98.cdf, through fwrite function, >> but when I pot the data through ctl file, >> it display for a single month, not for all the months, >> and the error which it showing is >> >> ga-> open taux.ctl >> Scanning description file: taux.ctl >> Data file tauxobs.grd is open as file 1 >> LON set to 64.5 160 >> LAT set to -15 10 >> LEV set to 1 1 >> Time values set: 1997:1:1:0 1997:1:1:0 >> ga-> set t 1 >> Time values set: 1997:1:1:0 1997:1:1:0 >> ga-> d taux >> Contouring: -0.08 to 0.08 interval 0.02 >> ga-> set t 2 >> Time values set: 1997:2:1:0 1997:2:1:0 >> ga-> d taux >> Low Level I/O Error: Read error on data file >> Data file name = tauxobs.grd >> Error reading 192 bytes at location 313344 >> Data Request Error: Error for variable 'taux' >> Error ocurred at column 1 >> DISPLAY error: Invalid expression >> Expression = taux >> >> >> same as for the 3 to 12 month, >> >> I am not understanding why it is so, >> Please tell me, how do I achieve this, >> >> with regards >> >> anushri >> >> >> -- >> This message (and any attachments) is for the recipient only. NERC >> is subject to the Freedom of Information Act 2000 and the contents >> of this email and any reply you make may be disclosed by NERC unless >> it is exempt from release under the Act. Any material supplied to >> NERC may be stored in an electronic records management system. >> > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/92cc04bd/attachment.html From anushridutta01 at GMAIL.COM Wed Feb 4 05:54:06 2009 From: anushridutta01 at GMAIL.COM (anushri dutta) Date: Wed, 4 Feb 2009 16:24:06 +0530 Subject: error in ctl file In-Reply-To: <726d83190902040154g45e87618s916780390786d28@mail.gmail.com> Message-ID: Respected sir, actually data is at 0.5 x 0.5 degree resolution, so nx=193 and ny=52 with regards anushri On 2/4/09, Madhusoodanan M.S. wrote: > > Anushri, > > your output (tauxobs.grd) should have 192*51*24*4 bytes (your grid size is > 1degx1deg). > > there could be errors in the .ctl file of tauxobs.grd. > > cheers, > Madhu > > 2009/2/4 anushri dutta > >> Respected Sir, >> the command which I am using is the >> >> 'reinit' >> 'sdfopen taux.cdf' >> 'set gxout fwrite' >> 'set fwrite tauxobs.grd' >> 'set x 129.5 320.5' >> 'set y 121.5 171.5' >> 'set z 1' >> 'set t 1 24' >> 'd taux' >> 'disable fwrite' >> >> and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, >> and t=24, >> so, where is the error in file i am not getting it, that is it display the >> variable for only time t=1, not for the other time, >> please help me. >> >> with redgards >> >> anushri >> >> >> On 2/4/09, Douglas Clark wrote: >>> >>> Anushri, >>> >>> It's difficult to know what the problem is unless you give more >>> information about how you created the file! >>> >>> There's potentially quite a lot of information that you need to provide, >>> but the following at least would be useful: >>> * the script or commands that you used to perform the fwrite, and what >>> GrADS output to the screen when you did this >>> * the ctl file you created for the output >>> * the results of the command "ncdump -h taux97-98.cdf" (if you have the >>> ncdump utility) >>> >>> In the meantime, is it possible that you only wrote out 1 time of the >>> netCDF file? How big is the file tauxobs.grd (e.g. ls -l in Unix)? If the >>> grid is of size nx*ny points, one year of monthly data should give a file of >>> size nx*ny*12*4 bytes. If your file is smaller, you haven't written all >>> times, so repeat the fwrite stage after 'set t 1 12' (or whatever is >>> appropriate). Try this first! >>> >>> Doug >>> >>> >>> >>> anushridutta01 at GMAIL.COM 04/02/2009 07:20 >>> >>> Dear Grads User, >>> Hi, >>> here I am facing a problem, that i want to write a binary data for >>> the 12 month by data taux97-98.cdf, through fwrite function, >>> but when I pot the data through ctl file, >>> it display for a single month, not for all the months, >>> and the error which it showing is >>> >>> ga-> open taux.ctl >>> Scanning description file: taux.ctl >>> Data file tauxobs.grd is open as file 1 >>> LON set to 64.5 160 >>> LAT set to -15 10 >>> LEV set to 1 1 >>> Time values set: 1997:1:1:0 1997:1:1:0 >>> ga-> set t 1 >>> Time values set: 1997:1:1:0 1997:1:1:0 >>> ga-> d taux >>> Contouring: -0.08 to 0.08 interval 0.02 >>> ga-> set t 2 >>> Time values set: 1997:2:1:0 1997:2:1:0 >>> ga-> d taux >>> Low Level I/O Error: Read error on data file >>> Data file name = tauxobs.grd >>> Error reading 192 bytes at location 313344 >>> Data Request Error: Error for variable 'taux' >>> Error ocurred at column 1 >>> DISPLAY error: Invalid expression >>> Expression = taux >>> >>> >>> same as for the 3 to 12 month, >>> >>> I am not understanding why it is so, >>> Please tell me, how do I achieve this, >>> >>> with regards >>> >>> anushri >>> >>> >>> -- >>> This message (and any attachments) is for the recipient only. NERC >>> is subject to the Freedom of Information Act 2000 and the contents >>> of this email and any reply you make may be disclosed by NERC unless >>> it is exempt from release under the Act. Any material supplied to >>> NERC may be stored in an electronic records management system. >>> >> >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/061268e5/attachment.html From pkarora66 at GMAIL.COM Wed Feb 4 07:02:15 2009 From: pkarora66 at GMAIL.COM (pradeep arora) Date: Wed, 4 Feb 2009 17:32:15 +0530 Subject: Findout SD In-Reply-To: <1233738215.12070.0.camel@eld313.desktop.frd.metoffice.com> Message-ID: Yes, The option given by Bernd is the best and easiest one to use. It has been tested also. PK Arora On Wed, Feb 4, 2009 at 2:33 PM, Bernd Becker wrote: > 'set gxout stat' > > play with the statistics output device. > > Good luck, > Bernd. > > > On Wed, 2009-02-04 at 00:55 -0800, sreenivas desamsetti wrote: > > > > > > Dear Sir/Madam, > > I have two different > > files with the same > > dimension. > > I wish to compute the > > Standard Deviation (SD) > > for each grid point > > basing on the time > > series of each grid > > point. > > Please let me know if > > any script available to > > compute the SD. > > Thanking you > > D Srinivas > > > > > > Desamsetti Srinivas > > Research Scholar > > Department of > > Meteorology and > > Oceanography > > Andhra University > > Visakhapatnam - 530 003 > > 91 9948530618 (M) > > sreeds82 at yahoo.com > > sreeds82 at rediffmail.com > > > > > > > > > > > -- > Bernd Becker The Monthly Outlook > Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.com - > http://www.metoffice.com > -- Wg Cdr PK Arora Joint Director (Met) AFCNWP, DIT Building (2nd Floor) Subroto Park, New Delhi 110010 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/166cba2d/attachment.html From madhusoodanan at GMAIL.COM Wed Feb 4 09:11:16 2009 From: madhusoodanan at GMAIL.COM (Madhusoodanan M.S.) Date: Wed, 4 Feb 2009 19:41:16 +0530 Subject: error in ctl file In-Reply-To: <8511ebae0902040254v2a80d7bfn13858cc86013e90f@mail.gmail.com> Message-ID: hi, your script appears to have some problems. check the header of taux.cdf (ncdump -h taux.cdf ==> it will give nx and ny.) 'set x 129.5 320.5' ==> ??? longitude will vary from 0 359.5 (if global) that means x will vary from 1 720 and y from 1 361 (if from 90S to 90N) .ctl file to read the tauxobs.grd should be correctly made then. cheers, Madhu > >>> >>> 'reinit' >>> 'sdfopen taux.cdf' >>> 'set gxout fwrite' >>> 'set fwrite tauxobs.grd' >>> 'set x 129.5 320.5' ** >>> 'set y 121.5 171.5' >>> 'set z 1' >>> 'set t 1 24' >>> 'd taux' >>> 'disable fwrite' >>> >>> and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, >>> and t=24, >>> so, where is the error in file i am not getting it, that is it display >>> the variable for only time t=1, not for the other time, >>> please help me. >>> >>> with redgards >>> >>> anushri >>> >>> >> > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/75bca173/attachment.html From antoyc at YAHOO.COM Wed Feb 4 20:04:15 2009 From: antoyc at YAHOO.COM (Antoy Chang) Date: Wed, 4 Feb 2009 17:04:15 -0800 Subject: Color legend in horizontal / vertical position (cbarn) In-Reply-To: <726d83190902040611l691bd3c2lf5a76173e31348ca@mail.gmail.com> Message-ID: Hi everybody, I have problem when displaying the 2?wave forecast using exactly the same script but?different color bar position (using cbarn) one is placed horizontally while the other vertically. How to solve this?, i want to place both vertically, the print area is exactly the same on both display.?Why this inconsistency happen?. Really appreciate your feedback.Thank You. ? p/s : Look at the attached files. ? --- On Wed, 2/4/09, Madhusoodanan M.S. wrote: From: Madhusoodanan M.S. Subject: Re: error in ctl file To: GRADSUSR at LIST.CINECA.IT Date: Wednesday, February 4, 2009, 10:11 PM hi, your script appears to have some problems. check the header of taux.cdf (ncdump -h taux.cdf ==> it will give nx and ny.) 'set x 129.5 320.5'?? ==> ??? longitude will vary from 0 359.5 (if global) that means x will vary from 1 720 and ???????????????? y????????????? from 1 361 (if from 90S to 90N) .ctl file to read the tauxobs.grd should be correctly made then. cheers, Madhu ? ? 'reinit' 'sdfopen taux.cdf' 'set gxout fwrite' 'set fwrite tauxobs.grd' 'set x 129.5 320.5' 'set y 121.5 171.5' 'set z 1' 'set t 1 24' 'd taux' 'disable fwrite' ? and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, and t=24, so, where is the error in file i am not getting it, that is it display the variable for only time t=1, not for the other time, please help me. ? with redgards ? anushri -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/dd33d41a/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: fng009.gif Type: image/gif Size: 53703 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090204/dd33d41a/attachment.gif -------------- next part -------------- A non-text attachment was scrubbed... Name: fng008_2.gif Type: image/gif Size: 44128 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090204/dd33d41a/attachment-0001.gif From antoyc at YAHOO.COM Thu Feb 5 00:07:21 2009 From: antoyc at YAHOO.COM (Antoy Chang) Date: Wed, 4 Feb 2009 21:07:21 -0800 Subject: Color legend in horizontal / vertical position (cbarn) In-Reply-To: <498A3EF9.3060500@hawaii.edu> Message-ID: Dear James, I've tried but still fail. Below is the part of gs for display the graphic : ? ? 'set grads off' ? 'set parea 0.5 7.7 0.5 9.5' ? 'set mpdset /usr/local/lib/grads/hires' ? 'set map 1 1 2' ? 'set annot 1' ? 'set gxout shaded' ? 'set clevs 'colndx ? 'set ccols 'colevs * 'd wht*1.25' ? 'd wht' ? 'run /usr/local/lib/grads/script/cbarn.gs 1.0 1 7.9 4.9' ? 'set gxout vector' ? 'set ccolor 1' ? 'set arrscl 0.15' ? 'define wamv=cos(0.0174532*wdir)' ? 'define wamu=sin(0.0174532*wdir)' ? skp=4 ? 'd skip(wamu,'skp','skp');wamv' ? 'set gxout contour' ? 'set cint 1' ? 'set clab forced' ? 'set ccolor 0' ? 'd mpr' ?? 'draw string 7.5 0.35 'metres * *'draw title MMS-WAM Significant Wave Heights(shaded contours), \ wave directions(vectors), mean wave period(seconds, white lines), *\? Valid Time '? FULL_DATE *Edited by Nursalleh on 8/8/08*************************************************************** *'draw title MMD-WAM Significant Wave Heights(metres,shaded contours),\ mean wave period(seconds, white lines), \? Valid Time ' the_day_LT ' '? FULL_DATE 'draw title MMD-WAM Significant Wave Heights(metres,shaded contours),\ mean wave period(seconds, white lines),\ Wave Direction(arrow) \? Valid Time ' the_day_LT ' '? FULL_DATE ******************************************************************************************** *'draw title 'MMS' 'WAM' 'Significant' 'Wave' 'Heights'\'Valid' 'Time' ' FULL_DATE *? 'draw title 'MMS' 'WAM' 'Significant' 'Wave' 'Heights'\'Valid' 'Time' ' vtime * *? for hard copy output * ? pfn='wam't'.pfn' ? 'enable print 'pfn ? 'print' ? 'disable print' --- On Thu, 2/5/09, James T. Potemra wrote: From: James T. Potemra Subject: Re: Color legend in horizontal / vertical position (cbarn) To: antoyc at yahoo.com Date: Thursday, February 5, 2009, 9:20 AM Antoy: I'm not sure why grads would change the color bar from vertical to horizontal if you are really not changing the page/graphics size. But, you can give the color bar function optional arguments for size, orientation and location. For example, "cbarn 0.5 0" will scale the colorbar by half (0.5) and make it horizontal (0). Here's the top of the cbarn script with more details: * The script will assume a colorbar is wanted even if there is * not room -- it will plot on the side or the bottom if there is * room in either place, otherwise it will plot along the bottom and * overlay labels there if any. This can be dealt with via * the 'set parea' command. In version 2 the default parea will * be changed, but we want to guarantee upward compatibility in * sub-releases. * * * modifications by mike fiorino 940614 * * - the extreme colors are plotted as triangles * - the colors are boxed in white * - input arguments in during a run execution: * * run cbarn sf vert xmid ymid * * sf - scale the whole bar 1.0 = original 0.5 half the size, etc. * vert - 0 FORCES a horizontal bar = 1 a vertical bar * xmid - the x position on the virtual page the center the bar * ymid - the x position on the virtual page the center the bar * * if vert,xmid,ymid are not specified, they are selected * as in the original algorithm * Hope this helps, Jim Antoy Chang wrote: > Hi everybody, > I have problem when displaying the 2 wave forecast using exactly the same script but different color bar position (using cbarn) one is placed horizontally while the other vertically. How to solve this?, i want to place both vertically, the print area is exactly the same on both display. Why this inconsistency happen?. Really appreciate your feedback.Thank You. > p/s : Look at the attached files. > --- On *Wed, 2/4/09, Madhusoodanan M.S. //* wrote: > > From: Madhusoodanan M.S. > Subject: Re: error in ctl file > To: GRADSUSR at LIST.CINECA.IT > Date: Wednesday, February 4, 2009, 10:11 PM > > hi, > > your script appears to have some problems. > > check the header of taux.cdf (ncdump -h taux.cdf ==> it will give > nx and ny.) > > 'set x 129.5 320.5' ==> ??? longitude will vary from 0 359.5 (if > global) > that means x will vary from 1 720 and > y from 1 361 (if from 90S to 90N) > > .ctl file to read the tauxobs.grd should be correctly made then. > > cheers, > Madhu > > > > > > 'reinit' > 'sdfopen taux.cdf' > 'set gxout fwrite' > 'set fwrite tauxobs.grd' > 'set x 129.5 320.5' > 'set y 121.5 171.5' > 'set z 1' > 'set t 1 24' > 'd taux' > 'disable fwrite' > and the data length in tauxobs.grd is 193*52*24*4, > where nx=193, ny=52, and t=24, > so, where is the error in file i am not getting it, > that is it display the variable for only time t=1, not > for the other time, > please help me. > with redgards > anushri > > > > > > > ------------------------------------------------------------------------ > > > ------------------------------------------------------------------------ > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090204/7148fb54/attachment.html From imspj at HAWAII.EDU Thu Feb 5 08:47:27 2009 From: imspj at HAWAII.EDU (Ilukpitiya Mudi S P Jayawardena) Date: Thu, 5 Feb 2009 18:47:27 +0500 Subject: how to plot a specific contour line with different thickness In-Reply-To: <8511ebae0902040254v2a80d7bfn13858cc86013e90f@mail.gmail.com> Message-ID: Hi Could you please tell me how to plot a specific contour line ( eg. 0 contour) with different thickness and different style. thanks Ilukpitiye -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/78a0cded/attachment.html From stiperski at CIRUS.DHZ.HR Thu Feb 5 08:59:49 2009 From: stiperski at CIRUS.DHZ.HR (Ivana Stiperski) Date: Thu, 5 Feb 2009 14:59:49 +0100 Subject: Script for arbitrary cross-sections with terrain In-Reply-To: Message-ID: Hi As far as I know a grads script that draws terrain in an arbitrary cross section has not been compiled yet. So I've made one :) and am attaching it for general use. hope it serves its purpose! the prerequisite is that you have either a terrain variable (in this script called topo) or something that can mimic the terrain in your grib file. cheers Ivana ------------------------------------------------ Ivana Stiperski Graduate Research Assistant Hydrological and Meteorological Service, Croatia Gric 3, HR-10000 Zagreb ph. + 385 1 4565 721 email: stiperski at cirus.dhz.hr -------------- next part -------------- A non-text attachment was scrubbed... Name: ArbitraryCSS_Terrain.gs Type: application/octet-stream Size: 3665 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090205/9960536c/attachment.obj From imspj at HAWAII.EDU Thu Feb 5 09:03:56 2009 From: imspj at HAWAII.EDU (Ilukpitiya Mudi S P Jayawardena) Date: Thu, 5 Feb 2009 19:03:56 +0500 Subject: plotting a specific contour Message-ID: Hi Could you please help me to plot a specific contour line ( eg. 0 contour) with different thickness and different style. thanks Ilukpitiye -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/ae58a531/attachment.html From ousmane at IRI.COLUMBIA.EDU Thu Feb 5 09:10:43 2009 From: ousmane at IRI.COLUMBIA.EDU (Ousmane Ndiaye) Date: Thu, 5 Feb 2009 09:10:43 -0500 Subject: plotting a specific contour In-Reply-To: Message-ID: use those three : set clevs lev (lev for you =0) set cstyle style (the style) set cthick thckns (the thickness) for more see http://grads.iges.org/grads/gadoc/gadocindex.html Ousmane ----- Original Message ----- From: "Ilukpitiya Mudi S P Jayawardena" To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, February 5, 2009 9:03:56 AM GMT -05:00 US/Canada Eastern Subject: plotting a specific contour Hi Could you please help me to plot a specific contour line ( eg. 0 contour) with different thickness and different style. thanks Ilukpitiye From teddyallen at YAHOO.COM Thu Feb 5 09:40:36 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Thu, 5 Feb 2009 06:40:36 -0800 Subject: how to plot a specific contour line with different thickness In-Reply-To: Message-ID: Dear Ilukpitiya Mudi S P Jayawardena, While this may not 100% answer your question, it may put you on the right track. You may be able to isolate your 0 line by "blacking out" the other contours. This would allow you to control your 0 line exclusively without affecting the other contours. ? 'set gxout contour' 'set ccolor 0'??? DEFINES THE COLOR OF YOUR CONTOUR 'set cthick 19'? DEFINES THE THICKNESS OF YOUR CONTOUR 'set black 1010 1016' DEFINES WHAT CONTOURS ARE "blacked-out" in this case I blacked out the 1010-1016 isobars ? It is somewhat of a cavalier approach, but it gets the job done. Probably a better way. Hope it helps. Teddy Allen ? ? web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Thu, 2/5/09, Ilukpitiya Mudi S P Jayawardena wrote: From: Ilukpitiya Mudi S P Jayawardena Subject: how to plot a specific contour line with different thickness To: GRADSUSR at LIST.CINECA.IT Date: Thursday, February 5, 2009, 8:47 AM Hi ? Could you please tell me how to plot a specific? contour? line ( eg. 0 contour) ?with different? thickness and different style. ? thanks ? Ilukpitiye -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/f2ed2159/attachment.html From jma at COLA.IGES.ORG Thu Feb 5 10:55:42 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 5 Feb 2009 10:55:42 -0500 Subject: New release GrADS version 2.0.a5 Message-ID: Hi, Everyone -- I have posted a new version of GrADS: 2.0.a5. The GrADS downloads web page (http://iges.org/grads/downloads.html) has been updated with links to the source code and pre-compiled binaries for a few OS flavors. This build has the initial implementation of a new capability -- graphics output that is compatible with GIS software. You can create files in the GeoTIFF format (a georeferenced raster image) as well as the KML format, which is actually two files: an image file in TIFF format and a text file in KML which contains the georeferencing information. The KML file is intended for use with Google Earth. There is some documentation on how to use these new graphics output options: http://iges.org/grads/gadoc/gradcomdsetgxout.html http://iges.org/grads/gadoc/gradcomdsetgeotiff.html http://iges.org/grads/gadoc/gradcomdsetkml.html There are other fixes and additions to this release, outlined below: Bug Fixes: * PDEF lcc now works properly in the Southern Hemisphere * sdfwrite output has correct time axis units Misc: * Added support for templating on T but not E for self-describing files * Added support for optional padding of sdfwrite output to include fixed dimensions * Added check in gribmap for non-supported grib2 product definition templates * Added support for jpeg image output format with printim If you build GrADS from source, please note there are some changes to the supplibs: I've added libtiff and libgeotiff and upgraded g2clib to version 1.1.8. Please see http://iges.org/grads/gadoc/supplibs.html for links to source tarballs. Enjoy! --Jennifer -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/504ba323/attachment.html From mmacleod at SCOTIAWEATHER.COM Thu Feb 5 11:40:05 2009 From: mmacleod at SCOTIAWEATHER.COM (mmacleod) Date: Thu, 5 Feb 2009 16:40:05 +0000 Subject: New release GrADS version 2.0.a5 In-Reply-To: <8647ECD8-97B2-4344-93CA-9F882288C049@cola.iges.org> Message-ID: Good afternoon Jennifer. The expanded capability of 2.0.a5 looks good. We have an application that would benefit from the KML and GeoTIFF formats. Our system operate on Ubuntu 8.4 O/S. How soon will the Ubuntu x8 - 64 binaries be available? Will you be sending an e-mail when they are posted to the web pages for download. Keep up the good work. The user community really appreciate it. Regards Mac MacLeod Jennifer Adams wrote: > Hi, Everyone -- > I have posted a new version of GrADS: 2.0.a5. The GrADS downloads web > page (_http://iges.org/grads/downloads.html_) has been updated with > links to the source code and pre-compiled binaries for a few OS flavors. > > This build has the initial implementation of a new capability -- > graphics output that is compatible with GIS software. You can create > files in the GeoTIFF format (a georeferenced raster image) as well as > the KML format, which is actually two files: an image file in TIFF > format and a text file in KML which contains the georeferencing > information. The KML file is intended for use with Google Earth. There > is some documentation on how to use these new graphics output options: > > http://iges.org/grads/gadoc/gradcomdsetgxout.html > http://iges.org/grads/gadoc/gradcomdsetgeotiff.html > http://iges.org/grads/gadoc/gradcomdsetkml.html > > There are other fixes and additions to this release, outlined below: > > Bug Fixes: > * PDEF lcc now works properly in the Southern Hemisphere > * sdfwrite output has correct time axis units > Misc: > * Added support for templating on T but not E for self-describing files > * Added support for optional padding of sdfwrite output to include fixed dimensions > * Added check in gribmap for non-supported grib2 product definition templates > * Added support for jpeg image output format with printim > If you build GrADS from source, please note there are some changes to the supplibs: I've added libtiff and libgeotiff and upgraded g2clib to version 1.1.8. Please see http://iges.org/grads/gadoc/supplibs.html for links to source tarballs. > Enjoy! > --Jennifer > > -- > Jennifer M. Adams > IGES/COLA > 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705 > jma at cola.iges.org > > > -- M.A. (Mac) MacLeod President and General Manager Scotia Weather Services Inc 192 Wyse Road, Suite 8, Dartmouth, N.S. B3A 1M9 Tele: 902-468-3866 Fax: 902-461-1768 E-mail: mmacleod at scotiaweather.com Visit us: www.scotiaweather.com From mmacleod at SCOTIAWEATHER.COM Thu Feb 5 11:47:40 2009 From: mmacleod at SCOTIAWEATHER.COM (mmacleod) Date: Thu, 5 Feb 2009 16:47:40 +0000 Subject: Question on new release of GrADS 2.0.a5 Message-ID: Good afternoon Jennifer. My original post came back as rejected by the GRADS-USER Server but you may have received it but here it si again. The expanded capability of 2.0.a5 looks good. We have an application that would benefit from the KML and GeoTIFF formats. Our system operate on Ubuntu 8.4 O/S. How soon will the Ubuntu x8 - 64 binaries be available? Will you be sending an e-mail when they are posted to the web pages for download. Keep up the good work. The user community really appreciate it. Regards Mac MacLeod -- M.A. (Mac) MacLeod President and General Manager Scotia Weather Services Inc 192 Wyse Road, Suite 8, Dartmouth, N.S. B3A 1M9 Tele: 902-468-3866 Fax: 902-461-1768 E-mail: mmacleod at scotiaweather.com Visit us: www.scotiaweather.com From jma at COLA.IGES.ORG Thu Feb 5 13:50:43 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 5 Feb 2009 13:50:43 -0500 Subject: Question on new release of GrADS 2.0.a5 In-Reply-To: <498B182C.3040709@scotiaweather.com> Message-ID: On Feb 5, 2009, at 11:47 AM, mmacleod wrote: > Good afternoon Jennifer. > > > My original post came back as rejected by the GRADS-USER Server but > you > may have received it but here it is again. Yes, that 'rejected post' message is a known problem -- I have asked the list owners to look into it. The messages do go through, no need to send them twice. > The expanded capability of 2.0.a5 looks good. We have an application > that would benefit from the KML and GeoTIFF formats. > Our system operate on Ubuntu 8.4 O/S. How soon will the Ubuntu x8 - > 64 > binaries be available? > Will you be sending an e-mail when they are posted to the web pages > for > download. I have just added an ubuntu x86_64 build. Windows, 32-bit linux, and other unix flavors are provided by volunteers in the community. Jennifer From mmacleod at SCOTIAWEATHER.COM Thu Feb 5 13:54:07 2009 From: mmacleod at SCOTIAWEATHER.COM (mmacleod) Date: Thu, 5 Feb 2009 18:54:07 +0000 Subject: Question on new release of GrADS 2.0.a5 In-Reply-To: <8DAFD7F4-1921-47B8-9A9A-F56A8217E7AB@cola.iges.org> Message-ID: Thanks Jennifer. We will be trying it within the next few days. Regards, Mac Jennifer Adams wrote: > On Feb 5, 2009, at 11:47 AM, mmacleod wrote: > >> Good afternoon Jennifer. >> >> >> My original post came back as rejected by the GRADS-USER Server but >> you >> may have received it but here it is again. > Yes, that 'rejected post' message is a known problem -- I have asked > the list owners to look into it. The messages do go through, no need > to send them twice. > >> The expanded capability of 2.0.a5 looks good. We have an application >> that would benefit from the KML and GeoTIFF formats. >> Our system operate on Ubuntu 8.4 O/S. How soon will the Ubuntu x8 - >> 64 >> binaries be available? >> Will you be sending an e-mail when they are posted to the web pages >> for >> download. > I have just added an ubuntu x86_64 build. Windows, 32-bit linux, and > other unix flavors are provided by volunteers in the community. > > Jennifer -- M.A. (Mac) MacLeod President and General Manager Scotia Weather Services Inc 192 Wyse Road, Suite 8, Dartmouth, N.S. B3A 1M9 Tele: 902-468-3866 Fax: 902-461-1768 E-mail: mmacleod at scotiaweather.com Visit us: www.scotiaweather.com From hrc at EAS.GATECH.EDU Thu Feb 5 17:11:42 2009 From: hrc at EAS.GATECH.EDU (Hai-Ru Chang) Date: Thu, 5 Feb 2009 17:11:42 -0500 Subject: Invalid QUERY argument: udct Message-ID: Hello, I installed grads v2.0.a4.oga.1 and tried to add libipc.gex extension into /usr/local/bin/gex. The GAUDXT is defined as GAUDXT=/ usr/local/bin/gex/udxt. But when I started grads and type 'query udct', I saw the message of Invalid QUERY argument: udct. Can anyone tell me what is happening? Thanks! Dr. Hai-Ru Chang School of Earth and Atmospheric Sciences Georgia Institute of Technology Atlanta, GA 30332-0340 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/85f2058b/attachment.html From antoyc at YAHOO.COM Thu Feb 5 20:51:27 2009 From: antoyc at YAHOO.COM (Antoy Chang) Date: Thu, 5 Feb 2009 17:51:27 -0800 Subject: Color legend in horizontal / vertical position (cbarn) In-Reply-To: Message-ID: Douglas/all, Thank you for your support, i manage to solve the problem. This is portion of my cbarn.gs used to plot my clour bar : ********************************************************************************************************* *? Script to plot a colorbar * *? The script will assume a colorbar is wanted even if there is *? not room -- it will plot on the side or the bottom if there is *? room in either place, otherwise it will plot along the bottom and *? overlay labels there if any.? This can be dealt with via *? the 'set parea' command.? In version 2 the default parea will *? be changed, but we want to guarantee upward compatibility in *? sub-releases. * * *?????? modifications by mike fiorino 940614 *?????? rvp added by diane stokes 20040729 *?????? cinc added by diane stokes 20050130 * *?????? - the extreme colors are plotted as triangles *?????? - the colors are boxed in white *?????? - input arguments in during a run execution: * *?????? run cbarn cinc sf vert xmid ymid rvp * *?????? cinc - increment for labels *?????? sf?? - scale the whole bar 1.0 = original 0.5 half the size, etc. *?????? vert - 0 FORCES a horizontal bar = 1 a vertical bar *?????? xmid - the x position on the virtual page the center the bar *?????? ymid - the x position on the virtual page the center the bar *?????? rvp? - 1 resets virtual page (one colorbar for several plots) * *?????? if cinc is not specified, every color interval is labeled *?????? if vert,xmid,ymid are not specified, they are selected *?????? as in the original algorithm * *?????? if rvp is not specified, it is 0 (don't reset virtual page) ******************************************************************************************************* if you can see, the cbarn arguments are cinc sf vert xmid ymid rvp, there is additional argument before the sf. in my script below i omit the cinc argument, i solve the problem by assign 1 for cinc argument. ? here is my script : ?'set grads off' ? 'set parea 0.5 7.7 0.5 9.5' ? 'set mpdset /usr/local/lib/grads/hires' ? 'set map 1 1 2' ? 'set annot 1' ? 'set gxout shaded' ? 'set clevs 'colndx ? 'set ccols 'colevs * 'd wht*1.25' ? 'd wht' ? 'run /usr/local/lib/grads/script/cbarn.gs 1.0 1 7.9 4.9'--->corrected cbarn.gs 1 1.0 1 7.9 4.9 ? 'set gxout vector' ? 'set ccolor 1' ? 'set arrscl 0.15' ? 'define wamv=cos(0.0174532*wdir)' ? 'define wamu=sin(0.0174532*wdir)' ? skp=4 ? 'd skip(wamu,'skp','skp');wamv' ? 'set gxout contour' ? 'set cint 1' ? 'set clab forced' ? 'set ccolor 0' ? 'd mpr' ?? 'draw string 7.5 0.35 'metres ? But i don't know why the i encounter the problem?since 2 days ago after using the same script for plotting since 2005..... ? Regards, Antoy --- On Thu, 2/5/09, Douglas Clark wrote: From: Douglas Clark Subject: Re: Color legend in horizontal / vertical position (cbarn) To: antoyc at yahoo.com Date: Thursday, February 5, 2009, 5:33 PM Sorry - just saw your later post in which you say the optional arguments didn't work. I've no idea. Doug >>> Douglas Clark 05/02/2009 09:28 >>> Antoy, No idea why the bar changes if everything else really is identical as you said. Could be some hard-to-find quirk. Can't you solve this using the optional arguments to cbarn? Form cbarn.gs: * run cbarn sf vert xmid ymid * * sf - scale the whole bar 1.0 = original 0.5 half the size, etc. * vert - 0 FORCES a horizontal bar = 1 a vertical bar * xmid - the x position on the virtual page the center the bar * ymid - the x position on the virtual page the center the bar So set vert=1 to ensure a vertical bar (and then xmid,ymid if it still moves about at all). Doug >>> Antoy Chang 05/02/2009 01:04 >>> Hi everybody, I have problem when displaying the 2 wave forecast using exactly the same script but different color bar position (using cbarn) one is placed horizontally while the other vertically. How to solve this?, i want to place both vertically, the print area is exactly the same on both display. Why this inconsistency happen?. Really appreciate your feedback.Thank You. p/s : Look at the attached files. --- On Wed, 2/4/09, Madhusoodanan M.S. wrote: From: Madhusoodanan M.S. Subject: Re: error in ctl file To: GRADSUSR at LIST.CINECA.IT Date: Wednesday, February 4, 2009, 10:11 PM hi, your script appears to have some problems. check the header of taux.cdf (ncdump -h taux.cdf ==> it will give nx and ny.) 'set x 129.5 320.5' ==> ??? longitude will vary from 0 359.5 (if global) that means x will vary from 1 720 and y from 1 361 (if from 90S to 90N) .ctl file to read the tauxobs.grd should be correctly made then. cheers, Madhu 'reinit' 'sdfopen taux.cdf' 'set gxout fwrite' 'set fwrite tauxobs.grd' 'set x 129.5 320.5' 'set y 121.5 171.5' 'set z 1' 'set t 1 24' 'd taux' 'disable fwrite' and the data length in tauxobs.grd is 193*52*24*4, where nx=193, ny=52, and t=24, so, where is the error in file i am not getting it, that is it display the variable for only time t=1, not for the other time, please help me. with redgards anushri -- This message (and any attachments) is for the recipient only. NERC is subject to the Freedom of Information Act 2000 and the contents of this email and any reply you make may be disclosed by NERC unless it is exempt from release under the Act. Any material supplied to NERC may be stored in an electronic records management system. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090205/8f753be8/attachment.html From arlindo.dasilva at GMAIL.COM Fri Feb 6 10:26:38 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Fri, 6 Feb 2009 10:26:38 -0500 Subject: Invalid QUERY argument: udct In-Reply-To: <7E9EB912-54B8-4A32-A30A-019A50314388@eas.gatech.edu> Message-ID: On Thu, Feb 5, 2009 at 5:11 PM, Hai-Ru Chang wrote: > Hello, I installed grads v2.0.a4.oga.1 and tried to add libipc.gex > extension into /usr/local/bin/gex. The GAUDXT is defined > as GAUDXT=/usr/local/bin/gex/udxt. But when I started grads and type 'query > udct', I saw the message of Invalid QUERY argument: udct. Can anyone tell > me what is happening? Thanks! > The OpenGrADS Extensions are not (yet) available with GrADS v2.0. Am I planning a "preview" release with v2.0.a5 in the next couple of weeks. However, you can use PyGrADS with GrADS v2.0 without having libipc. Arlindo -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090206/1317a538/attachment.html From cassar at PILOTBOOK.EU Fri Feb 6 10:36:13 2009 From: cassar at PILOTBOOK.EU (Pierre Cassar) Date: Fri, 6 Feb 2009 16:36:13 +0100 Subject: New release GrADS version 2.0.a5 Message-ID: > I have posted a new version of GrADS: 2.0.a5. > The GrADS downloads web page (http://iges.org/grads/downloads.html) > has been updated with links to the source code > and pre-compiled binaries for a few OS flavors. Do you plan to build a Freebsd 7 version in the near future? Thank you very much. Pierre Cassar From nieves at LTP.UPR.CLU.EDU Fri Feb 6 12:06:19 2009 From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves) Date: Fri, 6 Feb 2009 13:06:19 -0400 Subject: New release GrADS version 2.0.a5 In-Reply-To: Message-ID: > > I have posted a new version of GrADS: 2.0.a5. > > The GrADS downloads web page (http://iges.org/grads/downloads.html) > > has been updated with links to the source code > > and pre-compiled binaries for a few OS flavors. > > Do you plan to build a Freebsd 7 version in the near future? > Thank you very much. > > Pierre Cassar I plan to build the FreeBSD-7.0 and 7.1 versions as soon as the opengrads port is available. Jose From sreeds82 at YAHOO.COM Sun Feb 8 01:43:35 2009 From: sreeds82 at YAHOO.COM (sreenivas desamsetti) Date: Sat, 7 Feb 2009 22:43:35 -0800 Subject: Finding Correlation coefficient Message-ID: Dear Sir/Madam, I have two different files with the same dimension (x,y,t). I wish to compute the Standard Deviation (SD) and correlation coefficient. Lat time I asked the same question in the GrADS User List They suggested me to use 'set gxout stat' command. If I use the command it gives only one value of SD (it founds from all the available grid points (spatial) at one time. But I need to find the temporal correlation (temporal values have to be used for the computation of SD) and so finally we will get a spatial map of SD and Correlation coefficient. Please suggest me to get this. If any Mathematical functions are available in GrADS please let me know. Otherwise any GrADS script file will help me to find them. Thanking you D Srinivas Desamsetti Srinivas Research Scholar Department of Meteorology and Oceanography Andhra University Visakhapatnam - 530 003 91 9491686740 (M) sreeds82 at yahoo.com sreeds82 at rediffmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090207/8e3a7db8/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/jpeg Size: 11090 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090207/8e3a7db8/attachment.jpe From prjayakrishnan at GMAIL.COM Sun Feb 8 04:59:56 2009 From: prjayakrishnan at GMAIL.COM (Jayakrishnan PR) Date: Sun, 8 Feb 2009 15:29:56 +0530 Subject: Finding Correlation coefficient In-Reply-To: <49633.66927.qm@web36402.mail.mud.yahoo.com> Message-ID: Dear Srinivas, Actually I have not worked in statistical correlation methods till now. But I know that Prof.H.S.Ram Mohan and Prof.K.MohanKumar in my Department has worked a lot in this area. Their students Mr.Rajesh.J. and Lorna.P.Nayagam has published one paper regarding these techniques in an international journal. You can contact any of the professors for the clarification of this doubts. Both of them are friendly to me. On Sun, Feb 8, 2009 at 12:13 PM, sreenivas desamsetti wrote: > Dear Sir/Madam, > I have two different files with the same dimension (x,y,t). > I wish to compute the Standard Deviation (SD) and correlation coefficient. > Lat time I asked the same question in the GrADS User List > They suggested me to use 'set gxout stat' command. > If I use the command it gives only one value of SD (it founds from all the > available > grid points (spatial) at one time. > But I need to find the temporal correlation (temporal values have to be > used for the computation of SD) and so finally we will get > a spatial map of SD and Correlation coefficient. > Please suggest me to get this. > If any Mathematical functions are available in GrADS please let me know. > Otherwise any GrADS script file will help me to find them. > Thanking you > D Srinivas > Desamsetti Srinivas > Research Scholar > Department of Meteorology and Oceanography > Andhra University > Visakhapatnam - 530 003 > 91 9491686740 (M) > sreeds82 at yahoo.com > sreeds82 at rediffmail.com > > -- Sincerely *********************************************** Jayakrishnan.P.R CSIR-Research Fellow Department of Atmospheric Sciences Cochin University of Science and Technology (CUSAT), Cochin-682 016 Kerala, India. Mob: 09895417565 *********************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090208/7ed6acd9/attachment.html From am_kareem2001 at YAHOO.COM Mon Feb 9 01:14:47 2009 From: am_kareem2001 at YAHOO.COM (ahmed kareem) Date: Sun, 8 Feb 2009 22:14:47 -0800 Subject: correlation problem Message-ID: -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.ctl Type: application/octet-stream Size: 262 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: stn.ctl Type: application/octet-stream Size: 365 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment-0001.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: MAMZ.BIN Type: application/octet-stream Size: 720 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment-0002.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: sst6190.oct Type: application/octet-stream Size: 1922400 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment-0003.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: mam.GS Type: application/octet-stream Size: 1274 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090208/d10c24af/attachment-0004.obj From prjayakrishnan at GMAIL.COM Mon Feb 9 01:43:03 2009 From: prjayakrishnan at GMAIL.COM (Jayakrishnan PR) Date: Mon, 9 Feb 2009 12:13:03 +0530 Subject: correlation problem In-Reply-To: <295988.7243.qm@web34202.mail.mud.yahoo.com> Message-ID: Dear Kareem, It is nice to see that you yourself have solved the problem. I was about to solve your problem today and I started to work with that. Then only your mail came. Best wishes for your works. May god bless you and your guide. On Mon, Feb 9, 2009 at 11:44 AM, ahmed kareem wrote: > > -- Sincerely *********************************************** Jayakrishnan.P.R CSIR-Research Fellow Department of Atmospheric Sciences Cochin University of Science and Technology (CUSAT), Cochin-682 016 Kerala, India. Mob: 09895417565 *********************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090209/de817e2f/attachment.html From ucfasra at UCL.AC.UK Tue Feb 10 02:51:38 2009 From: ucfasra at UCL.AC.UK (Srivatsan Raghavan) Date: Tue, 10 Feb 2009 07:51:38 +0000 Subject: ocean mask Message-ID: Hi All, I have WRF generated precipitation fields over south east asia and would like to remove the data over the sea/ocean to make it more easily comparable to CRU, such as data only over land. Is there a way I can use GrADS for this? any scripts ? Thank you Srivatsan Raghavan Research Fellow Tropical Marine Science Institute National University of Singapore From prjayakrishnan at GMAIL.COM Tue Feb 10 03:44:59 2009 From: prjayakrishnan at GMAIL.COM (Jayakrishnan PR) Date: Tue, 10 Feb 2009 14:14:59 +0530 Subject: ocean mask In-Reply-To: Message-ID: Dear Srivatsan, Actually I am now in the learning stage of WRF and I dont know much about modelling. But I think the grads script basemap.gs will be suitable for the purpose. You can try with that. The following link will give the grads script. Hope you are doing well. Good luck. http://grads.iges.org/grads/gadoc/library.html On Tue, Feb 10, 2009 at 1:21 PM, Srivatsan Raghavan wrote: > Hi All, > > I have WRF generated precipitation fields over south east asia and would > like to remove the data over the sea/ocean to make it more easily > comparable to CRU, such as data only over land. > > Is there a way I can use GrADS for this? any scripts ? > > Thank you > > > > Srivatsan Raghavan > Research Fellow > Tropical Marine Science Institute > National University of Singapore > -- Sincerely *********************************************** Jayakrishnan.P.R CSIR-Research Fellow Department of Atmospheric Sciences Cochin University of Science and Technology (CUSAT), Cochin-682 016 Kerala, India. Mob: 09895417565 *********************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090210/d447240c/attachment.html From ucfasra at UCL.AC.UK Tue Feb 10 04:58:45 2009 From: ucfasra at UCL.AC.UK (Srivatsan Raghavan) Date: Tue, 10 Feb 2009 09:58:45 +0000 Subject: basemap files In-Reply-To: <5b5e360902100044s7419e85akaf8889b4d22347e7@mail.gmail.com> Message-ID: Hi All, I found that basemap.gs solves the problem of masking, as I wanted to mask waters in the southeast asian region, output derived from WRF. Despite the opoly and lpoly files in the directory that has fonts and udunits.dat, running it gives an error. ------------------ ga-> set mpdset lowres MPDSET file name = lowres ga-> d rainc+rainnc Contouring: 1000 to 14000 interval 1000 ga-> basemap O 7 2 L Error reading opoly_lowres.asc ga-> d rainc+rainnc ------------------------ I have set the GADDIR env. variable to the right path. Cannot figure out where I am wrong. Any help ? Many Thanks Srivatsan From anna.zacharioudaki at PLYMOUTH.AC.UK Tue Feb 10 12:57:17 2009 From: anna.zacharioudaki at PLYMOUTH.AC.UK (Anna Zacharioudaki) Date: Tue, 10 Feb 2009 17:57:17 +0000 Subject: have many grads files with same name, how can i give them different names and open them in grads? Message-ID: Dear users, I am running the WaveWatch III model for wave hindcasting. The model includes a gridded output post-processor for GRADS. The output files from the post-processor are: ww3.grads ww3.ctl I am running the model on a cluster and thus performing many runs simultaneously. I want then to compare/manipulate the output of several runs using GRADS. My problem is that all output files from the different runs share the same (see above). When I try to rename the file, e.g. ww3test.grads and ww3test.ctl and open it in grads it does not open because it does not recognise the new name. Is there a simple way to rename grads files? Or any other suggestions of how I could deal with this problem? Thank you, Anna Dr Anna Zacharioudaki Centre for Coastal Dynamics and Engineering, School of Engineering, University of Plymouth, Drake Circus, Plymouth, Devon PL4 8AA Te1: +44 (0)1752 586131 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090210/495d9e53/attachment.html From paul.earnshaw at METOFFICE.GOV.UK Tue Feb 10 13:04:55 2009 From: paul.earnshaw at METOFFICE.GOV.UK (Paul Earnshaw) Date: Tue, 10 Feb 2009 18:04:55 +0000 Subject: have many grads files with same name, how can i give them different names and open them in grads? In-Reply-To: <1038F6B357C5004FA619C32FB3DB99574956F128F6@ILS130.uopnet.plymouth.ac.uk> Message-ID: Anna, You can change the names of the files as you suggest, but there is a further change you need to make. There is a line in the .ctl file which references the .grads file. You will need to edit this line so it matches the new file name. The .ctl file is an ascii text file so it can be editted like any other text file, and the line will need to change should be the first line of the file. Paul On Tue, 2009-02-10 at 17:57 +0000, Anna Zacharioudaki wrote: > Dear users, > > > > I am running the WaveWatch III model for wave hindcasting. The model > includes a gridded output post-processor for GRADS. The output files > from the post-processor are: > > ww3.grads > > ww3.ctl > > I am running the model on a cluster and thus performing many runs > simultaneously. I want then to compare/manipulate the output of > several runs using GRADS. My problem is that all output files from the > different runs share the same (see above). > > When I try to rename the file, e.g. ww3test.grads and ww3test.ctl and > open it in grads it does not open because it does not recognise the > new name. > > > > Is there a simple way to rename grads files? Or any other suggestions > of how I could deal with this problem? > > > > Thank you, > > > > Anna > > > > > > Dr Anna Zacharioudaki > > Centre for Coastal Dynamics and Engineering, > > School of Engineering, > > University of Plymouth, > > Drake Circus, > Plymouth, Devon PL4 8AA > > Te1: +44 (0)1752 586131 > > > > -- Dr. Paul Earnshaw Global Model Development and Diagnostics Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel: +44 (0)1392 884525 Fax: +44 (0)1392 885681 email: paul.earnshaw at metoffice.gov.uk website: www.metoffice.gov.uk Met Office climate change predictions can now be viewed on Google Earth http://www.metoffice.gov.uk/research/hadleycentre/google/ From antoyc at YAHOO.COM Tue Feb 10 23:57:49 2009 From: antoyc at YAHOO.COM (Antoy Chang) Date: Tue, 10 Feb 2009 20:57:49 -0800 Subject: Finer grid wind interpolation Message-ID: Hi all, Sorry my my vague question. Just wondering if anyone here know hoe to do the wind interpolation to a finer grid (eg from 1 deg to 0.5 deg)?. Your response is highly appreciated ? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090210/6a2a7069/attachment.html From antoyc at YAHOO.COM Tue Feb 10 23:59:43 2009 From: antoyc at YAHOO.COM (Antoy Chang) Date: Tue, 10 Feb 2009 20:59:43 -0800 Subject: Finer grid wind interpolation Message-ID: Hi all., I just wondering if anyone here know how to do the wind field interpolation to a finer grid. Sorry for my vague question. Thank You very much for your response. ? Regards, Antoy -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090210/05e94ee4/attachment.html From stergios.misios at ZMAW.DE Wed Feb 11 10:06:14 2009 From: stergios.misios at ZMAW.DE (Stergios Misios) Date: Wed, 11 Feb 2009 16:06:14 +0100 Subject: Findout SD In-Reply-To: <905777.14810.qm@web36402.mail.mud.yahoo.com> Message-ID: You can do it easily with the following segment, mean=ave(variable,t=start_time,t=end_time) std=sqrt(ave(pow(variable-mean,2)),t=start_time,t=last_time)*(number_of_time_steps/(number_of_time_steps-1)) Good luck, Stergios From j.burn at GEOS.COM Wed Feb 11 11:00:51 2009 From: j.burn at GEOS.COM (Burn James) Date: Wed, 11 Feb 2009 17:00:51 +0100 Subject: Build GrADS 2.0.a5 on Opensuse 11.1 Message-ID: Hi I'm trying to build Grads 2.0 on a home PC running Opensuse 11.1 Linux. I've unpacked the latest source code to home/grads2.0 Unpacked the supp files from data.tar into usr/local/lib/grads Unpacked the supplib files from ftp://grads.iges.org/grads/Supplibs/2.0/bin/i686-pc-linux-gnu_supplibs2.0.tar.gz to /usr/local/lib/grads/supplibs When I run configure as below, it seems to find the the supplibs folder, (can see files such as readline.h), but doesn't find udunits.h and other files needed to enable an opendap install. I also had this problem when I was trying to build on Opensuse 10.2 but I "got around it" by finding some precompiled version of Grads which worked when I had the supplibs folder as above. However, now I'd like to be able to build Grads 2.0 myself now. I've started looking through the configure file, which helped me work out that I need to define the location of SUPPLIBS as per output below, but I can get no further. Help will be well appreciated, thanks James ******************* burn at linux-62ha:~> cd grads2.0 burn at linux-62ha:~/grads2.0> sh configure SUPPLIBS=/usr/local/lib/grads/supplibs_cola checking for a BSD-compatible install... /usr/bin/install -c checking whether build environment is sane... yes checking for a thread-safe mkdir -p... /bin/mkdir -p checking for gawk... gawk checking whether make sets $(MAKE)... yes checking whether to enable maintainer-specific portions of Makefiles... no Checking for programs ------------------ checking for gawk... (cached) gawk checking for gcc... gcc checking for C compiler default output file name... a.out checking whether the C compiler works... yes checking whether we are cross compiling... no checking for suffix of executables... checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking for g++... no checking for c++... no checking for gpp... no checking for aCC... no checking for CC... no checking for cxx... no checking for cc++... no checking for cl.exe... no checking for FCC... no checking for KCC... no checking for RCC... no checking for xlC_r... no checking for xlC... no checking whether we are using the GNU C++ compiler... no checking whether g++ accepts -g... no checking for a BSD-compatible install... /usr/bin/install -c checking whether ln -s works... yes Setting host specific options ------------------------- checking build system type... i686-pc-linux-gnu checking host system type... i686-pc-linux-gnu Host machine: i686-pc-linux-gnu Checking system libraries, headers, and compiler features ----------------------------------------- checking for cos in -lm... yes checking how to run the C preprocessor... gcc -E checking for X... no checking for grep that handles long lines and -e... /usr/bin/grep checking for egrep... /usr/bin/grep -E checking for ANSI C header files... yes checking for sys/types.h... yes checking for sys/stat.h... yes checking for stdlib.h... yes checking for string.h... yes checking for memory.h... yes checking for strings.h... yes checking for inttypes.h... yes checking for stdint.h... yes checking for unistd.h... yes checking malloc.h usability... yes checking malloc.h presence... yes checking for malloc.h... yes checking sys/file.h usability... yes checking sys/file.h presence... yes checking for sys/file.h... yes checking for unistd.h... (cached) yes checking for an ANSI C-conforming const... yes checking whether byte ordering is bigendian... no checking for strtod... yes checking for off_t... yes checking for special C compiler options needed for large files... no checking for _FILE_OFFSET_BITS value needed for large files... 64 checking for _LARGEFILE_SOURCE value needed for large files... no Checking for optional feature support ------------------------------------ checking for supplibs directory... /usr/local/lib/grads/supplibs_cola checking whether to use dynamic linking... yes (Warning: binaries may not be suitable for distribution) Checking for libraries to support GUI interface ... checking for main in -lXext... no checking for main in -lXt... no checking for pkg-config... /usr/bin/pkg-config checking pkg-config is at least version 0.9.0... yes checking for XAW7... no checking for main in -lXt... (cached) no - GUI disabled Checking in supplibs for libraries to support command line editing ... checking readline.h usability... no checking readline.h presence... no checking for readline.h... no Checking in system locations for libraries to support command line editing ... checking for a readline compatible library... no - Command line editing disabled Checking in supplibs for libraries to support geotiff ... checking tiff.h usability... no checking tiff.h presence... no checking for tiff.h... no checking tiffio.h usability... no checking tiffio.h presence... no checking for tiffio.h... no - geotiff disabled Checking in supplibs for libraries to support printim ... checking zlib.h usability... no checking zlib.h presence... no checking for zlib.h... no checking for main in -lpng12... no checking for main in -ljpeg... no Checking in system locations for libraries to support printim ... checking for GD... no checking for gdlib-config... no checking gd.h usability... no checking gd.h presence... no checking for gd.h... no - printim disabled (Warning: gxtran will not be built) Checking in supplibs for libraries to support grib2 ... checking for zlib.h... (cached) no checking for main in -lpng12... (cached) no Checking in system locations for grib2 libraries ... checking grib2.h usability... no checking grib2.h presence... no checking for grib2.h... no - grib2 disabled Checking in supplibs for libraries to support hdf4 ... checking for zlib.h... (cached) no checking udunits.h usability... no checking udunits.h presence... no checking for udunits.h... no checking for main in -ljpeg... (cached) no Checking in system locations for libraries to support hdf ... checking for udunits.h... (cached) no - hdf4 disabled Checking in supplibs for libraries to support netcdf ... checking for udunits.h... (cached) no Checking in system locations for libraries to support netcdf ... checking for udunits.h... (cached) no - netcdf disabled Checking in supplibs for libraries to support OPeNDAP access ... checking for udunits.h... (cached) no Checking in system locations for libraries to support OPeNDAP access ... checking for udunits.h... (cached) no checking gadap.h usability... no checking gadap.h presence... no checking for gadap.h... no - gradsdap (OPeNDAP/netCDF) build disabled +=========================================================================+ GrADS 2.0.a5 : Built Wed Feb 11 15:46:34 GMT 2009 for i686-pc-linux-gnu +=========================================================================+ +----------------------------------+ | | | Configuration Summary | | | | - readline disabled | | - printim disabled | | - netcdf disabled | | - hdf4 disabled | | - grib2 disabled | | - geotiff disabled | | - GUI disabled | | | | Build Summary | | | | + grads enabled | | - gradsdap (OPeNDAP) disabled | | | | + Dynamic linking enabled | | | +----------------------------------+ configure: creating src/VERSION configure: creating src/buildinfo.h configure: creating ./config.status config.status: creating Makefile config.status: creating src/Makefile config.status: creating src/config.h config.status: src/config.h is unchanged Run 'make && make install' to build GrADS. Executables will be installed to bin/ in current dir, unless prefix was specified. From stergios.misios at ZMAW.DE Wed Feb 11 11:17:48 2009 From: stergios.misios at ZMAW.DE (Stergios Misios) Date: Wed, 11 Feb 2009 17:17:48 +0100 Subject: How to set xlabels in a date/time axis In-Reply-To: <020e01c9808e$f6854410$d8867f0a@ad.arpa.piemonte.it> Message-ID: one not very clever solution (but working..) could be, 'set xaxis 1 41 4' 'set xlabs dat1|date2|date3|.....' The drawback is that you have to define the title of every time step.... hope it helps, stergios From jma at COLA.IGES.ORG Wed Feb 11 11:29:30 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Wed, 11 Feb 2009 11:29:30 -0500 Subject: Build GrADS 2.0.a5 on Opensuse 11.1 In-Reply-To: <20090211160151.B7B9D21288@mx2.cineca.it> Message-ID: On Feb 11, 2009, at 11:00 AM, Burn James wrote: > Hi > > I'm trying to build Grads 2.0 on a home PC running Opensuse 11.1 > Linux. > > I've unpacked the latest source code to home/grads2.0 > Unpacked the supp files from data.tar into usr/local/lib/grads > Unpacked the supplib files from > ftp://grads.iges.org/grads/Supplibs/2.0/bin/i686-pc-linux-gnu_supplibs2.0.tar.gz > to /usr/local/lib/grads/supplibs Down below it says you have the environment variable SUPPLIBS set to / usr/local/grads/supplibs_cola/ Which is it? Underneath the $SUPPLIBS directory, you should have to subdirs, ./ include and ./lib with all the libraries and include files. Are these present under /usr/local/grads/supplibs_cola? Jennifer > > > When I run configure as below, it seems to find the the supplibs > folder, > (can see files such as readline.h), but doesn't find udunits.h and > other > files needed to enable an opendap install. I also had this problem > when I > was trying to build on Opensuse 10.2 but I "got around it" by > finding some > precompiled version of Grads which worked when I had the supplibs > folder as > above. However, now I'd like to be able to build Grads 2.0 myself > now. > I've started looking through the configure file, which helped me > work out > that I need to define the location of SUPPLIBS as per output below, > but I > can get no further. > > Help will be well appreciated, thanks > > James > ******************* > > > > burn at linux-62ha:~> cd grads2.0 > burn at linux-62ha:~/grads2.0> sh configure > SUPPLIBS=/usr/local/lib/grads/supplibs_cola > checking for a BSD-compatible install... /usr/bin/install -c > > checking whether build environment is sane... yes > > checking for a thread-safe mkdir -p... /bin/mkdir -p > > checking for gawk... gawk > > checking whether make sets $(MAKE)... yes > > checking whether to enable maintainer-specific portions of > Makefiles... no > > > Checking for programs > ------------------ > checking for gawk... (cached) gawk > checking for gcc... gcc > checking for C compiler default output file name... a.out > checking whether the C compiler works... yes > checking whether we are cross compiling... no > checking for suffix of executables... > checking for suffix of object files... o > checking whether we are using the GNU C compiler... yes > checking whether gcc accepts -g... yes > checking for gcc option to accept ISO C89... none needed > checking for g++... no > checking for c++... no > checking for gpp... no > checking for aCC... no > checking for CC... no > checking for cxx... no > checking for cc++... no > checking for cl.exe... no > checking for FCC... no > checking for KCC... no > checking for RCC... no > checking for xlC_r... no > checking for xlC... no > checking whether we are using the GNU C++ compiler... no > checking whether g++ accepts -g... no > checking for a BSD-compatible install... /usr/bin/install -c > checking whether ln -s works... yes > > Setting host specific options > ------------------------- > checking build system type... i686-pc-linux-gnu > checking host system type... i686-pc-linux-gnu > Host machine: i686-pc-linux-gnu > > Checking system libraries, headers, and compiler features > ----------------------------------------- > checking for cos in -lm... yes > checking how to run the C preprocessor... gcc -E > checking for X... no > > checking for grep that handles long lines and -e... /usr/bin/grep > checking for egrep... /usr/bin/grep -E > checking for ANSI C header files... yes > checking for sys/types.h... yes > checking for sys/stat.h... yes > checking for stdlib.h... yes > checking for string.h... yes > checking for memory.h... yes > checking for strings.h... yes > checking for inttypes.h... yes > checking for stdint.h... yes > checking for unistd.h... yes > checking malloc.h usability... yes > checking malloc.h presence... yes > checking for malloc.h... yes > checking sys/file.h usability... yes > checking sys/file.h presence... yes > checking for sys/file.h... yes > checking for unistd.h... (cached) yes > > checking for an ANSI C-conforming const... yes > checking whether byte ordering is bigendian... no > checking for strtod... yes > checking for off_t... yes > checking for special C compiler options needed for large files... no > checking for _FILE_OFFSET_BITS value needed for large files... 64 > checking for _LARGEFILE_SOURCE value needed for large files... no > > Checking for optional feature support > ------------------------------------ > checking for supplibs directory... /usr/local/lib/grads/supplibs_cola > checking whether to use dynamic linking... yes (Warning: binaries > may not be > suitable for distribution) > > Checking for libraries to support GUI interface ... > checking for main in -lXext... no > checking for main in -lXt... no > checking for pkg-config... /usr/bin/pkg-config > checking pkg-config is at least version 0.9.0... yes > checking for XAW7... no > checking for main in -lXt... (cached) no > - GUI disabled > > Checking in supplibs for libraries to support command line editing ... > checking readline.h usability... no > checking readline.h presence... no > checking for readline.h... no > > Checking in system locations for libraries to support command line > editing ... > checking for a readline compatible library... no > - Command line editing disabled > > Checking in supplibs for libraries to support geotiff ... > checking tiff.h usability... no > checking tiff.h presence... no > checking for tiff.h... no > checking tiffio.h usability... no > checking tiffio.h presence... no > checking for tiffio.h... no > - geotiff disabled > > Checking in supplibs for libraries to support printim ... > checking zlib.h usability... no > checking zlib.h presence... no > checking for zlib.h... no > checking for main in -lpng12... no > checking for main in -ljpeg... no > > Checking in system locations for libraries to support printim ... > checking for GD... no > checking for gdlib-config... no > checking gd.h usability... no > checking gd.h presence... no > checking for gd.h... no > - printim disabled (Warning: gxtran will not be built) > > Checking in supplibs for libraries to support grib2 ... > checking for zlib.h... (cached) no > checking for main in -lpng12... (cached) no > > Checking in system locations for grib2 libraries ... > checking grib2.h usability... no > checking grib2.h presence... no > checking for grib2.h... no > - grib2 disabled > > Checking in supplibs for libraries to support hdf4 ... > checking for zlib.h... (cached) no > checking udunits.h usability... no > checking udunits.h presence... no > checking for udunits.h... no > checking for main in -ljpeg... (cached) no > > Checking in system locations for libraries to support hdf ... > checking for udunits.h... (cached) no > - hdf4 disabled > > Checking in supplibs for libraries to support netcdf ... > checking for udunits.h... (cached) no > > Checking in system locations for libraries to support netcdf ... > checking for udunits.h... (cached) no > - netcdf disabled > > Checking in supplibs for libraries to support OPeNDAP access ... > checking for udunits.h... (cached) no > > Checking in system locations for libraries to support OPeNDAP > access ... > checking for udunits.h... (cached) no > checking gadap.h usability... no > checking gadap.h presence... no > checking for gadap.h... no > - gradsdap (OPeNDAP/netCDF) build disabled > > > + > = > = > = > = > =====================================================================+ > GrADS 2.0.a5 : Built Wed Feb 11 15:46:34 GMT 2009 for i686-pc- > linux-gnu > > + > = > = > = > = > =====================================================================+ > > +----------------------------------+ > | | > | Configuration Summary | > | | > | - readline disabled | > | - printim disabled | > | - netcdf disabled | > | - hdf4 disabled | > | - grib2 disabled | > | - geotiff disabled | > | - GUI disabled | > | | > | Build Summary | > | | > | + grads enabled | > | - gradsdap (OPeNDAP) disabled | > | | > | + Dynamic linking enabled | > | | > +----------------------------------+ > configure: creating src/VERSION > configure: creating src/buildinfo.h > configure: creating ./config.status > config.status: creating Makefile > config.status: creating src/Makefile > config.status: creating src/config.h > config.status: src/config.h is unchanged > > Run 'make && make install' to build GrADS. > Executables will be installed to bin/ in current dir, unless prefix > was > specified. -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090211/4a7ccad5/attachment.html From jma at COLA.IGES.ORG Wed Feb 11 11:31:30 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Wed, 11 Feb 2009 11:31:30 -0500 Subject: New release GrADS version 2.0.a5 In-Reply-To: <8647ECD8-97B2-4344-93CA-9F882288C049@cola.iges.org> Message-ID: Dear All, Hoop has provided a new i686 build for 2.0.a5 linked to from the GrADS downloads page. Try this for your 32-bit linux needs. --Jennifer From j.burn at GEOS.COM Wed Feb 11 12:17:34 2009 From: j.burn at GEOS.COM (Burn James) Date: Wed, 11 Feb 2009 18:17:34 +0100 Subject: Build GrADS 2.0.a5 on Opensuse 11.1 Message-ID: Hi Sorry my mistake, I was experimenting and posted the wrong logfile. It is usr/local/grads/supplibs/ that I have setup and that contains ./lib and ./include That gave ... Checking in supplibs for libraries to support command line editing ... checking readline.h usability... yes checking readline.h presence... yes checking for readline.h... yes checking history.h usability... yes checking history.h presence... yes checking for history.h... yes ...but also... Checking in system locations for grib2 libraries ... checking grib2.h usability... no checking grib2.h presence... no checking for grib2.h... no - grib2 disabled Checking in supplibs for libraries to support hdf4 ... checking for zlib.h... (cached) no checking udunits.h usability... no checking udunits.h presence... no checking for udunits.h... no Checking in system locations for libraries to support hdf ... checking for udunits.h... (cached) no - hdf4 disabled which is odd because I can see the udunits.h and grib2.h files in /usr/local/lib/grads/supplibs/include Still puzzled. Thanks James On Wed, 11 Feb 2009 11:29:30 -0500, Jennifer Adams wrote: > >On Feb 11, 2009, at 11:00 AM, Burn James wrote: > >> Hi >> >> I'm trying to build Grads 2.0 on a home PC running Opensuse 11.1 >> Linux. >> >> I've unpacked the latest source code to home/grads2.0 >> Unpacked the supp files from data.tar into usr/local/lib/grads >> Unpacked the supplib files from >> ftp://grads.iges.org/grads/Supplibs/2.0/bin/i686-pc-linux-gnu_supplibs2.0.tar.gz >> to /usr/local/lib/grads/supplibs >Down below it says you have the environment variable SUPPLIBS set to / >usr/local/grads/supplibs_cola/ >Which is it? >Underneath the $SUPPLIBS directory, you should have to subdirs, ./ >include and ./lib with all the libraries and include files. Are these >present under /usr/local/grads/supplibs_cola? >Jennifer > >> >> >> When I run configure as below, it seems to find the the supplibs >> folder, >> (can see files such as readline.h), but doesn't find udunits.h and >> other >> files needed to enable an opendap install. I also had this problem >> when I >> was trying to build on Opensuse 10.2 but I "got around it" by >> finding some >> precompiled version of Grads which worked when I had the supplibs >> folder as >> above. However, now I'd like to be able to build Grads 2.0 myself >> now. >> I've started looking through the configure file, which helped me >> work out >> that I need to define the location of SUPPLIBS as per output below, >> but I >> can get no further. >> >> Help will be well appreciated, thanks >> >> James >> ******************* >> >> >> >> burn at linux-62ha:~> cd grads2.0 >> burn at linux-62ha:~/grads2.0> sh configure >> SUPPLIBS=/usr/local/lib/grads/supplibs_cola >> checking for a BSD-compatible install... /usr/bin/install -c >> >> checking whether build environment is sane... yes >> >> checking for a thread-safe mkdir -p... /bin/mkdir -p >> >> checking for gawk... gawk >> >> checking whether make sets $(MAKE)... yes >> >> checking whether to enable maintainer-specific portions of >> Makefiles... no >> >> >> Checking for programs >> ------------------ >> checking for gawk... (cached) gawk >> checking for gcc... gcc >> checking for C compiler default output file name... a.out >> checking whether the C compiler works... yes >> checking whether we are cross compiling... no >> checking for suffix of executables... >> checking for suffix of object files... o >> checking whether we are using the GNU C compiler... yes >> checking whether gcc accepts -g... yes >> checking for gcc option to accept ISO C89... none needed >> checking for g++... no >> checking for c++... no >> checking for gpp... no >> checking for aCC... no >> checking for CC... no >> checking for cxx... no >> checking for cc++... no >> checking for cl.exe... no >> checking for FCC... no >> checking for KCC... no >> checking for RCC... no >> checking for xlC_r... no >> checking for xlC... no >> checking whether we are using the GNU C++ compiler... no >> checking whether g++ accepts -g... no >> checking for a BSD-compatible install... /usr/bin/install -c >> checking whether ln -s works... yes >> >> Setting host specific options >> ------------------------- >> checking build system type... i686-pc-linux-gnu >> checking host system type... i686-pc-linux-gnu >> Host machine: i686-pc-linux-gnu >> >> Checking system libraries, headers, and compiler features >> ----------------------------------------- >> checking for cos in -lm... yes >> checking how to run the C preprocessor... gcc -E >> checking for X... no >> >> checking for grep that handles long lines and -e... /usr/bin/grep >> checking for egrep... /usr/bin/grep -E >> checking for ANSI C header files... yes >> checking for sys/types.h... yes >> checking for sys/stat.h... yes >> checking for stdlib.h... yes >> checking for string.h... yes >> checking for memory.h... yes >> checking for strings.h... yes >> checking for inttypes.h... yes >> checking for stdint.h... yes >> checking for unistd.h... yes >> checking malloc.h usability... yes >> checking malloc.h presence... yes >> checking for malloc.h... yes >> checking sys/file.h usability... yes >> checking sys/file.h presence... yes >> checking for sys/file.h... yes >> checking for unistd.h... (cached) yes >> >> checking for an ANSI C-conforming const... yes >> checking whether byte ordering is bigendian... no >> checking for strtod... yes >> checking for off_t... yes >> checking for special C compiler options needed for large files... no >> checking for _FILE_OFFSET_BITS value needed for large files... 64 >> checking for _LARGEFILE_SOURCE value needed for large files... no >> >> Checking for optional feature support >> ------------------------------------ >> checking for supplibs directory... /usr/local/lib/grads/supplibs_cola >> checking whether to use dynamic linking... yes (Warning: binaries >> may not be >> suitable for distribution) >> >> Checking for libraries to support GUI interface ... >> checking for main in -lXext... no >> checking for main in -lXt... no >> checking for pkg-config... /usr/bin/pkg-config >> checking pkg-config is at least version 0.9.0... yes >> checking for XAW7... no >> checking for main in -lXt... (cached) no >> - GUI disabled >> >> Checking in supplibs for libraries to support command line editing ... >> checking readline.h usability... no >> checking readline.h presence... no >> checking for readline.h... no >> >> Checking in system locations for libraries to support command line >> editing ... >> checking for a readline compatible library... no >> - Command line editing disabled >> >> Checking in supplibs for libraries to support geotiff ... >> checking tiff.h usability... no >> checking tiff.h presence... no >> checking for tiff.h... no >> checking tiffio.h usability... no >> checking tiffio.h presence... no >> checking for tiffio.h... no >> - geotiff disabled >> >> Checking in supplibs for libraries to support printim ... >> checking zlib.h usability... no >> checking zlib.h presence... no >> checking for zlib.h... no >> checking for main in -lpng12... no >> checking for main in -ljpeg... no >> >> Checking in system locations for libraries to support printim ... >> checking for GD... no >> checking for gdlib-config... no >> checking gd.h usability... no >> checking gd.h presence... no >> checking for gd.h... no >> - printim disabled (Warning: gxtran will not be built) >> >> Checking in supplibs for libraries to support grib2 ... >> checking for zlib.h... (cached) no >> checking for main in -lpng12... (cached) no >> >> Checking in system locations for grib2 libraries ... >> checking grib2.h usability... no >> checking grib2.h presence... no >> checking for grib2.h... no >> - grib2 disabled >> >> Checking in supplibs for libraries to support hdf4 ... >> checking for zlib.h... (cached) no >> checking udunits.h usability... no >> checking udunits.h presence... no >> checking for udunits.h... no >> checking for main in -ljpeg... (cached) no >> >> Checking in system locations for libraries to support hdf ... >> checking for udunits.h... (cached) no >> - hdf4 disabled >> >> Checking in supplibs for libraries to support netcdf ... >> checking for udunits.h... (cached) no >> >> Checking in system locations for libraries to support netcdf ... >> checking for udunits.h... (cached) no >> - netcdf disabled >> >> Checking in supplibs for libraries to support OPeNDAP access ... >> checking for udunits.h... (cached) no >> >> Checking in system locations for libraries to support OPeNDAP >> access ... >> checking for udunits.h... (cached) no >> checking gadap.h usability... no >> checking gadap.h presence... no >> checking for gadap.h... no >> - gradsdap (OPeNDAP/netCDF) build disabled >> >> >> + >> = >> = >> = >> = >> =====================================================================+ >> GrADS 2.0.a5 : Built Wed Feb 11 15:46:34 GMT 2009 for i686-pc- >> linux-gnu >> >> + >> = >> = >> = >> = >> =====================================================================+ >> >> +----------------------------------+ >> | | >> | Configuration Summary | >> | | >> | - readline disabled | >> | - printim disabled | >> | - netcdf disabled | >> | - hdf4 disabled | >> | - grib2 disabled | >> | - geotiff disabled | >> | - GUI disabled | >> | | >> | Build Summary | >> | | >> | + grads enabled | >> | - gradsdap (OPeNDAP) disabled | >> | | >> | + Dynamic linking enabled | >> | | >> +----------------------------------+ >> configure: creating src/VERSION >> configure: creating src/buildinfo.h >> configure: creating ./config.status >> config.status: creating Makefile >> config.status: creating src/Makefile >> config.status: creating src/config.h >> config.status: src/config.h is unchanged >> >> Run 'make && make install' to build GrADS. >> Executables will be installed to bin/ in current dir, unless prefix >> was >> specified. > >-- >Jennifer M. Adams >IGES/COLA >4041 Powder Mill Road, Suite 302 >Calverton, MD 20705 >jma at cola.iges.org > > > > From Frank_Colby at UML.EDU Wed Feb 11 13:05:57 2009 From: Frank_Colby at UML.EDU (Frank Colby) Date: Wed, 11 Feb 2009 13:05:57 -0500 Subject: script available? Message-ID: Hi Everyone, Does anyone have a script to compute frontogenesis? Thanks, Frank Colby UMass Lowell From arlindo.dasilva at GMAIL.COM Wed Feb 11 14:08:01 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Wed, 11 Feb 2009 14:08:01 -0500 Subject: Finer grid wind interpolation In-Reply-To: <320412.77412.qm@web33108.mail.mud.yahoo.com> Message-ID: On Tue, Feb 10, 2009 at 11:59 PM, Antoy Chang wrote: > Hi all., > I just wondering if anyone here know how to do the wind field interpolation > to a finer grid. Sorry for my vague question. Thank You very much for your > response. > > > The question for you is how accurate do you want to be? For a quick plot linear interpolation is fine; poles are always tricky, depending on whether your the poles are in the center or edge of the gridbox. Look for the built in linterp() function in grads v2, and the UDFs re()/regrid2() in grads v1.x. A proper wind interpolation should take as input both u and v components; currently grads functions can only return one argument, so you would need 2 functions: uinterp(u,v), and vinterp(u,v). I am not aware of any function implementing vector interpolation in grads, so you would need to write a user defined function for that. And as of this writing, only grads v1.x supports UDFs. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090211/4ad8316f/attachment.html From teddyallen at YAHOO.COM Thu Feb 12 11:45:20 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Thu, 12 Feb 2009 08:45:20 -0800 Subject: download troubles In-Reply-To: <4992E966.1080301@zmaw.de> Message-ID: Dear GrADS users, I had to uninstall GrADS v2.0? today and am now experiencing problems trying to install the v2.0.a3 for Windows XP. I remember in the past I was directed to an execution file that prompted me to simply "run" the program and then everything was automatically set up. Now when I extract the .zip files for version 2.0.a3, there is no indication of a "run" or execute file and the program is not visible on the programs start menu. Does anybody out there have any directions on how to download this version for the NON-COMPUTER expert? Thank you very much and I apologize for the amateur-ish posting. sincerely, teddy -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090212/188e6c5d/attachment.html From teddyallen at YAHOO.COM Thu Feb 12 13:54:40 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Thu, 12 Feb 2009 10:54:40 -0800 Subject: GrADS vanishing window In-Reply-To: <4994694F.8040008@uml.edu> Message-ID: Frank, Thank you for your support with this. I am running windows XP. Everything seems to download perfectly from the "superpack", but the window vanishes after the landscape / portrait prompt. Unfortunately, the window vanishes too quickly for me to read the related error message...It appears if there is some sort of disconnect after the landscape / portrait prompt for some mysterious reason....so, GrADS 2.0 is installed but is unable to function. Thank you again for your help, teddy web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Thu, 2/12/09, Frank Colby wrote: From: Frank Colby Subject: Re: download troubles To: teddyallen at yahoo.com Date: Thursday, February 12, 2009, 1:24 PM Terry, I had this problem with a computer with Vista -- had to run the installation as administrator (right click on the .exe file, and choose the "run as administrator" option). Are you running Vista? Frank -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090212/2dc49026/attachment.html From Ron.Goodson at EC.GC.CA Thu Feb 12 19:10:27 2009 From: Ron.Goodson at EC.GC.CA (Goodson,Ron [Edm]) Date: Thu, 12 Feb 2009 17:10:27 -0700 Subject: Anyone Have Script for ageostrophic wind components In-Reply-To: A<4992E966.1080301@zmaw.de> Message-ID: Hi Just wondering if anybody has a script for the ageostrophic wind already coded-up. Mostly just interested in isallobaric wind and inertial-advective wind. If not, I guess I have to do a bit of work to do. tks ron From Rick.Danielson at PHYS.OCEAN.DAL.CA Thu Feb 12 20:55:50 2009 From: Rick.Danielson at PHYS.OCEAN.DAL.CA (Rick Danielson) Date: Thu, 12 Feb 2009 21:55:50 -0400 Subject: Anyone Have Script for ageostrophic wind components In-Reply-To: <4A6AB38B55B49C44A22E021A83CBEDDB06937D4E@sr-pnr-exch3.prairie.int.ec.gc.ca> Message-ID: Hi Ron, The attached script should get you started, at least. Regards, Rick > Hi > > Just wondering if anybody has a script for the ageostrophic wind already > coded-up. Mostly just interested in isallobaric wind and > inertial-advective wind. If not, I guess I have to do a bit of work to > do. > > tks > > ron "I'd put my money on the sun During 1956-1965, "the world was and solar energy. What a consuming 4 billion barrels of oil per source of power! I hope we year and the average discovery was don't have to wait until oil around 30 billion. Today we consume and coal run out before we 30 billion barrels per year and the tackle that." - Edison to Ford discovery rate is now approaching 4 and Firestone in 1931 (from billion barrels of crude oil per year." "Uncommon Friends" by J.Newton) - Kjell Aleklett of Uppsala University -------------- next part -------------- * This script computes one of various diagnostic quantities using * basic NCEP Reanalysis grids stored in a template netcdf file. * An assumption is that the calculation domain of interest has been * specified - RD, Feb `99 function keb() * query file for dimensions of current viewing domain "q dims" rec = sublin(result,5) if (subwrd(rec,3) = "varying") tmin = subwrd(rec,11) tmax = subwrd(rec,13) tvar = "varying" else tval = subwrd(rec,9) tvar = "fixed" endif rec = sublin(result,4) if (subwrd(rec,3) = "varying") zmin = subwrd(rec,11) zmax = subwrd(rec,13) zvar = "varying" else zval = subwrd(rec,9) zvar = "fixed" endif rec = sublin(result,3) if (subwrd(rec,3) = "varying") ymin = subwrd(rec,11) ymax = subwrd(rec,13) yvar = "varying" else yval = subwrd(rec,9) yvar = "fixed" endif rec = sublin(result,2) if (subwrd(rec,3) = "varying") xmin = subwrd(rec,11) xmax = subwrd(rec,13) xvar = "varying" else xval = subwrd(rec,9) xvar = "fixed" endif * query file for dimensions of calculation domain "d timemin" rec = sublin(result,1) timemin = subwrd(rec,4) "d timemax" rec = sublin(result,1) timemax = subwrd(rec,4) "d levmin" rec = sublin(result,1) levmin = subwrd(rec,4) "d levmax" rec = sublin(result,1) levmax = subwrd(rec,4) "d moistmax" rec = sublin(result,1) moistmax = subwrd(rec,4) "d latmin" rec = sublin(result,1) latmin = subwrd(rec,4) "d latmax" rec = sublin(result,1) latmax = subwrd(rec,4) "d lonmin" rec = sublin(result,1) lonmin = subwrd(rec,4) "d lonmax" rec = sublin(result,1) lonmax = subwrd(rec,4) say "" say "calculating quasi geostrophic terms, pointwise" say "" say "using the following calculation domain:" say "time from "timemin" to "timemax say "level from "levmin" to "levmax say "lat from "latmin" to "latmax say "lon from "lonmin" to "lonmax say "" say "defining temporary variables" ; say "" defa = "EARTH = 6.371e6" say " "defa "define "defa defa = "D2R = 3.141592654 / 180.0" say " "defa "define "defa defa = "OMEG = 7.292e-5" say " "defa "define "defa defa = "GRAV = 9.8" say " "defa "define "defa say "" say "computing over lat and lon:" ; say "" "set t 1" "set z 1" "set lat "latmin" "latmax "set lon "lonmin" "lonmax defa = "delx = EARTH * cdiff(lon,x) * D2R * cos(lat*D2R)" say " "defa ; say "" "define "defa defa = "dely = EARTH * cdiff(lat,y) * D2R" say " "defa ; say "" "define "defa defa = "coriol = 2.0 * OMEG * sin(lat*D2R)" say " "defa ; say "" "define "defa say "computing over time, level, lat, and lon:" ; say "" "set t "timemin" "timemax "set lev "levmin" "levmax "set lat "latmin" "latmax "set lon "lonmin" "lonmax defa = "geou = -GRAV / coriol * cdiff(hgt,y) / dely" say " "defa ; say "" "define "defa defa = "geov = GRAV / coriol * cdiff(hgt,x) / delx" say " "defa ; say "" "define "defa defa = "newu = -GRAV / coriol * cdiff(hgt,y) / dely * (uwnd * -GRAV / coriol * cdiff(hgt,y) / dely + " defb = "vwnd * GRAV / coriol * cdiff(hgt,x) / delx) / " defc = "(-GRAV / coriol * cdiff(hgt,y) / dely * -GRAV / coriol * cdiff(hgt,y) / dely + " defd = "GRAV / coriol * cdiff(hgt,x) / delx * GRAV / coriol * cdiff(hgt,x) / delx)" say " "defa" "defb" "defc" "defd ; say "" "define "defa" "defb" "defc" "defd defa = "nevu = GRAV / coriol * cdiff(hgt,x) / delx * (uwnd * -GRAV / coriol * cdiff(hgt,y) / dely + " defb = "vwnd * GRAV / coriol * cdiff(hgt,x) / delx) / " defc = "(-GRAV / coriol * cdiff(hgt,y) / dely * -GRAV / coriol * cdiff(hgt,y) / dely + " defd = "GRAV / coriol * cdiff(hgt,x) / delx * GRAV / coriol * cdiff(hgt,x) / delx)" say " "defa" "defb" "defc" "defd ; say "" "define "defa" "defb" "defc" "defd say "undefining temporary variables" ; say "" "undefine EARTH" "undefine D2R" "undefine OMEG" "undefine GRAV" say "resetting the original viewing domain" ; say "" if (tvar = "varying") "set t "tmin" "tmax else "set t "tval endif if (zvar = "varying") "set z "zmin" "zmax else "set z "zval endif if (yvar = "varying") "set y "ymin" "ymax else "set y "yval endif if (xvar = "varying") "set x "xmin" "xmax else "set x "xval endif say "re-creating the display script" ; say "" "run gui_list" From fcheng at HRC-LAB.ORG Thu Feb 12 22:13:39 2009 From: fcheng at HRC-LAB.ORG (Fang-Yi Cheng) Date: Thu, 12 Feb 2009 19:13:39 -0800 Subject: does any one know how to plot wind vector in a timeseries form Message-ID: Dear all, Could you please let me know how to plot wind vector in a timeseries format? Thank you. Fang From luo at JAMSTEC.GO.JP Fri Feb 13 03:21:54 2009 From: luo at JAMSTEC.GO.JP (Jing-Jia Luo) Date: Fri, 13 Feb 2009 17:21:54 +0900 Subject: plot specific month or seasonal mean data from monthly time series Message-ID: Sorry to disturb you. Using grads v1.9b4, I have monthly time series data (say from Jan. 1950 to Dec. 2008) and want to plot only May or Jun-Aug mean data each year. I am wondering whether you know how to do so and would kindly give me your advice. Thanks a lot. Best regards, Jing-Jia Luo From Rodney.Teck at NUIM.IE Fri Feb 13 03:59:56 2009 From: Rodney.Teck at NUIM.IE (Rodney Teck) Date: Fri, 13 Feb 2009 08:59:56 +0000 Subject: script available? In-Reply-To: <49931385.6040602@uml.edu> Message-ID: Fank I've attached a script which I came across sometime back it may not do exactly as you want but see how you get on. Rod Teck+ On Wed, 2009-02-11 at 18:05, Frank Colby wrote: > Hi Everyone, > > Does anyone have a script to compute frontogenesis? > > Thanks, > > Frank Colby > UMass Lowell -- Dr.Rodney Teck ICARUS Climatic Research Unit St.Ann's NUI Maynooth Co Kildare Ireland tel: ++353 (1)7083769 email: rodney.teck at nuim.ie http://icarus.nuim.ie -------------- next part -------------- 'reinit' 'open lipr8710.ctl' 'open lisf8710.ctl' 'set lon 280 360' 'set lat 20 72.5' 'set t 5' 'set map 1' 'set gxout shaded' 'rhmax=max(rh, lev=850, lev=700)' 'n=100-(100-rhmax)*2' defcol_clouds_l() 'set csmooth linear' 'd n' 'run cbarc 1.8 7.6 0' 'set lev 700' 'set ccolor 16' 'set gxout contour' 'set cstyle 1' 'set cint 20' 'set clab on' 'set clskip 3' 'd z/10' 'set ccolor 2' 'set gxout contour' 'set cstyle 5' 'set cint 20' 'set clab on' 'set clskip 5' 'd 0.1*(z(lev=500)-z(lev=1000))' 'set csmooth linear' 'set gxout contour' 'set lev 850' 'dtx = cdiff(tmp,x) ' 'dty = cdiff(tmp,y) ' 'dt = sqrt(dtx*dtx+dty*dty)' 'dttx = cdiff(dt,x) ' 'dtty = cdiff(dt,y) ' 'dtt = sqrt(dttx*dttx+dtty*dtty)' 'rhmax=max(rh, lev=850, lev=500)' 'fr = dt*0.15*rhmax' 'dx = cdiff(lon,x)*3.1416/180 ' 'dy = cdiff(lat,y)*3.1416/180 ' 'adv = -1*((u*dtx)/(cos(lat*3.1416/180)*dx)+v*dty/dy)/6.37e6' 'zadv = abs(adv)/adv' defcol_fr_warm() 'd fr*zadv' defcol_fr_cold() 'd fr*zadv' pull dummy 'set clab forced' 'set cthick 6' 'set cstyle 3' 'set clskip 3' 'd mslp.2/100' *disp_pres() 'draw title Frontal Zones(AT-850), 1.5-3 km Clouds[%],H700,H500/1000' *'printim at112_8.png x800 y500 black' *'enable print f:/front110105.gmf' *'print' *'disable print' function defcol_clouds_l() 'set rgb 81 10 10 10' 'set rgb 82 20 20 20' 'set rgb 83 30 30 30' 'set rgb 84 40 40 40' 'set rgb 85 50 50 50' 'set rgb 86 60 60 60' 'set rgb 87 70 70 70' 'set rgb 88 80 80 80' 'set rgb 89 90 90 90' 'set rgb 90 100 100 100' 'set rgb 91 110 110 110' 'set rgb 92 120 120 120' 'set rgb 93 130 130 130' 'set rgb 94 140 140 140' 'set rgb 95 150 150 150' 'set rgb 96 160 160 160' 'set rgb 97 170 170 170' 'set rgb 98 180 180 180' 'set rgb 99 190 190 190' 'set clevs 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100' 'set ccols 0 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99' return function defcol_fr_warm() 'run define_colors' 'set clevs 32' 'set ccols 28' 'set cthick 10' 'set clab off' return function defcol_fr_cold() 'set clevs -29' 'set ccols 47' 'set clab off' return From teddyallen at YAHOO.COM Fri Feb 13 05:47:49 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Fri, 13 Feb 2009 02:47:49 -0800 Subject: plot specific month or seasonal mean data from monthly time series In-Reply-To: <20090213.172154.118628731.luo@jamstec.go.jp> Message-ID: Dear Jing-Jia, I have done the same in the past by following the below lines in a script: ? 'd ave (VARIABLE, t=FIRST MONTH OF INTEREST, t=LAST MONTH OF INTEREST, 12)' ? Essentially this says to average whatever the variable is from the first month of interest to the last month of interest, but averages at increments of 12 so that only the specific month is averaged and not all months. For example, if you had a time series where the first time step was Jan 1979 (t=1) and the final was?Nov 2016 (t=700)? [these are just hypothetical numbers] then you would write: 'd ave (variable, t=1, t=700, 12) This would average only the January months in your monthly dataset. I would guess for a Jun-Aug mean you would have to define a new variable that averages each month and then you could average them together? I am no expert, but?I hope this helps... cheers ? web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Fri, 2/13/09, Jing-Jia Luo wrote: From: Jing-Jia Luo Subject: plot specific month or seasonal mean data from monthly time series To: GRADSUSR at LIST.CINECA.IT Date: Friday, February 13, 2009, 3:21 AM Sorry to disturb you. Using grads v1.9b4, I have monthly time series data (say from Jan. 1950 to Dec. 2008) and want to plot only May or Jun-Aug mean data each year. I am wondering whether you know how to do so and would kindly give me your advice. Thanks a lot. Best regards, Jing-Jia Luo -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090213/0b034dce/attachment.html From pdjakow at GMAIL.COM Sat Feb 14 03:11:25 2009 From: pdjakow at GMAIL.COM (PiotrD) Date: Sat, 14 Feb 2009 09:11:25 +0100 Subject: script available? In-Reply-To: <1234515596.4304.3.camel@localhost.localdomain> Message-ID: Rodney Teck wrote: > Fank > > I've attached a script which I came across sometime back it may not do > exactly as you want but see how you get on. > > Rod Teck+ > Can you tell me what variable 'z' is?: 'set lev 700' 'set ccolor 16' 'set gxout contour' 'set cstyle 1' 'set cint 20' 'set clab on' 'set clskip 3' 'd z/10' Regards Piotr -- Prognoza modelu numerycznego GFS dla obszaru Polski: http://www.gfspl.rootnode.net From pemck at GMX.NET Sat Feb 14 07:42:09 2009 From: pemck at GMX.NET (Paul Emck) Date: Sat, 14 Feb 2009 13:42:09 +0100 Subject: Findout SD Message-ID: Hello, has somebody an elegant solution (in Stergios' style) for calculating the daily ltm SD from NCEP Reanalysis data? That is from 1968 - 1996. The skript would access the daily mean variable (e.g. air.19nn.nc) which in this case are 29 files. I tried using the concatenating Grads command for all these files and then used begin time, end time and increment but for some reasons, probably size, this crashes Grads. In addition, this would not do the job very well anyhow as a constant time increment (365 days) does not respect the leap year (and the calculation of SD jan 1 would use Nov 2/1996 e.g.). A solution that accesses the nc-files seperately and uses t=1 to t=365 would be fine (ignoring leap years), in my opinion. Cheers, Paul I use v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats printim On Wed, 11 Feb 2009 16:06:14 +0100, Stergios Misios wrote: >You can do it easily with the following segment, >mean=ave(variable,t=start_time,t=end_time) >std=sqrt(ave(pow(variable-mean,2)),t=start_time,t=last_time)* (number_of_time_steps/(number_of_time_steps-1)) > >Good luck, >Stergios From arlindo.dasilva at GMAIL.COM Sat Feb 14 10:24:13 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Sat, 14 Feb 2009 10:24:13 -0500 Subject: GrADS vanishing window In-Reply-To: <6429.79222.qm@web51511.mail.re2.yahoo.com> Message-ID: On Thu, Feb 12, 2009 at 1:54 PM, Teddy Allen wrote: > Frank, > Thank you for your support with this. I am running windows XP. Everything > seems to download perfectly from the "superpack", but the window vanishes > after the landscape / portrait prompt. Unfortunately, the window vanishes > too quickly for me to read the related error message...It appears if there > is some sort of disconnect after the landscape / portrait prompt for some > mysterious reason....so, GrADS 2.0 is installed but is unable to function. > Did you follow the troubleshooting instructions in the "Getting Started" document? In partcular, make sure your virus protection software is not preventing the X server from using local ports. Usually the first time you use it the virus software asks you if this OK for Xming to use such an such port. Arlindo > Thank you again for your help, > teddy > > *web: http://personal.ecu.edu/tla0923/HOME.html* > *"To move quickly, go alone. To move far, go* *together."* > *African Proverb ** > ** > ** > > > --- On *Thu, 2/12/09, Frank Colby * wrote: > > From: Frank Colby > Subject: Re: download troubles > To: teddyallen at yahoo.com > Date: Thursday, February 12, 2009, 1:24 PM > > Terry, > > I had this problem with a computer with Vista -- had to run the installation as > administrator (right click on the .exe file, and choose the "run as > administrator" option). Are you running Vista? > > Frank > > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090214/7d66aa91/attachment.html From teddyallen at YAHOO.COM Sat Feb 14 15:04:26 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Sat, 14 Feb 2009 12:04:26 -0800 Subject: GrADS vanishing window In-Reply-To: <77fcd6b20902140724p57a899adn5cabb495dbb5d96b@mail.gmail.com> Message-ID: Dear Arlindo, After further investigation, you are correct in your diagnosis: virus software! I made a very simple adjustment in the virus software settings and now everything works perfectly. Thank you very much! teddy web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Sat, 2/14/09, Arlindo da Silva wrote: From: Arlindo da Silva Subject: Re: GrADS vanishing window To: teddyallen at yahoo.com Cc: GRADSUSR at list.cineca.it Date: Saturday, February 14, 2009, 10:24 AM On Thu, Feb 12, 2009 at 1:54 PM, Teddy Allen wrote: Frank, Thank you for your support with this. I am running windows XP. Everything seems to download perfectly from the "superpack", but the window vanishes after the landscape / portrait prompt. Unfortunately, the window vanishes too quickly for me to read the related error message...It appears if there is some sort of disconnect after the landscape / portrait prompt for some mysterious reason....so, GrADS 2.0 is installed but is unable to function. Did you follow the troubleshooting instructions in the "Getting Started" document? In partcular, make sure your virus protection software is not preventing the X server from using local ports. Usually the first time you use it the virus software asks you if this OK for Xming to use such an such port. ?? ?Arlindo ? Thank you again for your help, teddy web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Thu, 2/12/09, Frank Colby wrote: From: Frank Colby Subject: Re: download troubles To: teddyallen at yahoo.com Date: Thursday, February 12, 2009, 1:24 PM Terry, I had this problem with a computer with Vista -- had to run the installation as administrator (right click on the .exe file, and choose the "run as administrator" option). Are you running Vista? Frank -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090214/419c5c5f/attachment.html From teddyallen at YAHOO.COM Sat Feb 14 15:04:26 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Sat, 14 Feb 2009 12:04:26 -0800 Subject: GrADS vanishing window In-Reply-To: <77fcd6b20902140724p57a899adn5cabb495dbb5d96b@mail.gmail.com> Message-ID: Dear Arlindo, After further investigation, you are correct in your diagnosis: virus software! I made a very simple adjustment in the virus software settings and now everything works perfectly. Thank you very much! teddy web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? --- On Sat, 2/14/09, Arlindo da Silva wrote: From: Arlindo da Silva Subject: Re: GrADS vanishing window To: teddyallen at yahoo.com Cc: GRADSUSR at list.cineca.it Date: Saturday, February 14, 2009, 10:24 AM On Thu, Feb 12, 2009 at 1:54 PM, Teddy Allen wrote: Frank, Thank you for your support with this. I am running windows XP. Everything seems to download perfectly from the "superpack", but the window vanishes after the landscape / portrait prompt. Unfortunately, the window vanishes too quickly for me to read the related error message...It appears if there is some sort of disconnect after the landscape / portrait prompt for some mysterious reason....so, GrADS 2.0 is installed but is unable to function. Did you follow the troubleshooting instructions in the "Getting Started" document? In partcular, make sure your virus protection software is not preventing the X server from using local ports. Usually the first time you use it the virus software asks you if this OK for Xming to use such an such port From sweiss at IAFRICA.COM Sun Feb 15 12:29:52 2009 From: sweiss at IAFRICA.COM (Steven) Date: Sun, 15 Feb 2009 18:29:52 +0100 Subject: Resolution Message-ID: Hi, I'm using a 0.5x0.5 resolution for GFS data. Ideally within a grid box, there could be a few towns and cities that I would like produce data for e.g. tmp2m at a fixed lat/lon. Is there anyway to get the resolution finer - say to 0.25x0.25 in order to get more values? Regards Steven From sweiss at IAFRICA.COM Sun Feb 15 12:42:00 2009 From: sweiss at IAFRICA.COM (Steven Weiss) Date: Sun, 15 Feb 2009 18:42:00 +0100 Subject: Resolution Message-ID: Hi, List server is rejecting messages - trying again I'm using a 0.5x0.5 resolution for GFS data. Ideally within a grid box, there could be a few towns and cities that I would like produce data for e.g. tmp2m at a fixed lat/lon. Is there anyway to get the resolution finer - say to 0.25x0.25 in order to get more values? Regards Steven From arlindo.dasilva at GMAIL.COM Sun Feb 15 13:35:59 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Sun, 15 Feb 2009 13:35:59 -0500 Subject: Resolution In-Reply-To: <20090215173053.4479F1FF05@mx2.cineca.it> Message-ID: On Sun, Feb 15, 2009 at 12:29 PM, Steven wrote: > Hi, > > I'm using a 0.5x0.5 resolution for GFS data. Ideally within a grid box, > there could be a few towns and cities that I would like produce data for > e.g. tmp2m at a fixed lat/lon. Is there anyway to get the resolution finer > - > say to 0.25x0.25 in order to get more values? > A simple, brute force way is use an interpolation function such as linterp() or re() to defined your variable on a higher resoution grid, then work from that defined variable. For this to work you have to create a fake.ctl that defines this high resolution grid, e.g., dset fake.dat options template xdef ... ydef ... zdef .. tdef ... vars 1 var 0 0 dymmy endvars (the data file fake.dat does not have to exist.) Then, ga-> open model ga-> open fake ga-> set dfile 2 (now you are in the fine resolution dimension environment) ga> define ts_fine = re(ts.1,0.25) then use ts_fine in your plots. Another possibility is to use PDEF to have this interpolation done on the fly, behind the scenes. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090215/083fe992/attachment.html From ela at COLA.IGES.ORG Mon Feb 16 13:16:43 2009 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Mon, 16 Feb 2009 13:16:43 -0500 Subject: does any one know how to plot wind vector in a timeseries form In-Reply-To: <33330188.100561234808190891.JavaMail.root@mail.iges.org> Message-ID: Hello Fang, Here is a trick from the grads documentation on graphics output types (http://grads.iges.org/grads/gadoc/gradcomdsetgxout.html). I'm not sure why it's done this way, but it seems to work: To draw a 1-D time series with graphics output types vector or barb, set the dimension environment so that only time is varying, then : display const(u,0); u ; v So you should have all dimensions fixed except time, with X, Y and Z set to the values you want, set gxout to 'vector', then execute the above command (assuming your wind variables are named u and v). Best regards, Eric L. Altshuler Assistant Research Scientist Center for Ocean-Land-Atmosphere Studies 4041 Powder Mill Road, Suite 302 Calverton, MD 20705-3106 USA E-mail: ela at cola.iges.org Phone: (301) 902-1257 Fax: (301) 595-9793 ----- Original Message ----- From: "Fang-Yi Cheng" To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, February 12, 2009 10:13:39 PM GMT -05:00 US/Canada Eastern Subject: does any one know how to plot wind vector in a timeseries form Dear all, Could you please let me know how to plot wind vector in a timeseries format? Thank you. Fang From ela at COLA.IGES.ORG Mon Feb 16 13:19:47 2009 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Mon, 16 Feb 2009 13:19:47 -0500 Subject: Resolution In-Reply-To: <77fcd6b20902151035p251534e5pd31bb9f287040792@mail.gmail.com> Message-ID: I think the file 'fake.dat' referred to below must exist, but it can be empty (size 0 bytes). ----- Original Message ----- From: "Arlindo da Silva" To: GRADSUSR at LIST.CINECA.IT Sent: Sunday, February 15, 2009 1:35:59 PM GMT -05:00 US/Canada Eastern Subject: Re: Resolution On Sun, Feb 15, 2009 at 12:29 PM, Steven < sweiss at iafrica.com > wrote: Hi, I'm using a 0.5x0.5 resolution for GFS data. Ideally within a grid box, there could be a few towns and cities that I would like produce data for e.g. tmp2m at a fixed lat/lon. Is there anyway to get the resolution finer - say to 0.25x0.25 in order to get more values? A simple, brute force way is use an interpolation function such as linterp() or re() to defined your variable on a higher resoution grid, then work from that defined variable. For this to work you have to create a fake.ctl that defines this high resolution grid, e.g., dset fake.dat options template xdef ... ydef ... zdef .. tdef ... vars 1 var 0 0 dymmy endvars (the data file fake.dat does not have to exist.) Then, ga-> open model ga-> open fake ga-> set dfile 2 (now you are in the fine resolution dimension environment) ga> define ts_fine = re(ts.1,0.25) then use ts_fine in your plots. Another possibility is to use PDEF to have this interpolation done on the fly, behind the scenes. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From mequitomz at YAHOO.COM.BR Tue Feb 17 05:47:26 2009 From: mequitomz at YAHOO.COM.BR (Arlindo Meque) Date: Tue, 17 Feb 2009 02:47:26 -0800 Subject: cluster analysis Message-ID: Dear GrADS users, I am working with GPCP pentad data and would like to group the grid points into climatologically homogenous regions. Is it possible to do cluster analysis on GrADS? Regards, Arlindo Veja quais s?o os assuntos do momento no Yahoo! +Buscados http://br.maisbuscados.yahoo.com From dasilva at OPENGRADS.ORG Tue Feb 17 14:56:46 2009 From: dasilva at OPENGRADS.ORG (Arlindo da Silva) Date: Tue, 17 Feb 2009 14:56:46 -0500 Subject: Fish() requires a *global* domain Message-ID: On Mon, Feb 16, 2009 at 10:08 PM, Isaac Tetteh wrote: > Hi Arlindo, > Please I'm currently using Win 32 Superpack/Grads v1.9.0-rc 1. > I'm using global NCEP/NCAR reanalysis data, which are at multiple pressure > levels. These are available at > > http://www.cdc.noaa.gov/data/gridded/data.ncep.reanalysis.derived.pressure.html > . > I have subset the data to my domain of interest eg -50 90 lon, -40 40 lat. > Even for the moisture budget, I need to further narrow down my domain > further to quantify the moisture in the area of interest to reflect the > focus of my analysis. > If you need further clarification, I will be happy to let you know you. That is the problem: in order to compute streamfunction fish() uses a global poisson solver; it does not work with anything that is not global. Although you can implement poisson solvers for a local domain, in this case you have to specify boundary conditions --- something that is not easy to properly specify. You can find in the literature variational algorithms for computing streamfunction from (u,v) on a regional domain, but we do not have this implemented yet. If you are incline to develop this, I'm sure it would be very welcome by the community. I understand the algorithm, I just do not have the time to code it right now. For your problem, though, it may be simpler to compute chi/psi in the global domain, and then subset to your domain of interest. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090217/06900d62/attachment.html From preuter at LABRI.FR Tue Feb 17 16:11:19 2009 From: preuter at LABRI.FR (Patrick Reuter) Date: Tue, 17 Feb 2009 22:11:19 +0100 Subject: gxout vector Message-ID: Dear Grads users! I am using 'set gxout vector' to draw arrows for vector arrows, with vector arrows of varying length. I observed that the end of the vector arrows always start at the grid points, no matter their length. Do you think it is possible to attach the center (middle) of *each* vector arrow to the grid points, instead of the end of the vector arrows ? Whenever I have to code in the source code of grads, I would be willing to do it and recompile. Could somebody tell me in which file and line this would be (grads 2.0) ? I think I only have to calculate the magnitude of the vector arrow, divide it by 2, and translate it into the opposite direction of the vector, and this for each vector separately. Or are there any user-defined functions ? Hope my question is clear. Thanks for any hint/experience ! Best regards Patrick From jma at COLA.IGES.ORG Wed Feb 18 12:22:22 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Wed, 18 Feb 2009 12:22:22 -0500 Subject: debugging print statement in gribmap (GrADS version 2.0.a5) In-Reply-To: <8647ECD8-97B2-4344-93CA-9F882288C049@cola.iges.org> Message-ID: Dear All, A colleague just pointed out to me that the gribmap utitlity in 2.0.a5 is printing out a debugging message that I forgot to remove. If you see "JMA(gagmap): gfld->ipdtnum" when running gribmap, please ignore -- it doesn't affect anything that gribmap is doing. I will get rid of it in the next release. Jennifer From cdouvis at YAHOO.COM Wed Feb 18 12:45:09 2009 From: cdouvis at YAHOO.COM (Costas Douvis) Date: Wed, 18 Feb 2009 09:45:09 -0800 Subject: display station data over grided data with the same color key Message-ID: Hi everyone Here's what I'm trying to do. I have 2 data sets. One is gridded (model output) and one is station data. So in order to compare them I want to display first the gridded data and then the station data as little circles whose color is determined by the color key of the shaded data. I thought that this would be very easy to do but so far it doesn't seem to work. I might be doing something wrong. These are the commands that I used: open ModelData.ctl open StationData.ctl set gxout shaded set clevs 9 12 15 18 21 24 27 30 33 d t2m set dfile 2 set gxout stnmark set cmark 3 set clevs 9 12 15 18 21 24 27 30 33 d temp This plots the circles in the right place but they are white. And if I press "d temp" again it does the same thing only with changing color and mark (square, cross, empty circle etc.) I'm using ubuntu 8.10 and I have tried with the precompiled binaries for grads 2.0.a2, 2.0.a3 and 2.0.a5 Any ideas? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090218/1a99915a/attachment.html From elisasiq at USP.BR Wed Feb 18 12:49:11 2009 From: elisasiq at USP.BR (MARIA ELISA SIQUEIRA SILVA) Date: Wed, 18 Feb 2009 14:49:11 -0300 Subject: template nc Message-ID: Hy, Does somebody know how to use template with netcdf data? I followed the example in Grads without success. Template with binary data woks well. Thanks, Elisa ------------------------------ DSET ^pr.%x30_%x39.nc OPTIONS template DTYPE netcdf TDEF 840 linear jan1880 1mo -- Maria Elisa Siqueira Silva Laborat?rio de Climatologia e Biogeografia - LCB Depto Geografia - FFLCH - USP Av. Prof. Lineu Prestes, 338 - Cidade Universit?ria S?o Paulo - SP 05508-900 (11) 3091-3769 From pavel at CYCLONE.TAU.AC.IL Wed Feb 18 14:47:48 2009 From: pavel at CYCLONE.TAU.AC.IL (Pavel Kishcha) Date: Wed, 18 Feb 2009 21:47:48 +0200 Subject: template nc In-Reply-To: <20090218144911.m98h47zw4cs8k00s@webmail.usp.br> Message-ID: Maria Elisa, As an example, below you can see my template ctl for multiple netcdf files. The netcdf files include monthly data: 2000-03.nc, 2000-04.nc, ?, 2008-01.nc, 2008-02.nc Use the command "open" in GRADS. Pavel DSET ^%y4-%m2.nc DTYPE netcdf OPTIONS template TITLE Variables UNDEF -9.99 XDEF 360 linear -179.500000 1.00000 YDEF 181 linear -89.5000000 1.000000 ZDEF 1 levels 0 TDEF 300 linear 01mar2000 1mo vars 8 aot 0 t,y,x total aerosol optical thickness aot_land 0 t,y,x aerosol optical thickness aot_water 0 t,y,x aerosol optical thickness aot_dust 0 t,y,x atmosphere content of dust aerosol angst 0 t,y,x angstrom exponent small 0 t,y,x aerosol optical thickness due to small particles medium 0 t,y,x aerosol optical thickness due to medium_particles large 0 t,y,x aerosol optical thickness due to large particles endvars On Wed, Feb 18, 2009 at 7:49 PM, MARIA ELISA SIQUEIRA SILVA wrote: > Hy, > > Does somebody know how to use template with netcdf data? > I followed the example in Grads without success. > Template with binary data woks well. > Thanks, > Elisa > > ------------------------------ > DSET ^pr.%x30_%x39.nc > OPTIONS template > DTYPE netcdf > TDEF 840 linear jan1880 1mo > > > -- > Maria Elisa Siqueira Silva > Laborat?rio de Climatologia e Biogeografia - LCB > Depto Geografia - FFLCH - USP > Av. Prof. Lineu Prestes, 338 - Cidade Universit?ria > S?o Paulo - SP > 05508-900 > (11) 3091-3769 > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090218/d34fb42c/attachment.html From pavel at CYCLONE.TAU.AC.IL Wed Feb 18 14:47:48 2009 From: pavel at CYCLONE.TAU.AC.IL (Pavel Kishcha) Date: Wed, 18 Feb 2009 21:47:48 +0200 Subject: template nc In-Reply-To: <20090218144911.m98h47zw4cs8k00s@webmail.usp.br> Message-ID: Maria Elisa, As an example, below you can see my template ctl for multiple netcdf files From jma at COLA.IGES.ORG Thu Feb 19 09:32:23 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 19 Feb 2009 09:32:23 -0500 Subject: kml question In-Reply-To: <77fcd6b20902181902w51ae0afdr9d5daf127854e74c@mail.gmail.com> Message-ID: On Feb 18, 2009, at 10:02 PM, Arlindo da Silva wrote: > Jennifer, Brian: > > The new kml option is very cool. Do you have any plans to support > PNG with "set gxout kml"? TIFF is proprietary, and have not had a > major upgrade since 1992. Do you know how well it supports > transparency? (Very important for combining layers on > GoogleEarth.) You could probably reuse most of the printim code > for a PNG option. Just curious. > > -- > Arlindo da Silva > dasilva at alum.mit.edu Hi, Arlindo -- I have cc'd gradsusr on this question, since I think it is of general interest. The image files created with 'gxout geotiff' and 'gxout kml' are roughly equivalent to the display inside the plot area of 'gxout grfill' -- a grid of pixels without anything else. So, if you issue the following commands: 'set parea 0 11 0 8.5' 'set grid off' 'set mproj scaled' 'set x 0.5 720.5' 'set y 0.5 381.5' 'set mpdraw off' 'set grads off' 'd variable' 'printim x720 y381' Then you would create a PNG that is roughly equivalent to the TIFF and you could substitute that file name in the KML and it would work just as well. In some way, this technique is more flexible because you can use shaded contours, overlay vectors, etc. and still draw them on Google Earth. The TIFF file created with 'gxout kml' has the same geolocation metadata embedded in it the way the 'gxout geotiff' output does. The difference is that the 'gxout geotiff' files have floating-point data values for each pixel and the 'gxout kml' files have color-numbered index values. I assume there are some GIS applications that would only need the image and not the data, so that's why I set up a way to create both. The KML provides a handy way to look at the images with Google Earth. For testing, I have been using a free program called Quantum GIS (qgis.org). Transparency is pretty easy to control within this application -- you can set any number of pixel values to be transparent (on a scale of 0-100%). I think Google Earth has some controls for transparency too. So far, I haven't seen a need to set up transparency in the GrADS geotiff output. And it's definitely not clear to me that if I put a transparency mask in the TIFF file, the GIS applications will know what to do with it. But I'm still a GIS neophyte... Please note that some problems with the 'gxout geotiff' output have been reported to me. I am doing some re-coding and testing and may have some patches and/or an updated version soon. Jennifer p.s. Hooray for TIFF for keeping the standard consistent for all these years. We need more standards like that. -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090219/f6bb0113/attachment.html From brian.walsh at AMEC.COM Thu Feb 19 10:05:20 2009 From: brian.walsh at AMEC.COM (Walsh, Brian J) Date: Thu, 19 Feb 2009 10:05:20 -0500 Subject: Help with NAM model Message-ID: Hi, I have an interesting problem with the NAM model (12 km resolution). When I download the grib files and process them, the values I get for the categorical precipitation types are rarely 0 or 1 (yes or no). I'm just wondering how other people that download and extract data for the NAM model deal with this issue; I've heard of doing linear interpolations or only extracting data for a sub-region, but I'm not sure if either will work. Thanks in advance. Brian The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. Its contents (including any attachments) may contain confidential and/or privileged information. If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090219/c75f8c72/attachment.html From hrc at EAS.GATECH.EDU Thu Feb 19 16:55:01 2009 From: hrc at EAS.GATECH.EDU (Hai-Ru Chang) Date: Thu, 19 Feb 2009 16:55:01 -0500 Subject: large grib files Message-ID: Hello, When I use grads to plot a field in a large grib file (~3.1GB), to the end of 1/3 of the file the output are always wrong. I doubt is is owing to the command gribmap can't handle such a big file. The following is the part of output from the gribmap command: !!!!! MATCH: 16462 2147303268 255 1 0 135 100 175 2147172943 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16463 2147433708 255 1 0 131 100 175 2147303383 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16464 2147564148 255 1 0 132 100 175 2147433823 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16465 2147694588 255 1 0 129 100 150 -2147403033 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16466 2147825028 255 1 0 130 100 150 -2147272593 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16467 2147955468 255 1 0 133 100 150 -2147142153 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 !!!!! MATCH: 16468 2148085908 255 1 0 135 100 150 -2147011713 0 16 btim: 1992012206:00 tau: 0 dtim: 1992012206:00 As you can see there are negative numbers after 150 in several lines. Dr. Hai-Ru Chang School of Earth and Atmospheric Sciences Georgia Institute of Technology Atlanta, GA 30332-0340 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090219/351e1628/attachment.html From sigmond at ATMOSP.PHYSICS.UTORONTO.CA Thu Feb 19 18:33:27 2009 From: sigmond at ATMOSP.PHYSICS.UTORONTO.CA (Michael Sigmond) Date: Fri, 20 Feb 2009 00:33:27 +0100 Subject: xdfopen, missing longitude dimension Message-ID: Hello, I'm trying to xdfopen a netcdf file that contains zonal-mean data. The problem is that the file does not contain the longitude dimension. Details below and suggestions are very welcome, thanks, Michael --------------------------------- ncdump -h test.nc: netcdf test { dimensions: time = UNLIMITED ; // (552 currently) plev = 31 ; lat = 96 ; variables: float lat(lat) ; lat:standard_name = "latitude" ; lat:axis = "Y" ; lat:units = "degrees_north" ; lat:long_name = "Latitude" ; lat:typeConversion_op_ncl = "double converted to float" ; float plev(plev) ; plev:long_name = "Pressure" ; plev:units = "hPa" ; plev:standard_name = "air_pressure" ; plev:positive = "down" ; plev:axis = "Z" ; .... --------------------------------------------------- test.xdf: DSET CCMVal2_ref_b1_2_CAM3.5_1_1960-2005_T2Mz_ua.nc XDEF lon 1 LINEAR 1 1 YDEF lat --------------------------------------- ga-> xdfopen test.xdf Scanning Descriptor File: test.xdf Alleged lat dimension (from XDEF) is not an SDF dimension. Cannot continue. Lat dimension name is: lon . SDF file test.nc was not successfully opened & parsed. From sigmond at ATMOSP.PHYSICS.UTORONTO.CA Thu Feb 19 18:33:27 2009 From: sigmond at ATMOSP.PHYSICS.UTORONTO.CA (Michael Sigmond) Date: Fri, 20 Feb 2009 00:33:27 +0100 Subject: xdfopen, missing longitude dimension Message-ID: Hello, I'm trying to xdfopen a netcdf file that contains zonal-mean data. The problem is that the file does not contain the longitude dimension. Details below and suggestions are very welcome, thanks, Michael --------------------------------- ncdump -h test.nc: netcdf test { dimensions: time = UNLIMITED ; // (552 currently) plev = 31 ; lat = 96 ; variables: float lat(lat) ; lat:standard_name = "latitude" ; lat:axis = "Y" ; lat:units = "degrees_north" ; lat:long_name = "Latitude" ; lat:typeConversion_op_ncl = "double converted to float" ; float plev(plev) ; plev:long_name = "Pressure" ; plev:units = "hPa" ; plev:standard_name = "air_pressure" ; plev:positive = "down" ; plev:axis = "Z" ; .... --------------------------------------------------- test.xdf: DSET CCMVal2_ref_b1_2_CAM3.5_1_1960-2005_T2Mz_ua.nc XDEF lon 1 LINEAR 1 1 YDEF lat --------------------------------------- ga-> xdfopen test.xdf Scanning Descriptor File: test.xdf Alleged lat dimension (from XDEF) is not an SDF dimension. Cannot continue. Lat dimension name is: lon From teddyallen at YAHOO.COM Tue Feb 24 09:50:05 2009 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Tue, 24 Feb 2009 06:50:05 -0800 Subject: area average script question Message-ID: Dear GrADS users, I have taken the many tips from previous postings and external templates to produce a script that displays a time series of rainfall anomalies taken from the monthly GPCP dataset. The script works great as long as the geographic coordinates are relayed as a point location (ex: lat=3N, lon=6N), but, I am more interested in determining rainfall anomalies for average area regions (ex: lat=6N,lat=11N,lon=5W,lon=2E). ? I have tried to incorporate the GrADS command AAVE to calculate the area average, but?I am not getting the desired results. Please see the below script for details: ? 'reinit' *open gpcp monthly data file 'sdfopen c:\testdata\netcdf\precip.mon.mean.nc' 'set grid off' 'set grads off' *create 12 x-axis labels; one for each month 'set t 1 12' ? *define new variable to create monthly climatology for rainfall starting from Jan 1997 'define rainclim = ave(precip, t+216, t=347, 12)' 'modify rainclim seasonal' 'set t 1 last' *define new variable for anomaly of monthly rainfall 'define anom = precip - rainclim' *set time period for anomaly time series 'set t 336 347' *define geographic point location to assess rainfall anomaly for the given time period 'set lat 8' 'set lon 1' 'd anom' ? I would like to replace the set lat and set lon commands with: 'define region = aave(anom,lat=6,lat=11,lon=355,lon=2) which should average rainfall anomaly for that area. I then change 'd anom' to 'd region', but nothing happens. Can anybody spot the error? Also, is there a problem when the x dimension varies from 0 to 360 and you are interested in the region from 357 to 1 (5W to 2E), obviously the order from 357 to 1 does not obey the proper value order? Has anybody else overcome issues with the AAVE command? ? Thank you for your time. sincerely, teddy allen ? web: http://personal.ecu.edu/tla0923/HOME.html "To move quickly, go alone. To move far, go together." African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? * ? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090224/990eeb35/attachment.html From bernd.becker at METOFFICE.GOV.UK Tue Feb 24 11:37:48 2009 From: bernd.becker at METOFFICE.GOV.UK (Bernd Becker) Date: Tue, 24 Feb 2009 16:37:48 +0000 Subject: kml question In-Reply-To: <6DFF8D37-4E9F-40A4-9C8E-F4406793568E@cola.iges.org> Message-ID: This is sooo nice! After creating a nice image as described below: On Thu, 2009-02-19 at 09:32 -0500, Jennifer Adams wrote: > 'set parea 0 11 0 8.5' > 'set grid off' > 'set mproj scaled' > 'set x 0.5 720.5' > 'set y 0.5 381.5' > 'set mpdraw off' > 'set grads off' > 'd variable' > 'printim wibble.tiff gif white x720 y381' how would one make the wibble.kml file? Thanks very much. Bernd. -- Bernd Becker The Monthly Outlook Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.com - http://www.metoffice.com From haley at UCAR.EDU Tue Feb 24 11:38:59 2009 From: haley at UCAR.EDU (Mary Haley) Date: Tue, 24 Feb 2009 09:38:59 -0700 Subject: kml question In-Reply-To: <1235493468.853.58.camel@eld313.desktop.frd.metoffice.com> Message-ID: This is an automated reply. I will be out of the office through February 29 and responding to email infrequently. If this email is a question pertaining to NCL, PyNGL, PyNIO, or NCAR Graphics, please post it to the appropriate email lists. You must subscribe first: NCL: http://mailman.ucar.edu/mailman/listinfo/ncl-talk NCL install: http://mailman.ucar.edu/mailman/listinfo/ncl-install PyNGL/PyNIO: http://mailman.ucar.edu/mailman/listinfo/pyngl-talk NCARG: http://mailman.ucar.edu/mailman/listinfo/ncarg-talk Thank you, --Mary Haley From jma at COLA.IGES.ORG Tue Feb 24 11:56:29 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Tue, 24 Feb 2009 11:56:29 -0500 Subject: kml question In-Reply-To: <1235493468.853.58.camel@eld313.desktop.frd.metoffice.com> Message-ID: On Feb 24, 2009, at 11:37 AM, Bernd Becker wrote: > This is sooo nice! > > > > After creating a nice image as described below: > On Thu, 2009-02-19 at 09:32 -0500, Jennifer Adams wrote: >> 'set parea 0 11 0 8.5' >> 'set grid off' >> 'set mproj scaled' >> 'set x 0.5 720.5' >> 'set y 0.5 381.5' >> 'set mpdraw off' >> 'set grads off' >> 'd variable' >> 'printim wibble.tiff gif white x720 y381' This is not good syntax. You'll have a GIF File with a ".tiff" extension. Use this instead: 'printim wibble.gif white x720 y381' From the printim documentation: If the filename ends with ".png" or ".PNG" then GrADS will automatically create the image in PNG format If the filename ends with ".gif" or ".GIF" then GrADS will automatically create the image in GIF format If the filename ends with ".jpg" or ".JPG" then GrADS will automatically create the image in JPEG format Remember that the printim image has to have the same aspect ratio as the data grid for it to display properly with the KML file. >> > > how would one make the wibble.kml file? 'set kml wibble' 'set gxout kml' 'd variable' This gives you two files: wibble.tif and wibble.kml. You can then edit wibble.kml and change the name of the image file in the tag to whatever you created with printim. Jennifer > > > Thanks very much. > Bernd. > > -- > Bernd Becker The Monthly Outlook > Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.com - http://www.metoffice.com -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090224/440a59b6/attachment.html From arlindo.dasilva at GMAIL.COM Tue Feb 24 12:25:33 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Tue, 24 Feb 2009 12:25:33 -0500 Subject: kml question In-Reply-To: <6DFF8D37-4E9F-40A4-9C8E-F4406793568E@cola.iges.org> Message-ID: On Thu, Feb 19, 2009 at 9:32 AM, Jennifer Adams wrote: > > Hi, Arlindo -- > I have cc'd gradsusr on this question, since I think it is of general > interest. > The image files created with 'gxout geotiff' and 'gxout kml' are roughly > equivalent to the display inside the plot area of 'gxout grfill' -- a grid > of pixels without anything else. So, if you issue the following commands: > 'set parea 0 11 0 8.5' > 'set grid off' > 'set mproj scaled' > 'set x 0.5 720.5' > 'set y 0.5 381.5' > 'set mpdraw off' > 'set grads off' > 'd variable' > 'printim x720 y381' > Then you would create a PNG that is roughly equivalent to the TIFF and you > could substitute that file name in the KML and it would work just as well. Would it make sense to add a convenience -kml option to printim? printim -kml z720 y381 image.png where the KML file image.kml would be based on the current dimension environment. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From dasilva at ALUM.MIT.EDU Tue Feb 24 12:50:27 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Tue, 24 Feb 2009 12:50:27 -0500 Subject: New supplibs-2.1.0 Message-ID: All, I have posted on sf.net a new version of the baselibs: http://sourceforge.net/project/showfiles.php?group_id=161773&package_id=241681&release_id=661716 This version adds several new packages * tiff: added v3.8.3 (for GrADS v2.0.a5) * geotiff: added v1.2.5 (for GrADS v2.0.a5) * cairo: added v1.8.6 (for gxyat and future use) * pixman: added v0.13.2 (needed by cairo) * fontconfig: added v2.6.0 (needed by cairo) * freetype: added v2.3.8 (needed by cairo) * shp: added v1.2.10 added The following packages have been updated: * netcdf: updated to v4.0.1beta3 * hdf5: updated to v1.8.2 * dap: updated to v3.7.10 for compatibility with NetCDF-4 * gadap: merged with COLA v2.0 and patched while some have been removed: * nc-dap: netcdf 4.0.1 includes this functionality * netcdf4: removed (main netcdf is now has v4.0.1) * neXtaw: removed as it does not work with Ubuntu and Xaw works just fine Besides the new libraries to support the new geotiff feature in GrADS v1.0.a5 we have introduced a major upgrade in NetCDF to v4.0.1 which has OPeNDAP support built in. This will allow us to have a single GrADS executable capable of handling all the supported formats: NetCDF-3 (classic NetCDF), NetCDF-4/HDF-5 (which supports compression among other things), GRIB-1/2, OPeNDAP URLs and sequential/stream binaries. The cairo and shape libraries are not yet being used to build the core grads package, but have been included to support upcoming features. VERY IMPORTANT: The current official GrADS v2.0.a5 sources from COLA do not yet build with these supplibs. So, at this point your best choice for building COLA's v2.0.a5 release are the supplibs-2.0.1. I'll announce in a separate message an OpenGrADS patch to COLA sources that build with these libraries. Please let me know of any problems and please contact me directly if you would like to contribute binaries for a platform not yet included. (I am currently working on the Windows port.) Cheers! Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From jma at COLA.IGES.ORG Tue Feb 24 12:57:30 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Tue, 24 Feb 2009 12:57:30 -0500 Subject: kml question In-Reply-To: <77fcd6b20902240925w628db8dfjedbf4dbc93ef957e@mail.gmail.com> Message-ID: On Feb 24, 2009, at 12:25 PM, Arlindo da Silva wrote: > > Would it make sense to add a convenience -kml option to printim? No. The printim code doesn't know anything about the grid structure or the dimension environment -- it only knows about what's in the graphics metafile. The geolocaiton information is not available at that point. This technique of substituting the .tif image file created with 'gxout kml' with another image created with printim is a hack -- an unofficial alternative if you want something besides a pixelated image of your data values. By the way, I have noticed that Google Earth's smoothing algorithm can make the GrADS images look really awful -- anyone know how to turn off that feature? Jennifer From dasilva at ALUM.MIT.EDU Tue Feb 24 15:19:57 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Tue, 24 Feb 2009 15:19:57 -0500 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available Message-ID: Dear GrADS Users, I have started uploading to sf.net sources and binaries for the OpenGrADS Bundle based on COLA's 2.0.a5 release: https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032&release_id=662661 The OpenGrADS Bundle is a new way of packaging GrADS builds that is somewhat inspired by my experience with the Win32 Superpacks. The idea is to have a single package that require no (minimum?) configuration and that is easily relocatable (say to a USB memory stick). The structure of an OpenGrADS Bundle is described here: http://opengrads.org/wiki/index.php?title=The_OpenGrADS_Bundle Once you download the tarball take a look at the INSTALL file. You can install it as a Bundle (highly recommended) or simply install the binaries the same way you did in the past. Note: The Windows build will be available soon, for now only Mac OS X, Linux and FreeBSD. New Features ------------ This version is based on COLA's 2.0.a5 release which includes support for GeoTIFF and KML, as well as bug fixes. See the ChangeLog for details. In addition, several OpenGrADS specific enhancements have been introduced: - The executable *grads* is now built with NetCDF v4.0.1beta3 which includes support for NetCDF-4/HDF-5 (similar to gradsnc4 in v1.9.0-rc1) and has built in OPeNDAP support. This binary can handle all supported formats: NetCDF-3 (classic NetCDF), NetCDF-4/HDF-5, HDF-4, GRIB-1/2 and grads binaries. - The executable *gradsdap* is no longer built as its functionality is now included in the single executable *grads*. In the OpenGrADS bundle *gradsdap* simply points to *grads*. - We have introduced option -C to enable colorized text (see attached). Try "-C 0", "-C 1" and "-C 2" for different color schemes. - User defined extensions, see below. Trying it out ------------- Assuming you did a "bundle instalation" you can start the *grads* executable at the top directory (under Contents/ in the tarball): % grads Alternatively you can start the executable *opengrads* which by default comes up with colorized text and a readline history that is preserved from session to session (grads standard option -H). You can also start the executable % merra which will automatically start the Athena GUI script described in this recipe: http://cookbooks.opengrads.org/index.php?title=Recipe-016:_Accessing_MERRA_data_with_a_Graphical_User_Interface You can try one of the extensions which are documented here: http://opengrads.org/doc/#udxt In particular the spherical harmonic filter: ga-> @ open $GADSET/model ga-> d sh_filt(ps,6) The output is attached. OpenGrADS Extensions -------------------- This version includes a *preview* release of the OpenGrADS Extensions. The same extensions previously released by the OpenGrADS project with GrADS 1.9.0-rc1 are now available with GrADS v2.0, with the addition of some new ones: libmf (from Mike Fiorino) and a new spherical harmonic filter, sh_filt(). Documentation for these user defined functions/commands can be found here: http://opengrads.org/doc/#udxt -------------------------------------------------------------- IMPORTANT CAVEAT As COLA has not yet published the official API for User Defined functions in GrADS v2.0 we have adopted here an API that is based on our work with v1.9.0-rc1. This is a very low-level API that is *not* endorsed by COLA. As such, it is *not* advisable that users adopt this API to write their own extensions. Furthermore, the extensions as implemented here are tied to a specific version of GrADS and should always be built in the context a given GrADS source tree. -------------------------------------------------------------- Disclaimer: Testers needed -------------------------- These extensions are still being fully tested and should be considered *experimental* at best. Please report any problem you encounter. Use them at your own risk. List of Extensions Included --------------------------- User Defined COMMAND Short Description Function at Library ---------- ----------------------------------- -------------------------- gsudf Initialize gs-function package c_gsudf@^gsudf.gex printenv Expand environment variables c_xenv@^env.gex runenv Expand env vars and run command c_env@^env.gex @ Expand env vars and run command c_env@^env.gex getenv Get value of environment variable c_getenv@^env.gex setenv Set value of environment variable c_setenv@^env.gex gxyat Save images in PNG/SVG/PDF/PS c_gxyat@^gxyat.gex hello Hello, World! sample command c_hello@^libhello.gex ipc_verb IPC verbose toggle c_Verb@^libipc.gex ipc_open Open stream for save/load c_Open@^libipc.gex ipc_close Close stream c_Close@^libipc.gex ipc_save Save expression to stream c_Save@^libipc.gex ipc_define Define variable (obsolete) c_Define@^libipc.gex ipc_error Print IPC error message c_Error@^libipc.gex mfhilo Find max/min or H/L in 2D field c_mfhilo@^libmf.gex cylprms Properties relative to lon/lat c_cylprms@^libmf.gex shp_lines Draw lines from shapefile c_lines@^shape.gex shp_polyf Draw polygons from shapefile c_polyf@^shape.gex ---------- ----------------------------------- -------------------------- User Defined FUNCTION Short Description Function at Library ---------- ----------------------------------- -------------------------- speed Wind-speed (sample gs-function) f_gsudf@^gsudf.gex lt Less than operator f_bjt@^libbjt.gex jd Julian day f_bjt@^libbjt.gex cosz Cosine solar zenith angle f_bjt@^libbjt.gex dayratio Daylight ratio f_bjt@^libbjt.gex if Conditional function f_bjt@^libbjt.gex maxv Maximum value f_bjt@^libbjt.gex minv Minimum value f_bjt@^libbjt.gex which Label gridpoints f_bjt@^libbjt.gex ftest F-test f_bjt@^libbjt.gex ttest T-test f_bjt@^libbjt.gex tfit Point linear regression f_bjt@^libbjt.gex fit Global linear regression f_bjt@^libbjt.gex tcorr2 Time correlation f_bjt@^libbjt.gex tregr2 Point linear regression f_bjt@^libbjt.gex tmave2 Time averaging w/masking f_bjt@^libbjt.gex madvu Calculates -d(u*EXPR)/dx f_bjt@^libbjt.gex madvv Calculates -d(V*EXPR)/dy f_bjt@^libbjt.gex madvw Calculates -d(W*EXPR)/dp f_bjt@^libbjt.gex muadv Zonal advection f_bjt@^libbjt.gex mvadv Meridional advection f_bjt@^libbjt.gex mwadv Vertical advection f_bjt@^libbjt.gex satvap Saturated vapor pressure f_bjt@^libbjt.gex dew Dew point temperature f_bjt@^libbjt.gex lw Thermal infrared fluxes f_bjt@^libbjt.gex lw2 Thermal infrared fluxes v2 f_bjt@^libbjt.gex pinterp Pressure interpolation f_bjt@^libbjt.gex zinterp Height interpolation f_bjt@^libbjt.gex line Draws a line f_bjt@^libbjt.gex vint2 Mass-weighted vertical integral f_bjt@^libbjt.gex fish Poisson solver f_fish@^fish.gex fish_psi Compute streamfunction f_psichi@^fish.gex fish_chi Compute velocity potential f_psichi@^fish.gex fish_vor Compute vorticity f_psichi@^fish.gex fish_div Compute divergence f_psichi@^fish.gex hello Hello, World! sample function f_hello@^libhello.gex ipc_save Save expression to stream f_Save@^libipc.gex ipc_load Load variable from file f_Load@^libipc.gex smth2d Shuman smoother/de-smoother f_smth2d@^libmf.gex uv2trw Find radial/tangential velocity f_uv2trw@^libmf.gex re General interpolator ffre@^re.gex sh_filt Spherical harmonic filter f_shfilt@^shfilt.gex sh_power Spherical harmonic spectra f_shpowr@^shfilt.gex ---------- ----------------------------------- -------------------------- Building it ------------ To build it, get pre-compiled (or build it yourself) supplibs-2.1.0 from https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=241681&release_id=661716 Untar and setup the supplibs as explained here http://opengrads.org/wiki/index.php?title=Building_GrADS_v1.9_from_Sources Then follow the usual procedure % tar xvfz grads-x.y.x-$platform.tar.gz % cd grads-x.y.z... % ./configure % make % make check (to run the automated tests) And then to create a bundle under ./opengras type % bundle/bundle_create.sh Give it a try % cd opengrads/ and take a look ar README and INSTALL. Try this: % cd opengrads/Contents % ./merra To make a distribution tarball with the bundle % make bundle-dist Keep checking the OpenGrADS WIki for updated documentation: http://opengrads.org/wiki/index.php?title=OpenGrADS_Documentation Let me know if you have problems. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- A non-text attachment was scrubbed... Name: colorized.png Type: image/png Size: 98769 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090224/52bff64b/attachment.png -------------- next part -------------- A non-text attachment was scrubbed... Name: sh_filt.png Type: image/png Size: 130560 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090224/52bff64b/attachment-0001.png From dasilva at OPENGRADS.ORG Tue Feb 24 15:32:39 2009 From: dasilva at OPENGRADS.ORG (Arlindo da Silva) Date: Tue, 24 Feb 2009 15:32:39 -0500 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available In-Reply-To: <77fcd6b20902241219o21fea758v582e018f73d1e0b8@mail.gmail.com> Message-ID: All, If you have used the OpenGrADS extensions in v1.9.0-rc1, a few things have changed: 1) Use the following commands to see what is available: ga-> q udf (querry user defined functions) ga-> q udc (query user defined commands) ga-> q udx (query both) 2) The environment variable used to specify your UDX table is now called GA2UDXT. However, if you use a bundle installation there is *no* environment variable of any kind to set, not even LD_LIBRARY_PATH! 3) You can also load additional UDX tables at run time: ga-> load udxt /path/to/myfunctions.udxt This will append the extensions in myfunctions.udxt to the ones already there (if any). 4) The format of the UDX has changed but this is an implementation detail. 5) Needless to say, previous extensions from v1.9.0-rc1 will *not* work with v2.0.a5.oga.1 Enjoy, Arlindo --- Arlindo da Silva dasilva at opengrads.org From meteo at SPORT123.IT Wed Feb 25 05:00:20 2009 From: meteo at SPORT123.IT (Valerio Guerra) Date: Wed, 25 Feb 2009 11:00:20 +0100 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available Message-ID: Arlindo da Silva wrote: > I have started uploading to sf.net sources and binaries for the > OpenGrADS Bundle based on COLA's 2.0.a5 release: Thank you very much, Arlindo! Do you plan to release also a version for Freebsd 7.x 64bit? > The idea is to have a single package that require > no (minimum?) configuration and that is easily relocatable (say to > a USB memory stick) This is really very interesting, a good idea. Anyway I placed the Linux 32bit version on an USB key and on SDHC memory stick and both didn't worked on my Asus EeePC-Xandros (strange utime problems with gzip/tar and permission problems if I copy the 'Contents' directory from a windows machine). BTW the standard Grads2 binaries for Linux work on my EeePC. Yes, of course EeePC is not the best solution to run Grads :-) Thank you again. Valerio Guerra From cassar at PILOTBOOK.EU Wed Feb 25 05:16:57 2009 From: cassar at PILOTBOOK.EU (Pierre Cassar) Date: Wed, 25 Feb 2009 11:16:57 +0100 Subject: Recent wgrib2 for Freebsd Message-ID: Hi All, I'm looking for recent Wgrib2 binaries for Freebsd (both i386 and x64) as I need the -ctl_ens option which was added in version 1.7.7g. On the official Wgrib2 site there are no binaries for freebsd and in Opengrads site (thank you) the binaries are quite old (version 1.7.2). Any suggestion? Thank you Pierre Cassar From dasilva at ALUM.MIT.EDU Wed Feb 25 22:35:08 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Wed, 25 Feb 2009 22:35:08 -0500 Subject: Wgrib2 v1.7.8e upload to sf.net Message-ID: All, After almost a year I have updated the sources and binaries of wgrib2 on sf.net to the latest v1.7.8e: https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=263736&release_id=664109 The sources are essentially the same sources fromWes Ebisuzaki's site at NOAA/CPC: http://www.cpc.noaa.gov/products/wesley/wgrib2/ with an autoconfiscated build. You can build it with the same supplibs you used to build grads. For now, binaries for Linux i686/x86_64, Mac OS X/Intel and FreeBSD i386. I'll post Win32 builds later. Enjoy, and let me know of any problems, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From slovoed123 at YAHOO.COM Thu Feb 26 05:09:27 2009 From: slovoed123 at YAHOO.COM (slovoed slovoed) Date: Thu, 26 Feb 2009 02:09:27 -0800 Subject: possible to improve Grads for some deficiencies? Message-ID: hi, I have been using Grads for many years. As I know, there are some deficiencies in Grads I can not slove it since I have no knowledge about the source code of Grads. 1. on the last line of script, one "Return" must be added, otherwise it fail to run the last line. 2.? * is used to comment some certain lines in script, but * must be placed in the first column, otherwide will report the error. It's not useful to?let the script nice-looking and easy to read. Does it like matlab, * can be placed anywhere for commenting. The points above might be easily solved by code developer, aren't they? slovoed -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/1b30825d/attachment.html From cassar at PILOTBOOK.EU Thu Feb 26 06:12:27 2009 From: cassar at PILOTBOOK.EU (Pierre Cassar) Date: Thu, 26 Feb 2009 12:12:27 +0100 Subject: Wgrib2 v1.7.8e upload to sf.net Message-ID: Arlindo da Silva wrote: > After almost a year I have updated the sources and binaries of > wgrib2 on sf.net to the latest v1.7.8e: I tested the Freebsd 6.x version on my Freebsd 6.1 machine (i386) with no problems and the -ctl_ens option is now ok. I also tested the 32bit build on a spare Freebsd 7.0 (x64) and... it worked :-) perhaps a little slower but it could be due to the "i386 emulation" Thank you very much! Pierre Cassar From info at PAOLOCIRACI.IT Thu Feb 26 06:41:15 2009 From: info at PAOLOCIRACI.IT (Paolo Ciraci) Date: Thu, 26 Feb 2009 12:41:15 +0100 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available In-Reply-To: Message-ID: Valerio Guerra wrote: > Yes, of course EeePC is not the best solution to run Grads :-) Well, perhaps that's not the best solution but it could be very useful to run Grads on a netbook or similar small devices. Weeks ago I successfully tested (batch mode only) the grads2.a3 binaries for Linux (i686) available on Cola site on my EeePC 900 equipped with the standard Xandros OS... no apparent problems. In the weekend I'll try this new interesting build and I'll let you know the results. Best regards Dr. Paolo Ciraci info at paolociraci.it From rajkumar_62 at REDIFFMAIL.COM Thu Feb 26 06:56:08 2009 From: rajkumar_62 at REDIFFMAIL.COM (Raj Kumar Sharma) Date: Thu, 26 Feb 2009 11:56:08 -0000 Subject: error in gs file Message-ID: I have one year of data and want to plot sst average of every 5 day. I tried to write following gs file. however it shows error, some can please tell that what is the mistake. when i trie dto write aa=(ave(s,t=1,5) it shows no error.ThanksRaj'reinit''open merged.ctl''set display color white'c'set gxout shaded'it = 1while it <385'set t 'it'set lon 30 150''set lat -70 25''set grid off'it1 = itit5 = it + 4'aa=mean(s,t=it1,t=it5)''set clevs .5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6''d aa'cbarn'q pos'cit = it + 5endwhilereturn -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/5137249f/attachment.html From cassar at PILOTBOOK.EU Thu Feb 26 06:59:29 2009 From: cassar at PILOTBOOK.EU (Pierre Cassar) Date: Thu, 26 Feb 2009 12:59:29 +0100 Subject: List closed? Message-ID: Hi All, are there known problems on the list subscription/management? A friend of mine attempted several times to subscribe to GRADSUSR but only obtained similar responses from listserver: > > SUBSCRIBE GRADSUSR ******** ********* > Sorry, the GRADSUSR list is closed. Contact the list owner > (GRADSUSR-request at LIST.CINECA.IT) for more information. All information request to the address above obtained no reply. Is it true? The Gradsusr list is really "closed"? Thank you. Pierre Cassar From rajkumar_62 at REDIFFMAIL.COM Thu Feb 26 07:17:25 2009 From: rajkumar_62 at REDIFFMAIL.COM (Raj Kumar Sharma) Date: Thu, 26 Feb 2009 12:17:25 -0000 Subject: problem in script, solved Message-ID: Thanks, The problem is solved, I was writing 'aa=ave(s,t=it1,t=it5)'wheras it should have been 'aa=ave(swh,t='it1',t='it5')'Raj -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/bdb5a947/attachment.html From info at PAOLOCIRACI.IT Thu Feb 26 08:14:35 2009 From: info at PAOLOCIRACI.IT (Paolo Ciraci) Date: Thu, 26 Feb 2009 14:14:35 +0100 Subject: List closed? In-Reply-To: <20090226125929.138E0A88.cassar@pilotbook.eu> Message-ID: Pierre Cassar wrote: > are there known problems on the list subscription/management? I don't know about subscriptions but sure the listserver at Cineca has problems with Gradsusr list (I do not know about the other lists) For instance when submitting a message to the Gradsusr list I receive an error response (Subject: Rejected posting to GRADSUSR at LIST.CINECA.IT) even if the post is successful. This problem has been reported several times on the list. Hope it could be fixed. Best regards. Dr. Paolo Ciraci info at paolociraci.it From jma at COLA.IGES.ORG Thu Feb 26 08:29:03 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 26 Feb 2009 08:29:03 -0500 Subject: List closed? In-Reply-To: <20090226125929.138E0A88.cassar@pilotbook.eu> Message-ID: From http://iges.org/grads/users.html: In order to keep out spammers and internet marketers, subscription is no longer automatic. Send an email to the address below giving your affiliation, electronic and postal addresses, phone, etc. and you will be added by the system administrator. To join, send email to: gradsusr-request at list.cineca.it To unsubscribe, send an email containing the message "SIGNOFF GRADSUSR" to listserv at list.cineca.it Jennifer On Feb 26, 2009, at 6:59 AM, Pierre Cassar wrote: > Hi All, > > are there known problems on the list subscription/management? > > A friend of mine attempted several times to subscribe to GRADSUSR but > only obtained similar responses from listserver: > >>> SUBSCRIBE GRADSUSR ******** ********* >> Sorry, the GRADSUSR list is closed. Contact the list owner >> (GRADSUSR-request at LIST.CINECA.IT) for more information. > > All information request to the address above obtained no reply. > > Is it true? > The Gradsusr list is really "closed"? > > Thank you. > > Pierre Cassar -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/df703029/attachment.html From cassar at PILOTBOOK.EU Thu Feb 26 09:39:50 2009 From: cassar at PILOTBOOK.EU (Pierre Cassar) Date: Thu, 26 Feb 2009 15:39:50 +0100 Subject: List closed? Message-ID: Jennifer Adams wrote: > To join, send email to: gradsusr-request at list.cineca.it Thank you Jennifer, I'll forward your note to my friend. Anyway he already asked info to that e-mail address without any response... Pierre Cassar From scchan at ATMOS.UMD.EDU Thu Feb 26 09:57:14 2009 From: scchan at ATMOS.UMD.EDU (Steven C Chan) Date: Thu, 26 Feb 2009 09:57:14 -0500 Subject: sdfwrite and sdfattr Message-ID: Hello everyone, The version of Grads involved is 2.0a3. I am trying out the new netCDF routines for Grads, and I am having some issues with how units of time is generated when sdfwrite and inability to manually corrected the problem using sdfattr. When the first sdfwrite command is used like: sdfwrite precipitation. It generates a new netCDF file with the variable precipitation with also the time defined with it, but the units of time is always defined like "minutes since dd-mm-yyyy hr:min". It seemed Grads itself using sdfopen does not like time to be defined that way; sdfopen likes "minutes since yyyy-mm-dd hr:min" (as correctly stated in COARDS standard). When the incorrectly time unit defined .nc file is opened with sdfopen, Grads will try to read it as yyyy-mm-dd. For example: "Minutes since 1-1-1948 00:00" as generated by sdfwrite will become: 01:01:19 00:00 with sdfopen (the 48 from 1948 is dropped as there is obviously no day 1948 in any calender month. Using sdfattr trying to override the units of time does not work; it just simply will not allow me to do it. Is this a known issue with sdfwrite / sdfattr? Thanks Steven Chan From jma at COLA.IGES.ORG Thu Feb 26 10:06:35 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 26 Feb 2009 10:06:35 -0500 Subject: sdfwrite and sdfattr In-Reply-To: <49A6ADCA.5000207@atmos.umd.edu> Message-ID: That's a bug that was fixed in 2.0.a5. The correct way to write the time axis units is: "minutes since yyyy-mm-dd hh:min" --Jennifer On Feb 26, 2009, at 9:57 AM, Steven C Chan wrote: > Hello everyone, > > The version of Grads involved is 2.0a3. > > I am trying out the new netCDF routines for Grads, and I am having > some > issues with how units of time is generated when sdfwrite and inability > to manually corrected the problem using sdfattr. > > When the first sdfwrite command is used like: sdfwrite > precipitation. It > generates a new netCDF file with the variable precipitation with also > the time defined with it, but the units of time is always defined like > "minutes since dd-mm-yyyy hr:min". > It seemed Grads itself using sdfopen does not like time to be defined > that way; sdfopen likes "minutes since yyyy-mm-dd hr:min" (as > correctly > stated in COARDS standard). When the incorrectly time unit defined .nc > file is opened with sdfopen, Grads will try to read it as yyyy-mm-dd. > For example: > > "Minutes since 1-1-1948 00:00" as generated by sdfwrite will become: > 01:01:19 00:00 with sdfopen (the 48 from 1948 is dropped as there is > obviously no day 1948 in any calender month. > > Using sdfattr trying to override the units of time does not work; it > just simply will not allow me to do it. > > Is this a known issue with sdfwrite / sdfattr? > > Thanks > Steven Chan -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/f38ae397/attachment.html From scchan at ATMOS.UMD.EDU Thu Feb 26 11:39:50 2009 From: scchan at ATMOS.UMD.EDU (Steven C Chan) Date: Thu, 26 Feb 2009 11:39:50 -0500 Subject: sdfwrite and sdfattr In-Reply-To: Message-ID: Thanks. The system administrator just installed the latest version, and the problem is resolved. Steven Chan Jennifer Adams wrote: > That's a bug that was fixed in 2.0.a5. > > The correct way to write the time axis units is: > "minutes since yyyy-mm-dd hh:min" > > --Jennifer > > On Feb 26, 2009, at 9:57 AM, Steven C Chan wrote: > >> Hello everyone, >> >> The version of Grads involved is 2.0a3. >> >> I am trying out the new netCDF routines for Grads, and I am having some >> issues with how units of time is generated when sdfwrite and inability >> to manually corrected the problem using sdfattr. >> >> When the first sdfwrite command is used like: sdfwrite precipitation. It >> generates a new netCDF file with the variable precipitation with also >> the time defined with it, but the units of time is always defined like >> "minutes since dd-mm-yyyy hr:min". > > >> It seemed Grads itself using sdfopen does not like time to be defined >> that way; sdfopen likes "minutes since yyyy-mm-dd hr:min" (as correctly >> stated in COARDS standard). When the incorrectly time unit defined .nc >> file is opened with sdfopen, Grads will try to read it as yyyy-mm-dd. >> For example: >> >> "Minutes since 1-1-1948 00:00" as generated by sdfwrite will become: >> 01:01:19 00:00 with sdfopen (the 48 from 1948 is dropped as there is >> obviously no day 1948 in any calender month. >> >> Using sdfattr trying to override the units of time does not work; it >> just simply will not allow me to do it. >> >> Is this a known issue with sdfwrite / sdfattr? >> >> Thanks >> Steven Chan > > -- > Jennifer M. Adams > IGES/COLA > 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705 > jma at cola.iges.org > > > From smcmillan at PLANALYTICS.COM Thu Feb 26 14:48:10 2009 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Thu, 26 Feb 2009 13:48:10 -0600 Subject: possible to improve Grads for some deficiencies? In-Reply-To: <63070.3441.qm@web63501.mail.re1.yahoo.com> Message-ID: Slovoed, 1. Do you mean the expression "return" or an end-of-line return (as when pressing the Enter key)? If the former, "return" required at end of function, otherwise you should include a "return" (carriage return) after last line of script. However, I am referring to Windows versions--don't know about others. 2. The * (or #) is not needed in the first column for commenting. You can also place it at the end of a line, for example: prompt 'What is your name? '* Last name only pull myname prompt 'What is your age? '# Age in whole years only pull myage Stephen slovoed slovoed Sent by: GRADSUSR at LIST.CINECA.IT 02/26/2009 04:09 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject possible to improve Grads for some deficiencies? hi, I have been using Grads for many years. As I know, there are some deficiencies in Grads I can not slove it since I have no knowledge about the source code of Grads. 1. on the last line of script, one "Return" must be added, otherwise it fail to run the last line. 2. * is used to comment some certain lines in script, but * must be placed in the first column, otherwide will report the error. It's not useful to let the script nice-looking and easy to read. Does it like matlab, * can be placed anywhere for commenting. The points above might be easily solved by code developer, aren't they? slovoed *************************************************** 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. 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. *************************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/23040481/attachment.html From slovoed123 at YAHOO.COM Thu Feb 26 15:01:56 2009 From: slovoed123 at YAHOO.COM (slovoed slovoed) Date: Thu, 26 Feb 2009 12:01:56 -0800 Subject: possible to improve Grads for some deficiencies? Message-ID: Stephen, thanks for reply. 1. what I mean is pressing the Enter key after last line, otherwise Grads will not run the last line 2. You are right, but * must be placed after the last symbol without any blanks between them. I thought these 2 points can be easily solved by code developer, or I hope ... Slovoed ________________________________ From: Stephen R McMillan To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, February 26, 2009 8:48:10 PM Subject: Re: possible to improve Grads for some deficiencies? Slovoed, 1. Do you mean the expression "return" or an end-of-line return (as when pressing the Enter key)? If the former, "return" required at end of function, otherwise you should include a "return" (carriage return) after last line of script. However, I am referring to Windows versions--don't know about others. 2. The * (or #) is not needed in the first column for commenting. You can also place it at the end of a line, for example: prompt 'What is your name? '* Last name only pull myname prompt 'What is your age? '# Age in whole years only pull myage Stephen slovoed slovoed Sent by: GRADSUSR at LIST.CINECA.IT 02/26/2009 04:09 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject possible to improve Grads for some deficiencies? hi, I have been using Grads for many years. As I know, there are some deficiencies in Grads I can not slove it since I have no knowledge about the source code of Grads. 1. on the last line of script, one "Return" must be added, otherwise it fail to run the last line. 2. * is used to comment some certain lines in script, but * must be placed in the first column, otherwide will report the error. It's not useful to let the script nice-looking and easy to read. Does it like matlab, * can be placed anywhere for commenting. The points above might be easily solved by code developer, aren't they? slovoed *************************************************** 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. 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. *************************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/70e7b795/attachment.html From romero619 at HOTMAIL.COM Thu Feb 26 15:04:37 2009 From: romero619 at HOTMAIL.COM (Pablo Romero) Date: Thu, 26 Feb 2009 12:04:37 -0800 Subject: possible to improve Grads for some deficiencies? In-Reply-To: Message-ID: I took a look at the grads source code, and from what I see/understand, commenting only works if the first character on a line is either '*' or '#', so attempting to comment anywhere else on a line won't work. I also tested this and confirmed that it doesnt work. Unfortunately, this is hard-coded into the 'gs script' parsing code within the grads source code; it would need to be re-coded by COLA. I would recommend adding it as a feature request within the tracker on the GRADS sourceforge site: http://sourceforge.net/tracker/?atid=821185&group_id=161773&func=browse Pablo ________________________________ > Date: Thu, 26 Feb 2009 13:48:10 -0600 > From: smcmillan at PLANALYTICS.COM > Subject: Re: possible to improve Grads for some deficiencies? > To: GRADSUSR at LIST.CINECA.IT > > Slovoed, > 1. Do you mean the expression > "return" or an end-of-line return (as when pressing the Enter > key)? If the former, "return" required at end of function, > otherwise you should include a "return" (carriage return) after > last line of script. However, I am referring to Windows versions--don't > know about others. > 2. The * (or #) is not needed > in the first column for commenting. You can also place it at the > end of a line, for example: > > prompt 'What is your name? '* > Last name only > pull myname > prompt 'What is your age? '# Age > in whole years only > pull myage > > Stephen > > > slovoed slovoed > Sent by: GRADSUSR at LIST.CINECA.IT > > 02/26/2009 04:09 AM > Please respond to > GRADSUSR at LIST.CINECA.IT > > To > GRADSUSR at LIST.CINECA.IT > cc > Subject > possible to improve Grads for some deficiencies? > > > > > hi, > > I have been using Grads for many years. As I know, there > are some deficiencies in Grads I can not slove it since I have no knowledge > about the source code of Grads. > > 1. on the last line of script, one "Return" > must be added, otherwise it fail to run the last line. > > 2. * is used to comment some certain lines in script, > but * must be placed in the first column, otherwide will report the error. > It's not useful to let the script nice-looking and easy to read. Does > it like matlab, * can be placed anywhere for commenting. > > The points above might be easily solved by code developer, > aren't they? > > > slovoed > > > > > > > *************************************************** > > 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. 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. > > *************************************************** _________________________________________________________________ Windows Live?: Discover 10 secrets about the new Windows Live. http://windowslive.com/connect/post/jamiethomson.spaces.live.com-Blog-cns!550F681DAD532637!7540.entry?ocid=TXT_TAGLM_WL_t2_ugc_post_022009 From romero619 at HOTMAIL.COM Thu Feb 26 15:14:06 2009 From: romero619 at HOTMAIL.COM (Pablo Romero) Date: Thu, 26 Feb 2009 12:14:06 -0800 Subject: possible to improve Grads for some deficiencies? In-Reply-To: <638943.79399.qm@web63507.mail.re1.yahoo.com> Message-ID: The 'gs script' parsing code relies on the final 'return/newline' character in order to acknowledge the 'end' of a line of commands in a script file; if it doesnt find the newline character at the end of a line, it simply ignores/discards any commands that it found on that specific line. This is why the last line of a script file requires the newline character; to 'complete' the line & allow the parser to continue... This should be fix-able, but it would require the developers at COLA to make the change, since it involves making changes to the core of the GRADS source code. if you're interested, this code is found in the gscrpt.c source code file within the GRADS source; its the *gsrtyp function that contains this code. Pablo ________________________________ > Date: Thu, 26 Feb 2009 12:01:56 -0800 > From: slovoed123 at YAHOO.COM > Subject: Re: possible to improve Grads for some deficiencies? > To: GRADSUSR at LIST.CINECA.IT > > > > > Stephen, > > thanks for reply. > > 1. what I mean is pressing the Enter > key after last line, otherwise Grads will not run the last line > 2. You are right, but * must be placed after the last symbol without any blanks between them. > > I thought these 2 points can be easily solved by code developer, or I hope ... > > > Slovoed > > ________________________________ > From: Stephen R McMillan > To: GRADSUSR at LIST.CINECA.IT > Sent: Thursday, February 26, 2009 8:48:10 PM > Subject: Re: possible to improve Grads for some deficiencies? > > > Slovoed, > 1. Do you mean the expression > "return" or an end-of-line return (as when pressing the Enter > key)? If the former, "return" required at end of function, > otherwise you should include a "return" (carriage return) after > last line of script. However, I am referring to Windows versions--don't > know about others. > 2. The * (or #) is not needed > in the first column for commenting. You can also place it at the > end of a line, for example: > > prompt 'What is your name? '* > Last name only > pull myname > prompt 'What is your age? '# Age > in whole years only > pull myage > > Stephen > > > slovoed slovoed > Sent by: GRADSUSR at LIST.CINECA.IT > > 02/26/2009 04:09 AM > Please respond to > GRADSUSR at LIST.CINECA.IT > > > > > To > GRADSUSR at LIST.CINECA.IT > cc > > > Subject > possible to improve Grads for some deficiencies? > > > > > > > > > hi, > > I have been using Grads for many years. As I know, there > are some deficiencies in Grads I can not slove it since I have no knowledge > about the source code of Grads. > > 1. on the last line of script, one "Return" > must be added, otherwise it fail to run the last line. > > 2. * is used to comment some certain lines in script, > but * must be placed in the first column, otherwide will report the error. > It's not useful to let the script nice-looking and easy to read. Does > it like matlab, * can be placed anywhere for commenting. > > The points above might be easily solved by code developer, > aren't they? > > > slovoed > > > > > > > *************************************************** > > 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. 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. > > *************************************************** > _________________________________________________________________ Windows Live? Hotmail??more than just e-mail. http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t2_hm_justgotbetter_howitworks_022009 From Charles.Seman at NOAA.GOV Thu Feb 26 15:33:28 2009 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Thu, 26 Feb 2009 15:33:28 -0500 Subject: possible to improve Grads for some deficiencies? In-Reply-To: <63070.3441.qm@web63501.mail.re1.yahoo.com> Message-ID: Slovoed, If you don't mind the appearance, try this: 'your GrADS command' ;* your comment and put no space between ; and * see for example http://grads.iges.org/grads/gadoc/script.html#intro: --- Alternatively, the three records could be on the same line separated by a semicolon: |if (i = 10) ; j = 20 ; endif| --- Chuck slovoed slovoed wrote: > hi, > > I have been using Grads for many years. As I know, there are some > deficiencies in Grads I can not slove it since I have no knowledge > about the source code of Grads. > > 1. on the last line of script, one "Return" must be added, otherwise > it fail to run the last line. > > 2. * is used to comment some certain lines in script, but * must be > placed in the first column, otherwide will report the error. It's not > useful to let the script nice-looking and easy to read. Does it like > matlab, * can be placed anywhere for commenting. > > The points above might be easily solved by code developer, aren't they? > > > slovoed > -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." From smcmillan at PLANALYTICS.COM Thu Feb 26 15:34:51 2009 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Thu, 26 Feb 2009 14:34:51 -0600 Subject: possible to improve Grads for some deficiencies? In-Reply-To: Message-ID: Pablo, Slavoed, Your responses are interesting, but I routinely use comments to the right of a line, i.e., no # or * as the first character. The comment portion must be separated from (and to the right of the active portion of the line) by ';' followed by a '*' or '#' then a comment (one or more spaces may separate the '*' or '#' and the comment). However, I don't know if this is true for all versions. I am currently using version 1.9.0-rc1 on a Windows XP. I have not tried on non-Windows or later versions. Stephen Pablo Romero Sent by: GRADSUSR at LIST.CINECA.IT 02/26/2009 02:04 PM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: possible to improve Grads for some deficiencies? I took a look at the grads source code, and from what I see/understand, commenting only works if the first character on a line is either '*' or '#', so attempting to comment anywhere else on a line won't work. I also tested this and confirmed that it doesnt work. Unfortunately, this is hard-coded into the 'gs script' parsing code within the grads source code; it would need to be re-coded by COLA. I would recommend adding it as a feature request within the tracker on the GRADS sourceforge site: http://sourceforge.net/tracker/?atid=821185&group_id=161773&func=browse Pablo ________________________________ > Date: Thu, 26 Feb 2009 13:48:10 -0600 > From: smcmillan at PLANALYTICS.COM > Subject: Re: possible to improve Grads for some deficiencies? > To: GRADSUSR at LIST.CINECA.IT > > Slovoed, > 1. Do you mean the expression > "return" or an end-of-line return (as when pressing the Enter > key)? If the former, "return" required at end of function, > otherwise you should include a "return" (carriage return) after > last line of script. However, I am referring to Windows versions--don't > know about others. > 2. The * (or #) is not needed > in the first column for commenting. You can also place it at the > end of a line, for example: > > prompt 'What is your name? '* > Last name only > pull myname > prompt 'What is your age? '# Age > in whole years only > pull myage > > Stephen > > > slovoed slovoed > Sent by: GRADSUSR at LIST.CINECA.IT > > 02/26/2009 04:09 AM > Please respond to > GRADSUSR at LIST.CINECA.IT > > To > GRADSUSR at LIST.CINECA.IT > cc > Subject > possible to improve Grads for some deficiencies? > > > > > hi, > > I have been using Grads for many years. As I know, there > are some deficiencies in Grads I can not slove it since I have no knowledge > about the source code of Grads. > > 1. on the last line of script, one "Return" > must be added, otherwise it fail to run the last line. > > 2. * is used to comment some certain lines in script, > but * must be placed in the first column, otherwide will report the error. > It's not useful to let the script nice-looking and easy to read. Does > it like matlab, * can be placed anywhere for commenting. > > The points above might be easily solved by code developer, > aren't they? > > > slovoed > > > > > > > *************************************************** > > 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. 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. > > *************************************************** _________________________________________________________________ Windows Live?: Discover 10 secrets about the new Windows Live. http://windowslive.com/connect/post/jamiethomson.spaces.live.com-Blog-cns!550F681DAD532637!7540.entry?ocid=TXT_TAGLM_WL_t2_ugc_post_022009 *************************************************** 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. 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. *************************************************** -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/63e68e15/attachment.html From romero619 at HOTMAIL.COM Thu Feb 26 16:02:44 2009 From: romero619 at HOTMAIL.COM (Pablo Romero) Date: Thu, 26 Feb 2009 13:02:44 -0800 Subject: possible to improve Grads for some deficiencies? In-Reply-To: Message-ID: Stephen, your method WILL work if it includes the semicolon, as Charles had mentioned in his response(I didnt see it in your examples, Im guessing it was just a typo and you forgot to add the semicolons). The grads script parsing code treats semicolons in the same way it treats newline characters, EXCEPT for the final newline character that it expects in order to complete a line of commands in a script file; so, the ';' will still not solve the issue with GRADS ignoring the last line of a script file *if* the final newline character is missing... i.e., the following example wont work 'set lon -120 -100' 'set lat 10 30' 'd tmp' say 'this is the last line of script code';{EOF: no newline character here} the last line won't be processed, even if it ends with a semi-colon. I *believe* this might be more of a windows issue, since I believe that most *nix text editors automatically will insert the last newline character, so its a non-issue in unix/linux (however, I might be wrong about this; need confirmation). Pablo ________________________________ > Date: Thu, 26 Feb 2009 14:34:51 -0600 > From: smcmillan at PLANALYTICS.COM > Subject: Re: possible to improve Grads for some deficiencies? > To: GRADSUSR at LIST.CINECA.IT > > Pablo, Slavoed, > Your responses are interesting, but > I routinely use comments to the right of a line, i.e., no # or * as the > first character. The comment portion must be separated from (and > to the right of the active portion of the line) by ';' followed by a '*' > or '#' then a comment (one or more spaces may separate the '*' or '#' and > the comment). However, I don't know if this is true for all versions. > I am currently using version 1.9.0-rc1 on a Windows XP. I have > not tried on non-Windows or later versions. > Stephen > > > Pablo Romero > Sent by: GRADSUSR at LIST.CINECA.IT > > 02/26/2009 02:04 PM > Please respond to > GRADSUSR at LIST.CINECA.IT > > To > GRADSUSR at LIST.CINECA.IT > cc > Subject > Re: possible to improve Grads for some > deficiencies? > > > > > I took a look at the grads source code, and from what > I see/understand, commenting only works if the first character on a line > is either '*' or '#', so attempting to comment anywhere else on a line > won't work. I also tested this and confirmed that it doesnt work. > > Unfortunately, this is hard-coded into the 'gs script' > parsing code within the grads source code; it would need to be re-coded > by COLA. > > I would recommend adding it as a feature request within > the tracker on the GRADS sourceforge site: > > http://sourceforge.net/tracker/?atid=821185&group_id=161773&func=browse > > Pablo > > ________________________________ >> Date: Thu, 26 Feb 2009 13:48:10 -0600 >> From: smcmillan at PLANALYTICS.COM >> Subject: Re: possible to improve Grads for some deficiencies? >> To: GRADSUSR at LIST.CINECA.IT >> >> Slovoed, >> 1. Do you mean the expression >> "return" or an end-of-line return (as when pressing the > Enter >> key)? If the former, "return" required at end of function, >> otherwise you should include a "return" (carriage return) > after >> last line of script. However, I am referring to Windows versions--don't >> know about others. >> 2. The * (or #) is not needed >> in the first column for commenting. You can also place it at the >> end of a line, for example: >> >> prompt 'What is your name? '* >> Last name only >> pull myname >> prompt 'What is your age? '# Age >> in whole years only >> pull myage >> >> Stephen >> >> >> slovoed slovoed >> Sent by: GRADSUSR at LIST.CINECA.IT >> >> 02/26/2009 04:09 AM >> Please respond to >> GRADSUSR at LIST.CINECA.IT >> >> To >> GRADSUSR at LIST.CINECA.IT >> cc >> Subject >> possible to improve Grads for some deficiencies? >> >> >> >> >> hi, >> >> I have been using Grads for many years. As I know, there >> are some deficiencies in Grads I can not slove it since I have no > knowledge >> about the source code of Grads. >> >> 1. on the last line of script, one "Return" >> must be added, otherwise it fail to run the last line. >> >> 2. * is used to comment some certain lines in script, >> but * must be placed in the first column, otherwide will report the > error. >> It's not useful to let the script nice-looking and easy to read. Does >> it like matlab, * can be placed anywhere for commenting. >> >> The points above might be easily solved by code developer, >> aren't they? >> >> >> slovoed >> >> >> >> >> >> >> *************************************************** >> >> 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. 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. >> >> *************************************************** > _________________________________________________________________ > Windows Live?: Discover 10 secrets about the new Windows Live. > http://windowslive.com/connect/post/jamiethomson.spaces.live.com-Blog-cns!550F681DAD532637!7540.entry?ocid=TXT_TAGLM_WL_t2_ugc_post_022009 > > > > > > *************************************************** > > 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. 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. > > *************************************************** _________________________________________________________________ Windows Live? Hotmail??more than just e-mail. http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t2_hm_justgotbetter_howitworks_022009 From jma at COLA.IGES.ORG Thu Feb 26 16:22:30 2009 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 26 Feb 2009 16:22:30 -0500 Subject: possible to improve Grads for some deficiencies? In-Reply-To: <63070.3441.qm@web63501.mail.re1.yahoo.com> Message-ID: On Feb 26, 2009, at 5:09 AM, slovoed slovoed wrote: > hi, > > I have been using Grads for many years. As I know, there are some > deficiencies in Grads I can not slove it since I have no knowledge > about the source code of Grads. > > 1. on the last line of script, one "Return" must be added, otherwise > it fail to run the last line. > > 2. * is used to comment some certain lines in script, but * must be > placed in the first column, otherwide will report the error. It's > not useful to let the script nice-looking and easy to read. Does it > like matlab, * can be placed anywhere for commenting. > > The points above might be easily solved by code developer, aren't > they? Perhaps, but GrADS development resources are limited, and both problems are solved much more easily by 1) putting a carriage return at the end of your script 2) putting a ; right before the * in your comment Jennifer -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/b60784a4/attachment.html From Dan.Schwarzkopf at NOAA.GOV Thu Feb 26 17:43:38 2009 From: Dan.Schwarzkopf at NOAA.GOV (Dan Schwarzkopf) Date: Thu, 26 Feb 2009 23:43:38 +0100 Subject: running integral Message-ID: does anybody know how to do a "running integral", ie r(n) = integral (from 1 to n) of f(n)*wt(n) ? (the wt may be 1) From slovoed123 at YAHOO.COM Thu Feb 26 18:38:30 2009 From: slovoed123 at YAHOO.COM (slovoed slovoed) Date: Thu, 26 Feb 2009 15:38:30 -0800 Subject: possible to improve Grads for some deficiencies? Message-ID: You are right, that's what I can do now. I may ask for too much. I expect grads like matlab, which has no this kind of problems. slovoed ________________________________ From: Jennifer Adams To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, February 26, 2009 10:22:30 PM Subject: Re: possible to improve Grads for some deficiencies? On Feb 26, 2009, at 5:09 AM, slovoed slovoed wrote: hi, I have been using Grads for many years. As I know, there are some deficiencies in Grads I can not slove it since I have no knowledge about the source code of Grads. 1. on the last line of script, one "Return" must be added, otherwise it fail to run the last line. 2. * is used to comment some certain lines in script, but * must be placed in the first column, otherwide will report the error. It's not useful to let the script nice-looking and easy to read. Does it like matlab, * can be placed anywhere for commenting. The points above might be easily solved by code developer, aren't they? Perhaps, but GrADS development resources are limited, and both problems are solved much more easily by 1) putting a carriage return at the end of your script 2) putting a ; right before the * in your comment Jennifer -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/b6e295a1/attachment.html From nieves at LTP.UPR.CLU.EDU Thu Feb 26 20:56:00 2009 From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves) Date: Thu, 26 Feb 2009 21:56:00 -0400 Subject: freebsd-7.0 build Message-ID: Hello I have built the latest supplibs, wgrib2 and the opengrads bundle for FreeBSD-7.0 on i386 and amd64 platforms. The corresponding package files have been uploaded to the opengrads download area. In a few days I will be updating my machiens to FreeBSD-7.1. I plan to build whatever is needed again for 7.1. I will keep my 7.0 working environment for a couple of weeks in case there any problems and I need to reproduce them and rebuild the 7.0 versions. Jose From sreeds82 at YAHOO.COM Fri Feb 27 02:29:31 2009 From: sreeds82 at YAHOO.COM (sreenivas desamsetti) Date: Thu, 26 Feb 2009 23:29:31 -0800 Subject: Walker Hadley Circulatioin Message-ID: Dear Sir/ Madam, I have u, v, and omega values over the globe. Could you please tell me how to plot the Walker and Hadley circulations. If any script will be available please share with me. Thanking you DS Desamsetti Srinivas Research Scholar Department of Meteorology and Oceanography Andhra University Visakhapatnam - 530 003 91 9491686740 (M) sreeds82 at yahoo.com sreeds82 at rediffmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090226/e88cc8b1/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/jpeg Size: 11090 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090226/e88cc8b1/attachment.jpe From meteo at SPORT123.IT Fri Feb 27 05:47:08 2009 From: meteo at SPORT123.IT (Valerio Guerra) Date: Fri, 27 Feb 2009 11:47:08 +0100 Subject: freebsd-7.0 build Message-ID: Jose F. Nieves wrote: > I have built the latest supplibs, wgrib2 and the opengrads bundle for > FreeBSD-7.0 on i386 and amd64 platforms. The corresponding package > files have been uploaded to the opengrads download area. Thank you very much for the Freebsd 7 build. I placed the package in my home dir and I tried a simple test in batch mode (I do not use X) with the sample 'model.ctl' The package works but gives strange messages when launching Grads. The following is the commands I used to install the package and to run Grads: cd /usr/home/guerra gzip -d grads-2.0.a5.oga.1-bundle-amd64-unknown-freebsd7.0.tar.gz tar -xf grads-2.0.a5.oga.1-bundle-amd64-unknown-freebsd7.0.tar cd grads-2.0.a5.oga.1 cd Contents ./grads -bl And the following are the "strange" init messages: ---cut here--- Grid Analysis and Display System (GrADS) Version 2.0.a5.oga.1 Copyright (c) 1988-2008 by Brian Doty and the Institute for Global Environment and Society (IGES) GrADS comes with ABSOLUTELY NO WARRANTY See file COPYRIGHT for more information Config: v2.0.a5.oga.1 little-endian readline printim grib2 netcdf hdf4-sds opendap-grids,stn athena geotiff Issue 'q config' command for more information. Loading User Defined Extensions table ... ok. GX Package Initialization: Size = 11 8.5 Running in Batch mode "", line 407, terminal 'dtc': exit_ca_mode but no enter_ca_mode "", line 472, terminal '5410': enter_dim_mode but no exit_attribute_mode "", line 472, terminal '5410': enter_reverse_mode but no exit_attribute_mode "", line 503, terminal '5420': enter_dim_mode but no exit_attribute_mode "", line 503, terminal '5420': enter_reverse_mode but no exit_attribute_mode "", line 515, terminal '5425': enter_dim_mode but no exit_attribute_mode "", line 515, terminal '5425': enter_reverse_mode but no exit_attribute_mode "", line 525, terminal 'pc6300plus': enter_dim_mode but no exit_attribute_mode "", line 525, terminal 'pc6300plus': enter_reverse_mode but no exit_attribute_mode "", line 560, terminal 's4': enter_dim_mode but no exit_attribute_mode "", line 560, terminal 's4': enter_secure_mode but no exit_attribute_mode "", line 560, terminal 's4': enter_reverse_mode but no exit_attribute_mode "", line 563, terminal 'b610': enter_dim_mode but no exit_attribute_mode "", line 563, terminal 'b610': enter_reverse_mode but no exit_attribute_mode "", line 585, terminal 'dm2500': exit_standout_mode but no enter_standout_mode "", line 822, terminal 'scoansi': remove_clock but no display_clock "", line 1239, terminal 'dmchat': exit_standout_mode but no enter_standout_mode "", line 1267, terminal 'wsiris': enter_dim_mode but no exit_attribute_mode "", line 2156, terminal '4105': enter_ca_mode but no exit_ca_mode "", line 2156, terminal '4105': enter_insert_mode but no exit_insert_mode "", line 2181, terminal 'tek4125': enter_insert_mode but no exit_insert_mode "", line 2841, terminal 'vt520': restore_cursor but no save_cursor "", line 2861, terminal 'vt520nam': restore_cursor but no save_cursor "", line 3016, terminal 'xterm-sco': enter_alt_charset_mode but no acs_chars "", line 3016, terminal 'xterm-sco': exit_alt_charset_mode but no acs_chars "", line 4012, terminal 'wy99': from_status_line but no to_status_line "", line 4285, terminal 'ibmx': remove_clock but no display_clock "", line 4301, terminal 'ibmcx': remove_clock but no display_clock "", line 4313, terminal 'zen50': display_clock but no remove_clock "", line 4446, terminal 'fos': remove_clock but no display_clock ---cut here--- Anyway I can use Grads as usual. Perhaps those messages are debug information? Or something related to X, which I do not use? Thank you very much. Valerio Guerra From nieves at LTP.UPR.CLU.EDU Fri Feb 27 09:02:39 2009 From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves) Date: Fri, 27 Feb 2009 10:02:39 -0400 Subject: freebsd-7.0 build In-Reply-To: Message-ID: Valerio As I mentioned those messages refer to entries in the termcap file. I did what I mentioned, installed the bundle in another computer and I cannot reproduce that problem Starting "/home/users/nieves/grads-2.0.a5.oga.1/Contents/FreeBSD/amd64/grads -bl " ... Grid Analysis and Display System (GrADS) Version 2.0.a5.oga.1 Copyright (c) 1988-2008 by Brian Doty and the Institute for Global Environment and Society (IGES) GrADS comes with ABSOLUTELY NO WARRANTY See file COPYRIGHT for more information Config: v2.0.a5.oga.1 little-endian readline printim grib2 netcdf hdf4-sds opendap-grids,stn athena geotiff Issue 'q config' command for more information. Loading User Defined Extensions table ... ok. GX Package Initialization: Size = 11 8.5 Running in Batch mode ga-> I will retry after disabling X to see once more is there is any differemce. Jose From meteo at SPORT123.IT Fri Feb 27 09:49:44 2009 From: meteo at SPORT123.IT (Valerio Guerra) Date: Fri, 27 Feb 2009 15:49:44 +0100 Subject: freebsd-7.0 build Message-ID: Thank you Jose for your note > As I mentioned those messages refer to entries in the termcap > file. I did what I mentioned, installed the bundle in another > computer and I cannot reproduce that problem but unfortunately I have no idea on how to solve this problem. Should I modify termcap by hand? And if so, how? Sorry but I'm not a computer programmer, above all on Freebsd... Thank you very much. Valerio Guerra From info at PAOLOCIRACI.IT Fri Feb 27 12:13:05 2009 From: info at PAOLOCIRACI.IT (Paolo Ciraci) Date: Fri, 27 Feb 2009 18:13:05 +0100 Subject: freebsd-7.0 build Message-ID: Hi, > And the following are the "strange" init messages: we are switching our systems from Debian 4 + Freebsd 6.1 (i386) to Freebsd 7.0 (x64) so I'm coming to be interested again in Grads for Freebsd. Now Grads runs without problems on an Debian machine but in the new hardware+software configuration we want only one S.O. (Freebsd 7) and one hardware platform (Dell, Intel Xeon). So I downloaded the Freebsd 7 "bundle" package and I gave it I try... Well I had nearly the same initialization errors described by Valerio Guerra but Grads seems to work in spite of them. I run both the "Contents" version and the "Classic" even if they seem to my the same binaries with symlinks. I also gave a try to 'merra' but in this case I obtain the usual errors (line 'xxxx', terminal 'xxxx': etc...) on a "colored" display ended by an [EOF] and then grads exits to the shell prompt. Like Valerio I do not use X-windows and I run Grads in batch mode only (grads -b) with printim to generate graphical output. The machine is a Dell equipped with Intel Xeon 3.0 GHz amd 4Gb ram + Freebsd 7.0 stable. Thank you very much. Dr. Paolo Ciraci info at paolociraci.it From arlindo.dasilva at GMAIL.COM Fri Feb 27 12:42:05 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Fri, 27 Feb 2009 12:42:05 -0500 Subject: freebsd-7.0 build In-Reply-To: <200902271402.n1RE2dbQ041568@ltp.upr.clu.edu> Message-ID: On Fri, Feb 27, 2009 at 9:02 AM, Jose F. Nieves wrote: > Valerio > > As I mentioned those messages refer to entries in the termcap > file. I did what I mentioned, installed the bundle in another > computer and I cannot reproduce that problem > Jose/Valerio: Could you both send me the output of: % ldd grads We need to make sure it links statically against ncurses in the supplibs, not termcap. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From nieves at LTP.UPR.CLU.EDU Fri Feb 27 14:00:00 2009 From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves) Date: Fri, 27 Feb 2009 15:00:00 -0400 Subject: freebsd-7.0 build Message-ID: [nieves at ltp1 FreeBSD]% ldd amd64/grads amd64/grads: libXext.so.6 => /usr/local/lib/libXext.so.6 (0x800e49000) libXaw8.so.8 => /usr/local/lib/libXaw8.so.8 (0x800f59000) libXpm.so.4 => /usr/local/lib/libXpm.so.4 (0x8010c6000) libXmu.so.6 => /usr/local/lib/libXmu.so.6 (0x8011d7000) libXt.so.6 => /usr/local/lib/libXt.so.6 (0x8012ef000) libSM.so.6 => /usr/local/lib/libSM.so.6 (0x80144e000) libICE.so.6 => /usr/local/lib/libICE.so.6 (0x801556000) libX11.so.6 => /usr/local/lib/libX11.so.6 (0x801670000) libthr.so.3 => /lib/libthr.so.3 (0x801874000) libstdc++.so.6 => /usr/lib/libstdc++.so.6 (0x80198a000) libm.so.5 => /lib/libm.so.5 (0x801b96000) libgcc_s.so.1 => /lib/libgcc_s.so.1 (0x801cb0000) libc.so.7 => /lib/libc.so.7 (0x801dbd000) libXau.so.6 => /usr/local/lib/libXau.so.6 (0x801fda000) libXp.so.6 => /usr/local/lib/libXp.so.6 (0x8020dd000) libXdmcp.so.6 => /usr/local/lib/libXdmcp.so.6 (0x8021e5000) librpcsvc.so.4 => /usr/lib/librpcsvc.so.4 (0x8022ea000) Jose From meteo at SPORT123.IT Fri Feb 27 14:08:54 2009 From: meteo at SPORT123.IT (Valerio Guerra) Date: Fri, 27 Feb 2009 20:08:54 +0100 Subject: freebsd-7.0 build Message-ID: Arlindo da Silva wrote: > Jose/Valerio: Could you both send me the output of: > % ldd grads Hi Arlindo, the following is the output of ldd grads issued both in 'Contents' dir and in 'Classic' dir. Hope it could be useful. ---cut here--- $ cd /usr/home/guerra/grads-2.0.a5.oga.1 $ cd Contents $ ldd grads ldd: grads: not a dynamic executable $ cd ../Classic $ ldd bin/grads bin/grads: libXext.so.6 => /usr/local/lib/libXext.so.6 (0x800e49000) libXaw8.so.8 => /usr/local/lib/libXaw8.so.8 (0x800f59000) libXpm.so.4 => /usr/local/lib/libXpm.so.4 (0x8010c6000) libXmu.so.6 => /usr/local/lib/libXmu.so.6 (0x8011d7000) libXt.so.6 => /usr/local/lib/libXt.so.6 (0x8012ef000) libSM.so.6 => /usr/local/lib/libSM.so.6 (0x80144e000) libICE.so.6 => /usr/local/lib/libICE.so.6 (0x801556000) libX11.so.6 => /usr/local/lib/libX11.so.6 (0x801670000) libthr.so.3 => /lib/libthr.so.3 (0x801874000) libstdc++.so.6 => /usr/lib/libstdc++.so.6 (0x80198a000) libm.so.5 => /usr/lib/libm.so.5 (0x801b96000) libgcc_s.so.1 => /usr/local/lib/gcc-4.2.5/libgcc_s.so.1 (0x801cb0000) libc.so.7 => /lib/libc.so.7 (0x801dbd000) libXau.so.6 => /usr/local/lib/libXau.so.6 (0x801fe0000) libXp.so.6 => /usr/local/lib/libXp.so.6 (0x8020e3000) libXdmcp.so.6 => /usr/local/lib/libXdmcp.so.6 (0x8021eb000) librpcsvc.so.4 => /usr/lib/librpcsvc.so.4 (0x8022f0000) ---cut here--- Valerio Guerra From arlindo.dasilva at GMAIL.COM Sat Feb 28 09:20:57 2009 From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva) Date: Sat, 28 Feb 2009 09:20:57 -0500 Subject: freebsd-7.0 build In-Reply-To: Message-ID: On Fri, Feb 27, 2009 at 2:08 PM, Valerio Guerra wrote: > Arlindo da Silva wrote: >> Jose/Valerio: Could you both send me the output of: >> % ldd grads > > Hi Arlindo, > > the following is the output of ldd grads issued both in 'Contents' dir and in 'Classic' dir. > Hope it could be useful. > ... Thank you for sending me the ldd outtput. Unfortunately, it was not too conclusive. A couple of points, 1) If this is indeed a termcap/ncurses issue then it should go away if you start GrADS with command line editing off: % ./grads -E This is fine for off-line, batch work, but command-line editing is essential for productive interactive work. 2) Did you experience similar problems in Jose's previous build for FreeBSD 7 for amd64? 3) Now, if removing readline as in 1) works and the answer to 2) was no (no such problem in the previous build), then the solution would be to rebuld readline in the supplibs-2.1.0 the very same way it was built in supplibs-2.0.1: without ncurses. Jose: in the supplibs-2.1.0 build environment % gmake ncurses.purge readline.purge % make readline.install and then rebuilt grads with this revised supplbs-2.1.0. I hope it helps. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From dasilva at ALUM.MIT.EDU Sat Feb 28 10:42:00 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 28 Feb 2009 10:42:00 -0500 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available In-Reply-To: <77fcd6b20902241219o21fea758v582e018f73d1e0b8@mail.gmail.com> Message-ID: All, I have refresh the sources and binaries for the OpenGrADS Bundle on sourceforge: https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032&release_id=662661 The main changes are: 1) A bug has been fixed in the UDX interface which was causing "exec script" to fail. This is an important one, please update. 2) I have revised the internals of the Bundle implementation so that it should work now with USB flash drives. As of this writing, most binaries have been updated except for FreeBSD 7 (FreeBSD 6.3 has been updated). Make sure the files you download start with "grads-2.0.a5.oga.2" (notice the "2" at the end.) Let me know of any problems. Enjoy, Arlindo On Tue, Feb 24, 2009 at 3:19 PM, Arlindo da Silva wrote: > Dear GrADS Users, > > I have started uploading to sf.net sources and binaries for the > OpenGrADS Bundle based on COLA's 2.0.a5 release: > > https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032&release_id=662661 > > The OpenGrADS Bundle is a new way of packaging GrADS builds that > is somewhat inspired by my experience with the Win32 > Superpacks. The idea is to have a single package that require > no (minimum?) configuration and that is easily relocatable (say to > a USB memory stick). The structure of an OpenGrADS Bundle is > described here: > > http://opengrads.org/wiki/index.php?title=The_OpenGrADS_Bundle > > Once you download the tarball take a look at the INSTALL > file. You can install it as a Bundle (highly recommended) or > simply install the binaries the same way you did in the past. > > Note: The Windows build will be available soon, for now only Mac OS X, > Linux and FreeBSD. > > New Features > ------------ > > This version is based on COLA's 2.0.a5 release which includes > support for GeoTIFF and KML, as well as bug fixes. See the > ChangeLog for details. In addition, several OpenGrADS specific > enhancements have been introduced: > > - The executable *grads* is now built with NetCDF v4.0.1beta3 > which includes support for NetCDF-4/HDF-5 (similar to gradsnc4 > in v1.9.0-rc1) and has built in OPeNDAP support. This binary > can handle all supported formats: NetCDF-3 (classic NetCDF), > NetCDF-4/HDF-5, HDF-4, GRIB-1/2 and grads binaries. > > - The executable *gradsdap* is no longer built as its > functionality is now included in the single executable > *grads*. In the OpenGrADS bundle *gradsdap* simply points > to *grads*. > > - We have introduced option -C to enable colorized text (see attached). > Try "-C 0", "-C 1" and "-C 2" for different color schemes. > > - User defined extensions, see below. > > > Trying it out > ------------- > > Assuming you did a "bundle instalation" you can start the *grads* > executable at the top directory (under Contents/ in the tarball): > > % grads > > Alternatively you can start the executable *opengrads* which by > default comes up with colorized text and a readline history that > is preserved from session to session (grads standard option -H). > > You can also start the executable > > % merra > > which will automatically start the Athena GUI script described in > this recipe: > > http://cookbooks.opengrads.org/index.php?title=Recipe-016:_Accessing_MERRA_data_with_a_Graphical_User_Interface > > You can try one of the extensions which are documented here: > > http://opengrads.org/doc/#udxt > > In particular the spherical harmonic filter: > > ga-> @ open $GADSET/model > ga-> d sh_filt(ps,6) > > The output is attached. > > > OpenGrADS Extensions > -------------------- > > This version includes a *preview* release of the OpenGrADS > Extensions. The same extensions previously released by the > OpenGrADS project with GrADS 1.9.0-rc1 are now available with > GrADS v2.0, with the addition of some new ones: libmf (from Mike > Fiorino) and a new spherical harmonic filter, sh_filt(). > Documentation for these user defined functions/commands can be > found here: > > http://opengrads.org/doc/#udxt > > -------------------------------------------------------------- > IMPORTANT CAVEAT > > As COLA has not yet published the official API for User Defined > functions in GrADS v2.0 we have adopted here an API that is > based on our work with v1.9.0-rc1. This is a very low-level > API that is *not* endorsed by COLA. As such, it is *not* > advisable that users adopt this API to write their own > extensions. Furthermore, the extensions as implemented here > are tied to a specific version of GrADS and should always be > built in the context a given GrADS source tree. > -------------------------------------------------------------- > > > Disclaimer: Testers needed > -------------------------- > > These extensions are still being fully tested and should be > considered *experimental* at best. Please report any problem you > encounter. Use them at your own risk. > > List of Extensions Included > --------------------------- > > User > Defined > COMMAND Short Description Function at Library > ---------- ----------------------------------- -------------------------- > gsudf Initialize gs-function package c_gsudf@^gsudf.gex > printenv Expand environment variables c_xenv@^env.gex > runenv Expand env vars and run command c_env@^env.gex > @ Expand env vars and run command c_env@^env.gex > getenv Get value of environment variable c_getenv@^env.gex > setenv Set value of environment variable c_setenv@^env.gex > gxyat Save images in PNG/SVG/PDF/PS c_gxyat@^gxyat.gex > hello Hello, World! sample command c_hello@^libhello.gex > ipc_verb IPC verbose toggle c_Verb@^libipc.gex > ipc_open Open stream for save/load c_Open@^libipc.gex > ipc_close Close stream c_Close@^libipc.gex > ipc_save Save expression to stream c_Save@^libipc.gex > ipc_define Define variable (obsolete) c_Define@^libipc.gex > ipc_error Print IPC error message c_Error@^libipc.gex > mfhilo Find max/min or H/L in 2D field c_mfhilo@^libmf.gex > cylprms Properties relative to lon/lat c_cylprms@^libmf.gex > shp_lines Draw lines from shapefile c_lines@^shape.gex > shp_polyf Draw polygons from shapefile c_polyf@^shape.gex > ---------- ----------------------------------- -------------------------- > > User > Defined > FUNCTION Short Description Function at Library > ---------- ----------------------------------- -------------------------- > speed Wind-speed (sample gs-function) f_gsudf@^gsudf.gex > lt Less than operator f_bjt@^libbjt.gex > jd Julian day f_bjt@^libbjt.gex > cosz Cosine solar zenith angle f_bjt@^libbjt.gex > dayratio Daylight ratio f_bjt@^libbjt.gex > if Conditional function f_bjt@^libbjt.gex > maxv Maximum value f_bjt@^libbjt.gex > minv Minimum value f_bjt@^libbjt.gex > which Label gridpoints f_bjt@^libbjt.gex > ftest F-test f_bjt@^libbjt.gex > ttest T-test f_bjt@^libbjt.gex > tfit Point linear regression f_bjt@^libbjt.gex > fit Global linear regression f_bjt@^libbjt.gex > tcorr2 Time correlation f_bjt@^libbjt.gex > tregr2 Point linear regression f_bjt@^libbjt.gex > tmave2 Time averaging w/masking f_bjt@^libbjt.gex > madvu Calculates -d(u*EXPR)/dx f_bjt@^libbjt.gex > madvv Calculates -d(V*EXPR)/dy f_bjt@^libbjt.gex > madvw Calculates -d(W*EXPR)/dp f_bjt@^libbjt.gex > muadv Zonal advection f_bjt@^libbjt.gex > mvadv Meridional advection f_bjt@^libbjt.gex > mwadv Vertical advection f_bjt@^libbjt.gex > satvap Saturated vapor pressure f_bjt@^libbjt.gex > dew Dew point temperature f_bjt@^libbjt.gex > lw Thermal infrared fluxes f_bjt@^libbjt.gex > lw2 Thermal infrared fluxes v2 f_bjt@^libbjt.gex > pinterp Pressure interpolation f_bjt@^libbjt.gex > zinterp Height interpolation f_bjt@^libbjt.gex > line Draws a line f_bjt@^libbjt.gex > vint2 Mass-weighted vertical integral f_bjt@^libbjt.gex > fish Poisson solver f_fish@^fish.gex > fish_psi Compute streamfunction f_psichi@^fish.gex > fish_chi Compute velocity potential f_psichi@^fish.gex > fish_vor Compute vorticity f_psichi@^fish.gex > fish_div Compute divergence f_psichi@^fish.gex > hello Hello, World! sample function f_hello@^libhello.gex > ipc_save Save expression to stream f_Save@^libipc.gex > ipc_load Load variable from file f_Load@^libipc.gex > smth2d Shuman smoother/de-smoother f_smth2d@^libmf.gex > uv2trw Find radial/tangential velocity f_uv2trw@^libmf.gex > re General interpolator ffre@^re.gex > sh_filt Spherical harmonic filter f_shfilt@^shfilt.gex > sh_power Spherical harmonic spectra f_shpowr@^shfilt.gex > ---------- ----------------------------------- -------------------------- > > > Building it > ------------ > > To build it, get pre-compiled (or build it yourself) supplibs-2.1.0 from > > https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=241681&release_id=661716 > > Untar and setup the supplibs as explained here > > http://opengrads.org/wiki/index.php?title=Building_GrADS_v1.9_from_Sources > > Then follow the usual procedure > > % tar xvfz grads-x.y.x-$platform.tar.gz > % cd grads-x.y.z... > % ./configure > % make > % make check (to run the automated tests) > > And then to create a bundle under ./opengras type > > % bundle/bundle_create.sh > > Give it a try > > % cd opengrads/ > > and take a look ar README and INSTALL. Try this: > > % cd opengrads/Contents > % ./merra > > To make a distribution tarball with the bundle > > % make bundle-dist > > Keep checking the OpenGrADS WIki for updated documentation: > > http://opengrads.org/wiki/index.php?title=OpenGrADS_Documentation > > Let me know if you have problems. > > Arlindo > > > -- > Arlindo da Silva > dasilva at alum.mit.edu > -- Arlindo da Silva dasilva at alum.mit.edu From dasilva at ALUM.MIT.EDU Sat Feb 28 10:53:25 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 28 Feb 2009 10:53:25 -0500 Subject: OpenGrADS Bundle v2.0.a5.oga.1 available In-Reply-To: <1.3.200902262321.67546@sport123.it> Message-ID: On Thu, Feb 26, 2009 at 5:21 PM, wrote: > Hi Arlindo > >> 1) Can you do the untar on the Eee rather than on the windows machine? >> Just copy the tarball over. > > I have no problem to untar the package if I put it on the EeePC hardisk but if I put > it on the SDHC memory stick or on a USB key I can't untar it with strange utime > problems. If I try to copy the Contents directory from the hardisk to the SDHC or USB key > I obtain two errors on creating symbolic the links Linux and Resources. > No diffence if I'm user or root. > I hope to have fixed this problem in the 2.0.a5.oga.2 release I just uploaded to source forge. The main problem was that I was using traditional Unix/Linux symbolic links (symlinks) in the bundle directory tree. Even under Linux, USB flash drives tend to use a MS-DOS (FAT) filesystem which does not understand symlinks (therefore the errors when you untar'ed to the flash drive directly). If you are curious: I now implement the "symbolic link" functionality internally using simple text files ending with the character '@' (the end-user should not have to be concerned about these implementation details.) Could you please let me know if it works for you on the Eee now? Mine hasn't arrived yet, I only tested it with a 4GB flash drive on my MacBook. Thank you, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From dasilva at ALUM.MIT.EDU Sat Feb 28 11:50:58 2009 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 28 Feb 2009 11:50:58 -0500 Subject: Wgrib2 v1.7.8e upload to sf.net In-Reply-To: <77fcd6b20902251935i76a93748x95f6d20b8c0288c6@mail.gmail.com> Message-ID: A note for Windows users of wgrib2, I have uploaded 2 type of wgrib2.exe binaries for windows: 1) wgrib2-v0.1.7.8e-bin-i686-pc-cygwin.tar.gz This i what you want to use with Win32 GrADS; just drop the wgrib2.exe under the win32/ directory and you should be set to go (the dependency on cygwin1.dll will be taken care of this way). By itself, if you do not have cygwin install it will not work. 2) wgrib2-v0.1.7.8e-bin-i686-pc-win32.tar.gz This is a self-contained native win32 build which does not depend on cygwin (POSIX layer for Windows). However, the ASCII output follow the Windows/MS-DOS convention for text files and it may or may not work with Win32 GrADS --- use the cygwin build in 1) for grads work. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu