From horvath at CIRUS.DHZ.HR Thu Feb 1 03:57:11 2007 From: horvath at CIRUS.DHZ.HR (Kristian Horvath) Date: Thu, 1 Feb 2007 09:57:11 +0100 Subject: polynomial interpolation In-Reply-To: <45BE57A0.8080205@noaa.gov> Message-ID: Hello, Is there a way to draw a higher order polynomial interpolation through the desired points (other then linear with 'draw line') in GrADS or I would need to write an UDF? Thanks, Kristian _________________________________ Kristian Horvath email: horvath at cirus.dhz.hr Graduate Research Associate tel:+385-1-4565752 Meteorological and Hydrological Service www: http://www.dhmz.htnet.hr/index.php Gric 3, 10000 Zagreb, Croatia personal www: http://radar.dhz.hr/~horvath From arubawx at SETARNET.AW Thu Feb 1 11:01:30 2007 From: arubawx at SETARNET.AW (Marck Oduber) Date: Thu, 1 Feb 2007 17:01:30 +0100 Subject: wave watch data on dods stuck Message-ID: either i've missed it or not...but is there a reason why the wave watch data is stuck on jan 24 2007? on nomad5? on jan 25 2007? on nomad6? marck From John.Papineau at NOAA.GOV Thu Feb 1 12:42:03 2007 From: John.Papineau at NOAA.GOV (John Papineau) Date: Thu, 1 Feb 2007 08:42:03 -0900 Subject: time loop In-Reply-To: <7.0.1.0.0.20070201091958.01fe6c58@cirus.dhz.hr> Message-ID: Greetings Grad's World: I have 37 output files (grb, ctl, idx) from a model and I want to use *tloop* or *looping* and have no luck so far. It is unclear to me from the documentation whether I need to combine all of these files into one megafile or its it possible to use them individually? My initial *.ctl* file looks like: *tdef 1 linear 06z11oct2006 1mo* Actually, I have output every two hours instead of once per month. I did manually change the 1mo to 2hr and it screwed things up, I had to re-run post-processing to generate new ones. Any thoughts? Many Thanks- John John Papineau NWS Anchorage, Alaska "climb high, dive deep" From mulero at WISC.EDU Thu Feb 1 18:44:06 2007 From: mulero at WISC.EDU (PEDRO J MULERO) Date: Thu, 1 Feb 2007 17:44:06 -0600 Subject: Azimuthally-averaged fields Message-ID: Hello GrADS Users, Has anyone ever calculated azimuthally-averaged fields in GrADS? If so, may you share to code that does so? I'd greatly appreciate any help. Best regards, Pedro **************************************************** Pedro J. Mulero Research Assistant University of Wisconsin Space Science and Engineering Center email: mulero at wisc.edu office: AOSS 1439, (608) 262-8954 mobile: (787) 667-8653 **************************************************** From karin.meier-fleischer at ZMAW.DE Fri Feb 2 03:45:48 2007 From: karin.meier-fleischer at ZMAW.DE (Karin Meier-Fleischer) Date: Fri, 2 Feb 2007 09:45:48 +0100 Subject: leap year Message-ID: Hi, I've got a data file containing a leap year (1964). If I try to set the time to 23:59Z29Feb1964 GrADS always changes the time to 23:59Z02MAR1964. How could I prevent this? Thanks, Karin ---------------------------------\\--------------------------------------- Karin Meier-Fleischer \\ Max-Planck-Institut fuer Meteorologie Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten FAX: +49-40-41173-476 \\ Bundesstr. 53 Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg http://www.mad.zmaw.de/ \\ Germany ---------------------------------------\\--------------------------------- From smcmillan at PLANALYTICS.COM Fri Feb 2 13:01:50 2007 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Fri, 2 Feb 2007 13:01:50 -0500 Subject: Fwrite Gr2stn Output Message-ID: Hello, Is it possible to fwrite time series output, for a point using the gr2stn interpolation function? I can do it for a specified gridpoint, but not for an interpolated point using gr2stn. For example: I can do following without a problem (say, output from the GFS model): set t 1 last set lat 40 set lon -76 set gxout fwrite set fwrite d:\temp\test.dat d tmp2m-273.16 disable fwrite However, using same file and settings but for a point interpolated to -75.83 and 40.05: d gr2stn(tmp2m-273.16, -75.83, 40.05) Above will display fine on screen, but will not fwrite to specified file. Any help/suggestions would be appreciated (including "not possible"). I am aware there are other methods such as using the 'set gxout print' and write functions, for example (and have used successfully), but it would be much simpler and more efficient to output to a binary file using fwrite. Thanks-- Stephen McMillan Planalytics *************************************************** 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/20070202/e6959314/attachment.html From Charles.Seman at NOAA.GOV Fri Feb 2 16:22:44 2007 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Fri, 2 Feb 2007 16:22:44 -0500 Subject: Fwrite Gr2stn Output In-Reply-To: Message-ID: Stephen, Have you tried defining a variable, like: define tmp2mc = gr2stn(tmp2m-273.16, -75.83, 40.05) and then tried writing this to an fwrite file? Hope this helps, Chuck Stephen R McMillan wrote: > > Hello, > > Is it possible to fwrite time series output, for a point using the > gr2stn interpolation function? I can do it for a specified gridpoint, > but not for an interpolated point using gr2stn. For example: > > I can do following without a problem (say, output from the GFS model): > > set t 1 last > set lat 40 > set lon -76 > set gxout fwrite > set fwrite d:\temp\test.dat > d tmp2m-273.16 > disable fwrite > > However, using same file and settings but for a point interpolated to > -75.83 and 40.05: > > d gr2stn(tmp2m-273.16, -75.83, 40.05) > > Above will display fine on screen, but will not fwrite to specified > file. Any help/suggestions would be appreciated (including "not > possible"). I am aware there are other methods such as using the 'set > gxout print' and write functions, for example (and have used > successfully), but it would be much simpler and more efficient to > output to a binary file using fwrite. > > Thanks-- > Stephen McMillan > Planalytics > > *************************************************** > 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. > *************************************************** -- 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 reflect any position of the Government or NOAA." From jma at COLA.IGES.ORG Fri Feb 2 17:33:51 2007 From: jma at COLA.IGES.ORG (Jennifer M. Adams) Date: Fri, 2 Feb 2007 17:33:51 -0500 Subject: Fwrite Gr2stn Output In-Reply-To: Message-ID: The 'set gxout fwrite' option doesn't work with station expressions, which is what gr2stn returns. You need to wrap your expression inside a s2g1d ("station to grid 1D") call, to change your time series at an interpolated location back into a grid type. Then fwrite will work. Jennifer On Feb 2, 2007, at 1:01 PM, Stephen R McMillan wrote: > > Hello, > > Is it possible to fwrite time series output, for a point using the > gr2stn interpolation function? I can do it for a specified > gridpoint, but not for an interpolated point using gr2stn. For > example: > > I can do following without a problem (say, output from the GFS model): > > set t 1 last > set lat 40 > set lon -76 > set gxout fwrite > set fwrite d:\temp\test.dat > d tmp2m-273.16 > disable fwrite > > However, using same file and settings but for a point interpolated > to -75.83 and 40.05: > > d gr2stn(tmp2m-273.16, -75.83, 40.05) > > Above will display fine on screen, but will not fwrite to specified > file. Any help/suggestions would be appreciated (including "not > possible"). I am aware there are other methods such as using the > 'set gxout print' and write functions, for example (and have used > successfully), but it would be much simpler and more efficient to > output to a binary file using fwrite. > > Thanks-- > Stephen McMillan > Planalytics > > *************************************************** > 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. > *************************************************** Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Beltsville, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070202/9ca306dc/attachment.html From jweather at SDSC.EDU Fri Feb 2 17:44:35 2007 From: jweather at SDSC.EDU (John Weatherford) Date: Fri, 2 Feb 2007 14:44:35 -0800 Subject: Library Issues Message-ID: Hello all, I am working to get a full installation of grads-1.9b4 up and running including all the goodies (netCDF, OPeNDAP, etc.). I believe the best way to get all this is to get the Supplibs package. I am running on a solaris machine and had to complie the source to get a basic installation running. I have tried to download the binaries for the Supplibs (ftp://grads.iges.org/grads/Supplibs/1.9) but the solaris binaries will not uncompress, and the source files have an error in the makefile (Unexpected EOF). The only binaries I found that seemed to work is under the "old" directory. However, as I should have expected, these give a compile error due to a different function definition. So my question is, does anyone know of an alternate place to download the supplibs, and/or the source that I might try out. Thanks for any help! I really appreciate it. -John Weatherford San Diego Supercomputer Center From herb.philbrick at MITTALSTEEL.COM Fri Feb 2 20:23:14 2007 From: herb.philbrick at MITTALSTEEL.COM (Herb Philbrick) Date: Fri, 2 Feb 2007 19:23:14 -0600 Subject: Regrid2 error on XP - array bounds exceeded Message-ID: Hello all, I'm trying to use Regrid2 on a Windows XP machine (with SP2). Grads writes the input file to the udf and my compiled Regrid2.exe reads the contents fine. Regrid2 encounters an "array bounds exceeded" error while in the "BOX_AVE" subroutine. The exact line is area_box(icnt)=dx*dy*area_in(ii0,jj). The problem is that ii0 exceeds the dimension of the area_in array by 1, i.e. ii0=46 (1st subscript) and dim area_in(45,47). It seems like the condition is invited when a related upper limit (niip1) is set to 1 + array size earlier in the program. The source code of Regrid2 is dated 11/30/2001. My compiler is Digital Fortran Version 6.0a I sure would appreciate the some guidance or tips. I'd be much obliged Thank you, Herb Philbrick herb.philbrick at mittalsteel.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070202/70af0710/attachment.html From smcmillan at PLANALYTICS.COM Fri Feb 2 22:02:47 2007 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Fri, 2 Feb 2007 22:02:47 -0500 Subject: Fwrite Gr2stn Output Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070202/bd3059e8/attachment.html From smcmillan at PLANALYTICS.COM Fri Feb 2 22:25:46 2007 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Fri, 2 Feb 2007 22:25:46 -0500 Subject: Fwrite Gr2stn Output Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070202/0b94f9ed/attachment.html From arabonis at EGS.UCT.AC.ZA Mon Feb 5 01:10:41 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Mon, 5 Feb 2007 08:10:41 +0200 Subject: GFS files corrupted on dod server Message-ID: Hi Jennifer Seems the first few days of the gfs files on monsoondata are corrupted. Thanks Jean Pierre -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free From jma at COLA.IGES.ORG Mon Feb 5 07:13:18 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Mon, 5 Feb 2007 07:13:18 -0500 Subject: GFS files corrupted on dod server In-Reply-To: <45C6CA61.5070109@egs.uct.ac.za> Message-ID: Jean Pierre -- You'll have to be a bit more specific. I don't see anything wrong with the 00z forecast for today, and everything looks normal. Which files are you having trouble with? On Feb 5, 2007, at 1:10 AM, Jean Pierre Arabonis wrote: > Hi Jennifer > Seems the first few days of the gfs files on monsoondata are > corrupted. > Thanks > Jean Pierre > > > -- > Jean Pierre Arabonis > arabonis at egs.uct.ac.za > Tel 021 780 1021 cell 084 401 1365 > This email is 100% Microsoft Free > -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Beltsville, MD 20705 jma at cola.iges.org -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070205/fcfad182/attachment.html From f.piani at IBIMET.CNR.IT Mon Feb 5 08:53:51 2007 From: f.piani at IBIMET.CNR.IT (Francesco Piani) Date: Mon, 5 Feb 2007 14:53:51 +0100 Subject: Eliassen-Palm Flux Message-ID: Hi GradsUsers. Does anyone have a script to calculate Eliassen-Palm Flux? Thanks in advance, Francesco -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070205/af95ecbb/attachment.html From pbehling at WISC.EDU Mon Feb 5 14:09:45 2007 From: pbehling at WISC.EDU (Pat Behling) Date: Mon, 5 Feb 2007 13:09:45 -0600 Subject: websites that use grads Message-ID: http://ccr.aos.wisc.edu/model/grads/grads1.html Pat Behling Senior Information Processing Consultant Center for Climatic Research Gaylord Nelson Institute for Environmental Studies University of Wisconsin 1225 W. Dayton St. Madison, WI 53706-1695 phone: 608-262-8765 fax: 608-263-4190 email: pbehling at wisc.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070205/a97333a5/attachment.html From munnich at ATMOS.UCLA.EDU Mon Feb 5 16:48:00 2007 From: munnich at ATMOS.UCLA.EDU (Matthias Munnich) Date: Mon, 5 Feb 2007 13:48:00 -0800 Subject: leap year In-Reply-To: Message-ID: -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Karin, without a look at the descriptor file and a complete log of an example session it's hard to figure out what exactly is going on. Can you use "set t" to set time to 29Feb? ... Matt Karin Meier-Fleischer wrote: > Hi, > > I've got a data file containing a leap year (1964). If I try to set the > time to 23:59Z29Feb1964 GrADS always > changes the time to 23:59Z02MAR1964. > > How could I prevent this? > > Thanks, Karin > ---------------------------------\\--------------------------------------- > > Karin Meier-Fleischer \\ Max-Planck-Institut fuer > Meteorologie > Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten > FAX: +49-40-41173-476 \\ Bundesstr. 53 > Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg > http://www.mad.zmaw.de/ \\ Germany > ---------------------------------------\\--------------------------------- > - -- - -------------------------------- Matthias Munnich Univ. of California, Los Angeles Dept. of Atmos. and Oceanic Sc. and Inst. of Geophysics and Planetary Physics 3845 Slichter Hall Los Angeles, CA 90095-1567 Phone: +1-310-794 5899 Fax: +1-310-206 3051 Email: munnich at atmos.ucla.edu - -------------------------------- -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFFx6YPidaDTh/cfhYRAqmdAKCDrKfeM5v5YBdiS1DaNnQg2ceS4ACgtMgV ZDPrtiY2YUNsOTSCIhYh1nE= =uP3h -----END PGP SIGNATURE----- From karin.meier-fleischer at ZMAW.DE Tue Feb 6 04:31:38 2007 From: karin.meier-fleischer at ZMAW.DE (Karin Meier-Fleischer) Date: Tue, 6 Feb 2007 10:31:38 +0100 Subject: leap year In-Reply-To: <45C7A610.5030306@atmos.ucla.edu> Message-ID: Hallo Matthias, schoen von Dir zu hoeren. Here is my problem with the leap year: the data file contains 12 monthly time steps (from 23:59Z31JAN1964 to 23:59Z31DEC1964) -1 : Date Time Code Level Size Miss : Minimum Mean Maximum 1 : 1964-01-31 23:59 -1 0 42126 0 : 0.0000 25.063 557.66 2 : 1964-02-29 23:59 -1 0 42126 0 : 0.0000 25.040 1086.3 3 : 1964-03-31 23:59 -1 0 42126 0 : 0.0000 20.098 539.31 4 : 1964-04-30 23:59 -1 0 42126 0 : 0.0000 20.345 265.12 5 : 1964-05-31 23:59 -1 0 42126 0 : 0.0000 22.864 315.59 6 : 1964-06-30 23:59 -1 0 42126 0 : 0.0000 29.650 388.09 7 : 1964-07-31 23:59 -1 0 42126 0 : 0.0000 27.702 287.31 8 : 1964-08-31 23:59 -1 0 42126 0 : 0.0000 26.029 243.57 9 : 1964-09-30 23:59 -1 0 42126 0 : 0.0000 22.945 355.16 10 : 1964-10-31 23:59 -1 0 42126 0 : 0.0000 24.714 359.15 11 : 1964-11-30 23:59 -1 0 42126 0 : 0.0000 21.487 424.59 12 : 1964-12-31 23:59 -1 0 42126 0 : 0.0000 17.004 294.62 cdo info: Processed 1 variable 12 timesteps. ( 0.02s ) The GrADS script leap_year.gs has the following output: ------------------------------------------------------- Scanning self-describing file: /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc Found 1 data variables in SDF file. Found displayable variable rain_con with 0 levels in SDF file. dset /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc title CLM simulation of the 20th century climate with anthropogenic forcing (C20_1) for the SGA project (CLM_C20_1) undef 9.96921e+36 xdef 238 linear -10.6 0.200001 ydef 177 linear 34.6 0.200001 zdef 1 linear 0 0 tdef 12 linear 23:59Z31JAN1964 1mo vars 1 rain_con 0 -999 convective rainfall -- t=1 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is fixed Time = 23:59Z31JAN1964 T = 1 -- t=2 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is fixed Time = 23:59Z02MAR1964 T = 3 -- t=3 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is fixed Time = 23:59Z31MAR1964 T = 3 -- t=12 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is fixed Time = 23:59Z31DEC1964 T = 12 -- time=23:59Z29FEB1964 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is fixed Time = 23:59Z02MAR1964 T = 3 -- t=1-12 -- Default file number is: 1 X is varying Lon = -10.6 to 36.8002 X = 1 to 238 Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 Z is fixed Lev = 0 Z = 1 T is varying Time = 23:59Z31JAN1964 to 23:59Z31DEC1964 T = 1 to 12 -------------------------------------------------------------------------- The plot shows a wrong value for february (it ist the same as march): instead of value 1086.3 it plots 539.31. Any idea what is going wrong? Thanks in advance, Karin ---------------------------------\\--------------------------------------- Karin Meier-Fleischer \\ Max-Planck-Institut fuer Meteorologie Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten FAX: +49-40-41173-476 \\ Bundesstr. 53 Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg http://www.mad.zmaw.de/ \\ Germany ---------------------------------------\\--------------------------------- On Mon, 5 Feb 2007, Matthias Munnich wrote: > -----BEGIN PGP SIGNED MESSAGE----- > Hash: SHA1 > > Hi Karin, > > without a look at the descriptor file and a complete log of an example > session > it's hard to figure out what exactly is going on. Can you use "set t" > to set time > to 29Feb? > > ... Matt > > > Karin Meier-Fleischer wrote: >> Hi, >> >> I've got a data file containing a leap year (1964). If I try to set the >> time to 23:59Z29Feb1964 GrADS always >> changes the time to 23:59Z02MAR1964. >> >> How could I prevent this? >> >> Thanks, Karin >> ---------------------------------\\--------------------------------------- >> >> Karin Meier-Fleischer \\ Max-Planck-Institut fuer >> Meteorologie >> Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten >> FAX: +49-40-41173-476 \\ Bundesstr. 53 >> Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg >> http://www.mad.zmaw.de/ \\ Germany >> ---------------------------------------\\--------------------------------- >> > > > - -- > - -------------------------------- > Matthias Munnich > Univ. of California, Los Angeles > Dept. of Atmos. and Oceanic Sc. and > Inst. of Geophysics and Planetary Physics > 3845 Slichter Hall > Los Angeles, CA 90095-1567 > Phone: +1-310-794 5899 > Fax: +1-310-206 3051 > Email: munnich at atmos.ucla.edu > - -------------------------------- > -----BEGIN PGP SIGNATURE----- > Version: GnuPG v1.4.2.2 (GNU/Linux) > Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org > > iD8DBQFFx6YPidaDTh/cfhYRAqmdAKCDrKfeM5v5YBdiS1DaNnQg2ceS4ACgtMgV > ZDPrtiY2YUNsOTSCIhYh1nE= > =uP3h > -----END PGP SIGNATURE----- > -------------- next part -------------- function LeapYear(args) _DATA='/pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc' 'sdfopen '_DATA 'q ctlinfo'; say result say '-- t=1 --' 'set t 1' 'q dims'; say result say '-- t=2 --' 'set t 2' 'q dims'; say result say '-- t=3 --' 'set t 3' 'q dims'; say result say '-- t=12 --' 'set t 12' 'q dims'; say result say '-- time=23:59Z29FEB1964 --' 'set time 23:59Z29FEB1964' 'q dims'; say result say '-- t=1-12 --' 'set t 1 12' 'q dims'; say result 'set x 1' 'set y 1' 'd tloop(rain_con)' 'printim leap_year.gif gif white' return -------------- next part -------------- A non-text attachment was scrubbed... Name: leap_year.gif Type: image/gif Size: 32041 bytes Desc: GIF plot output file of leap_year.gs Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070206/1e4c9039/attachment.gif From f.piani at IBIMET.CNR.IT Tue Feb 6 10:29:14 2007 From: f.piani at IBIMET.CNR.IT (Francesco Piani) Date: Tue, 6 Feb 2007 16:29:14 +0100 Subject: Eddy Heat Flux Message-ID: Hi all, Grads Users. Does anyone have a script to calculate Meridional Eddy Heat Flux? Thanks in advance, Francesco -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/8bf4cd00/attachment.html From mcroke at AIRDAT.COM Tue Feb 6 11:35:02 2007 From: mcroke at AIRDAT.COM (Meredith Croke) Date: Tue, 6 Feb 2007 11:35:02 -0500 Subject: Using fwrite and controlling decimal places Message-ID: Hi Everyone, I am new to GrADS and am trying to write a file that has 10 columns of temperature readings for over 600 times. The problem I am running into is controlling how far the number extends beyond the decimal point. I would like there to only be two digits. I tried using 'set dignum 2', but this didn't work, and from what I can tell it is mainly used for plotting. My write line is as follows: f=write (test.dat, ds' 'hs' 'df' 'hf' 0'fff' 'TF' 'TF1' 'TF2' 'TF3' 'TF4' 'TF5' 'TF6' 'TF7' 'TF8' 'TF9,append) Where all the TF's are the temperatures that I would like to only have 2 numbers after the decimal point. Is there some format statement I can put in that would fix this problem? Thanks for your help! Sincerely, Meredith -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/a0716478/attachment.html From smcmillan at PLANALYTICS.COM Tue Feb 6 12:03:44 2007 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Tue, 6 Feb 2007 12:03:44 -0500 Subject: Using fwrite and controlling decimal places In-Reply-To: <090531DCD6745B4E8B49300999500668F24D70@rdu-w2k3ex01.airdat.com> Message-ID: Try this Meredith: tmp=math_format('%0.2f',tmp) where tmp is your temperature variable. Good luck-- Stephen McMillan Meredith Croke Sent by: GRADSUSR at LIST.CINECA.IT 02/06/2007 11:35 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Using fwrite and controlling decimal places Hi Everyone, I am new to GrADS and am trying to write a file that has 10 columns of temperature readings for over 600 times. The problem I am running into is controlling how far the number extends beyond the decimal point. I would like there to only be two digits. I tried using 'set dignum 2', but this didn't work, and from what I can tell it is mainly used for plotting. My write line is as follows: f=write (test.dat, ds' 'hs' 'df' 'hf' 0'fff' 'TF' 'TF1' 'TF2' 'TF3' 'TF4' 'TF5' 'TF6' 'TF7' 'TF8' 'TF9,append) Where all the TF's are the temperatures that I would like to only have 2 numbers after the decimal point. Is there some format statement I can put in that would fix this problem? Thanks for your help! Sincerely, Meredith *************************************************** 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/20070206/8b0e55f4/attachment.html From mcroke at AIRDAT.COM Tue Feb 6 14:02:02 2007 From: mcroke at AIRDAT.COM (Meredith Croke) Date: Tue, 6 Feb 2007 20:02:02 +0100 Subject: fwrite, decimal places Message-ID: Hi, I am trying to write a grads script that outputs a .dat file, the problem I am running into is I can't seem to figure out how to tell grads how many decimal places I want to use. I tried 'set dignum 2' but this didn't do anything. My write line is as follows: f=write(test.dat,ds' 'hs' 'df' 'hf' 0'fff' 'TF' 'TF1' 'TF2' 'TF3' 'TF4' 'TF5' 'TF6' 'TF7' 'TF8' 'TF9,append) Is there anyway I can set the T's to REAL and declare how many numbers beyond the decimal point I want? Thanks, Meredith From jsg229 at PSU.EDU Tue Feb 6 11:31:23 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Tue, 6 Feb 2007 11:31:23 -0500 Subject: grads question Message-ID: I am attempting to create a .dat file and .ctl file for use in creating radius-height (r-z) plots in Grads. I have 4 different times and the dimensions of the data are 32 points in the radial direction and 30 points in the vertical direction. Suppose that I want to create a r-z plot of temperature. This is what I have done: DO T=1,4 !T=time counter WRITE(541,REC=C) (R(I,T),I=1,IR) !R=radius, IR=32 WRITE(541,REC=C) (Z(K,T),K=1,KX) !Z=height, KX=30 WRITE(541,REC=C) ((T_GRADS(K,I,C),K=1,KX),I=1,IR) !T_GRADS = temperature END DO Now, my .ctl file looks like:DSET 0.3sf_grads.dat byteswapped TITLE Sample UNDEF -9999 XDEF 30 linear 1 1 YDEF 32 linear 1 1 ZDEF 1 linear 1 1 TDEF 4 linear 12:00Z15SEP1999 180MN vars 3 z 0 99 height r 0 99 radius t_az 0 99 TEMP AZ. AVG. endvars When I attempt to create my r-z plot, I get a field is constant, value=0 error. What am I doing wrong here. Any help/advice in fixing this problem would be most appreciated. Thanks, Jeff -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/6ab55c53/attachment.html From kieucq at YAHOO.CO.UK Tue Feb 6 14:48:38 2007 From: kieucq at YAHOO.CO.UK (Lemmy) Date: Tue, 6 Feb 2007 19:48:38 +0000 Subject: grads question In-Reply-To: <1170779482l.1114114l.0l@psu.edu> Message-ID: The record length of R, Z (1D) and T_grads (2D) are not consistent JEFFREY S GALL wrote: I am attempting to create a .dat file and .ctl file for use in creating radius-height (r-z) plots in Grads. I have 4 different times and the dimensions of the data are 32 points in the radial direction and 30 points in the vertical direction. Suppose that I want to create a r-z plot of temperature. This is what I have done: DO T=1,4 !T=time counter WRITE(541,REC=C) (R(I,T),I=1,IR) !R=radius, IR=32 WRITE(541,REC=C) (Z(K,T),K=1,KX) !Z=height, KX=30 WRITE(541,REC=C) ((T_GRADS(K,I,C),K=1,KX),I=1,IR) !T_GRADS = temperature END DO Now, my .ctl file looks like: DSET 0.3sf_grads.dat byteswapped TITLE Sample UNDEF -9999 XDEF 30 linear 1 1 YDEF 32 linear 1 1 ZDEF 1 linear 1 1 TDEF 4 linear 12:00Z15SEP1999 180MN vars 3 z 0 99 height r 0 99 radius t_az 0 99 TEMP AZ. AVG. endvars When I attempt to create my r-z plot, I get a field is constant, value=0 error. What am I doing wrong here. Any help/advice in fixing this problem would be most appreciated. Thanks, Jeff --------------------------------- Now you can scan emails quickly with a reading pane. Get the new Yahoo! Mail. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070206/e16ff26e/attachment.html From bnarapus at GMU.EDU Tue Feb 6 14:09:20 2007 From: bnarapus at GMU.EDU (Balachandrudu Narapusetty) Date: Tue, 6 Feb 2007 14:09:20 -0500 Subject: fwrite, decimal places In-Reply-To: <20070206190245.EA16A227AA@mx2.cineca.it> Message-ID: Hi, You might want to use "prnopts" command. You can find a quick usage of that command at http://grads.iges.org/grads/gadoc/gadocindex.html hope that helps, bala ----- Original Message ----- From: Meredith Croke Date: Tuesday, February 6, 2007 2:02 pm Subject: fwrite, decimal places > Hi, > > I am trying to write a grads script that outputs a .dat file, the > problem I > am running into is I can't seem to figure out how to tell grads how > manydecimal places I want to use. I tried 'set dignum 2' but this > didn't do > anything. My write line is as follows: > f=write(test.dat,ds' 'hs' 'df' 'hf' 0'fff' 'TF' 'TF1' 'TF2' 'TF3' > 'TF4''TF5' 'TF6' 'TF7' 'TF8' 'TF9,append) > > Is there anyway I can set the T's to REAL and declare how many numbers > beyond the decimal point I want? > > Thanks, > > Meredith > From munnich at ATMOS.UCLA.EDU Tue Feb 6 15:27:15 2007 From: munnich at ATMOS.UCLA.EDU (Matthias Munnich) Date: Tue, 6 Feb 2007 12:27:15 -0800 Subject: leap year Jennifer, Brian please note. In-Reply-To: Message-ID: -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hallo Karin, you are struck by a time conversion bug for monthly increments. I usually use mid-month for monthly averaged data. If this is OK with you, you can circumvent this bug by using a descriptor file and xdfopen. - -----------8<--------- dset /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc tdef 12 linear 15JAN1964 1mo - ----------->8--------- If the data point has to be at near the end of month use - -----------8<--------- dset /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc tdef 12 linear 28JAN1964 1mo - ----------->8--------- ** Jennifer, Brian ** I was able to reproduce this bug. Any file with - -----------8<--------- tdef 12 linear NNjan1964 1mo - ----------->8--------- with NN>29 gives me this error. It has nothing to do with leap years even tdef 12 linear 29jan1963 1mo give the same behavior: ga-> set t 1 Time values set: 1963:1:29:0 1963:1:29:0 ga-> set t 2 Time values set: 1963:3:1:0 1963:3:1:0 ga-> q dim Default file number is: 1 X is varying Lon = 299.252 to 420.754 X = 1.00001 to 109 Y is varying Lat = 9.533 to 80.185 Y = 1 to 64 Z is fixed Lev = 0 Z = 1 T is fixed Time = 00Z01MAR1963 T = 3 ga-> set t 3 Time values set: 1963:3:29:0 1963:3:29:0 ga-> q dim Default file number is: 1 X is varying Lon = 299.252 to 420.754 X = 1.00001 to 109 Y is varying Lat = 9.533 to 80.185 Y = 1 to 64 Z is fixed Lev = 0 Z = 1 T is fixed Time = 00Z29MAR1963 T = 3 ga-> q config Config: v1.9b4 32-bit little-endian readline sdf/xdf netcdf lats printim Note that this is not just labeling. Data at t=2 cannot be accessed at all. Instead data at t=3 are silently displayed. This is bad! ... Matt PS: Karin, ich war erstaunt Dich wieder auf der Gradsusr-liste zu finden. Ich hatte gedacht, Du hast beim DRKZ aufgehoert. Gruesse nach HH! Karin Meier-Fleischer wrote: > Hallo Matthias, > > schoen von Dir zu hoeren. > > > Here is my problem with the leap year: > > the data file contains 12 monthly time steps (from 23:59Z31JAN1964 to > 23:59Z31DEC1964) > > -1 : Date Time Code Level Size Miss : Minimum Mean Maximum > 1 : 1964-01-31 23:59 -1 0 42126 0 : 0.0000 25.063 557.66 > 2 : 1964-02-29 23:59 -1 0 42126 0 : 0.0000 25.040 1086.3 > 3 : 1964-03-31 23:59 -1 0 42126 0 : 0.0000 20.098 539.31 > 4 : 1964-04-30 23:59 -1 0 42126 0 : 0.0000 20.345 265.12 > 5 : 1964-05-31 23:59 -1 0 42126 0 : 0.0000 22.864 315.59 > 6 : 1964-06-30 23:59 -1 0 42126 0 : 0.0000 29.650 388.09 > 7 : 1964-07-31 23:59 -1 0 42126 0 : 0.0000 27.702 287.31 > 8 : 1964-08-31 23:59 -1 0 42126 0 : 0.0000 26.029 243.57 > 9 : 1964-09-30 23:59 -1 0 42126 0 : 0.0000 22.945 355.16 > 10 : 1964-10-31 23:59 -1 0 42126 0 : 0.0000 24.714 359.15 > 11 : 1964-11-30 23:59 -1 0 42126 0 : 0.0000 21.487 424.59 > 12 : 1964-12-31 23:59 -1 0 42126 0 : 0.0000 17.004 294.62 > cdo info: Processed 1 variable 12 timesteps. ( 0.02s ) > > > The GrADS script leap_year.gs has the following output: > ------------------------------------------------------- > > Scanning self-describing file: > /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc > Found 1 data variables in SDF file. > Found displayable variable rain_con with 0 levels in SDF file. > dset /pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc > title CLM simulation of the 20th century climate with anthropogenic > forcing (C20_1) for the SGA project (CLM_C20_1) > undef 9.96921e+36 > xdef 238 linear -10.6 0.200001 > ydef 177 linear 34.6 0.200001 > zdef 1 linear 0 0 > tdef 12 linear 23:59Z31JAN1964 1mo > vars 1 > rain_con 0 -999 convective rainfall > > -- t=1 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is fixed Time = 23:59Z31JAN1964 T = 1 > > -- t=2 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is fixed Time = 23:59Z02MAR1964 T = 3 > > -- t=3 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is fixed Time = 23:59Z31MAR1964 T = 3 > > -- t=12 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is fixed Time = 23:59Z31DEC1964 T = 12 > > -- time=23:59Z29FEB1964 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is fixed Time = 23:59Z02MAR1964 T = 3 > > -- t=1-12 -- > Default file number is: 1 > X is varying Lon = -10.6 to 36.8002 X = 1 to 238 > Y is varying Lat = 34.6 to 69.8001 Y = 1 to 177 > Z is fixed Lev = 0 Z = 1 > T is varying Time = 23:59Z31JAN1964 to 23:59Z31DEC1964 T = 1 to 12 > -------------------------------------------------------------------------- > > The plot shows a wrong value for february (it ist the same as march): > instead of value 1086.3 it plots 539.31. > > Any idea what is going wrong? > > Thanks in advance, > Karin > ---------------------------------\\--------------------------------------- > Karin Meier-Fleischer \\ Max-Planck-Institut fuer Meteorologie > Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten > FAX: +49-40-41173-476 \\ Bundesstr. 53 > Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg > http://www.mad.zmaw.de/ \\ Germany > ---------------------------------------\\--------------------------------- > > On Mon, 5 Feb 2007, Matthias Munnich wrote: > > Hi Karin, > > without a look at the descriptor file and a complete log of an example > session > it's hard to figure out what exactly is going on. Can you use "set t" > to set time > to 29Feb? > > ... Matt > > > Karin Meier-Fleischer wrote: > >>> Hi, > >>> > >>> I've got a data file containing a leap year (1964). If I try to > set the > >>> time to 23:59Z29Feb1964 GrADS always > >>> changes the time to 23:59Z02MAR1964. > >>> > >>> How could I prevent this? > >>> > >>> Thanks, Karin > >>> > ---------------------------------\\--------------------------------------- > >>> > >>> Karin Meier-Fleischer \\ Max-Planck-Institut fuer > >>> Meteorologie > >>> Phone: +49-40-41173-431 \\ Abteilung Modelle & Daten > >>> FAX: +49-40-41173-476 \\ Bundesstr. 53 > >>> Email: karin.meier-fleischer at zmaw.de \\ D-20146 Hamburg > >>> http://www.mad.zmaw.de/ \\ Germany > >>> > ---------------------------------------\\--------------------------------- > >>> > > >> > ------------------------- > function LeapYear(args) > _DATA='/pf/k/k204045/data/CLM/CLM_C20_1_D3_ms_RAIN_CON_49-60.nc' > 'sdfopen '_DATA > 'q ctlinfo'; say result > say '-- t=1 --' > 'set t 1' > 'q dims'; say result > say '-- t=2 --' > 'set t 2' > 'q dims'; say result > say '-- t=3 --' > 'set t 3' > 'q dims'; say result > say '-- t=12 --' > 'set t 12' > 'q dims'; say result > say '-- time=23:59Z29FEB1964 --' > 'set time 23:59Z29FEB1964' > 'q dims'; say result > say '-- t=1-12 --' > 'set t 1 12' > 'q dims'; say result > 'set x 1' > 'set y 1' > 'd tloop(rain_con)' > 'printim leap_year.gif gif white' > return > ------------------------- - -- - -------------------------------- Matthias Munnich Univ. of California, Los Angeles Dept. of Atmos. and Oceanic Sc. and Inst. of Geophysics and Planetary Physics 3845 Slichter Hall Los Angeles, CA 90095-1567 Phone: +1-310-794 5899 Fax: +1-310-206 3051 Email: munnich at atmos.ucla.edu - -------------------------------- -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFFyOSiidaDTh/cfhYRAnnMAJ9oAJH5MRGLPxXzoqDQenf2uGnuCQCgoyt0 TCKiIWgHYErksY5mVUsipnw= =xai3 -----END PGP SIGNATURE----- From mulero at WISC.EDU Tue Feb 6 18:38:51 2007 From: mulero at WISC.EDU (PEDRO J MULERO) Date: Tue, 6 Feb 2007 17:38:51 -0600 Subject: grads2nc question Message-ID: Hello GrADS Users, Has anybody ever used the grads2nc application? I'm trying to execute the grads2nc application using control files for gradsv1.9, not 1.5x as preferred. This is the error message I get: zvals (53) = 0.010000 0.050000 0.100000 0.250000 0.500000 0.750000 1.000000 1.250000 1.500000 1.700000 2.000000 2.250000 2.500000 2.750000 3.000000 3.250000 3.500000 3.750000 4.000000 4.250000 4.500000 4.750000 5.000000 5.250000 5.500000 5.750000 6.000000 6.250000 6.500000 6.750000 7.000000 7.250000 7.500000 7.750000 8.000000 8.250000 8.500000 8.750000 9.000000 9.250000 9.500000 9.750000 10.000000 10.500000 11.000000 11.500000 12.000000 12.500000 13.000000 13.500000 14.000000 14.500000 15.000000 Invalid TDEF starting time "12z04aug2004"! Cannot read the IEEE data file! Parsing aborted! ---------------- Do you know what changes do I need to do to the control file in order for the program to run effectively? Does any of you have a grads1.5x control file that I can use to make the proper changes? I'd greatly appreciate any help on this matter. Best regards, Pedro **************************************************** Pedro J. Mulero Research Assistant University of Wisconsin Space Science and Engineering Center email: mulero at wisc.edu office: AOSS 1439, (608) 262-8954 mobile: (787) 667-8653 **************************************************** From Charles.Seman at NOAA.GOV Tue Feb 6 18:46:57 2007 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Tue, 6 Feb 2007 18:46:57 -0500 Subject: Fwrite Gr2stn Output In-Reply-To: Message-ID: Stephen, Thanks, I wasn't aware of this... Chuck Stephen R McMillan wrote: > Chuck, > I tried your suggestion but got the error "define does not yet support > station data...." I forgot to mention I'm using V1.8SL11 on a Windows > XP. Perhaps it works on the newer version? > > I'll pursue Jennifer Adams' subsequent response to my question. Thanks > anyway! > Stephen Mc > -----GRADSUSR at LIST.CINECA.IT wrote: ----- > > To: GRADSUSR at LIST.CINECA.IT > From: Charles Seman > Sent by: GRADSUSR at LIST.CINECA.IT > Date: 02/02/2007 04:22PM > Subject: Re: Fwrite Gr2stn Output > > Stephen, > > Have you tried defining a variable, like: > > define tmp2mc = gr2stn(tmp2m-273.16, -75.83, 40.05) > > and then tried writing this to an fwrite file? > > Hope this helps, > Chuck > > Stephen R McMillan wrote: > > > > Hello, > > > > Is it possible to fwrite time series output, for a point using the > > gr2stn interpolation function? I can do it for a specified > gridpoint, > > but not for an interpolated point using gr2stn. For example: > > > > I can do following without a problem (say, output from the GFS > model): > > > > set t 1 last > > set lat 40 > > set lon -76 > > set gxout fwrite > > set fwrite d:\temp\test.dat > > d tmp2m-273.16 > > disable fwrite > > > > However, using same file and settings but for a point > interpolated to > > -75.83 and 40.05: > > > > d gr2stn(tmp2m-273.16, -75.83, 40.05) > > > > Above will display fine on screen, but will not fwrite to specified > > file. Any help/suggestions would be appreciated (including "not > > possible"). I am aware there are other methods such as using the > 'set > > gxout print' and write functions, for example (and have used > > successfully), but it would be much simpler and more efficient to > > output to a binary file using fwrite. > > > > Thanks-- > > Stephen McMillan > > Planalytics > > > > *************************************************** > > 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. > > *************************************************** > > -- > > 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 reflect > any position of the Government or NOAA." > > > > > *************************************************** > 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. > *************************************************** -- 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 reflect any position of the Government or NOAA." From arabonis at EGS.UCT.AC.ZA Wed Feb 7 01:10:39 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Wed, 7 Feb 2007 08:10:39 +0200 Subject: GFS files corrupted on monsoondata Message-ID: Hi I posted this a few days ago and we still seem to be missing the first few days on the gfs dataset on monsoondata. Thanks Jean Pierre Hi Jennifer Seems the first few days of the gfs files on monsoondata are corrupted. Thanks Jean Pierre -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free From mahesh_cgmd at REDIFFMAIL.COM Wed Feb 7 08:22:26 2007 From: mahesh_cgmd at REDIFFMAIL.COM (Mahesh Shinde) Date: Wed, 7 Feb 2007 13:22:26 -0000 Subject: zonal average Message-ID: ?Dear frds, I wanted to do zonal average of a parameter between 20-40N:0-360.My code is as follows ... ... set lon -90 set lat 20 40 d ave(z,lon=0,lon=360) ... ... but i am not getting desireable results.... can anybody help me to calculate zonally average, if somebody has done it already. thank you -mahesh Mr.Mahesh A Shinde Climate & Global Modeling Division Indian Institute of Tropical Meteorology Pashan, Pune:411008 Maharashtra INDIA -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070207/26585f65/attachment.html From jsg229 at PSU.EDU Wed Feb 7 09:56:30 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Wed, 7 Feb 2007 09:56:30 -0500 Subject: height-radius plot Message-ID: I was hoping someone could help me with this problem that I have been stuck on for a while. I want to make a height-radius contour plot of temperature at 1 given time, i.e. have radius from 0 to 150km on the x-axis and height z=39m to 18000m on the y-axis. I open the file using: OPEN(541,FILE='0.3sf_grads.dat',FORM='UNFORMATTED', & ACCESS='DIRECT',RECL=KX*IX*4) where IX=32 and KX=30, that is 32 points in the x-dir and 30 points in the y-dir. I then write the radius, height and temperature to the file: IREC=1 WRITE(541,REC=IREC) ((R(K,I),K=1,KX),I=1,IR) IREC=IREC+1 WRITE(541,REC=IREC) ((Z(K,I),K=1,KX),I=1,IR) IREC=IREC+1 WRITE(541,REC=IREC) ((T(K,I),K=1,KX),I=1,IR) IREC=IREC+1 It should be noted that R(K,1)=0.0 and R(K,IR)=155.0km. Also, Z(1,I)=39.0m and Z(KX,I)=18000m. Now, I want to create the corresponding .ctl file (0.3sf_grads.ctl) such that temperature is contoured on the r-z place. I have checked my tempearture, radius, and height values and I know they are correct. I would like the distance (radius) in km to be plotted on the x-axis and height in m on the y-axis. Any help in fixing this problem would be most appreciated, as I havent been able to find a solution to the problem. Thanks. Jeff -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070207/e63ac1c2/attachment.html From Mary-Jo.Nath at NOAA.GOV Wed Feb 7 10:13:13 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Wed, 7 Feb 2007 10:13:13 -0500 Subject: height-radius plot In-Reply-To: <1170860188l.798856l.0l@psu.edu> Message-ID: Jeffrey, Assuming that R is not a function of Z, instead of writing R and Z to the file as data, put them into the CTL file as the dimension levels, something like this: DSET 0.3sf_grads.dat byteswapped TITLE Sample UNDEF -9999 XDEF 30 levels 0 (add other radius values here) 155 YDEF 1 linear 1 1 ZDEF 32 levels 39 (add other heights here) 18000 TDEF 4 linear 12:00Z15SEP1999 180MN vars 1 t_az 32 99 TEMP AZ. AVG. endvars (NOTE: If X is equally spaced, then you can use XDEF 30 LINEAR startingX incrementX Same with the ZDEF.) Write only the t_az data to the file. When you plot, "set z 1 32", set your time, and display. -MJ JEFFREY S GALL wrote: > I was hoping someone could help me with this problem that I have been > stuck on for a while. > I want to make a height-radius contour plot of temperature at 1 given > time, i.e. have > radius from 0 to 150km on the x-axis and height z=39m to 18000m on the > y-axis. > > I open the file using: > > OPEN(541,FILE='0.3sf_grads.dat',FORM='UNFORMATTED', > & ACCESS='DIRECT',RECL=KX*IX*4) > > where IX=32 and KX=30, that is 32 points in the x-dir and 30 points in > the y-dir. > > I then write the radius, height and temperature to the file: > > IREC=1 > WRITE(541,REC=IREC) ((R(K,I),K=1,KX),I=1,IR) > IREC=IREC+1 > WRITE(541,REC=IREC) ((Z(K,I),K=1,KX),I=1,IR) > IREC=IREC+1 > WRITE(541,REC=IREC) ((T(K,I),K=1,KX),I=1,IR) > IREC=IREC+1 > > It should be noted that R(K,1)=0.0 and R(K,IR)=155.0km. Also, > Z(1,I)=39.0m and > Z(KX,I)=18000m. > > Now, I want to create the corresponding .ctl file (0.3sf_grads.ctl) > such that temperature > is contoured on the r-z place. I have checked my tempearture, radius, > and height values > and I know they are correct. I would like the distance (radius) in km > to be plotted on > the x-axis and height in m on the y-axis. > > Any help in fixing this problem would be most appreciated, as I havent > been able to find > a solution to the problem. Thanks. > > Jeff > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070207/b414deeb/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: Mary-Jo.Nath.vcf Type: text/x-vcard Size: 217 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070207/b414deeb/attachment.vcf From horvath at CIRUS.DHZ.HR Wed Feb 7 10:55:48 2007 From: horvath at CIRUS.DHZ.HR (Kristian Horvath) Date: Wed, 7 Feb 2007 16:55:48 +0100 Subject: ncdump In-Reply-To: Message-ID: Dear GrADS-ers, I am wandering how to dump .nc files in a table-like (row-column) ascii format e.g. lat-lon for a 2D variable in one time instant? I tried to use ncdump, but its output looks somewhat different. Of course, I could write a program to change the format of the ncdump output, but is it the only way to go? Many thanks. Kristian _________________________________ Kristian Horvath email: horvath at cirus.dhz.hr Graduate Research Associate tel:+385-1-4565752 Meteorological and Hydrological Service www: http://www.dhmz.htnet.hr/index.php Gric 3, 10000 Zagreb, Croatia personal www: http://radar.dhz.hr/~horvath From heiner at MISU.SU.SE Wed Feb 7 11:28:38 2007 From: heiner at MISU.SU.SE (=?ISO-8859-1?Q?Heiner_K=F6rnich?=) Date: Wed, 7 Feb 2007 17:28:38 +0100 Subject: ncdump In-Reply-To: <7.0.1.0.0.20070207164537.02056688@cirus.dhz.hr> Message-ID: Hi, you probably could use ncl_filedump. You find it here: http://www.ncl.ucar.edu/Document/Tools/ncl_filedump.shtml But you have to register to use it. Heiner Kristian Horvath wrote: > Dear GrADS-ers, > > I am wandering how to dump .nc files in a table-like (row-column) > ascii format e.g. lat-lon for a 2D variable in one time instant? I > tried to use ncdump, but its output looks somewhat different. Of > course, I could write a program to change the format of the ncdump > output, but is it the only way to go? > > Many thanks. > > Kristian > > _________________________________ > Kristian Horvath email: > horvath at cirus.dhz.hr > Graduate Research Associate tel:+385-1-4565752 > Meteorological and Hydrological Service www: > http://www.dhmz.htnet.hr/index.php > Gric 3, 10000 Zagreb, Croatia personal www: > http://radar.dhz.hr/~horvath -- Heiner K?rnich Dept. of Meteorology Tel: +46 8 164333 Stockholms University, SE-106 91 Stockholm, Sweden Email: heiner at misu.su.se www.misu.su.se/~heiner/ From labyerle at MET.UTAH.EDU Wed Feb 7 11:25:07 2007 From: labyerle at MET.UTAH.EDU (Lee Byerle) Date: Wed, 7 Feb 2007 09:25:07 -0700 Subject: ncdump In-Reply-To: <7.0.1.0.0.20070207164537.02056688@cirus.dhz.hr> Message-ID: Kristian, You may use unix to get netcdf data to ascii (w/o ",", ";" etc): for example: ncdump -v u10 myfile.nc > u10_din_dump #remove extraneous characters sed 's/,/ /g' u10_din_dump > temp1.out sed 's/;/ /g' temp1.out > temp2.out sed 's/}/ /g' temp2.out > temp3.out #Strip top lines of header tail +94 temp3.out > new_u10.out new_u10.out is ascii file. Lee > Dear GrADS-ers, > > I am wandering how to dump .nc files in a table-like (row-column) > ascii format e.g. lat-lon for a 2D variable in one time instant? I > tried to use ncdump, but its output looks somewhat different. Of > course, I could write a program to change the format of the ncdump > output, but is it the only way to go? > > Many thanks. > > Kristian > > _________________________________ > Kristian Horvath email: > horvath at cirus.dhz.hr > Graduate Research Associate tel:+385-1-4565752 > Meteorological and Hydrological Service www: > http://www.dhmz.htnet.hr/index.php > Gric 3, 10000 Zagreb, Croatia personal www: > http://radar.dhz.hr/~horvath > From doty at COLA.IGES.ORG Wed Feb 7 12:46:54 2007 From: doty at COLA.IGES.ORG (Brian E. Doty) Date: Wed, 7 Feb 2007 12:46:54 -0500 Subject: leap year Jennifer, Brian please note. Message-ID: Matthias, thanks, we will look into it! Around here, we usually use the first of the month for monthly means (eg, 1may1998) ...Brian From arabonis at EGS.UCT.AC.ZA Thu Feb 8 08:57:00 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 8 Feb 2007 15:57:00 +0200 Subject: Contour Labels Message-ID: Hi All I'm trying to figure out how to make contour labels without a white background, ie I would like them to "float" on the background image. Many Thanks Jean Pierre -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free From jsg229 at PSU.EDU Thu Feb 8 10:43:28 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Thu, 8 Feb 2007 10:43:28 -0500 Subject: change x-axis label Message-ID: When I create a z-x contour plot, I have radius plotted on the x-axis. When I display my contour plot, everything is correct except the label. It goes something like 0 20E 40E 60E.... Rather, the units of the x-axis are 0 20km 40km 60km.... Is there anyway I can turn off the plotting of the 'E' or add the 'km' instead of the 'E'? Thanks, Jeff -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/dc54f0ae/attachment.html From Mary-Jo.Nath at NOAA.GOV Thu Feb 8 10:52:25 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Thu, 8 Feb 2007 10:52:25 -0500 Subject: change x-axis label In-Reply-To: <1170949405l.1249324l.0l@psu.edu> Message-ID: |Jeffrey, | |Try: | |set mproj off set xlab %0fkm | | You can look up these commands and their options in the index at http://grads.iges.org/grads/gadoc/gadocindex.html | |MJ | JEFFREY S GALL wrote: > When I create a z-x contour plot, I have radius plotted on the > x-axis. When I display my contour plot, everything is correct except > the label. It goes something like > > 0 20E 40E 60E.... > > Rather, the units of the x-axis are > 0 20km 40km 60km.... > > Is there anyway I can turn off the plotting of the 'E' or add the 'km' > instead of the 'E'? > > Thanks, > > Jeff > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/92c71e8b/attachment.html From Mary-Jo.Nath at NOAA.GOV Thu Feb 8 10:57:16 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Thu, 8 Feb 2007 10:57:16 -0500 Subject: change x-axis label In-Reply-To: <1170949405l.1249324l.0l@psu.edu> Message-ID: Jeffrey, Oops! Try this instead: set mproj off set xlab %.0km The decimal point in the "set xlab" is important. MJ > > |Jeffrey, > | > > |Try: > | > > |set mproj off > set xlab %0fkm > | > > | > You can look up these commands and their options in the index at > http://grads.iges.org/grads/gadoc/gadocindex.html > | > > |MJ > | > JEFFREY S GALL wrote: > When I create a z-x contour plot, I have radius plotted on the > x-axis. When I display my contour plot, everything is correct except > the label. It goes something like > > 0 20E 40E 60E.... > > Rather, the units of the x-axis are > 0 20km 40km 60km.... > > Is there anyway I can turn off the plotting of the 'E' or add the 'km' > instead of the 'E'? > > Thanks, > > Jeff > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/79e6ac12/attachment.html From robock at ENVSCI.RUTGERS.EDU Thu Feb 8 11:03:24 2007 From: robock at ENVSCI.RUTGERS.EDU (Alan Robock) Date: Thu, 8 Feb 2007 11:03:24 -0500 Subject: stripes in eps output Message-ID: Dear Friends, Has anyone found a solution to the age-old problem of horizontal lines between bands of color in shaded grads output? We have eps files to submit for an accepted JGR paper and are concerned about how the final version will look. Is there anything we can do to the files before we submit them to fix the problem? Alan Alan Robock, Professor II Department of Environmental Sciences Phone: +1-732-932-9478 Rutgers University Fax: +1-732-932-8644 14 College Farm Road E-mail: robock at envsci.rutgers.edu New Brunswick, NJ 08901-8551 USA http://envsci.rutgers.edu/~robock From doty at COLA.IGES.ORG Thu Feb 8 11:37:11 2007 From: doty at COLA.IGES.ORG (Brian E. Doty) Date: Thu, 8 Feb 2007 11:37:11 -0500 Subject: stripes in eps output Message-ID: Working on it... Brian From jsg229 at PSU.EDU Thu Feb 8 12:05:50 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Thu, 8 Feb 2007 12:05:50 -0500 Subject: gxeps error Message-ID: I use the gxeps command to convert from grads meta file to a postscript file. When I use the command gxeps -R -c -i 'filename'.gm -o ' filename'.eps I get a segmentation fault. Any idea as to why this is happening? I am running on a 64-bit linux system. Thanks, Jeff -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/c91ceeb4/attachment.html From jma at COLA.IGES.ORG Thu Feb 8 12:21:51 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 8 Feb 2007 12:21:51 -0500 Subject: stripes in eps output In-Reply-To: Message-ID: Alan et al., Those horizontal stripes in the ps/eps output are a real nuisance. Potential sources of the problem start with the contouring algorithm in the GrADS graphics layer, but the presence of the stripes also seem to be influenced by gxps, gxeps, any external or internal utility that converts the .eps or .ps files to something else as they get incorporated into written documents, and the printers used by scientists and publishers also seem to be a factor. My recent experience is that it's nearly impossible to predict whether you'll get stripes or not. It's seems more likely that you'll have stripes than no stripes, and the stripes also seem to be most garish in plots with shaded contours and white space. Here at COLA, we've begun using image format files (.png and .gif) instead of vector graphic formats, simply because they're more reliable even though they don't always look as good. COLA users are being weaned from using PCs running MS Windows, but my sense is that the gv32 program allowed for conversion of GrADS metafiles into something else that did not have stripes. I can't verify that, perhaps the windows users out there can confirm or deny. Brian is planning to rework the contouring algorithm, so there's hope for the future; unfortunately, that hope doesn't help in the short term. I would see if you can get the publisher to provide a test print of your figures, and if that looks bad, see if you can use an image format instead. A .png at 850x1100 resolution, shrunk down to journal size, might not look so bad. Jennifer On Feb 8, 2007, at 11:03 AM, Alan Robock wrote: > Dear Friends, > > Has anyone found a solution to the age-old problem of horizontal lines > between bands of color in shaded grads output? We have eps files to > submit for an accepted JGR paper and are concerned about how the final > version will look. > > Is there anything we can do to the files before we submit them to fix > the problem? > > Alan > > Alan Robock, Professor II > Department of Environmental Sciences Phone: > +1-732-932-9478 > Rutgers University Fax: > +1-732-932-8644 > 14 College Farm Road E-mail: > robock at envsci.rutgers.edu > New Brunswick, NJ 08901-8551 USA http://envsci.rutgers.edu/ > ~robock > -- 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/20070208/a9022ece/attachment.html From jsg229 at PSU.EDU Thu Feb 8 12:44:02 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Thu, 8 Feb 2007 12:44:02 -0500 Subject: fix for stripes Message-ID: to fix stripes try: * print file to .png format 'enable print 'filename'.gm' 'print' 'disable print' *'!gxeps -R -c -i 'filename'.gm -o ' filename'.eps' '!gxps -c -i 'filename'.gm -o ' filename'.eps' '!convert +antialias -density 300 -geometry 1024x768 'filename'.eps 'filename'.png' -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/c0ef5e94/attachment.html From herb.philbrick at MITTALSTEEL.COM Thu Feb 8 14:02:43 2007 From: herb.philbrick at MITTALSTEEL.COM (Herb Philbrick) Date: Thu, 8 Feb 2007 13:02:43 -0600 Subject: Regrid2 on MS XP Message-ID: Hello all, I'm trying to use Regrid2 on a MS Windows XP machine. Grads writes the input file to the udf and my compiled Regrid2.exe reads the contents fine. The trouble is Regrid2 encounters an "array bounds exceeded" error. I sure would appreciate hearing from a user of Regrid2 on an MS Windows PC. Thanks kindly, Herb Philbrick -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070208/95967d9c/attachment.html From munnich at ATMOS.UCLA.EDU Thu Feb 8 14:13:05 2007 From: munnich at ATMOS.UCLA.EDU (Matthias Munnich) Date: Thu, 8 Feb 2007 11:13:05 -0800 Subject: gxeps error In-Reply-To: <1170954349l.1441834l.0l@psu.edu> Message-ID: -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi Jeff, I am getting this as well. When I find time, I will investigate and get back to you. I didn't realized this before because I nowadays rarely, if ever use gxeps in standalone mode. Gxeps has been incorporated into the print command of v1.9b4. If "enable print" hasn't been issued, you can use ga-> print myplot.eps to create an EPS of the current plot. See ga -> print -h for more options(*). Somehow this doesn't segfault on 64bit Linux. Hope this helps, ... Matt (*) Jennifer, this appears not yet to be documented. JEFFREY S GALL wrote: > I use the gxeps command to convert from grads meta file to a > postscript file. When I use the command > > gxeps -R -c -i 'filename'.gm -o ' filename'.eps > > I get a segmentation fault. Any idea as to why this is happening? > I am running on a 64-bit linux system. > > > > Thanks, > > > > Jeff > - -- - -------------------------------- Matthias Munnich Univ. of California, Los Angeles Dept. of Atmos. and Oceanic Sc. and Inst. of Geophysics and Planetary Physics 3845 Slichter Hall Los Angeles, CA 90095-1567 Phone: +1-310-794 5899 Fax: +1-310-206 3051 Email: munnich at atmos.ucla.edu - -------------------------------- -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFFy3ZBidaDTh/cfhYRAuH8AKC1v8fx/5AXQGfwNrAWSK/JzFijQQCgjLPN N8TtpifNj0qMzzSF9KWD7ZQ= =PWqd -----END PGP SIGNATURE----- From munnich at ATMOS.UCLA.EDU Thu Feb 8 15:27:39 2007 From: munnich at ATMOS.UCLA.EDU (Matthias Munnich) Date: Thu, 8 Feb 2007 12:27:39 -0800 Subject: stripes in eps output In-Reply-To: Message-ID: -----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 The trouble is that one of the most common PostScript renderers, Ghostscript (gs), is nowadays doing antialiasing by default to soften ragged lines and corners. It blends neighboring bit even for purely horizontal lines. This leads to the pixel-width horizontal lines due to the way Grads is filling contours. Antialiasing can be turned of with "-dGraphicsAlphaBits=1" in gs or "-noantialias" in ghostview I still don't believe that hardcopy publications are likely to show lines. These lines would be close to invisible anyway on a 1200dpi printer . I suspect people are rendering their graphics themselves with tools like ImageMagick. Almost any free software is using Ghostscript as the rendering engine. These tools usually use default Ghostscript settings and you get the lines. Well, Jennifer's email suggests otherwise. Unfortunately you do not control how the graphic is rendered in print and submitting a rendered figure as suggested by Jennifer and Jeff is one solution. However, you loose print quality. My approach would be to look through a recent issues of JGR and see if the figures show these lines. I always recognize grads plots -- I guess I stared at too many by now. I am pretty sure it's fine to submit gxeps created EPS files to JGR. There were no lines in the figures of my recent GRL paper. ... Matt Jennifer Adams wrote: > Alan et al., Those horizontal stripes in the ps/eps output are a > real nuisance. Potential sources of the problem start with the > contouring algorithm in the GrADS graphics layer, but the presence > of the stripes also seem to be influenced by gxps, gxeps, any > external or internal utility that converts the .eps or .ps files to > something else as they get incorporated into written documents, and > the printers used by scientists and publishers also seem to be a > factor. My recent experience is that it's nearly impossible to > predict whether you'll get stripes or not. It's seems more likely > that you'll have stripes than no stripes, and the stripes also seem > to be most garish in plots with shaded contours and white space. > > Here at COLA, we've begun using image format files (.png and .gif) > instead of vector graphic formats, simply because they're more > reliable even though they don't always look as good. > > COLA users are being weaned from using PCs running MS Windows, but > my sense is that the gv32 program allowed for conversion of GrADS > metafiles into something else that did not have stripes. I can't > verify that, perhaps the windows users out there can confirm or > deny. > > Brian is planning to rework the contouring algorithm, so there's > hope for the future; unfortunately, that hope doesn't help in the > short term. I would see if you can get the publisher to provide a > test print of your figures, and if that looks bad, see if you can > use an image format instead. A .png at 850x1100 resolution, shrunk > down to journal size, might not look so bad. > > Jennifer > > > > On Feb 8, 2007, at 11:03 AM, Alan Robock wrote: > >> Dear Friends, >> >> Has anyone found a solution to the age-old problem of horizontal >> lines between bands of color in shaded grads output? We have eps >> files to submit for an accepted JGR paper and are concerned about >> how the final version will look. >> >> Is there anything we can do to the files before we submit them to >> fix the problem? >> >> Alan >> >> Alan Robock, Professor II Department of Environmental Sciences >> Phone: +1-732-932-9478 Rutgers University >> Fax: +1-732-932-8644 14 College Farm Road >> E-mail: robock at envsci.rutgers.edu >> New Brunswick, NJ 08901-8551 >> USA http://envsci.rutgers.edu/~robock >> > > -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705 jma at cola.iges.org > > > - -- - -------------------------------- Matthias Munnich Univ. of California, Los Angeles Dept. of Atmos. and Oceanic Sc. and Inst. of Geophysics and Planetary Physics 3845 Slichter Hall Los Angeles, CA 90095-1567 Phone: +1-310-794 5899 Fax: +1-310-206 3051 Email: munnich at atmos.ucla.edu - -------------------------------- -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.2.2 (GNU/Linux) Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org iD8DBQFFy4e7idaDTh/cfhYRAkhyAKDYrNmXqc5O4m5L+DWrypw/kvsk4QCgn2rl Oo2eoszDJQ4TTyuFviUL6g0= =/obd -----END PGP SIGNATURE----- From Mary-Jo.Nath at NOAA.GOV Thu Feb 8 15:49:54 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Thu, 8 Feb 2007 15:49:54 -0500 Subject: stripes in eps output In-Reply-To: Message-ID: I create figures for AMS journals all the time using GrADs. I use "gxps" to convert the metafile to postscript, run them through Adobe Illustrator to convert from RGB to CMYK and tweak the colors, labels, contour labels, etc.. I've never had problems with those funky shading lines showing up in the final figures in the publication. MJ Alan Robock wrote: > Dear Friends, > > Has anyone found a solution to the age-old problem of horizontal lines > between bands of color in shaded grads output? We have eps files to > submit for an accepted JGR paper and are concerned about how the final > version will look. > > Is there anything we can do to the files before we submit them to fix > the problem? > > Alan > > Alan Robock, Professor II > Department of Environmental Sciences Phone: +1-732-932-9478 > Rutgers University Fax: +1-732-932-8644 > 14 College Farm Road E-mail: robock at envsci.rutgers.edu > New Brunswick, NJ 08901-8551 USA http://envsci.rutgers.edu/~robock > From salathe at WASHINGTON.EDU Thu Feb 8 16:13:15 2007 From: salathe at WASHINGTON.EDU (=?ISO-8859-1?Q?Eric_Salath=E9?=) Date: Thu, 8 Feb 2007 13:13:15 -0800 Subject: stripes in eps output In-Reply-To: <45CB87BB.7080403@atmos.ucla.edu> Message-ID: My experience has also been that EPS files generated from grads using gxps look just fine when printed in journals -- most of my paper have GrADs images and I would not be using Grads if the images looked bad in print. I usually tweak and reformat graphics in Illustrator, which deals with the the lines just fine and makes very nice raster files if needed. The problem comes with looking at the EPS file (eg in ghostview or Preview) or when converting to a raster image. Since most images end up on the web as PNG, and not in print, this is important. Grad's printim does a decent job and is quick, but without a lot of control over image size and format. The solution is to turn of anti aliasing. Using +antialias in ImageMagick convert as Jeffery Gall suggests works reasonably well (convert actually calls gs), but then lines are rendered too thick and pixelated. If you specify a fine dpi, then scale down, you can get better line art while eliminating the striping. Compare these results on an eps generated with gxps (note I've added -trim to remove blank margins): convert -trim +antialias -density 150 test.eps test.png convert -trim +antialias -scale 25% -density 600 test.eps test_scale.png To do this, I use a grads script that calls a csh script that can handle any output format: Grads script saveimg in $GASCRP: function saveimg(file) * print screen to eps file specified as argument * first print to gx file 'enable print tmp.gx' 'print' 'disable print' * now translate to image format '! gxconvert tmp.gx ' file '! \rm -f tmp.gx' csh script in $PATH: #!/bin/csh if (! $#argv == 1) then echo "usage: gxconvert file.gx file.ext where ext indicates output format" exit endif gxps -c -i $1 -o tmp.ps > /dev/null if ($2:e == 'ps') then mv tmp.ps $2 else if ($2:e == 'eps') then ps2epsi tmp.ps $2 \rm -f tmp.ps else convert +antialias -trim -scale 50% -density 300 tmp.ps $2 \rm -f tmp.ps endif exit -Eric -- Eric Salath? CSES Climate Impacts Group University of Washington 206-616-5351 On Feb 8, 2007, at 12:27 PM, Matthias Munnich wrote: > -----BEGIN PGP SIGNED MESSAGE----- > Hash: SHA1 > > The trouble is that one of the most common PostScript renderers, > Ghostscript (gs), is nowadays doing antialiasing by default to soften > ragged lines and corners. It blends neighboring bit even for purely > horizontal lines. This leads to the pixel-width horizontal lines due > to the way Grads is filling contours. Antialiasing can be turned of > with "-dGraphicsAlphaBits=1" in gs or "-noantialias" in ghostview > > I still don't believe that hardcopy publications are likely to show > lines. These lines would be close to invisible anyway on a 1200dpi > printer . I suspect people are rendering their graphics themselves > with tools like ImageMagick. Almost any free software is using > Ghostscript as the rendering engine. These tools usually use default > Ghostscript settings and you get the lines. > > Well, Jennifer's email suggests otherwise. Unfortunately you do not > control how the graphic is rendered in print and submitting a rendered > figure as suggested by Jennifer and Jeff is one solution. However, > you loose print quality. > > My approach would be to look through a recent issues of JGR and see if > the figures show these lines. I always recognize grads plots -- I > guess I stared at too many by now. I am pretty sure it's fine to > submit gxeps created EPS files to JGR. There were no lines in the > figures of my recent GRL paper. > > > ... Matt > > > > Jennifer Adams wrote: >> Alan et al., Those horizontal stripes in the ps/eps output are a >> real nuisance. Potential sources of the problem start with the >> contouring algorithm in the GrADS graphics layer, but the presence >> of the stripes also seem to be influenced by gxps, gxeps, any >> external or internal utility that converts the .eps or .ps files to >> something else as they get incorporated into written documents, and >> the printers used by scientists and publishers also seem to be a >> factor. My recent experience is that it's nearly impossible to >> predict whether you'll get stripes or not. It's seems more likely >> that you'll have stripes than no stripes, and the stripes also seem >> to be most garish in plots with shaded contours and white space. >> >> Here at COLA, we've begun using image format files (.png and .gif) >> instead of vector graphic formats, simply because they're more >> reliable even though they don't always look as good. >> >> COLA users are being weaned from using PCs running MS Windows, but >> my sense is that the gv32 program allowed for conversion of GrADS >> metafiles into something else that did not have stripes. I can't >> verify that, perhaps the windows users out there can confirm or >> deny. >> >> Brian is planning to rework the contouring algorithm, so there's >> hope for the future; unfortunately, that hope doesn't help in the >> short term. I would see if you can get the publisher to provide a >> test print of your figures, and if that looks bad, see if you can >> use an image format instead. A .png at 850x1100 resolution, shrunk >> down to journal size, might not look so bad. >> >> Jennifer >> >> >> >> On Feb 8, 2007, at 11:03 AM, Alan Robock wrote: >> >>> Dear Friends, >>> >>> Has anyone found a solution to the age-old problem of horizontal >>> lines between bands of color in shaded grads output? We have eps >>> files to submit for an accepted JGR paper and are concerned about >>> how the final version will look. >>> >>> Is there anything we can do to the files before we submit them to >>> fix the problem? >>> >>> Alan >>> >>> Alan Robock, Professor II Department of Environmental Sciences >>> Phone: +1-732-932-9478 Rutgers University >>> Fax: +1-732-932-8644 14 College Farm Road >>> E-mail: robock at envsci.rutgers.edu >>> New Brunswick, NJ 08901-8551 >>> USA http://envsci.rutgers.edu/~robock >>> >> >> -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 >> Calverton, MD 20705 jma at cola.iges.org >> >> >> > > > - -- > - -------------------------------- > Matthias Munnich > Univ. of California, Los Angeles > Dept. of Atmos. and Oceanic Sc. and > Inst. of Geophysics and Planetary Physics > 3845 Slichter Hall > Los Angeles, CA 90095-1567 > Phone: +1-310-794 5899 > Fax: +1-310-206 3051 > Email: munnich at atmos.ucla.edu > - -------------------------------- > -----BEGIN PGP SIGNATURE----- > Version: GnuPG v1.4.2.2 (GNU/Linux) > Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org > > iD8DBQFFy4e7idaDTh/cfhYRAkhyAKDYrNmXqc5O4m5L+DWrypw/kvsk4QCgn2rl > Oo2eoszDJQ4TTyuFviUL6g0= > =/obd > -----END PGP SIGNATURE----- From gsclarke at EMIRATES.NET.AE Fri Feb 9 02:05:32 2007 From: gsclarke at EMIRATES.NET.AE (Gregory Clarke) Date: Fri, 9 Feb 2007 08:05:32 +0100 Subject: Shading Undefined Values Message-ID: Good Day All, I would like to know if it is possible to shade the areas of gridded data that are undefined. I have the the following line in my control file... undef -99999 If I try to set the contour values in the script to work with this value, GrADS does not shade the area; e.g. 'set cmin -100000' 'set cmax -99990' 'set cint 1' 'd myvar' I have gone through the GrADS manual and archive, and at the moment the only possible way is to "define" the undefined area with some loop testing the value at each grid point, and then plotting the redefined variable. Any suggestions will be greatly appreciated. Kind regards Greg Clarke Support: OMAA From davide.sacchetti at ARPAL.ORG Fri Feb 9 03:13:45 2007 From: davide.sacchetti at ARPAL.ORG (Davide Sacchetti) Date: Fri, 9 Feb 2007 09:13:45 +0100 Subject: Shading Undefined Values In-Reply-To: <20070209070614.000EE21B7F@mx2.cineca.it> Message-ID: grads can't display undef values by definition of "undefined value" ... if you wish to display undefined values a chance is: suppose you know that your data are all positive or undefined (ex temperature in K), 'set clevs -5' 'set ccols 2 0' 'display const(yourdata,-10,-u) it should work bye bye Davide On Fri, 2007-02-09 at 08:05 +0100, Gregory Clarke wrote: > Good Day All, > > I would like to know if it is possible to shade the areas of gridded data > that are undefined. > > I have the the following line in my control file... > > undef -99999 > > If I try to set the contour values in the script to work with this value, > GrADS does not shade the area; e.g. > > 'set cmin -100000' > 'set cmax -99990' > 'set cint 1' > 'd myvar' > > I have gone through the GrADS manual and archive, and at the moment the > only possible way is to "define" the undefined area with some loop testing > the value at each grid point, and then plotting the redefined variable. > > Any suggestions will be greatly appreciated. > > Kind regards > Greg Clarke > Support: OMAA -- Sacchetti Davide ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica V.le Brigare Partigiane 2 16121 Genova (I) tel: +39 010 6437535 mail: davide.sacchetti at arpal.org web: www.meteoliguria.it From ruthmota at CICESE.MX Fri Feb 9 05:14:19 2007 From: ruthmota at CICESE.MX (Ruth Cerezo Mota) Date: Fri, 9 Feb 2007 02:14:19 -0800 Subject: stripes in eps output In-Reply-To: <32FECA65-B6D4-4151-B1A0-3282EDB11140@washington.edu> Message-ID: Hi all! I use the Jasc Paint Shop Pro software to open and edit my GrADs graphics, and with this software I did not have any probles with stripes... cheers Ruth Cerezo-Mota ********************** University of Oxford Atmospheric, Oceanic & Planetary Physics Clarendon Laboratory Parks Road Oxford, OX1 3PU U.K. Tel (44)1865 272920 email: cerezo at atm.ox.ac.uk *************************** > My experience has also been that EPS files > generated from grads using > gxps look just fine when printed in journals -- > most of my paper have > GrADs images and I would not be using Grads if > the images looked bad > in print. I usually tweak and reformat graphics > in Illustrator, which > deals with the the lines just fine and makes > very nice raster files > if needed. > > The problem comes with looking at the EPS file > (eg in ghostview or > Preview) or when converting to a raster image. > Since most images end > up on the web as PNG, and not in print, this is > important. Grad's > printim does a decent job and is quick, but > without a lot of control > over image size and format. The solution is to > turn of anti aliasing. > Using +antialias in ImageMagick convert as > Jeffery Gall suggests > works reasonably well (convert actually calls > gs), but then lines are > rendered too thick and pixelated. If you specify > a fine dpi, then > scale down, you can get better line art while > eliminating the > striping. Compare these results on an eps > generated with gxps (note > I've added -trim to remove blank margins): > > convert -trim +antialias -density 150 test.eps > test.png > convert -trim +antialias -scale 25% -density 600 > test.eps test_scale.png > > > > > To do this, I use a grads script that calls a > csh script that can > handle any output format: > > Grads script saveimg in $GASCRP: > > function saveimg(file) > > * print screen to eps file specified as argument > > > * first print to gx file > > 'enable print tmp.gx' > 'print' > 'disable print' > > * now translate to image format > > '! gxconvert tmp.gx ' file > '! \rm -f tmp.gx' > > csh script in $PATH: > > #!/bin/csh > > if (! $#argv == 1) then > echo "usage: gxconvert file.gx file.ext where > ext indicates output > format" > exit > endif > > gxps -c -i $1 -o tmp.ps > /dev/null > > if ($2:e == 'ps') then > mv tmp.ps $2 > else if ($2:e == 'eps') then > ps2epsi tmp.ps $2 > \rm -f tmp.ps > else > convert +antialias -trim -scale 50% -density 300 > tmp.ps $2 > \rm -f tmp.ps > endif > > exit > > > -Eric > > -- > Eric Salath? > CSES Climate Impacts Group > > University of Washington > ~salathe> > 206-616-5351 > > On Feb 8, 2007, at 12:27 PM, Matthias Munnich > wrote: > >> -----BEGIN PGP SIGNED MESSAGE----- >> Hash: SHA1 >> >> The trouble is that one of the most common >> PostScript renderers, >> Ghostscript (gs), is nowadays doing >> antialiasing by default to soften >> ragged lines and corners. It blends >> neighboring bit even for purely >> horizontal lines. This leads to the >> pixel-width horizontal lines due >> to the way Grads is filling contours. >> Antialiasing can be turned of >> with "-dGraphicsAlphaBits=1" in gs or >> "-noantialias" in ghostview >> >> I still don't believe that hardcopy >> publications are likely to show >> lines. These lines would be close to >> invisible anyway on a 1200dpi >> printer . I suspect people are rendering >> their graphics themselves >> with tools like ImageMagick. Almost any free >> software is using >> Ghostscript as the rendering engine. These >> tools usually use default >> Ghostscript settings and you get the lines. >> >> Well, Jennifer's email suggests otherwise. >> Unfortunately you do not >> control how the graphic is rendered in print >> and submitting a rendered >> figure as suggested by Jennifer and Jeff is >> one solution. However, >> you loose print quality. >> >> My approach would be to look through a recent >> issues of JGR and see if >> the figures show these lines. I always >> recognize grads plots -- I >> guess I stared at too many by now. I am >> pretty sure it's fine to >> submit gxeps created EPS files to JGR. There >> were no lines in the >> figures of my recent GRL paper. >> >> >> ... Matt >> >> >> >> Jennifer Adams wrote: >>> Alan et al., Those horizontal stripes in the >>> ps/eps output are a >>> real nuisance. Potential sources of the >>> problem start with the >>> contouring algorithm in the GrADS graphics >>> layer, but the presence >>> of the stripes also seem to be influenced by >>> gxps, gxeps, any >>> external or internal utility that converts >>> the .eps or .ps files to >>> something else as they get incorporated into >>> written documents, and >>> the printers used by scientists and >>> publishers also seem to be a >>> factor. My recent experience is that it's >>> nearly impossible to >>> predict whether you'll get stripes or not. >>> It's seems more likely >>> that you'll have stripes than no stripes, and >>> the stripes also seem >>> to be most garish in plots with shaded >>> contours and white space. >>> >>> Here at COLA, we've begun using image format >>> files (.png and .gif) >>> instead of vector graphic formats, simply >>> because they're more >>> reliable even though they don't always look >>> as good. >>> >>> COLA users are being weaned from using PCs >>> running MS Windows, but >>> my sense is that the gv32 program allowed for >>> conversion of GrADS >>> metafiles into something else that did not >>> have stripes. I can't >>> verify that, perhaps the windows users out >>> there can confirm or >>> deny. >>> >>> Brian is planning to rework the contouring >>> algorithm, so there's >>> hope for the future; unfortunately, that hope >>> doesn't help in the >>> short term. I would see if you can get the >>> publisher to provide a >>> test print of your figures, and if that looks >>> bad, see if you can >>> use an image format instead. A .png at >>> 850x1100 resolution, shrunk >>> down to journal size, might not look so bad. >>> >>> Jennifer >>> >>> >>> >>> On Feb 8, 2007, at 11:03 AM, Alan Robock >>> wrote: >>> >>>> Dear Friends, >>>> >>>> Has anyone found a solution to the age-old >>>> problem of horizontal >>>> lines between bands of color in shaded grads >>>> output? We have eps >>>> files to submit for an accepted JGR paper >>>> and are concerned about >>>> how the final version will look. >>>> >>>> Is there anything we can do to the files >>>> before we submit them to >>>> fix the problem? >>>> >>>> Alan >>>> >>>> Alan Robock, Professor II Department of >>>> Environmental Sciences >>>> Phone: +1-732-932-9478 Rutgers University >>>> Fax: +1-732-932-8644 14 College Farm Road >>>> E-mail: robock at envsci.rutgers.edu >>>> New >>>> Brunswick, NJ 08901-8551 >>>> USA http://envsci.rutgers.edu/~robock >>>> >>> >>> -- Jennifer M. Adams IGES/COLA 4041 Powder >>> Mill Road, Suite 302 >>> Calverton, MD 20705 jma at cola.iges.org >>> >>> >>> >>> >> >> >> - -- >> - -------------------------------- >> Matthias Munnich >> Univ. of California, Los Angeles >> Dept. of Atmos. and Oceanic Sc. and >> Inst. of Geophysics and Planetary Physics >> 3845 Slichter Hall >> Los Angeles, CA 90095-1567 >> Phone: +1-310-794 5899 >> Fax: +1-310-206 3051 >> Email: munnich at atmos.ucla.edu >> - -------------------------------- >> -----BEGIN PGP SIGNATURE----- >> Version: GnuPG v1.4.2.2 (GNU/Linux) >> Comment: Using GnuPG with Mozilla - >> http://enigmail.mozdev.org >> >> iD8DBQFFy4e7idaDTh/cfhYRAkhyAKDYrNmXqc5O4m5L+DWrypw/kvsk4QCgn2rl >> Oo2eoszDJQ4TTyuFviUL6g0= >> =/obd >> -----END PGP SIGNATURE----- > From oluseun_idowu at YAHOO.COM Sat Feb 10 04:30:10 2007 From: oluseun_idowu at YAHOO.COM (oluseun idowu) Date: Sat, 10 Feb 2007 09:30:10 +0000 Subject: Help with Data Analysis In-Reply-To: <4579DE64.4000806@noaa.gov> Message-ID: Dear Users, Is there anyone who know what I should do for a solution to these questions: 1) I have a grib data with an array of 278 columns and 111 rows for 8 time steps of 3hrs interval (00,03,06, ..., 21) 2) I need to automatically extract the Maximum and Minimum data between (00 hr and 21hr) at 36 selected grid points from this data set. 3) The data consist of 2 years data. Your help will be appreciated. Thanks, Idowu --------------------------------- The all-new Yahoo! Mail goes wherever you go - free your email address from your Internet provider. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070210/720f574e/attachment.html From Chi.Nguyen at SCI.MONASH.EDU.AU Sat Feb 10 08:40:10 2007 From: Chi.Nguyen at SCI.MONASH.EDU.AU (Chi Nguyen) Date: Sun, 11 Feb 2007 00:40:10 +1100 Subject: Help with Data Analysis Message-ID: Hi Idowu, I have an idea which may help to solve your problem. 1) Make a mask field, with has value =1 at your 36 selected grid points, and =0 elsewhere. Say, its name is mask 2) Suppose the variable you need to analyse is F. Then define: F_new=mask*F so F_new will be equal F at that 36 points, and 0 elsewhere. 3) It will be easy to find max and min of F_new. In the case that F is all positive or negative, you will need to do two times: i.e. for finding max do as above, but for finding min F_new=-1*mask*F Hope that this makes sense to you. Good luck! Cheers, Mai oluseun idowu wrote:> > Dear Users, > > Is there anyone who know what I should do for a solution to these > questions: > > 1) I have a grib data with an array of 278 columns and 111 rows for 8 > time steps of 3hrs interval (00,03,06, ..., 21) > > 2) I need to automatically extract the Maximum and Minimum data between > (00 hr and 21hr) at 36 selected grid points from this data set. > 3) The data consist of 2 years data. > > Your help will be appreciated. > > Thanks, > Idowu > > > > > --------------------------------- > The all-new Yahoo! Mail goes wherever you go - free your email address > from your Internet provider. -- ********************************* Nguyen Chi Mai School of Mathematical Sciences R. 462, Bld. 28, Clayton Campus Monash University Work Ph. +61 3 9905 4424 From ivtoman at INET.HR Sat Feb 10 16:20:42 2007 From: ivtoman at INET.HR (Ivan) Date: Sat, 10 Feb 2007 22:20:42 +0100 Subject: Updated background map? Message-ID: Hi! Maybe I will ask something that exist on GrADS users list, but I cannot browse it, and that is my first question - how to browse it? If I go to http://list.cineca.it/cgi-bin/wa?A1=ind07&L=gradsusr there is required login, and I have not any password. How to get one?? OK, now real question... how to edit background map (lores, mres or hires are binary files)? Does exist an updated version? I need new political borders on Balkan region. Original maps have old Yugoslavia on them. Thanks! And last. Lot of meteo-related webistes has gfs-master (0.5?) charts *before* data become available on nomad3 dods server http://nomad3.ncep.noaa.gov:9090/dods/gfs_master/ How is that possible? Does that server recieves delayed data, that sites takes data from somewhere else or I missed something? Thanks again! Ivan From zenica3 at GMAIL.COM Sun Feb 11 08:12:30 2007 From: zenica3 at GMAIL.COM (ibrahim Hadzismailovic) Date: Sun, 11 Feb 2007 05:12:30 -0800 Subject: Updated background map? In-Reply-To: <45CE372A.3050900@inet.hr> Message-ID: Zdravo Ivane, Imam novu mapu, sa novim granicama. Javite mi se na moj mail pa cu Vam poslati. zenica3 at gmail.com Ibrahim, Sarajevo -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070211/c2cb0c59/attachment.html From jsg229 at PSU.EDU Mon Feb 12 12:53:33 2007 From: jsg229 at PSU.EDU (JEFFREY S GALL) Date: Mon, 12 Feb 2007 12:53:33 -0500 Subject: increase font size Message-ID: I want to increase the y-axis and x-axis tick label font size. How do I go about doing this? I would also like to increase the fontsize of contour labels. Any help in doing this would be much appreciated. Thanks, Jeff -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070212/1a49da44/attachment.html From Mary-Jo.Nath at NOAA.GOV Mon Feb 12 15:14:31 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Mon, 12 Feb 2007 15:14:31 -0500 Subject: increase font size In-Reply-To: <1171302811l.1286354l.0l@psu.edu> Message-ID: See the following in the documentation, Jeff: "set xlopts", "set ylopts", "set clopts". MJ JEFFREY S GALL wrote: > I want to increase the y-axis and x-axis tick label font size. How > do I go about doing this? > > I would also like to increase the fontsize of contour labels. Any > help in doing this would be much appreciated. > > Thanks, > > Jeff > > > > From camilo.uribe at GMAIL.COM Mon Feb 12 18:49:25 2007 From: camilo.uribe at GMAIL.COM (Camilo) Date: Mon, 12 Feb 2007 18:49:25 -0500 Subject: How Create Data Files? Message-ID: I've read the documentation about Creating Data Files[1] and try to compile the examples in gfortran but still can't figure out how to create the data files. How should I format the data? i.e. how much spacing (if any) is there between each column? A sample data file that grads can read and its explanation would be a huge help. Thanks for any help in advance. [1] http://www.iges.org/grads/gadoc/aboutgriddeddata.html From massimo.aceti at METEOGIORNALE.IT Tue Feb 13 04:53:43 2007 From: massimo.aceti at METEOGIORNALE.IT (Massimo Aceti) Date: Tue, 13 Feb 2007 10:53:43 +0100 Subject: Source for GENS datasets In-Reply-To: <45BA1663.2030905@noaa.gov> Message-ID: Wesley Ebisuzaki wrote: > NCO has http access to their FTP server. You can find the GENS data at > > http://www.ftp.ncep.noaa.gov/data/nccf/com/gens/prod/ > > > (( For those who aren't familiar with .inv files, they are wgrib > inventories of grib files. With the proper > software, one can randomly download selected records from a http > server. That can save huge amounts > of download time. See > http://www.cpc.ncep.noaa.gov/products/wesley/fast_downloading_grib.html)) > > Wesley > fast downloading grib is an excellent way to download grib files from ncep server in a few time, but, there aren't inv (or idx) files in the gens grib2 directories (http://www.ftp.ncep.noaa.gov/data/nccf/com/gens/prod/gefs.20070213/00/pgrb2a/) and in the grib1 directory there are only 10 ens perturbations How can I download selected records with get_inv.pl and get_grib.pl from gens grib2 dataset? Thanks Max From Wesley.Ebisuzaki at NOAA.GOV Tue Feb 13 10:01:17 2007 From: Wesley.Ebisuzaki at NOAA.GOV (Wesley Ebisuzaki) Date: Tue, 13 Feb 2007 10:01:17 -0500 Subject: Source for GENS datasets In-Reply-To: <45D18AA7.6040509@meteogiornale.it> Message-ID: Max, NCO's policy is to generate the .idx files upon request and I've heard that they have been very good to helping out. (Bandwidth costs money and .idx files can save a lot of bandwidth.) I've sent the point of contact by email because people move around at NCO. Wesley Massimo Aceti wrote: > Wesley Ebisuzaki wrote: > >> NCO has http access to their FTP server. You can find the GENS >> data at >> >> http://www.ftp.ncep.noaa.gov/data/nccf/com/gens/prod/ >> > >> >> (( For those who aren't familiar with .inv files, they are wgrib >> inventories of grib files. With the proper >> software, one can randomly download selected records from a http >> server. That can save huge amounts >> of download time. See >> http://www.cpc.ncep.noaa.gov/products/wesley/fast_downloading_grib.html)) >> >> >> Wesley >> > > fast downloading grib is an excellent way to download grib files from > ncep server in a few time, but, there aren't inv (or idx) files in the > gens grib2 directories > (http://www.ftp.ncep.noaa.gov/data/nccf/com/gens/prod/gefs.20070213/00/pgrb2a/) > > and in the grib1 directory there are only 10 ens perturbations > > How can I download selected records with get_inv.pl and get_grib.pl from > gens grib2 dataset? > > Thanks > > Max From kahill at NCSU.EDU Tue Feb 13 16:35:20 2007 From: kahill at NCSU.EDU (Kevin Hill) Date: Tue, 13 Feb 2007 16:35:20 -0500 Subject: Making plots with separate y-axes Message-ID: Hello everyone, I am trying to plot time-series of minimum sea level pressure for two different simulations, along with a line showing the difference between the two. I would like the minimum pressure values to correspond to one set of y-axis values (say 900 - 1000 hPa, shown on the left side of the plot), and the line showing the difference in sea level pressure to correspond to another set of y-axis values (say -10 to 10 hPa, shown on the right side of the plot). So far I have not encountered any way to do this type of plot with grads; can any of you provide assistance in how to do this? I know I can do it using other programs (even Excel), but would like to do this with grads if possible. Thank you, Kevin From fengwh at ENV.LEEDS.AC.UK Wed Feb 14 04:17:04 2007 From: fengwh at ENV.LEEDS.AC.UK (Wuhu Feng) Date: Wed, 14 Feb 2007 09:17:04 +0000 Subject: Making plots with separate y-axes In-Reply-To: <45D22F18.7040509@ncsu.edu> Message-ID: Hi Kevin, Try to use 'set ylpos 0 r' and 'set ylpos 0 l' Cheers. Wuhu -------------------------------------- Dr. Wuhu Feng Institute for Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT UK Tel:(44) 0113 343 6766 Fax:(44) 0113 343 6716 E-mail:fengwh at env.leeds.ac.uk Web: http://www.env.leeds.ac.uk/~fengwh -------------------------------------- On Tue, 13 Feb 2007, Kevin Hill wrote: > Hello everyone, > > I am trying to plot time-series of minimum sea level pressure for two > different simulations, along with a line showing the difference between > the two. I would like the minimum pressure values to correspond to one > set of y-axis values (say 900 - 1000 hPa, shown on the left side of the > plot), and the line showing the difference in sea level pressure to > correspond to another set of y-axis values (say -10 to 10 hPa, shown on > the right side of the plot). So far I have not encountered any way to do > this type of plot with grads; can any of you provide assistance in how > to do this? I know I can do it using other programs (even Excel), but > would like to do this with grads if possible. > > Thank you, > > Kevin > From preuter at LABRI.FR Wed Feb 14 14:14:48 2007 From: preuter at LABRI.FR (Patrick Reuter) Date: Wed, 14 Feb 2007 20:14:48 +0100 Subject: Updated background map? In-Reply-To: <45CE372A.3050900@inet.hr> Message-ID: Hi there, I cannot browse the archive neither. > And last. Lot of meteo-related webistes has gfs-master (0.5?) charts > *before* data become available on nomad3 dods server > http://nomad3.ncep.noaa.gov:9090/dods/gfs_master/ How is that possible? > Does that server recieves delayed data, that sites takes data from > somewhere else or I missed something? Did you get any insight or answer in this ? I would be really interested ! Cheers Patrick > > Thanks again! > > Ivan From ivtoman at INET.HR Wed Feb 14 14:56:17 2007 From: ivtoman at INET.HR (Ivan Toman) Date: Wed, 14 Feb 2007 20:56:17 +0100 Subject: Updated background map? In-Reply-To: <20070214201448.v90i8scz4okw0og4@webmail.labri.fr> Message-ID: Hi Patrick! I managed how to get password. Go to http://list.cineca.it/cgi-bin/wa?GETPW1=A1=ind07%26L=gradsusr and you need to register your pass. Type in e-mail and choose password, after that you will recieve e-mail with instructions how to proceed. I cannot believe I didn't see that before asking here (sorry people). About data... Yes. Data od nomads servers is delayed, and nomads generally are *very* overloaded machines, so they become unreliable more and more with time. ftp2u protocol is somewhat better than dods, it seems. But, if you need really reliable solution, you will have to download whole data (in GRIB2 format) from other locations than nomads. See for example: http://www.cpc.ncep.noaa.gov/products/wesley/get_gfs.html I don't know better solution (*), but also looking for it! (*) = small download, only for limited area of interest, without much delay and missing data Ivan Patrick Reuter wrote: > Hi there, > > I cannot browse the archive neither. > >> And last. Lot of meteo-related webistes has gfs-master (0.5?) charts >> *before* data become available on nomad3 dods server >> http://nomad3.ncep.noaa.gov:9090/dods/gfs_master/ How is that possible? >> Does that server recieves delayed data, that sites takes data from >> somewhere else or I missed something? > > Did you get any insight or answer in this ? I would be really > interested ! > > Cheers > Patrick > > >> >> Thanks again! >> >> Ivan > > From preuter at LABRI.FR Wed Feb 14 17:44:24 2007 From: preuter at LABRI.FR (Patrick Reuter) Date: Wed, 14 Feb 2007 23:44:24 +0100 Subject: Updated background map? In-Reply-To: <45D36961.3030700@inet.hr> Message-ID: Ok thanx! Get_gfs uses the following adress : # $URL='http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.$YYYY$MM$DD$HH/gfs.t${HH}z.pgrb2f${FHR}'; # # grib1 files from www.ftp.ncep.noaa.gov $URL='http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.$YYYY$MM$DD$HH/gfs.t${HH}z.pgrbf${FHR}'; But this seems to be 1x1 resolution data, and not 0.5 0.5. It is hard to get limited area data, I always download the whole world and throw away what I don t need ... Patrick Ivan Toman a ????crit??? : > Hi Patrick! > > I managed how to get password. Go to > http://list.cineca.it/cgi-bin/wa?GETPW1=A1=ind07%26L=gradsusr > > and you need to register your pass. Type in e-mail and choose password, > after that you will recieve e-mail with instructions how to proceed. I > cannot believe I didn't see that before asking here (sorry people). > > About data... Yes. Data od nomads servers is delayed, and nomads > generally are *very* overloaded machines, so they become unreliable more > and more with time. ftp2u protocol is somewhat better than dods, it > seems. But, if you need really reliable solution, you will have to > download whole data (in GRIB2 format) from other locations than nomads. > See for example: > http://www.cpc.ncep.noaa.gov/products/wesley/get_gfs.html > > I don't know better solution (*), but also looking for it! > > > (*) = small download, only for limited area of interest, without much > delay and missing data > > Ivan > > > Patrick Reuter wrote: >> Hi there, >> >> I cannot browse the archive neither. >> >>> And last. Lot of meteo-related webistes has gfs-master (0.5??) charts >>> *before* data become available on nomad3 dods server >>> http://nomad3.ncep.noaa.gov:9090/dods/gfs_master/ How is that possible? >>> Does that server recieves delayed data, that sites takes data from >>> somewhere else or I missed something? >> >> Did you get any insight or answer in this ? I would be really >> interested ! >> >> Cheers >> Patrick >> >> >>> >>> Thanks again! >>> >>> Ivan >> >> From John.Papineau at NOAA.GOV Wed Feb 14 18:03:23 2007 From: John.Papineau at NOAA.GOV (John Papineau) Date: Wed, 14 Feb 2007 14:03:23 -0900 Subject: defval ?? In-Reply-To: <20070214201448.v90i8scz4okw0og4@webmail.labri.fr> Message-ID: Greetings GRADS warriors- 1) I am attempting to write a script using 'defval'. Specifically, I want to _retrieve a value like U wind_ (ugrdprs in my nomenclature) at a point. If I run: 'q defval ugrdprs 'i' 'j the following error message results: "defined name not found" I have tried all possible combinations of ' around the variables, etc Have also tried specific grid values like 'q defval ugrdprs '30' '25 same error message. 2) It is unclear from the documentation if I have ## number of files (each at a different time) that I want to query for specific values, such as noted above, do I need to use the "read" command? or just have them all open? 3) What does "rc" stand for? record code, if so what does it mean? Many thanks- John From davide.sacchetti at ARPAL.ORG Thu Feb 15 03:22:47 2007 From: davide.sacchetti at ARPAL.ORG (Davide Sacchetti) Date: Thu, 15 Feb 2007 09:22:47 +0100 Subject: defval ?? In-Reply-To: <45D3953B.5020401@noaa.gov> Message-ID: warriors answer ... 1) q defval works for DEFINED variables, NOT for variables described in ctl file ... try: define uuu=ugrdprs and apply your syntax to uuu, it should work 2) on my opinion the documentation is clear enough. Anyway if you have more than one data file you could describe all your data by a single ctl file using the template option. 3) The "read" command is not used for data files but for ascii files: it returns a return code (you may wish calling it "rc") and the content of what you have read. Hoping to help a bit ... Bye bye Davide On Wed, 2007-02-14 at 14:03 -0900, John Papineau wrote: > Greetings GRADS warriors- > > 1) I am attempting to write a script using 'defval'. > Specifically, I want to _retrieve a value like U wind_ (ugrdprs in > my nomenclature) at a point. If I run: > > 'q defval ugrdprs 'i' 'j > > the following error message results: > > "defined name not found" > I > have tried all possible combinations of ' around the variables, etc > Have also tried specific grid values like 'q defval ugrdprs '30' '25 > same error message. > > 2) It is unclear from the documentation if I have ## number of files > (each at a > different time) that I want to query for specific values, such as noted > above, do I need to use the "read" command? or just have them all open? > > 3) What does "rc" stand for? record code, if so what does it mean? > > Many thanks- John -- Sacchetti Davide ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica V.le Brigare Partigiane 2 16121 Genova (I) tel: +39 010 6437535 mail: davide.sacchetti at arpal.org web: www.meteoliguria.it From bernd.becker at METOFFICE.GOV.UK Thu Feb 15 04:28:29 2007 From: bernd.becker at METOFFICE.GOV.UK (bernd.becker at METOFFICE.GOV.UK) Date: Thu, 15 Feb 2007 09:28:29 +0000 Subject: gribmap binary from 1.9 distribution AIX Message-ID: Hello, Who has a gribmap binary from the Grads 1.9 distribution that runs on AIX 5 PowerPc ? Can I have a copy, please? Thanks! -- 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.gov.uk - http://www.metoffice.gov.uk From mcroke at AIRDAT.COM Thu Feb 15 11:27:29 2007 From: mcroke at AIRDAT.COM (Meredith Croke) Date: Thu, 15 Feb 2007 17:27:29 +0100 Subject: dat file does not exit Message-ID: Hello Everyone, I am trying to produce a .dat file of 15 days worth of data, each with forecast hours out to 384. The problem I am running into is, I am missing a few of these files due to server problems, when the script can't find the data file I get a segmentation fault and grads closes. Is there a way I can write in my script after the open command that if the file does not exist, skip and move onto the next forecast hour? I tried creating fake data files with nothing in them in the hopes I would get either some kindof text output or 9999 in my tab delimited file but I just got an I/O error and the script stopped running. Thanks in advance for your help. Meredith From davide.sacchetti at ARPAL.ORG Thu Feb 15 11:39:43 2007 From: davide.sacchetti at ARPAL.ORG (Davide Sacchetti) Date: Thu, 15 Feb 2007 17:39:43 +0100 Subject: dat file does not exit In-Reply-To: <20070215162813.E7E6F22F56@mx2.cineca.it> Message-ID: you could use what "open" returns (variable "result", available only in a script context) Bye bye Davide On Thu, 2007-02-15 at 17:27 +0100, Meredith Croke wrote: > Hello Everyone, > > I am trying to produce a .dat file of 15 days worth of data, each with > forecast hours out to 384. The problem I am running into is, I am missing a > few of these files due to server problems, when the script can't find the > data file I get a segmentation fault and grads closes. Is there a way I can > write in my script after the open command that if the file does not exist, > skip and move onto the next forecast hour? I tried creating fake data files > with nothing in them in the hopes I would get either some kindof text output > or 9999 in my tab delimited file but I just got an I/O error and the script > stopped running. > > Thanks in advance for your help. > > Meredith -- Sacchetti Davide ARPAL UO3 Centro Meteo Idrologico Regione Liguria - Dir. Scientifica V.le Brigare Partigiane 2 16121 Genova (I) tel: +39 010 6437535 mail: davide.sacchetti at arpal.org web: www.meteoliguria.it From mlsilver at NCSU.EDU Thu Feb 15 12:13:52 2007 From: mlsilver at NCSU.EDU (Morgan Silverman) Date: Thu, 15 Feb 2007 18:13:52 +0100 Subject: 4panel plot from separate .ctl files Message-ID: Hello, I am trying to create a 4 panel plot where each panel shows data from a different .ctl file. Does anyone know if this is possible and how? I was thinking there maybe a way to create the first panel and then append that file with the second, third, and fourth panels, but I'm not sure if that would work. Thanks. -Morgan Silverman From gwest at MET.UTAH.EDU Thu Feb 15 12:35:38 2007 From: gwest at MET.UTAH.EDU (Greg West) Date: Thu, 15 Feb 2007 10:35:38 -0700 Subject: 4panel plot from separate .ctl files In-Reply-To: <20070215171436.63F2820A99@mx2.cineca.it> Message-ID: Morgan, Opening multiple control files at once might solve your problem. To do this, simply open another ctl file without closing the other ones first. To refer to variables in the 2nd ctl file you open, you'd refer to the variable as variable_name.2, ..and for the 3rd ctl file you open, variable_name.3, etc. Greg Morgan Silverman wrote: > Hello, > I am trying to create a 4 panel plot where each panel shows data from a > different .ctl file. Does anyone know if this is possible and how? I was > thinking there maybe a way to create the first panel and then append that > file with the second, third, and fourth panels, but I'm not sure if that > would work. Thanks. > -Morgan Silverman -- ****************************************************** Gregory West Graduate Research Assistant University of Utah - Dept of Meteorology Office: INSCC 480-10 135 S 1460 E WBB Rm. 819 Salt Lake City, UT 84112-0110 801-585-1410 gwest at met.utah.edu ****************************************************** From eduardo.agosta at GMAIL.COM Thu Feb 15 15:01:38 2007 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Thu, 15 Feb 2007 17:01:38 -0300 Subject: READING BINARY FILES Message-ID: Hi you all! I am having problems reading binary files using a *.ctl file. I have created the binary file using Matlab after preparing the data on a mx1 matrix (in fact a one-column vector containing the equally spaced lon-lat data). Is there any specific flag that I may use in the OPTIONS comand within the ctl file? Below I copy the ctl file I am using: DSET ^ozojun TITLE TOC anomaly UNDEF -99999 OPTIONS SEQUENTIAL XDEF 288 linear 0 1.25 YDEF 31 linear -60 1.0 ZDEF 1 levels 50 TDEF 1 linear 1jan1950 1yr VARS 1 o3 0 99 TOC anomaly ENDVARS The answer is "Open error: Can`t open binary data file". Any help? Thanks in advance. Eduardo. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070215/e12a3e25/attachment.html From Mary-Jo.Nath at NOAA.GOV Thu Feb 15 16:33:48 2007 From: Mary-Jo.Nath at NOAA.GOV (Mary Jo Nath) Date: Thu, 15 Feb 2007 16:33:48 -0500 Subject: READING BINARY FILES In-Reply-To: <3367a2e40702151201t5fb91bf9tfca3949ef8afac74@mail.gmail.com> Message-ID: Eduardo, It appears that Grads isn't finding your data file. Make sure the file named "ozojun" is in the same directory/folder as the *.ctl file. Otherwise put the file's path name on the "DSET" line to indicate where "ozojun" resides. MJ Eduardo Agosta Scarel wrote: > Hi you all! > > I am having problems reading binary files using a *.ctl file. I have > created the binary file using Matlab after preparing the data on a mx1 > matrix (in fact a one-column vector containing the equally spaced > lon-lat data). Is there any specific flag that I may use in the > OPTIONS comand within the ctl file? > Below I copy the ctl file I am using: > > DSET ^ozojun > TITLE TOC anomaly > UNDEF -99999 > OPTIONS SEQUENTIAL > XDEF 288 linear 0 1.25 > YDEF 31 linear -60 1.0 > ZDEF 1 levels 50 > TDEF 1 linear 1jan1950 1yr > VARS 1 > o3 0 99 TOC anomaly > ENDVARS > > The answer is "Open error: Can`t open binary data file". > > Any help? > Thanks in advance. > > Eduardo. > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070215/6ca1661a/attachment.html From eduardo.agosta at GMAIL.COM Thu Feb 15 17:20:10 2007 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Thu, 15 Feb 2007 19:20:10 -0300 Subject: READING BINARY FILES In-Reply-To: <45D4D1BC.8050004@noaa.gov> Message-ID: Mary Jo, thanks for your suggestion, but even it doesn?t work. I will try to move all the folder to a new one and try. Eduardo. 2007/2/15, Mary Jo Nath : > > Eduardo, > > It appears that Grads isn't finding your data file. Make sure the file > named "ozojun" is in the same directory/folder as the *.ctl file. > Otherwise put the file's path name on the "DSET" line to indicate where > "ozojun" resides. > > MJ > > Eduardo Agosta Scarel wrote: > > Hi you all! > > I am having problems reading binary files using a *.ctl file. I have > created the binary file using Matlab after preparing the data on a mx1 > matrix (in fact a one-column vector containing the equally spaced lon-lat > data). Is there any specific flag that I may use in the OPTIONS comand > within the ctl file? > Below I copy the ctl file I am using: > > DSET ^ozojun > TITLE TOC anomaly > UNDEF -99999 > OPTIONS SEQUENTIAL > XDEF 288 linear 0 1.25 > YDEF 31 linear -60 1.0 > ZDEF 1 levels 50 > TDEF 1 linear 1jan1950 1yr > VARS 1 > o3 0 99 TOC anomaly > ENDVARS > > The answer is "Open error: Can`t open binary data file". > > Any help? > Thanks in advance. > > Eduardo. > > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070215/a02ac921/attachment.html From polarlow at GMAIL.COM Thu Feb 15 19:09:06 2007 From: polarlow at GMAIL.COM (Robert Reed) Date: Thu, 15 Feb 2007 20:09:06 -0400 Subject: PC Grads path issue Message-ID: Hi, I`ve just installed PCGrads to my windowsXP pro edition, with a view to learning the grads basic commands and scripting language. My first problem, and somewhat embarassing is regarding the install structure. For example I want to place my grads scripts in a `scripts` folder within the main directory such as: C:\Program Files\PCGrADS\scripts How would I execute the scripts within that folder? I`ve been trying for the last twenty minutes with various path formats, but with no success. Thanks. From axl419 at YAHOO.COM Thu Feb 15 19:56:51 2007 From: axl419 at YAHOO.COM (See Hai Ooi) Date: Fri, 16 Feb 2007 00:56:51 +0000 Subject: PC Grads path issue In-Reply-To: Message-ID: Hi, Your problem may be due to the spacing between "Program" and "Files". If this is the case, suggest you do the following:- 1. create a subdirectory such as "work" for you to keep all your working files,data etc. 2. then try this within GrADS: run c:\work\scripts\xyz.gs Hope this helps. Ooi --- Robert Reed wrote: > Hi, > > I`ve just installed PCGrads to my windowsXP pro > edition, with a view > to learning the grads basic commands and scripting > language. > > My first problem, and somewhat embarassing is > regarding the install > structure. For example I want to place my grads > scripts in a `scripts` > folder within the main directory such as: > > C:\Program Files\PCGrADS\scripts > > How would I execute the scripts within that folder? > > I`ve been trying for the last twenty minutes with > various path > formats, but with no success. > > Thanks. > Send instant messages to your online friends http://uk.messenger.yahoo.com From polarlow at GMAIL.COM Thu Feb 15 20:43:07 2007 From: polarlow at GMAIL.COM (Robert Reed) Date: Thu, 15 Feb 2007 21:43:07 -0400 Subject: PC Grads path issue In-Reply-To: <513862.76281.qm@web31003.mail.mud.yahoo.com> Message-ID: Thank you, that worked :) On 2/15/07, See Hai Ooi wrote: > Hi, > > Your problem may be due to the spacing between > "Program" and "Files". If this is the case, suggest > you do the following:- > 1. create a subdirectory such as "work" for you to > keep all your working files,data etc. > 2. then try this within GrADS: > > run c:\work\scripts\xyz.gs > > Hope this helps. > > Ooi > > --- Robert Reed wrote: > > > Hi, > > > > I`ve just installed PCGrads to my windowsXP pro > > edition, with a view > > to learning the grads basic commands and scripting > > language. > > > > My first problem, and somewhat embarassing is > > regarding the install > > structure. For example I want to place my grads > > scripts in a `scripts` > > folder within the main directory such as: > > > > C:\Program Files\PCGrADS\scripts > > > > How would I execute the scripts within that folder? > > > > I`ve been trying for the last twenty minutes with > > various path > > formats, but with no success. > > > > Thanks. > > > > > Send instant messages to your online friends http://uk.messenger.yahoo.com > From paul.earnshaw at METOFFICE.GOV.UK Tue Feb 20 05:35:53 2007 From: paul.earnshaw at METOFFICE.GOV.UK (Earnshaw Paul) Date: Tue, 20 Feb 2007 10:35:53 +0000 Subject: Problems reading NetCDF files In-Reply-To: <49124.163.1.242.169.1171016059.squirrel@arrecife.cicese.mx> Message-ID: Hi all, Just a quick query related to reading NetCDF files. I have attached two (small) files to this message to illustrate the point I am making. When I open the two files I get following response from GrADS. ------------------------------------------------------------------------- Grid Analysis and Display System (GrADS) Version 1.9b4 Copyright (c) 1988-2005 by Brian Doty and IGES Center for Ocean-Land-Atmosphere Studies (COLA) Institute for Global Environment and Society (IGES) GrADS comes with ABSOLUTELY NO WARRANTY See file COPYRIGHT for more information Config: v1.9b4 32-bit little-endian readline sdf/xdf netcdf lats athena printim Issue 'q config' command for more information. GX Package Initialization: Size = 11 8.5 ga-> sdfopen test_1.nc Scanning self-describing file: test_1.nc SDF file test_1.nc is open as file 1 LON set to 1 1 LAT set to 1 1 LEV set to 97231.4 97231.4 Time values set: 1997:6:23:0 1997:6:23:0 ga-> sdfopen test_2.nc Scanning self-describing file: test_2.nc Error parsing time units in SDF file. ga-> ------------------------------------------------------------------------- Now, as far as I can tell there is only one difference between the two files. The units attribute for the dimension "pres" in test_1.nc is "Pa" and in test_2.nc is "pa". Should this difference really prevent GrADS from being able to read the file and get the error message above? Cheers, Paul -------------- next part -------------- A non-text attachment was scrubbed... Name: test_1.nc Type: application/x-netcdf Size: 1412 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070220/f0f1e846/attachment.nc -------------- next part -------------- A non-text attachment was scrubbed... Name: test_2.nc Type: application/x-netcdf Size: 1412 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070220/f0f1e846/attachment-0001.nc From lalvarez at IDEMOC.COM Wed Feb 21 07:47:44 2007 From: lalvarez at IDEMOC.COM (Leopoldo Alvarez) Date: Wed, 21 Feb 2007 13:47:44 +0100 Subject: Smooth the variables in model output Message-ID: Dear users GrADS, I need to smooth the variables that i obtain with MM5 model output. Due to the presence of terrain, etc... I have much noise in the temperature, humidity, geopotencial variables... I have used ?set csmooth? without success in grads. regards, Leopoldo Alvarez - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Leopoldo Alvarez Adviser - Tecnico Investigator Meteorology and Oceanographic Studies Institute Canary Islands - Spain E-mail: lalvarez at idemoc.com Tlfno: (+34) 600 253 237 - - - - - - - - - - - - - - - - - - - - - - - - - - - From hallak at MODEL.IAG.USP.BR Wed Feb 21 08:46:17 2007 From: hallak at MODEL.IAG.USP.BR (Ricardo Hallak) Date: Wed, 21 Feb 2007 10:46:17 -0300 Subject: Smooth the variables in model output In-Reply-To: <20070221124829.8D602224E7@mx2.cineca.it> Message-ID: Dear Leopoldo, try: "d smth9(var)" or add more power in this function: "d smth9(smth9(var))" or yet: "d smth9(smth9(...(smth9(var))...))" Regards, Ricardo On Wed, 21 Feb 2007 13:47:44 +0100, Leopoldo Alvarez wrote > Dear users GrADS, > > I need to smooth the variables that i obtain with MM5 model output. > Due to the presence of terrain, etc... I have much noise in the > temperature, humidity, geopotencial variables... > > I have used “set csmooth” without success in grads. > > regards, > Leopoldo Alvarez > - - > - - - - - - - - - - - - - - - - - - - - - - - - - - - > Leopoldo Alvarez > Adviser - Tecnico Investigator > Meteorology and Oceanographic Studies Institute > Canary Islands - Spain > E-mail: lalvarez at idemoc.com > Tlfno: (+34) 600 253 237 > - - - - - - - - - - - - - - - - - - - - - - - - - - - Hallak Laborat?rio de Hidrometeorologia (LABHIDRO) Departamento de Ci?ncias Atmosf?ricas (DCA) Instituto de Astronomia, Geof?sica e Ci?ncias Atmosf?ricas (IAG) Universidade de S?o Paulo (USP) Rua do Mat?o, 1226 - Cidade Universit?ria -S?o Paulo SP -05508-900 Fone.: +55 (0-XX-11) 3091-4735 Fax: +55 (0-XX-11) 3091-4714 From fuentes at POSTIT.GSFC.NASA.GOV Wed Feb 21 11:40:00 2007 From: fuentes at POSTIT.GSFC.NASA.GOV (M Fuentes) Date: Wed, 21 Feb 2007 11:40:00 -0500 Subject: Create a 3x3 box around a tropical cyclone Message-ID: Dear GrADS Users: I need to create a 3x3 box around the center of a tropical cyclone to make some kinetic energy calculations. I will like to know if someone has some ideas in how to create that box. Thanks -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean. From arabonis at EGS.UCT.AC.ZA Wed Feb 21 22:58:55 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 22 Feb 2007 05:58:55 +0200 Subject: Monsoondata Message-ID: Hi Seems the whole of Monsoondata has gone down, anyone else experiancing the problem? The gds cannot be accesses and the website gives the error below Many Thanks Jean Pierre Network Error (tcp_error) A communication error occurred: "Operation timed out" The Web Server may be down, too busy, or experiencing other problems preventing it from responding to requests. You may wish to try again at a later time. For assistance, contact your network support team. -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070222/d0aeecf5/attachment.html From jma at COLA.IGES.ORG Thu Feb 22 07:48:37 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 22 Feb 2007 07:48:37 -0500 Subject: Monsoondata In-Reply-To: <45DD14FF.9090300@egs.uct.ac.za> Message-ID: Yes, we know there's a problem -- this time it's a network issue that we'll fix as soon as we can. --Jennifer On Feb 21, 2007, at 10:58 PM, Jean Pierre Arabonis wrote: > Hi > Seems the whole of Monsoondata has gone down, anyone else > experiancing the problem? The gds cannot be accesses and the > website gives the error below > > Many Thanks > Jean Pierre > > > Network Error (tcp_error) > > A communication error occurred: "Operation timed out" The Web > Server may be down, too busy, or experiencing other problems > preventing it from responding to requests. You may wish to try > again at a later time. > For assistance, contact your network support team. > > > > > -- > Jean Pierre Arabonis > arabonis at egs.uct.ac.za > Tel 021 780 1021 cell 084 401 1365 > This email is 100% Microsoft Free -- 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/20070222/e0a6b3e3/attachment.html From Charles.Seman at NOAA.GOV Thu Feb 22 10:28:37 2007 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Thu, 22 Feb 2007 10:28:37 -0500 Subject: Create a 3x3 box around a tropical cyclone In-Reply-To: <45DC75E0.9040102@postit.gsfc.nasa.gov> Message-ID: M Fuentes, Please find attached a script plot_hurricane_center2c.gs which defines a box around the minimum surface pressure associated with a "hurricane"... I hope this helps. Let me know if you have any questions about it, OK? Thanks, Chuck M Fuentes wrote: > Dear GrADS Users: > I need to create a 3x3 box around the center of a tropical cyclone to > make some kinetic energy calculations. I will like to know if someone > has some ideas in how to create that box. > Thanks > > -- > This message has been scanned for viruses and > dangerous content by MailScanner, and is > believed to be clean. -- 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 reflect any position of the Government or NOAA." -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: plot_hurricane_center2c.gs Url: http://gradsusr.org/pipermail/gradsusr/attachments/20070222/ce2af4bd/attachment.pl From jma at COLA.IGES.ORG Thu Feb 22 10:28:37 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 22 Feb 2007 10:28:37 -0500 Subject: Monsoondata In-Reply-To: Message-ID: monsoondata.org is up again. --Jennifer From arabonis at EGS.UCT.AC.ZA Thu Feb 22 14:10:11 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 22 Feb 2007 21:10:11 +0200 Subject: Monsoondata In-Reply-To: <77E4FC97-0DC5-49D7-8876-EB4A24E72E9D@cola.iges.org> Message-ID: Thanks Jennifer, much appreciated. JP Jennifer Adams wrote: > monsoondata.org is up again. --Jennifer > -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free From arabonis at EGS.UCT.AC.ZA Thu Feb 22 14:19:49 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 22 Feb 2007 21:19:49 +0200 Subject: Monsoondata - GDS In-Reply-To: <77E4FC97-0DC5-49D7-8876-EB4A24E72E9D@cola.iges.org> Message-ID: Sorry to bring bad news ..... seems the GDS on monsoondata is still not working. Cheers JP Jennifer Adams wrote: > monsoondata.org is up again. --Jennifer > -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free From jma at COLA.IGES.ORG Thu Feb 22 14:23:17 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 22 Feb 2007 14:23:17 -0500 Subject: Monsoondata - GDS In-Reply-To: <45DDECD5.9010903@egs.uct.ac.za> Message-ID: ok, now try... On Feb 22, 2007, at 2:19 PM, Jean Pierre Arabonis wrote: > Sorry to bring bad news ..... seems the GDS on monsoondata is still > not > working. > Cheers > JP > > Jennifer Adams wrote: > >> monsoondata.org is up again. --Jennifer >> > > -- > Jean Pierre Arabonis > arabonis at egs.uct.ac.za > Tel 021 780 1021 cell 084 401 1365 > This email is 100% Microsoft Free > -- 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/20070222/570ae0b9/attachment.html From lalvarez at IDEMOC.COM Thu Feb 22 14:52:00 2007 From: lalvarez at IDEMOC.COM (Leopoldo Alvarez) Date: Thu, 22 Feb 2007 20:52:00 +0100 Subject: Smooth the variables in model output Message-ID: Dear Ricardo, Thanks for this method. I need another form to smooth much more, for example, when i attempt to smooth the geopotencial 850 hPa that crosses a strong relief. regards, Leopoldo From arabonis at EGS.UCT.AC.ZA Thu Feb 22 16:02:28 2007 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 22 Feb 2007 23:02:28 +0200 Subject: Monsoondata - GDS - good news In-Reply-To: Message-ID: All great now, again lots of thanks. Cheers JP Jennifer Adams wrote: > ok, now try... > > On Feb 22, 2007, at 2:19 PM, Jean Pierre Arabonis wrote: > >> Sorry to bring bad news ..... seems the GDS on monsoondata is still not >> working. >> Cheers >> JP >> >> Jennifer Adams wrote: >> >>> monsoondata.org is up again. --Jennifer >>> >> >> -- >> Jean Pierre Arabonis >> arabonis at egs.uct.ac.za >> Tel 021 780 1021 cell 084 401 1365 >> This email is 100% Microsoft Free >> > > -- > Jennifer M. Adams > IGES/COLA > 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705 > jma at cola.iges.org > > > -- Jean Pierre Arabonis arabonis at egs.uct.ac.za Tel 021 780 1021 cell 084 401 1365 This email is 100% Microsoft Free -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070222/a6763641/attachment.html From yitbek at YAHOO.COM Thu Feb 22 16:34:29 2007 From: yitbek at YAHOO.COM (Yitbarek Kifle) Date: Thu, 22 Feb 2007 22:34:29 +0100 Subject: Binary to Ascii conversion Message-ID: Dear all, As a new and inexperienced user of Grads,i have come with a problem of converting the Grads time series binary output to Ascii format.I have read almost all such problems on the users group list but that doesn't seem to help me. Here is the discription of what i did, 'set gxout fwrite' 'set fwrite yiki.dat' set t 613 697 d tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5)) disable fwrite i have the file Yiki.dat. Would somebody please tell me what to do next in step by step basis,as I do not have a programming skill? thanks in advance. From Diane.Stokes at NOAA.GOV Thu Feb 22 17:05:58 2007 From: Diane.Stokes at NOAA.GOV (Diane Stokes) Date: Thu, 22 Feb 2007 17:05:58 -0500 Subject: Binary to Ascii conversion In-Reply-To: <20070222213513.8486A20F9E@mx2.cineca.it> Message-ID: If you don't need the data in binary (which is what you have in file yiki.dat), you can skip the "gxout fwrite" lines and use "gxout print" instead. Eg, in a grads script: 'set gxout print' 'set t 613 697' 'display tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' rc=write(yiki.txt,result) say 'return code from write is: 'rc rc=close(yiki.txt) say 'return code from close is: 'rc In the above, the output from the display command gets stored in variable "result" and you write the contents of "result" to ascii file yiki.txt. See the doc for 'set prnopts' at http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html if you want to change the format of the text output. Diane Yitbarek Kifle wrote: > Dear all, > > As a new and inexperienced user of Grads,i have come with a problem of > converting the Grads time series binary output to Ascii format.I have read > almost all such problems on the users group list but that doesn't seem to > help me. Here is the discription of what i did, > > 'set gxout fwrite' > 'set fwrite yiki.dat' > set t 613 697 > d tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5)) > disable fwrite > > i have the file Yiki.dat. > > Would somebody please tell me what to do next in step by step basis,as I do > not have a programming skill? > > thanks in advance. > From nakazawa at MRI-JMA.GO.JP Fri Feb 23 00:13:24 2007 From: nakazawa at MRI-JMA.GO.JP (Tetsuo Nakazawa) Date: Fri, 23 Feb 2007 14:13:24 +0900 Subject: display global data Message-ID: Dear everyone, My friend asked me the following questions, but I could not answer. I hope someone helps me. Question: There is a following ctl file. The data is global, starting from -179.75 with 0.5 deg resolution. When the display command is issued, then GrADS says that the display area is set to like, LON set to 0 360 LAT set to .... LEV set to .... Time .......... The question is that the original ctl shows that the western boundary of the data is from -179.75, however, the display says GrADS changed the longitude, starting from 0. Do you know why? ******************************************** Tetsuo Nakazawa nakazawa at mri-jma.go.jp Typhoon Research Department Meteorological Research Institute Tel: 029-853-8671 Fax: 029-853-8735 ******************************************** From horvath at CIRUS.DHZ.HR Fri Feb 23 05:05:22 2007 From: horvath at CIRUS.DHZ.HR (Kristian Horvath) Date: Fri, 23 Feb 2007 11:05:22 +0100 Subject: Binary to Ascii conversion In-Reply-To: <45DE13C6.6050704@noaa.gov> Message-ID: Hello, I also benefited from the discussion and tried the following script, which resulted in a peculiar result: 'reinit' 'sdfopen /doma/horvath/uwnd.sig995.1961.nc' 'set t 1' 'set lon 15 25' 'set lat 38 48' 'set gxout print' 'set prnopts %10.3e 5 1' 'd uwnd' rc=write(ascii.txt,result) rc=close(ascii.txt) 'quit' return The script works and the result is: Printing Grid -- 30 Values -- Undef = 32766 3.100e+00 1.900e+00 1.400e+00 1.300e+00 -1.200e+00 1.400e+00 1.100e+00 2.000e-01 -6.000e-01 -2.100e+00 2.000e-01 -2.000e-01 -7.000e-01 -5.000e-01 -8.000e-01 -1.900e+00 -2.100e+00 -2.500e+00 -1.800e+00 -8.000e-01 1.600e+00 1.400e+00 -1.000e-01 -8.000e-01 1.987e-06 5.700e+00 3.900e+00 1.700e+00 6.000e-01 7.000e-01 Now if I modify the script: 'set lon 15 25' 'set lat 38' ; result is Printing Grid -- 5 Values -- Undef = 32766 3.100e+00 1.900e+00 1.400e+00 1.300e+00 -1.200e+00 ; which is ok. and the second modification: 'set lon 15 25' 'set lat 48' Printing Grid -- 5 Values -- Undef = 32766 1.600e+00 1.400e+00 -1.000e-01 -8.000e-01 1.987e-06 ; which is not in accordance with the "area" example What is the reason for the above discrepancy? Thanks and regards, Kristian At 23:05 22.2.2007, you wrote: >If you don't need the data in binary (which is what you have in file >yiki.dat), you can skip the "gxout fwrite" lines and use "gxout print" >instead. Eg, in a grads script: > >'set gxout print' >'set t 613 697' >'display tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' >rc=write(yiki.txt,result) >say 'return code from write is: 'rc >rc=close(yiki.txt) >say 'return code from close is: 'rc > >In the above, the output from the display command gets stored in >variable "result" and you write the contents of "result" to ascii file >yiki.txt. > >See the doc for 'set prnopts' at >http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html >if you want to change the format of the text output. > > Diane > > >Yitbarek Kifle wrote: >>Dear all, >> >>As a new and inexperienced user of Grads,i have come with a problem of >>converting the Grads time series binary output to Ascii format.I have read >>almost all such problems on the users group list but that doesn't seem to >>help me. Here is the discription of what i did, >> >>'set gxout fwrite' >>'set fwrite yiki.dat' >>set t 613 697 >>d tloop(mean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5)) >>disable fwrite >> >>i have the file Yiki.dat. >> >>Would somebody please tell me what to do next in step by step basis,as I do >>not have a programming skill? >> >>thanks in advance. > > >________ Information from NOD32 ________ >This message was checked by NOD32 Antivirus System for Linux Mail Server. >http://www.nod32.com > >__________ NOD32 2072 (20070220) Information __________ > >This message was checked by NOD32 antivirus system. >http://www.eset.com > _________________________________ Kristian Horvath email: horvath at cirus.dhz.hr Graduate Research Associate tel:+385-1-4565752 Meteorological and Hydrological Service www: http://www.dhmz.htnet.hr/index.php Gric 3, 10000 Zagreb, Croatia personal www: http://radar.dhz.hr/~horvath From vello at EANS.EE Fri Feb 23 05:13:23 2007 From: vello at EANS.EE (=?ISO-8859-1?Q?Vello_M=FC=FCrsepp?=) Date: Fri, 23 Feb 2007 12:13:23 +0200 Subject: Removing from userlist In-Reply-To: A<7.0.1.0.0.20070223104047.02074ed0@cirus.dhz.hr> Message-ID: Hi, Please remove my e-mail address from this list. Thank you, Vello M??rsepp Estonian Air Navigation Services From anavfranco at YAHOO.COM.BR Fri Feb 23 08:06:33 2007 From: anavfranco at YAHOO.COM.BR (Ana Carol) Date: Fri, 23 Feb 2007 10:06:33 -0300 Subject: time average Message-ID: Hello, I would like to know how to make time (sazonal) average in grADS. Thanks, Ana Franco LECO/UDESC Laborat?rio de Estudos Clim?ticos e Oce?nicos Rua Visconde de Ouro Preto, 457 - Centro Florian?polis/SC Fone:(48)3222-9168 ramal 34 __________________________________________________ Fale com seus amigos de gra?a com o novo Yahoo! Messenger http://br.messenger.yahoo.com/ -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/a77c0039/attachment.html From yitbek at YAHOO.COM Fri Feb 23 09:34:18 2007 From: yitbek at YAHOO.COM (Yitbarek Kifle) Date: Fri, 23 Feb 2007 15:34:18 +0100 Subject: Binary to Ascii conversion Message-ID: I would like to thank Diane Stokes and Kristian Horvath for the Script that you sent me.I am still having a minor problem while excuting the script file.The problem is,when i run the script it turns out an error on line 9 of the scipt as follows; 'reinit' 'sdfopen D:\reanalysis\soilw.mon.mean.nc' 'set t 613 697' 'set lon 95.5 112.5' 'set lat 31.5 42.5' 'set gxout print' 'set prnopts %10.3e 1 1' 'd tloop(amean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' rc=write(D:\reanalysis\ascii.txt,result) rc=close(D:\reanalysis\ascii.txt) 'quit' return and the other script i have is; 'reinit' 'sdfopen D:\reanalysis\soilw.mon.mean.nc' 'set gxout print' 'set prnopts %8.5e 1' 'set t 613 697' 'display tloop(amean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' say result write ('D:\reanalysis\ascii.txt',result,append) 'reinit' and the later gives me an ascii file as follows;(What i didn't like is that the the text "Printing Grid -- 1 Values -- Undef = 32766"is repeated after every numeric output.) Printing Grid -- 1 Values -- Undef = 32766 1.89156e+02 Printing Grid -- 1 Values -- Undef = 32766 1.73351e+02 Printing Grid -- 1 Values -- Undef = 32766 1.64099e+02 Printing Grid -- 1 Values -- Undef = 32766 1.60951e+02 Printing Grid -- 1 Values -- Undef = 32766 1.70652e+02 Printing Grid -- 1 Values -- Undef = 32766 1.85685e+02 Printing Grid -- 1 Values -- Undef = 32766 2.07550e+02 Printing Grid -- 1 Values -- Undef = 32766 2.19629e+02 Printing Grid -- 1 Values -- Undef = 32766 2.19556e+02 Printing Grid -- 1 Values -- Undef = 32766 2.24317e+02 . . . Would somebody please make a correction of these scripts or send me another one? Thanks in advance. From lfmontanar at UNAL.EDU.CO Fri Feb 23 09:45:59 2007 From: lfmontanar at UNAL.EDU.CO (Luis Fernando Montana Roa) Date: Fri, 23 Feb 2007 09:45:59 -0500 Subject: Removing from userlist Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/04a33560/attachment.html From Charles.Seman at NOAA.GOV Fri Feb 23 10:27:50 2007 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Fri, 23 Feb 2007 10:27:50 -0500 Subject: time average In-Reply-To: <566317.60647.qm@web50405.mail.yahoo.com> Message-ID: Ana, Please see the "modify" command: http://grads.iges.org/grads/gadoc/gradcomdmodify.html I hope this helps, Chuck Ana Carol wrote: > Hello, > > I would like to know how to make time (sazonal) average in grADS. > > Thanks, > > > Ana Franco > > LECO/UDESC > Laborat?rio de Estudos Clim?ticos e Oce?nicos > Rua Visconde de Ouro Preto, 457 - Centro > Florian?polis/SC > Fone:(48)3222-9168 ramal 34 > > __________________________________________________ > Fale com seus amigos de gra?a com o novo Yahoo! Messenger > http://br.messenger.yahoo.com/ > -- 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 reflect any position of the Government or NOAA." From hrc at EAS.GATECH.EDU Fri Feb 23 10:06:19 2007 From: hrc at EAS.GATECH.EDU (Hai-Ru Chang) Date: Fri, 23 Feb 2007 10:06:19 -0500 Subject: q dims problem Message-ID: Hello everybody! I have a opteron workstation with the following architecture: Linux master 2.6.17-ck1-suse101-osmp #264 SMP Tue Oct 3 16:44:26 EDT 2006 x86_64 x86_64 x86_64 GNU/Linux. When I issued the 'q dims' command inside grads, I saw the following output Default file number is: 1 X is varying Lon = 147.258 to 232.741 X = 1.00002 to 89 Y is varying Lat = -39.2537 to 39.2511 Y = 1 to 89 Z is fixed Lev = 1000 Z = 1 T is fixed Time = 12Z01JUL30064773072 T = 1 There are two errors in the output. First, the lower limit of x is 1.00002 instead of 1. Second, the year on the time line is totally wrong. I used Portland group compiler to build up the software from source. Does any one know how to fix the problem? Thanks in advance! Dr. Hai-Ru Chang School of Earth and Atmospheric Sciences Georgia Institute of Technology Atlanta, GA 30332-0340 tel: 404-385-4404 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/c20050e3/attachment.html From radojevic.milka at OURANOS.CA Fri Feb 23 10:36:04 2007 From: radojevic.milka at OURANOS.CA (Milka Radojevic) Date: Fri, 23 Feb 2007 10:36:04 -0500 Subject: Removing from userlist Message-ID: Hi, Please remove my e-mail address from this list. Milka Radojevic (MSc) Research Assistant Global Environmental and Climate Change Centre (GEC3) & Adaptation and Impacts Research Division (AIRD) Environment Canada at Ouranos ==================== Tel: 514-282-6464 ext. 342 Fax: 514-282-7131 milka.radojevic at ec.gc.ca radojevic.milka at ouranos.ca ==================== Address: 550 Sherbrooke West 19th floor West Tower Montreal QC Canada H3A 1B9 Web links : http://www.mcgill.ca/gec3/ and http:// www.ccsn.ca -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/2a60e61b/attachment.html From jma at COLA.IGES.ORG Fri Feb 23 11:16:09 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Fri, 23 Feb 2007 11:16:09 -0500 Subject: q dims problem In-Reply-To: <448FA341-672B-4EE4-97C3-7C6632565474@eas.gatech.edu> Message-ID: On Feb 23, 2007, at 10:06 AM, Hai-Ru Chang wrote: > Hello everybody! > I have a opteron workstation with the > following architecture: > > > Linux master 2.6.17-ck1-suse101-osmp #264 SMP > Tue Oct 3 16:44:26 EDT 2006 x86_64 x86_64 x86_64 GNU/Linux. > > > When I issued the 'q dims' command inside grads, I saw the > following output > > Default file number is: 1 > X is varying Lon = 147.258 to 232.741 X = 1.00002 to 89 > Y is varying Lat = -39.2537 to 39.2511 Y = 1 to 89 > Z is fixed Lev = 1000 Z = 1 > T is fixed Time = 12Z01JUL30064773072 T = 1 > > > There are two errors in the output. First, the lower limit of x is > 1.00002 instead of 1. This is a precision issue because your lat/lon grid increment is small. This is mostly harmless -- doesn't affect displays, is usually only a problem when using set gxout fwrite. > Second, the year on the time line is totally wrong. I used > Portland group compiler to build up the software from source. > Does any one know how to fix the problem? This is a bug that appears on 64-bit linux systems. It has been fixed. Since you are building from source, you can edit the subroutine gat2ch in gautil.c. There were substitution strings for long integers ("%l") in the sprintf statements that should be regular integers ("%i") int gat2ch (struct dt *dtim, int tinc, char *ch) { int mn1,mn2,hr1,hr2,dy1,dy2,len,mnth; mnth = dtim->mo - 1L; mn1 = dtim->mn/10L; mn2 = dtim->mn - (mn1*10); hr1 = dtim->hr/10L; hr2 = dtim->hr - (hr1*10); dy1 = dtim->dy/10L; dy2 = dtim->dy - (dy1*10); if (tinc==1) { sprintf(ch,"%04i",dtim->yr); } else if (tinc==2) { if (dtim->yr==9999L) { sprintf(ch,"%s",monc[mnth]); } else { sprintf(ch,"%s%04i",monc[mnth],dtim->yr); } } else if (tinc==3) { sprintf(ch,"%i%i%s%04i",dy1,dy2,monc[mnth],dtim->yr); } else if (tinc==4) { sprintf(ch,"%i%iZ%i%i%s%04i",hr1,hr2,dy1,dy2, monc[mnth],dtim->yr); } else if (tinc==5) { sprintf(ch,"%i%i:%i%iZ%i%i%s%04i",hr1,hr2,mn1,mn2,dy1,dy2, monc[mnth],dtim->yr); } else sprintf(ch,"???"); len=0; while (ch[len]) len++; return (len); } -- 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/20070223/53b9b412/attachment.html From jma at COLA.IGES.ORG Fri Feb 23 11:40:56 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Fri, 23 Feb 2007 11:40:56 -0500 Subject: display global data In-Reply-To: <27500ACE-6FA2-43B8-83DF-EC856963A865@mri-jma.go.jp> Message-ID: On Feb 23, 2007, at 12:13 AM, Tetsuo Nakazawa wrote: > Dear everyone, > > My friend asked me the following questions, but I could not answer. > I hope someone helps me. > > Question: > There is a following ctl file. The data is global, starting from > -179.75 with 0.5 deg resolution. > When the display command is issued, then GrADS says that the > display area is set to like, > > LON set to 0 360 > LAT set to .... > LEV set to .... > Time .......... > > The question is that the original ctl shows that the western > boundary of the data is from -179.75, > however, the display says GrADS changed the longitude, starting from > 0. Do you know why? When GrADS detects that a grid 'wraps' the globe, it sets the LON dimension environment to be 0 to 360. If it doesn't wrap, it sets the X dimension environment to be 1 to XSIZE. The test for whether a grid wraps goes like this: In your descriptor file you have: XDEF n linear v1 v2 temp = (v1 + (n*v2) - 360) if (fabs(temp-v1) < 0.01) then the grid wraps. Non-linear X axes are never tested to see if they wrap. 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/20070223/e3066c77/attachment.html From Diane.Stokes at NOAA.GOV Fri Feb 23 11:55:19 2007 From: Diane.Stokes at NOAA.GOV (Diane Stokes) Date: Fri, 23 Feb 2007 11:55:19 -0500 Subject: Binary to Ascii conversion In-Reply-To: <20070223143503.415612331D@mx2.cineca.it> Message-ID: Since you are doing an area average in the call to amean, you don't need to set varying lat and lon ranges before your display. Instead set lon and lat (or x and y) to a fixed point. Eg: 'set x 1' 'set y 1' If you do so, that "Printing grid" header message should show up only at the top of the output. Diane Yitbarek Kifle wrote: > I would like to thank Diane Stokes and Kristian Horvath for the Script that > you sent me.I am still having a minor problem while excuting the script > file.The problem is,when i run the script it turns out an error on line 9 > of the scipt as follows; > > 'reinit' > 'sdfopen D:\reanalysis\soilw.mon.mean.nc' > 'set t 613 697' > 'set lon 95.5 112.5' > 'set lat 31.5 42.5' > 'set gxout print' > 'set prnopts %10.3e 1 1' > 'd tloop(amean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' > rc=write(D:\reanalysis\ascii.txt,result) > rc=close(D:\reanalysis\ascii.txt) > 'quit' > return > > and the other script i have is; > > 'reinit' > 'sdfopen D:\reanalysis\soilw.mon.mean.nc' > 'set gxout print' > 'set prnopts %8.5e 1' > 'set t 613 697' > 'display tloop(amean(soilw,lon=95.5,lon=112.5,lat=31.5,lat=42.5))' > say result > write ('D:\reanalysis\ascii.txt',result,append) > 'reinit' > > and the later gives me an ascii file as follows;(What i didn't like is that > the the text "Printing Grid -- 1 Values -- Undef = 32766"is repeated after > every numeric output.) > > Printing Grid -- 1 Values -- Undef = 32766 > 1.89156e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 1.73351e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 1.64099e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 1.60951e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 1.70652e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 1.85685e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 2.07550e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 2.19629e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 2.19556e+02 > Printing Grid -- 1 Values -- Undef = 32766 > 2.24317e+02 > . > . > . > Would somebody please make a correction of these scripts or send me another > one? > > Thanks in advance. > From sansaad at AJATO.COM.BR Fri Feb 23 11:50:46 2007 From: sansaad at AJATO.COM.BR (sansaad) Date: Fri, 23 Feb 2007 16:50:46 +0000 Subject: problems with aave In-Reply-To: <45DF07F6.9090002@noaa.gov> Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/1427357d/attachment.html From Diane.Stokes at NOAA.GOV Fri Feb 23 13:24:14 2007 From: Diane.Stokes at NOAA.GOV (Diane Stokes) Date: Fri, 23 Feb 2007 13:24:14 -0500 Subject: problems with aave In-Reply-To: <21418Y1172249446Y24657Ykon2@velop3> Message-ID: I see two (maybe three) contributing factors to the difference: First, the result from aave is latitude weighted. From the documentation: "The aave function assumes that the world coordinates are longitude in the X dimension and latitude in the Y dimension, and does weighting in the latitude dimension by the difference between the sines of the latitude at the northern and southern edges of the grid box. For areal averaging without latitude weighting, use the amean function." The impact of using aave vs amean depends on the latitude range of your region. Second, in your fortran code, the denominator is too small when you divide mh(i) by the number of gridsquares. You have: mh(i)=mh(i)/((35-26)*(50-43)) Try: mh(i)=mh(i)/((35-26+1)*(50-43+1)) A third (possible) factor is that you have missing value(s) that would get excluded from the aave (or amean) calculation but are not excluded in the code. The impact from that would depend on the value you use for "missing". Diane sansaad wrote: > Hi everybody > > I'd done some calculations using both GRADS and Fortran, and I realised > that the spatial average using GRADS (aave) differs by a ratio of > approximately 0.787 from the spatial average using Fortran. By this I > mean that if I multiply the Fortran's values by 0.787, they ALWAYS match > with the GRADS' values (with small errors). > > Can anyone explain the reason for this "magic number" (0.787)? > > A part of the program using GRADS and Fortran respectively are given > bellow. > > I would appreciate any kind of commentaries. > > Sandra Saad > > > ### Spatial Average using GRADS ######## > > set t 1 last > d aave(h,x=26,x=35,y=43,y=50) > > > ### Spatial Average using Fortran (90) ######## > > do i=1,it > do ny=43,50 > do nx=26,35 > mh(i)=mh(i)+h(i,nx,ny) > enddo > enddo > mh(i)=mh(i)/((35-26)*(50-43)) > enddo > > do i=1,it > write(3,*)i,mh(i) > enddo > From renato at WEB.HIDRO.UFCG.EDU.BR Fri Feb 23 13:48:13 2007 From: renato at WEB.HIDRO.UFCG.EDU.BR (Renato Oliveira) Date: Fri, 23 Feb 2007 15:48:13 -0300 Subject: userlist In-Reply-To: Message-ID: Please, remove me from userlist. **Renato Oliveira** Laborat?rio de Hidr?ulica II Universidade Federal de Campina Grande - UFCG From sansaad at AJATO.COM.BR Fri Feb 23 14:48:36 2007 From: sansaad at AJATO.COM.BR (sansaad) Date: Fri, 23 Feb 2007 19:48:36 +0000 Subject: problems with aave In-Reply-To: <45DF314E.8040707@noaa.gov> Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070223/b7099604/attachment.html From hrc at EAS.GATECH.EDU Fri Feb 23 14:56:34 2007 From: hrc at EAS.GATECH.EDU (Hai-Ru Chang) Date: Fri, 23 Feb 2007 14:56:34 -0500 Subject: gxeps on opteron Message-ID: Hello, I compiled gradsnc with pgi compiler on my opteron workstation (Suse 10.1). The gxeps created a file of 30 bytes, but gxps created a file of 8084139 bytes. Dose anyone know how to fix the problem. Thanks! Hai-Ru From adhikary at CGRER.UIOWA.EDU Fri Feb 23 16:38:17 2007 From: adhikary at CGRER.UIOWA.EDU (Bhupesh Adhikary) Date: Fri, 23 Feb 2007 15:38:17 -0600 Subject: Change email in the userlist Message-ID: Dear Grads users: Please change my GRADS userlist email address from this to bhupesh74 at hotmail.com Thanks bhupesh From nakazawa at MRI-JMA.GO.JP Sun Feb 25 22:06:34 2007 From: nakazawa at MRI-JMA.GO.JP (Tetsuo Nakazawa) Date: Mon, 26 Feb 2007 12:06:34 +0900 Subject: display global data Message-ID: Dear Jennifer-san, >> My friend asked me the following questions, but I could not answer. >> I hope someone helps me. >> >> Question: >> There is a following ctl file. The data is global, starting from >> -179.75 with 0.5 deg resolution. >> When the display command is issued, then GrADS says that the >> display area is set to like, >> >> LON set to 0 360 >> LAT set to .... >> LEV set to .... >> Time .......... >> >> The question is that the original ctl shows that the western >> boundary of the data is from -179.75, >> however, the display says GrADS changed the longitude, starting from >> 0. Do you know why? >> > When GrADS detects that a grid 'wraps' the globe, it sets the LON > dimension environment to be 0 to 360. If it doesn't wrap, it sets > the X dimension environment to be 1 to XSIZE. The test for whether > a grid wraps goes like this: > In your descriptor file you have: XDEF n linear v1 v2 > temp = (v1 + (n*v2) - 360) > if (fabs(temp-v1) < 0.01) then the grid wraps. > Non-linear X axes are never tested to see if they wrap. > Thanks. Sometimes one would like to plot globally as the data is stored. And also, one suspects that the data is wrong or GrADS is wrong if GrADS changes the starting longitude from 0 deg. ******************************************** Tetsuo Nakazawa nakazawa at mri-jma.go.jp Typhoon Research Department Meteorological Research Institute Tel: 029-853-8671 Fax: 029-853-8735 ******************************************** From johnguhin at GMAIL.COM Mon Feb 26 02:14:13 2007 From: johnguhin at GMAIL.COM (John Guhin) Date: Mon, 26 Feb 2007 00:14:13 -0700 Subject: Template problem Message-ID: Hi, I have three files like these.... 199701_01-09_2D.ctl 199701_10-19_2D.ctl 199701_20-31_2D.ctl I wrote a control file like this... dset ^199701_%d2-%d2_2D.dat options template ................................................ I can open the control file, but while plotting it says Cannot contour grid - all undefined values But I can plot from individual files. Where am I doing wrong? -- John Guhin -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070226/a1638787/attachment.html From mmahakur at TROPMET.RES.IN Mon Feb 26 03:22:58 2007 From: mmahakur at TROPMET.RES.IN (M.Mahakur) Date: Mon, 26 Feb 2007 13:52:58 +0530 Subject: Template problem In-Reply-To: <37b4e6a50702252314k7805f748s9488b0007df10b37@mail.gmail.com> Message-ID: Hello, Use string substitution for your non standard template: http://www.iges.org/grads/gadoc/descriptorfile.html#CHSUB Regards, Mahakur John Guhin wrote: > > Hi, > I have three files like these.... > 199701_01-09_2D.ctl > 199701_10-19_2D.ctl > 199701_20-31_2D.ctl > > I wrote a control file like this... > dset ^199701_%d2-%d2_2D.dat > options template > ................................................ > > I can open the control file, but while plotting it says > Cannot contour grid - all undefined values > > But I can plot from individual files. Where am I doing wrong? > > -- > John Guhin From dbcl at CEH.AC.UK Mon Feb 26 04:02:50 2007 From: dbcl at CEH.AC.UK (Douglas Clark) Date: Mon, 26 Feb 2007 09:02:50 +0000 Subject: Template problem Message-ID: John, Assuming that your files have 9 or 10 days' of data in each, I think you'll need to use the CHSUB option - see GrADS index http://grads.iges.org/grads/gadoc/gadocindex.html under control file > options > template (or the likes). With your current template, the way the template works is such that every day's data is assumed to be in a different file - because you have %d2 in the template. So when you set the time to 1997/01/01, GrADS looks for data in a file called 199701_01-01_2D.dat, which doesn't exist - so it says all values are undefined (arguably it might be more useful if GrADS said "can't find file" at this point, but it doesn't). With time 1997/01/02, GrADS looks for data in a file called 199701_02-02_2D.dat, and so on. I think CHSUB will come to your rescue (although I have never used it myself). You need a ctl file with something along the lines of DSET ^199701_%ch_2D.dat CHSUB 1 9 01-09 CHSUB 10 19 10-19 CHSUB 20 31 20-31 OPTIONS template The first two integers on each CHSUB line will depend on the time step given in yout TDEF line - I have assumed daily data. See the GrADS documentation for details. Doug >>> johnguhin at GMAIL.COM 26/02/2007 07:14 >>> Hi, I have three files like these.... 199701_01-09_2D.ctl 199701_10-19_2D.ctl 199701_20-31_2D.ctl I wrote a control file like this... dset ^199701_%d2-%d2_2D.dat options template ................................................ I can open the control file, but while plotting it says Cannot contour grid - all undefined values But I can plot from individual files. Where am I doing wrong? -- John Guhin -- 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 jma at COLA.IGES.ORG Mon Feb 26 14:34:44 2007 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Mon, 26 Feb 2007 14:34:44 -0500 Subject: Other GrADS distribution sites Message-ID: Dear All, We are wondering how many other sites besides COLA are distributing a verison of GrADS. Please don't be shy; we're not angry or jealous, just curious -- we'd like to know where else GrADS is available. If you are distributing a version of GrADS, please send me the URL. Thanks! 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/20070226/bd7da3e7/attachment.html From Wesley.Ebisuzaki at NOAA.GOV Mon Feb 26 15:09:48 2007 From: Wesley.Ebisuzaki at NOAA.GOV (Wesley Ebisuzaki) Date: Mon, 26 Feb 2007 15:09:48 -0500 Subject: Other GrADS distribution sites In-Reply-To: Message-ID: from Jennifer Adams wrote: > Dear All, > We are wondering how many other sites besides COLA are distributing a > verison of GrADS. Please don't be shy; we're not angry or jealous, > just curious -- we'd like to know where else GrADS is available. If > you are distributing a version of GrADS, please send me the URL. > Thanks! > Jennifer > > > -- > Jennifer M. Adams > IGES/COLA > 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705 > jma at cola.iges.org > > > Jennifer, Not a URL but have been including GrADS on various CD-ROMs and DVDs. Annual NCEP/NCAR Reanalysis (1950-2005), distributed by NCAR Kalnay et al., 1996: The NCEP/NCAR 40 Year Reanalysis Project. Bull. Amer. Meteor. Soc., 77, 437-471. Kistler et al. 2001: The NCEP-NCAR 50 Year Reanalysis: Monthly Means CD-ROM and Documentation, Bull. Amer. Meteor. Soc., 82, 247-267. P. J. Koch and L. W. Uccellini, 2004: "Northeast Snowstorms, Volume II: The Cases", Meteorological Monographs, vol. 32, num. 54, Amer. Meteor. Soc. F. Mesinger, et al., 2006: "North American Regional Reanalysis", Bulletin of the Amer. Meter. Soc., 87, 343-360. Wesley Ebisuzaki CPC/NCEP From mezenga100 at IG.COM.BR Mon Feb 26 16:45:02 2007 From: mezenga100 at IG.COM.BR (mezenga100) Date: Mon, 26 Feb 2007 18:45:02 -0300 Subject: display global data Message-ID: Hi, Please remove my e-mail address from this list. Thank you, Sergio Buarque UFAL - Brasil/Alagoas/Maceio From wt at ATMOS.COLOSTATE.EDU Mon Feb 26 16:55:54 2007 From: wt at ATMOS.COLOSTATE.EDU (Warren Turkal) Date: Mon, 26 Feb 2007 14:55:54 -0700 Subject: compiling GrADS Message-ID: Does anyone have any information with regard to compiling GrADS with full NetCDF functionality. I have the 3.6.2-beta6 netcdf libs on my system and GrADS doesn't seem to want to recognize them. Does anyone have any ideas? wt -- Warren Turkal From pertusus at FREE.FR Mon Feb 26 16:26:59 2007 From: pertusus at FREE.FR (Patrice Dumas) Date: Mon, 26 Feb 2007 22:26:59 +0100 Subject: Other GrADS distribution sites In-Reply-To: Message-ID: On Mon, Feb 26, 2007 at 02:34:44PM -0500, Jennifer Adams wrote: > Dear All, > We are wondering how many other sites besides COLA are distributing a > verison of GrADS. Please don't be shy; we're not angry or jealous, > just curious -- we'd like to know where else GrADS is available. If > you are distributing a version of GrADS, please send me the URL. I have one here (with changes in the autotools stuff and some portability changes), you certainly already know: http://www.environnement.ens.fr/perso/dumas/fc-srpms/grads-1.9b4.1.tar.gz There is also the package in fedora that I maintain. It is grads-1.9b4 without non GPL compatible code. For example (for FC-5): http://download.fedora.redhat.com/pub/fedora/linux/extras/5/i386/repodata/repoview/grads-0-1.9b4-19.fc5.html This package is also in fedora extras rebuild for CentOS: http://centos.karan.org/el4/extras/stable/i386/RPMS/repodata/repoview/grads-0-1.9b4-11.el4.kb.html The corresponding SRPMS have the grads tarball within them. The source and binary packages are available on all the mirrors, I can't provide a list ;-). -- Pat From pertusus at FREE.FR Mon Feb 26 17:06:11 2007 From: pertusus at FREE.FR (Patrice Dumas) Date: Mon, 26 Feb 2007 23:06:11 +0100 Subject: compiling GrADS In-Reply-To: <200702261455.54775.wt@atmos.colostate.edu> Message-ID: On Mon, Feb 26, 2007 at 02:55:54PM -0700, Warren Turkal wrote: > Does anyone have any information with regard to compiling GrADS with full > NetCDF functionality. I have the 3.6.2-beta6 netcdf libs on my system and > GrADS doesn't seem to want to recognize them. Does anyone have any ideas? You may want to try what I have here: http://www.environnement.ens.fr/perso/dumas/fc-srpms/grads-1.9b4.1.tar.gz It will search the libs on the system if you pass --enable-dyn-supplibs, and there are other ./configure switches to specify where netcdf, hdf, or libnc-dap are. With the original grads-1.9b4 from COLA, you can still do something similar with a call along: ./configure LDFLAGS="-L/some/where/netcdf/" CPPFLAGS="-I/some/where/include/netcd/" However to use a system OPeNDAP I suggest using the tarball I provide. -- Pat From portis at ATMOS.UIUC.EDU Mon Feb 26 16:11:49 2007 From: portis at ATMOS.UIUC.EDU (Diane Portis) Date: Mon, 26 Feb 2007 15:11:49 -0600 Subject: Remote analysis with GrADSDODS Message-ID: Hello, I have been trying out the remote analysis capabilities of gds. I was able to do it when the remote dataset was not in a subdirectory (e.g. the ssta dataset that is located at http://cola8.iges.org:9090/dods), however when the dataset was located in a subdirectory, I received an error message that the dataset did not exist or was not an netcdf file. The first paste below shows that the file is readable with the sdfopen command, the second paste lists my script to do remote analysis on that dataset and the error message when I run the script. I suspect that it is a syntax error. Thanks for any insights, Diane ---Opening the file with sdfopen - okay----------------------------------------------------------------------- > ga-> sdfopen > http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete > Scanning self-describing file: > http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete > SDF file > http://nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete > is open as file 1 > LON set to -220 -0.625 > LAT set to 0 89.625 > LEV set to 1000 1000 > Time values set: 1979:1:1:0 1979:1:1:0 --Trying to do analysis on server side------------------------------------- script: climt1.gs > baseurl = 'http://nomads.ncdc.noaa.gov:9091/dods/_expr_' > datasets = 'NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete' > expr='tloop(aave(tcdc,lon=-220,lon=0.625,lat=70,lat=89.625))' > xlim='1:1' > ylim='1:1' > zlim='1:1' > tlim='00Z01jan1979:00Z01may1979' > 'sdfopen 'baseurl'{'datasets'}{'expr'}{'xlim','ylim','zlim','tlim'}' > 'define var = result.1' > 'd var' > ------------------------------------------------------------------ > ga-> run climt1.gs > nc_open failure: > Unknown error 1001 > http://nomads.ncdc.noaa.gov:9091/dods/_expr_{NCEP_NARR_MONTHLY_AGGREGATIONS/narrmon-a_221_complete} (broke this line up when pasting, so it could be read) > {tloop(aave(tcdc,lon=-220,lon=0.625,lat=70,lat=89.625))}{1:1,1:1,1:1,00Z01jan1979:00Z01may1979} > does not exist or is not a netCDF file. > > Couldn't ingest SDF metadata. > If this was an HDF-SDS file, try gradshdf. > DEFINE error: no file open yet > DISPLAY error: no file open yet > ga-> From ronny.meyer at ENERGYMETEO.DE Tue Feb 27 10:46:27 2007 From: ronny.meyer at ENERGYMETEO.DE (Ronny Meyer) Date: Tue, 27 Feb 2007 16:46:27 +0100 Subject: wgrib --> wgrib2 Message-ID: Dear GRADS-USERS, I have a question concering the output of wgrib2. Is the output of wgrib2 the same as the output of wgrib, in that sense, that the formatting of extracted points is the same? Thanks in advance. (and on this way: Thanks for all that great programs) Ronny Meyer -- Dipl.-Phys. Ronny Meyer Projektleiter Forschung und Entwicklung / Project Manager Research and Development energy & meteo systems GmbH Marie-Curie-Str. 1 26129 Oldenburg Germany fon +49(0) 441 36116 474 fax +49(0) 441 36116 479 ronny.meyer at ergymeteo.de www.energymeteo.de energy & meteo systems GmbH, Oldenburg Gesch?ftsf?hrer: Dr. Ulrich Focken, Dr. Matthias Lange Amtsgericht Oldenburg HRB Nr. 5596 Der Inhalt dieser E-Mail ist vertraulich und ausschlie?lich f?r den beabsichtigten Empf?nger bestimmt. Falls Sie nicht der beabsichtigte Empf?nger sind oder falls diese E-Mail irrt?mlich an Sie versandt wurde, verst?ndigen Sie bitte umgehend den Absender und l?schen Sie anschlie?end die E-Mail einschlie?lich aller Anlagen von Ihrem System. Jeder(r) unberechtigte Lekt?re, Gebrauch, Ver?ffentlichung oder Weitergabe ist untersagt. Die Sicherheit von E-Mail-?bermittlungen kann nicht garantiert werden. This e-mail contains confidential information and is for the intended recipient only. If you are not the intended recipient, or if you have received this email in error, please notify the sender immediately, and delete the e-mail (including any attachment/s) subsequently. Please do not read, copy, use, distribute or disclose the contents of this communication to other. The confidentiality of email cannot be guaranteed. -------------- next part -------------- A non-text attachment was scrubbed... Name: ronny.meyer.vcf Type: text/x-vcard Size: 532 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070227/18eed5f9/attachment.vcf From mcroke at AIRDAT.COM Tue Feb 27 13:02:54 2007 From: mcroke at AIRDAT.COM (Meredith Croke) Date: Tue, 27 Feb 2007 13:02:54 -0500 Subject: Skipping files and continue Message-ID: Hi All, I'm new to GrADS and am trying to write a script that opens one file at a time, performs a funtion, writes the data to a .dat file, then closes and opens the next ctl file. I have several hundred files that I am opening. The problem I am running into is this; all the CTL files exist, however there are a few data files that do not exist. When it reaches these files I get a segmentation fault. Is there anyway I can write a statement that says if the file does not exist to print some value (either 999 or NA), add 3 to the hour and continue to the next file? Thank you very much for your help! Meredith From Michael.Fiorino at NOAA.GOV Tue Feb 27 13:15:03 2007 From: Michael.Fiorino at NOAA.GOV (Michael Fiorino) Date: Tue, 27 Feb 2007 13:15:03 -0500 Subject: display global data In-Reply-To: <2EB2E83C-4B4E-4E72-B20A-A513FDB63E2E@mri-jma.go.jp> Message-ID: Nakazawa-san, the best way to do that would be to 'set x 1 N points' this forces to display the exact grid. best /r Mike Tetsuo Nakazawa wrote: > Dear Jennifer-san, > >>> My friend asked me the following questions, but I could not answer. >>> I hope someone helps me. >>> >>> Question: >>> There is a following ctl file. The data is global, starting from >>> -179.75 with 0.5 deg resolution. >>> When the display command is issued, then GrADS says that the >>> display area is set to like, >>> >>> LON set to 0 360 >>> LAT set to .... >>> LEV set to .... >>> Time .......... >>> >>> The question is that the original ctl shows that the western >>> boundary of the data is from -179.75, >>> however, the display says GrADS changed the longitude, starting from >>> 0. Do you know why? >>> >> When GrADS detects that a grid 'wraps' the globe, it sets the LON >> dimension environment to be 0 to 360. If it doesn't wrap, it sets >> the X dimension environment to be 1 to XSIZE. The test for whether >> a grid wraps goes like this: >> In your descriptor file you have: XDEF n linear v1 v2 >> temp = (v1 + (n*v2) - 360) >> if (fabs(temp-v1) < 0.01) then the grid wraps. >> Non-linear X axes are never tested to see if they wrap. >> > > Thanks. > > Sometimes one would like to plot globally as the data is stored. > And also, one suspects > that the data is wrong or GrADS is wrong if GrADS changes the > starting longitude from 0 deg. > > ******************************************** > Tetsuo Nakazawa > nakazawa at mri-jma.go.jp > Typhoon Research Department > Meteorological Research Institute > Tel: 029-853-8671 Fax: 029-853-8735 > ******************************************** From mjin at ATMOS.UMD.EDU Tue Feb 27 15:08:39 2007 From: mjin at ATMOS.UMD.EDU (Menglin Jin) Date: Tue, 27 Feb 2007 15:08:39 -0500 Subject: read multi data files Message-ID: Hi, I use the following .ctl file to read July data from 1982-1999, and the data files' names are "82.07.***", "83.07.**8"... dset ../data/%y2.07.tau.grads options template undef -9.99E+33 xdef 360 linear -180.000 1.00 ydef 180 linear -90.000 1.00 zdef 1 levels 1000 tdef 15 linear 00:00Z01jul1982 1mo vars 1 aod 0 99 MODIS product endvars But I run into error message when I read it in grads: Low Level I/O Error: Read error on data file Data file name = ../data/%y2.07.tau.grads Error reading 360 bytes at location 64800 Each data file itself is good, as I can plot it by reading in each file separately. Caould anyone give me some suggestion? Many thanks. ML From singh at SATELLITE.HYARC.NAGOYA-U.AC.JP Tue Feb 27 21:42:14 2007 From: singh at SATELLITE.HYARC.NAGOYA-U.AC.JP (Prasamsa Singh) Date: Wed, 28 Feb 2007 11:42:14 +0900 Subject: plotting river map Message-ID: Dear GRADS-USERS, I would like to plot permanent rivers map in my figures. Is it possible to plot rivers on map using grads ? Thanking you in advance, Prasamsa Singh --------- Prasamsa Singh Doctoral student, Hydrospheric & Atmospheric Research Center (HyARC) Nagoya University From matias at TROPOSFERA.ATMOSFCU.UNAM.MX Wed Feb 28 01:39:59 2007 From: matias at TROPOSFERA.ATMOSFCU.UNAM.MX (J. Matias Mendez Perez) Date: Wed, 28 Feb 2007 00:39:59 -0600 Subject: plotting river map In-Reply-To: <20070228114214.2c907386.singh@satellite.hyarc.nagoya-u.ac.jp> Message-ID: Dear Prasamsa, I've created hidrography (rivers) from Mexico on GrADS. I used gradsmap.f, a fortran program developed by Dr. Kumar from Indian Institute of Tropical Meteorology. His program prepares backgroud map for GrADS from a text file containig line segments in coordinates. I'm sending you an attachment with the following files: gradsmap.f fortran code hidrografia.asc input data hidro_mex output background map (hidrography) plot_rivers.gs script to plot rivers from Mexico map.ctl control file dummy.bin binary file hidro_mex.ps output graph You will have to prepare your own data according with instructions inside fortran code. I hope this can help you. Best regards from Mexico, Matias > Dear GRADS-USERS, > > I would like to plot permanent rivers map in my figures. > Is it possible to plot rivers on map using grads ? > > Thanking you in advance, > Prasamsa Singh > > > --------- > Prasamsa Singh > Doctoral student, > Hydrospheric & Atmospheric Research Center (HyARC) > Nagoya University > -------------- next part -------------- A non-text attachment was scrubbed... Name: hidro_mex.tar Type: application/x-tar Size: 1699840 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20070228/fcb50e67/attachment.tar From singh at SATELLITE.HYARC.NAGOYA-U.AC.JP Wed Feb 28 02:05:48 2007 From: singh at SATELLITE.HYARC.NAGOYA-U.AC.JP (Prasamsa Singh) Date: Wed, 28 Feb 2007 16:05:48 +0900 Subject: plotting river map In-Reply-To: <56695.189.136.150.17.1172644799.squirrel@atmosfera.unam.mx> Message-ID: Dear Matias, Thank you very much for your suggestion. I will try it. With regards, Prasamsa Singh On Wed, 28 Feb 2007 00:39:59 -0600 "J. Matias Mendez Perez" wrote: > Matias --------- Prasamsa Singh Doctoral student, Hydrospheric & Atmospheric Research Center (HyARC) Nagoya University From giuliani at LAMMA.RETE.TOSCANA.IT Wed Feb 28 03:26:22 2007 From: giuliani at LAMMA.RETE.TOSCANA.IT (Graziano Giuliani) Date: Wed, 28 Feb 2007 09:26:22 +0100 Subject: Other GrADS distribution sites In-Reply-To: Message-ID: On Mon, 2007-02-26 at 14:34 -0500, Jennifer Adams wrote: > Dear All, > We are wondering how many other sites besides COLA are distributing a > verison of GrADS. Please don't be shy; we're not angry or jealous, > just curious -- we'd like to know where else GrADS is available. If > you are distributing a version of GrADS, please send me the URL. [...] I am currently the mantainer of the archlinux grads package at http://aur.archlinux.org/packages.php?do_Details=1&ID=7868 It is not a binary package, You have to recompile it, but it includes a patch which: 1) Solves some warnings in compilation 2) Use Netcdf-3 interface instead of Netcdf-2 for netcdf sdf 3) Add shapefile format capabilities for mpdset data sources (links against shapelib (http://shapelib.maptools.org/). 4) Correct a small bug for the linking of imagemagick for the deprecated wi capability. 5) Link against new libnc-dap and dap++ 6) Revamps autotools layer I had another GRADS "custom" version, with capability to open stdout as file target for images (useful in CGI environment where I did not want for megabytes of scratch images lying around), but I did distribute it more than five years ago and the server where it was hosted has long gone. Was not hard to implement, but now with cheap-as-chips storage I think it should be not nedded anymore. Ciao, Graziano. -- Graziano Giuliani LaMMA - CNR Ibimet / Regiona Toscana From eduardo.agosta at GMAIL.COM Wed Feb 28 07:02:56 2007 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Wed, 28 Feb 2007 09:02:56 -0300 Subject: filtering Message-ID: Hi you all! Does anybody have a script or certain instructions to perform time filtering (gaussian, binomial, or whatever) on the gridded fields? Thanks in advance. Eduardo -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070228/88897844/attachment.html From remm69 at GMAIL.COM Wed Feb 28 08:34:26 2007 From: remm69 at GMAIL.COM (Rafael Mundaray) Date: Wed, 28 Feb 2007 09:34:26 -0400 Subject: Declaration of Variable Message-ID: Hi to all! I have a problem to open of my output file of my WRF model. The output file is in CDF format and the GrADS gives a referred error me to that it does not find variable X. I used the xdfopen command and the GrADS says that it needs a description file. Then, I need to know, HOW DO BUILD THIS DESCRIPTOR FILE? The file contains 50 atmospheric variables approximately. Thank you...!!! -- Ing. Rafael E. Mundaray M. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20070228/fb4a5def/attachment.html From chintalu at TROPMET.RES.IN Wed Feb 28 09:24:11 2007 From: chintalu at TROPMET.RES.IN (chintalu) Date: Wed, 28 Feb 2007 19:54:11 +0530 Subject: How to draw wind rose using Bob Harts Plotskew In-Reply-To: <43331319.80309@free.fr> Message-ID: > Dear grads users I need some help regarding how to plot wind rose using Bob Harts Plotskew if any body has already done it . Regarding Bob Harts Plotskew I have few questions i) I have radio sonde data for several stations I need to compute thermodynamic parameters for e.g theta thetae , u, v, q, qs, es, ep Do I need to make the undesirable part of the script as comment ii ) How to make input file which can be readable by the Bob Harts plotsqew script Some examples of input and output files and modified script as an example could be very helpful. Thanks in advance for your help G.R.Chinthalu *************************************************8 > > From bruyerec at UCAR.EDU Wed Feb 28 11:14:36 2007 From: bruyerec at UCAR.EDU (Cindy Bruyere) Date: Wed, 28 Feb 2007 09:14:36 -0700 Subject: Declaration of Variable In-Reply-To: <1f59b3870702280534p6bac5d7dx921dd743ade4bb1f@mail.gmail.com> Message-ID: For WRF model output you may want to read this: http://www.mmm.ucar.edu/wrf/users/graphics/ARWpost/ARWpost.htm The WRF model has a staggered grid, so in general you way want to use the above converter to convert the data to unstaggered grids, which are easier to deal with. GrADS 1.9 has a "gradsnc" script that can display the staggered grid of WRF ARW. But you need 4 .ctl files to display the data correctly. Inside the ARWpost code there is a utility (under the directory util/) that you can use to generate the 4 .ctl files. Cindy On Feb 28, 2007, at 6:34 AM, Rafael Mundaray wrote: > Hi to all! > > I have a problem to open of my output file of my WRF model. The > output file is in CDF format and the GrADS gives a referred error > me to that it does not find variable X. I used the xdfopen command > and the GrADS says that it needs a description file. Then, I need > to know, HOW DO BUILD THIS DESCRIPTOR FILE? The file contains 50 > atmospheric variables approximately. > > Thank you...!!! > > -- > Ing. Rafael E. Mundaray M. From lalvarez at IDEMOC.COM Wed Feb 21 14:17:27 2007 From: lalvarez at IDEMOC.COM (Leopoldo Alvarez - IDEMOC -) Date: Wed, 21 Feb 2007 20:17:27 +0100 Subject: Smooth the variables in model output In-Reply-To: <20070221134225.M17423@model.iag.usp.br> Message-ID: Dear Ricardo, Thanks for this method. I need another form to smooth much more, for example, when i attempt to smooth the geopotencial 850 hPa that crosses a strong relief. regards, Leopoldo Quoting Ricardo Hallak : > Dear Leopoldo, > > try: > > "d smth9(var)" > > or add more power in this function: > > "d smth9(smth9(var))" > > or yet: > > "d smth9(smth9(...(smth9(var))...))" > > Regards, > Ricardo > > On Wed, 21 Feb 2007 13:47:44 +0100, Leopoldo Alvarez wrote >> Dear users GrADS, >> >> I need to smooth the variables that i obtain with MM5 model output. >> Due to the presence of terrain, etc... I have much noise in the >> temperature, humidity, geopotencial variables... >> >> I have used “set csmooth” without success in grads. >> >> regards, >> Leopoldo Alvarez >> - - >> - - - - - - - - - - - - - - - - - - - - - - - - - - - >> Leopoldo Alvarez >> Adviser - Tecnico Investigator >> Meteorology and Oceanographic Studies Institute >> Canary Islands - Spain >> E-mail: lalvarez at idemoc.com >> Tlfno: (+34) 600 253 237 >> - - - - - - - - - - - - - - - - - - - - - - - - - - - > > > Hallak > > Laborat?rio de Hidrometeorologia (LABHIDRO) > Departamento de Ci?ncias Atmosf?ricas (DCA) > Instituto de Astronomia, Geof?sica e Ci?ncias Atmosf?ricas (IAG) > Universidade de S?o Paulo (USP) > Rua do Mat?o, 1226 - Cidade Universit?ria -S?o Paulo SP -05508-900 > Fone.: +55 (0-XX-11) 3091-4735 > Fax: +55 (0-XX-11) 3091-4714 >