From medha.deshpande at REDIFFMAIL.COM Tue Jul 1 03:10:52 2008 From: medha.deshpande at REDIFFMAIL.COM (medha deshpande) Date: Tue, 1 Jul 2008 07:10:52 -0000 Subject: Vertical Cross Section of hurricane Message-ID: ? Dear GrADS users I wish to draw the line for radius of maximum wind in the vertical cross section of a hurricane. Can anyone help me in this regard? Thanks in advance Medha Deshpande -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080701/7d4922da/attachment.html From brian.walsh at AMEC.COM Tue Jul 1 09:20:56 2008 From: brian.walsh at AMEC.COM (Walsh, Brian J) Date: Tue, 1 Jul 2008 09:20:56 -0400 Subject: Creating Shaded Plots from Text File Message-ID: Hello, I am looking for a way to use GrADS to plot up some climate data from a text file. For example, I have some average snowfall amounts for a region and I would like to create a shaded contoured image using this data. How hard is this to do and could someone send me some sample scripts/data files, etc. to get me started. Thanks in advance, Brian Walsh The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. Its contents (including any attachments) may contain confidential and/or privileged information. If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080701/ca7d511e/attachment.html From dasilva at ALUM.MIT.EDU Tue Jul 1 10:04:38 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Tue, 1 Jul 2008 10:04:38 -0400 Subject: Creating Shaded Plots from Text File In-Reply-To: Message-ID: On Tue, Jul 1, 2008 at 9:20 AM, Walsh, Brian J wrote: > Hello, > > I am looking for a way to use GrADS to plot up some climate data from a > text file. For example, I have some average snowfall amounts for a region > and I would like to create a shaded contoured image using this data. How > hard is this to do and could someone send me some sample scripts/data files, > etc. to get me started. > Send me a sample file and I'll see what I can do. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080701/0feaf149/attachment.html From vpk1 at LEICESTER.AC.UK Tue Jul 1 18:06:30 2008 From: vpk1 at LEICESTER.AC.UK (Kanawade, V.P.) Date: Tue, 1 Jul 2008 23:06:30 +0100 Subject: plotting error Message-ID: Dear All, I am running GrADS on Windows vista. I am trying to plot 5 deg. gridded CO data at 150hPa level. I have created binary file (attached) and ctl file (attached). I am able to open ctl file but when try to plot data it gives me erro. I am unable to figure out whether there is problem with binary file or running GrADS on windows. Any suggestion in this regard would be much appreciated. thank you in adv. Please see below for error message - Landscape mode? (no for portrait): reinit GX Package Initialization: Size = 11 8.5 ga-> set display color white ga-> c ga-> open C:\GrADS19\data\CO150.ctl Scanning description file: C:\GrADS19\data\CO150.ctl Data file grads.dat is open as file 1 LON set to 0 360 LAT set to -60 55 LEV set to 1 1 Time values set: 1948:1:1:0 1948:1:1:0 ga-> query file 1 File 1 : land portion at 5 deg resolution Descriptor: C:\GrADS19\data\CO150.ctl Binary: grads.dat Type = Gridded Xsize = 72 Ysize = 24 Zsize = 1 Tsize = 1 Number of Variables = 1 co 1 99 elevation in m ga-> d co Low Level I/O Error: Read error on data file Data file name = grads.dat Error reading 72 bytes at location 432 Data Request Error: Error for variable 'co' Error ocurred at column 1 DISPLAY error: Invalid expression Expression = co ga-> -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080701/f1818667/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: CO150.ctl Type: application/octet-stream Size: 225 bytes Desc: CO150.ctl Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080701/f1818667/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: grads.dat Type: application/octet-stream Size: 1728 bytes Desc: grads.dat Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080701/f1818667/attachment-0001.obj From maldmm5 at GMAIL.COM Tue Jul 1 22:34:34 2008 From: maldmm5 at GMAIL.COM (Mald MM5) Date: Wed, 2 Jul 2008 14:34:34 +1200 Subject: How to use station data in grads Message-ID: Hi there, I have station data with me. lat lon and the meteorological fields. I want to use this data in grads. Its only one level (surface) and daily data for about 30 years. How can I use this data in grads. Anybody has any programs to grid such kind of data to make grads understandable format and how would the ctl file be? Thank you Mald > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/d17566cf/attachment.html From dasilva at ALUM.MIT.EDU Wed Jul 2 01:31:46 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Wed, 2 Jul 2008 01:31:46 -0400 Subject: plotting error In-Reply-To: <22BD51C71909654C8B94797767D468B2317468E836@EXC-MBX1.cfs.le.ac.uk> Message-ID: On Tue, Jul 1, 2008 at 6:06 PM, Kanawade, V.P. wrote: > Dear All, > > I am running GrADS on Windows vista. I am trying to plot 5 deg. gridded CO > data at 150hPa level. I have created binary file (attached) and ctl file > (attached). I am able to open ctl file but when try to plot data it gives me > erro. I am unable to figure out whether there is problem with binary file or > running GrADS on windows. > > Any suggestion in this regard would be much appreciated. thank you in adv. > Your file does not have the correct size. It should have 72*24*4 = 6,912 bytes and not the 2K it has. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/b3ae094c/attachment.html From eduardo.agosta at GMAIL.COM Wed Jul 2 05:38:25 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Wed, 2 Jul 2008 06:38:25 -0300 Subject: linear trend with monthly series Message-ID: Dear Folks, I am trying to extimate linear trend in an interannual series of one-monthly data from a multy-year monthly data file. I can estimate properly the linear trend using the script below, if the data file has the one-monthly data series previously separately. The problem is when using the full monthly multy-year data file. How can I just select one given month from the data , for instance June, and estimate the linear trend for the June-series alone, having a multy-year data file with all the months of the year? I have tried using the function monmask (with argument 6), but I cannot solve it properly. Any suggestion I will much appreciate. Eduardo. reinit' 'sdfopen d:\anterior\era\mes.nc' 'set clopts -1 -1 0.16' 'set mproj sps' k=70 'set grads off' 'set lev 'k'' 'set mproj sps' 'set lat -90 0' 'set lon -180 180' 'set lat -60' 'set lon -180' 'set t 1 132' * contruyo variable lineal 1, 2, ..., 276 to para calcular m y c de y = mx + c 'xvar = 100*t*0 + 1' 'xvar = const(xvar,1,-u)' 'xvar = 133 - sum(xvar,t+0,t=132)' 'set lat -90 0' 'set lon -180 180' 'set t 1' 'define coeff = tregr(xvar,100*t,t=1, t=132)' 'define tave = ave(100*t, t=1, t=132)' 'define xvarave = ave(xvar, t=1, t=132)' 'define trend = coeff*(xvar - xvarave) + tave' -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/fcd4162c/attachment.html From mau_imh at VKTTV.EDU.VN Wed Jul 2 01:48:35 2008 From: mau_imh at VKTTV.EDU.VN (Nguyen Dang Mau) Date: Wed, 2 Jul 2008 12:48:35 +0700 Subject: Create text file from binary file In-Reply-To: <77fcd6b20807012231pe37b080m8c67df109e3b8f5e@mail.gmail.com> Message-ID: Dear Sir/Madam, I am trying to create a series data (arranged by years: 1971-2000) from a binary file. At here, I have some average temperature and precipitation files; they are output of RegCM3 model. So I hope to use GRADS to create text file using this file. So How I do and if someone have some scripts, ect. Please give me some advices or send your scripts to me. Thank you very much, All the best Nguyen Dang Mau Center for Meteorology and Climatology Vietnam Institute of Meteorology, Hydrology and Environment 5/62 Nguyen Chi Thanh Rd., Dongda, Hanoi, Vietnam -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/b4f32568/attachment.html From theedge981 at GMAIL.COM Wed Jul 2 12:04:41 2008 From: theedge981 at GMAIL.COM (Dan Leins) Date: Wed, 2 Jul 2008 11:04:41 -0500 Subject: Cross section using maskout for terrain Message-ID: All, I am attempting to draw a simple cross section showing temperature across a model domain and I would like to maskout the area below ground. I have seen some traffic on this issue. I see that I should be able to define the terrain by subtracting pressfc/100 from lev, with all negative values being areas that are underground. Then using maskout, I can do the following: d maskout(tmpprs,t1pm1-lev) This gives me something in the ballpark of what I'm looking for (image1.png), but the shaded area looks very gridded in nature. When I simply plot a contoured plot of 't1pm1-lev', it doesn't look gridded at all and appears nice and smooth (image2.png). The model grid spacing is 4km and I'm covering about 500km in my cross section so I'm guessing things should appear relatively smooth. Why does maskout produce results that look so gridded? Is there any other way I can maskout the underground areas but still have a nice crisp outline of the terrain? Thanks! Dan Leins -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/b6c0f42f/attachment.html -------------- next part -------------- A non-text attachment was scrubbed... Name: image1.png Type: image/png Size: 9931 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080702/b6c0f42f/attachment.png -------------- next part -------------- A non-text attachment was scrubbed... Name: image2.png Type: image/png Size: 2961 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080702/b6c0f42f/attachment-0001.png From Bill.Bua at NOAA.GOV Wed Jul 2 12:18:50 2008 From: Bill.Bua at NOAA.GOV (Bill Bua) Date: Wed, 2 Jul 2008 12:18:50 -0400 Subject: Cross section using maskout for terrain In-Reply-To: <46e821750807020904y4e850ffcvb7009b4922acd4c7@mail.gmail.com> Message-ID: Dan -- I suspect it's because you're using lev, which is discretized into layers, and shading, resulting in 'boxy' plots to the nearest level set up in your control file and "masking out" anything that has a value of your variable t1pm1-lev < 0. The plot of t1pm1-lev is smooth because of the plotting routine used for contours. I'm not sure if this will work better, but perhaps you can use the model topographic height if it's available and subtract the height of the 'lev' pressure surfaces as an alternative. Your mask might look more like the topography than the levels if you do so. Hope this helps. Bill Bua Dan Leins wrote: > All, > > I am attempting to draw a simple cross section showing temperature > across a model domain and I would like to maskout the area below > ground. I have seen some traffic on this issue. I see that I should be > able to define the terrain by subtracting pressfc/100 from lev, with > all negative values being areas that are underground. Then using > maskout, I can do the following: > > d maskout(tmpprs,t1pm1-lev) > > This gives me something in the ballpark of what I'm looking for > (image1.png), but the shaded area looks very gridded in nature. When I > simply plot a contoured plot of 't1pm1-lev', it doesn't look gridded > at all and appears nice and smooth (image2.png). The model grid > spacing is 4km and I'm covering about 500km in my cross section so I'm > guessing things should appear relatively smooth. > > Why does maskout produce results that look so gridded? Is there any > other way I can maskout the underground areas but still have a nice > crisp outline of the terrain? > > Thanks! > Dan Leins > > > ------------------------------------------------------------------------ > > > ------------------------------------------------------------------------ > From eduardo.agosta at GMAIL.COM Wed Jul 2 13:26:02 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Wed, 2 Jul 2008 14:26:02 -0300 Subject: Help with linear trend, please! Message-ID: Dear Grads' Users, I am trying to extimate linear trend in an interannual series of one-monthly data from a multy-year monthly data file. I can estimate properly the linear trend using the script below, if the data file has the one-monthly data series previously separately. The problem is when using the full monthly multy-year data file. How can I just select one given month from the data , for instance June, and estimate the linear trend for the June-series alone, having a multy-year data file with all the months of the year? I have tried using the function monmask (with argument 6), but I cannot solve it properly. Any suggestion I will much appreciate. Eduardo. reinit' 'sdfopen d:\anterior\era\mes.nc' 'set clopts -1 -1 0.16' 'set mproj sps' k=70 'set grads off' 'set lev 'k'' 'set mproj sps' 'set lat -90 0' 'set lon -180 180' 'set lat -60' 'set lon -180' 'set t 1 132' * contruyo variable lineal 1, 2, ..., 276 to para calcular m y c de y = mx + c 'xvar = 100*t*0 + 1' 'xvar = const(xvar,1,-u)' 'xvar = 133 - sum(xvar,t+0,t=132)' 'set lat -90 0' 'set lon -180 180' 'set t 1' 'define coeff = tregr(xvar,100*t,t=1, t=132)' 'define tave = ave(100*t, t=1, t=132)' 'define xvarave = ave(xvar, t=1, t=132)' 'define trend = coeff*(xvar - xvarave) + tave' -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/9cb7c19c/attachment.html From theedge981 at GMAIL.COM Wed Jul 2 15:06:35 2008 From: theedge981 at GMAIL.COM (Dan Leins) Date: Wed, 2 Jul 2008 14:06:35 -0500 Subject: Cross section using maskout for terrain In-Reply-To: <486BAA6A.5030309@noaa.gov> Message-ID: Bill, Thanks for the suggestion. I attempted to set up a mask using the difference between hgtprs and hgtsfc. Unfortunately this produced similar "boxy" results like the images I attached initially. What routine does maskout use to "mask out" the negative areas? Is there a way for it to use the same smooth routine used when contouring? Dan On Wed, Jul 2, 2008 at 11:18 AM, Bill Bua wrote: > Dan -- > > I suspect it's because you're using lev, which is discretized into > layers, and shading, resulting in 'boxy' plots to the nearest level set > up in your control file and "masking out" anything that has a value of > your variable t1pm1-lev < 0. The plot of t1pm1-lev is smooth because of > the plotting routine used for contours. > > I'm not sure if this will work better, but perhaps you can use the model > topographic height if it's available and subtract the height of the > 'lev' pressure surfaces as an alternative. Your mask might look more > like the topography than the levels if you do so. > > Hope this helps. > > Bill Bua > > Dan Leins wrote: > >> All, >> >> I am attempting to draw a simple cross section showing temperature >> across a model domain and I would like to maskout the area below >> ground. I have seen some traffic on this issue. I see that I should be >> able to define the terrain by subtracting pressfc/100 from lev, with >> all negative values being areas that are underground. Then using >> maskout, I can do the following: >> >> d maskout(tmpprs,t1pm1-lev) >> >> This gives me something in the ballpark of what I'm looking for >> (image1.png), but the shaded area looks very gridded in nature. When I >> simply plot a contoured plot of 't1pm1-lev', it doesn't look gridded >> at all and appears nice and smooth (image2.png). The model grid >> spacing is 4km and I'm covering about 500km in my cross section so I'm >> guessing things should appear relatively smooth. >> >> Why does maskout produce results that look so gridded? Is there any >> other way I can maskout the underground areas but still have a nice >> crisp outline of the terrain? >> >> Thanks! >> Dan Leins >> >> >> ------------------------------------------------------------------------ >> >> >> ------------------------------------------------------------------------ >> >> -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/5fab1c99/attachment.html From eduardo.agosta at GMAIL.COM Wed Jul 2 17:06:07 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Wed, 2 Jul 2008 18:06:07 -0300 Subject: Help with linear trend, please! In-Reply-To: <3367a2e40807021026j6140ebd5i7f35303bd7a1cbfa@mail.gmail.com> Message-ID: Hi you all! I have found myself the answer to my inquiry. For this may help someone else, the way to estimate a linear trend (or regression) on one-monthly series from a multi-year multi-monthly data file, wihtout having to seperate the chosen month interannual series into a single file, is just to perform the regression calculus for that month! I help myself from other previous prompts from the GRADSUSRS Archive. Here it is the script that can do that: 'sdfopen file.nc' 'set lev 'k'' 'set lat -90 0' 'set lon -180 180' 'set t 1 276' 'define yvar = z' *this is the variable on which to estimate the linear trend * the lineal variable 1, 2, ..., 276 to estimate m and c of y = mx + c * y = mx + c * y = amp*x + ave 'set lat -60' 'set lon -180' 'set t 1 276' 'xvar = yvar*0 + 1' 'xvar = const(xvar,1,-u)' 'xvar = 277 - sum(xvar,t+0,t=276)' 'set lat -90 0' 'set lon -180 180' 'set t 1' 'define p = sum(xvar,t=10,t=276,12)' * starttime=10, supposing you want to estimate October linear trend 'define q = sum(yvar,t=10,t=276,12)' 'define r = sum(xvar*yvar,t=10,t=276,12)' 'define s = sum(xvar*xvar,t=10,t=276,12)' 'define rn = ave(xvar*yvar,t=10,t=276,12)' 'define n = r/rn' 'define amp = (n*r - p*q)/(n*s - p*p)' *this is m 'define ave = (q - amp*p)/n' 'set t 1 276' 'monmask 10' * this is for month October, you need the monmas.gs script from the Archive 'define trend=amp*xvar*mmask + ave' * this is the final linear trend for Oct alone That`s all!! Bye. Eduardo. 2008/7/2 Eduardo Agosta Scarel : > Dear Grads' Users, > > I am trying to extimate linear trend in an interannual series of > one-monthly data from a multy-year monthly data file. > I can estimate properly the linear trend using the script below, if > the data file has the one-monthly data series previously separately. The > problem is when using the full monthly multy-year data file. > > How can I just select one given month from the data , for instance June, > and estimate the linear trend for the June-series alone, having a multy-year > data file with all the months of the year? > I have tried using the function monmask (with argument 6), but I cannot > solve it properly. > > Any suggestion I will much appreciate. > > Eduardo. > > reinit' > 'sdfopen d:\anterior\era\mes.nc' > 'set clopts -1 -1 0.16' > 'set mproj sps' > k=70 > > 'set grads off' > 'set lev 'k'' > 'set mproj sps' > 'set lat -90 0' > 'set lon -180 180' > > 'set lat -60' > 'set lon -180' > 'set t 1 132' > * contruyo variable lineal 1, 2, ..., 276 to para calcular m y c de y = mx > + c > 'xvar = 100*t*0 + 1' > 'xvar = const(xvar,1,-u)' > 'xvar = 133 - sum(xvar,t+0,t=132)' > > 'set lat -90 0' > 'set lon -180 180' > 'set t 1' > > 'define coeff = tregr(xvar,100*t,t=1, t=132)' > 'define tave = ave(100*t, t=1, t=132)' > 'define xvarave = ave(xvar, t=1, t=132)' > 'define trend = coeff*(xvar - xvarave) + tave' > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/642e3810/attachment.html From Charles.Seman at NOAA.GOV Wed Jul 2 17:19:56 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Wed, 2 Jul 2008 17:19:56 -0400 Subject: Cross section using maskout for terrain In-Reply-To: <46e821750807021206n2735ee6cl93d618d03241c31f@mail.gmail.com> Message-ID: Dan, Found a file containing a script "orog.gs" posted to the Listserv by Wesley Ebisuzaki (attached). I hope this helps, Chuck Dan Leins wrote: > Bill, > > Thanks for the suggestion. I attempted to set up a mask using the > difference between hgtprs and hgtsfc. Unfortunately this produced > similar "boxy" results like the images I attached initially. > > What routine does maskout use to "mask out" the negative areas? Is > there a way for it to use the same smooth routine used when contouring? > > Dan > > On Wed, Jul 2, 2008 at 11:18 AM, Bill Bua > wrote: > > Dan -- > > I suspect it's because you're using lev, which is discretized into > layers, and shading, resulting in 'boxy' plots to the nearest > level set > up in your control file and "masking out" anything that has a value of > your variable t1pm1-lev < 0. The plot of t1pm1-lev is smooth > because of > the plotting routine used for contours. > > I'm not sure if this will work better, but perhaps you can use the > model > topographic height if it's available and subtract the height of the > 'lev' pressure surfaces as an alternative. Your mask might look more > like the topography than the levels if you do so. > > Hope this helps. > > Bill Bua > > Dan Leins wrote: > > All, > > I am attempting to draw a simple cross section showing temperature > across a model domain and I would like to maskout the area below > ground. I have seen some traffic on this issue. I see that I > should be > able to define the terrain by subtracting pressfc/100 from > lev, with > all negative values being areas that are underground. Then using > maskout, I can do the following: > > d maskout(tmpprs,t1pm1-lev) > > This gives me something in the ballpark of what I'm looking for > (image1.png), but the shaded area looks very gridded in > nature. When I > simply plot a contoured plot of 't1pm1-lev', it doesn't look > gridded > at all and appears nice and smooth (image2.png). The model grid > spacing is 4km and I'm covering about 500km in my cross > section so I'm > guessing things should appear relatively smooth. > > Why does maskout produce results that look so gridded? Is > there any > other way I can maskout the underground areas but still have a > nice > crisp outline of the terrain? > > Thanks! > Dan Leins > > > ------------------------------------------------------------------------ > > > ------------------------------------------------------------------------ > > -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: orography_on_vertical_cross_section Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/a1c9025e/attachment.pl From vovanhoa80 at YAHOO.COM Wed Jul 2 22:14:06 2008 From: vovanhoa80 at YAHOO.COM (Van Hoa Vo) Date: Wed, 2 Jul 2008 19:14:06 -0700 Subject: Create text file from binary file In-Reply-To: <20080702111751.017EE4F8@as2.cineca.com> Message-ID: Hi Hay viet 1 chuong trinh Fortran nho ma doc cai file binary do, neu da biet duoc cai control file cua fbinary file thi co nghia la da hieu duoc cau truc file do roi, viet mot chuong trinh Fortran mo cai file do ra voi cai tuy chon nhu luc tao file do (i.e ACCESS='DIRECT', FORM='BINARY', RECL=4*NX*NY), sau do doc tung record ra la duoc, can gi phai dung den Grads cho no mat cong Bye Hoa Van Vo Researcher, Research and Application Division National Center for Hydro-Meteorological Forecasting (NCHMF) Tel (+084) 9330942 Email: vovanhoa80 at yahoo.com --- On Wed, 7/2/08, Nguyen Dang Mau wrote: From: Nguyen Dang Mau Subject: Create text file from binary file To: GRADSUSR at LIST.CINECA.IT Date: Wednesday, July 2, 2008, 12:48 PM Dear Sir/Madam, I am trying to create a series data (arranged by years: 1971-2000) from a binary file. At here, I have some average temperature and precipitation files; they are output of RegCM3 model. So I hope to use GRADS to create text file using this file. So How ?I do and if someone have some scripts, ect? Please give me some advices or send your scripts to me. Thank you very much, All the best Nguyen Dang Mau Center for Meteorology and Climatology Vietnam Institute of Meteorology, Hydrology and Environment 5/62 Nguyen Chi Thanh Rd., Dongda, Hanoi, Vietnam ? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080702/42fa89d2/attachment.html From Herbert.termaat at WUR.NL Thu Jul 3 03:28:10 2008 From: Herbert.termaat at WUR.NL (Maat, Herbert ter) Date: Thu, 3 Jul 2008 09:28:10 +0200 Subject: skipping time Message-ID: Dear GrADS-users, I have a very simple question which is probably easy to answer as well. I have temperature data in a timeseries of, let's say, 120 timesteps (30 days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to display a selection of this data, namely only temperature at 00Z and skipping the other timesteps, thus display times 1, 4, 9, and so on. Is there a simple statement to do this in GrADS or do I need to end up writing a script for this? Kind regards, Herbert ter Maat -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080703/b991d2cb/attachment.html From d_akritidis at YAHOO.COM Thu Jul 3 06:22:27 2008 From: d_akritidis at YAHOO.COM (Akritidis Dimitrios) Date: Thu, 3 Jul 2008 03:22:27 -0700 Subject: greyscale output Message-ID: Dear all, I want to get my output (gif,png,...) on a greyscale. I set the display on greyscale but using the ''printim'' command my output is colored. Thanks in advance! -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080703/360d7605/attachment.html From sachinghude at TROPMET.RES.IN Thu Jul 3 07:12:30 2008 From: sachinghude at TROPMET.RES.IN (Sachin Ghude) Date: Thu, 3 Jul 2008 17:42:30 +0630 Subject: greyscale output In-Reply-To: <578093.26201.qm@web43134.mail.sp1.yahoo.com> Message-ID: Hi! I hope follwing commands will help you, 'set display color white' 'set display grey white' 'enable print C:filename.gm' 'print' 'disable print' '!gxeps C:filename.gm' 'printim C:filename.jpg' Sachin D Ghude ---------- Original Message ----------- From: Akritidis Dimitrios To: GRADSUSR at LIST.CINECA.IT Sent: Thu, 3 Jul 2008 03:22:27 -0700 Subject: greyscale output > Dear all, > > I want to get my output (gif,png,...) on a greyscale. I set the display on greyscale but using the ''printim'' command my output is colored. > > Thanks in advance! > > ------- End of Original Message ------- From rjoseph at ATMOS.UMD.EDU Thu Jul 3 11:43:34 2008 From: rjoseph at ATMOS.UMD.EDU (Renu R. Joseph) Date: Thu, 3 Jul 2008 11:43:34 -0400 Subject: Template option for files with Julian Day names In-Reply-To: <3367a2e40807021406x3408978dn63106722c58d59bf@mail.gmail.com> Message-ID: Hello All, I need suggestions/solutions on dealing with the problem of opening files with julian day included in the names. I have files that named with the following name format: m6s4_3hYYDDDHH.bin, where YY represents the last two digits of the years (01 is 2001); DDD is julian day; HH is the hour of the start of the measurement If grads allowed a julian day formats it could have been opened with "m6s4_3h%y2%jd%h200.bin" in the dset file location as follows. dset ^m6s4_3h%y2%jd3%h2.bin OPTIONS template OPTIONS big_endian OPTIONS yrev UNDEF -9999 TITLE Precipitation estimates XDEF 1440 LINEAR 0.125 0.25 YDEF 480 LINEAR -59.875 0.25 ZDEF 01 LEVELS 1 TDEF 1 LINEAR 00z01jan2001 3 hr VARS 1 per 1 99 (mm/hr) ENDVARS Since I cannot do that, I know that one of the options is to use DSET ^m6s4_3h%ch.bin CHSUB 1 1 m6s4_3h0100100.bin CHSUB 2 2 m6s4_3h0100100.bin .................................................. CHSUB 20448 20448 m6s4_3h0736521.bin OPTIONS template but that includes so many lines. Is there a way around this, unless I physically rename files to be consistent with grads format names? Thanks, -Renu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080703/df71cf35/attachment.html From mike.bosilovich at GMAIL.COM Thu Jul 3 12:24:01 2008 From: mike.bosilovich at GMAIL.COM (Mike Bosilovich) Date: Thu, 3 Jul 2008 12:24:01 -0400 Subject: skipping time In-Reply-To: <1FA9A5649C33644B85AE596DBAF16D7CDCCDD9@scomp0040.wurnet.nl> Message-ID: I came up with a way to do this. It is not elegant, but works, and I would be happy to see another option. The idea is to maskout all other times except the one you want to plot. However, grads won't let you take a single time variable, and expand that into a time series so that it can mask the other times. So, here goes, for 00Z: set time 18Z1Jul2004 00Z2Jul2004 define m1=1 set time 00Z2Jul2004 06Z2Jul2004 define m2=1 set time 00Z2Jul2004 18Z2jul2004 define m3=m2*m1 *m3=1 at 00Z2Jul, undef all others modify m3 diurnal * now m3=1 for all 00Z set time 00Z1Jul2004 18Z31jul2004 display var*m3 The display should have a time series of 00Z points. If you want those points connected then prior to display, set missconn on. Hope it helps, Mike On Thu, Jul 3, 2008 at 3:28 AM, Maat, Herbert ter wrote: > Dear GrADS-users, > > I have a very simple question which is probably easy to answer as well. > > I have temperature data in a timeseries of, let's say, 120 timesteps (30 > days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to > display a selection of this data, namely only temperature at 00Z and > skipping the other timesteps, thus display times 1, 4, 9, and so on. Is > there a simple statement to do this in GrADS or do I need to end up writing > a script for this? > > Kind regards, > Herbert ter Maat > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080703/80a2b10a/attachment.html From smcmillan at PLANALYTICS.COM Thu Jul 3 13:06:32 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Thu, 3 Jul 2008 12:06:32 -0500 Subject: skipping time In-Reply-To: <79aecae70807030924r158fa5c8ycad864bd57e1f141@mail.gmail.com> Message-ID: Herbert, You could also try using the "skip" command. Perhaps something like this will also work for you: set time 00z01jun2007 00z01jul2007 set missconn on (if connecting lines desired, as Mike suggested) d skip(air,4) Above should display data only at 00z timesteps (assuming 4x daily as in your sample). Stephen McMillan Mike Bosilovich Sent by: GRADSUSR at LIST.CINECA.IT 07/03/2008 11:24 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: skipping time I came up with a way to do this. It is not elegant, but works, and I would be happy to see another option. The idea is to maskout all other times except the one you want to plot. However, grads won't let you take a single time variable, and expand that into a time series so that it can mask the other times. So, here goes, for 00Z: set time 18Z1Jul2004 00Z2Jul2004 define m1=1 set time 00Z2Jul2004 06Z2Jul2004 define m2=1 set time 00Z2Jul2004 18Z2jul2004 define m3=m2*m1 *m3=1 at 00Z2Jul, undef all others modify m3 diurnal * now m3=1 for all 00Z set time 00Z1Jul2004 18Z31jul2004 display var*m3 The display should have a time series of 00Z points. If you want those points connected then prior to display, set missconn on. Hope it helps, Mike On Thu, Jul 3, 2008 at 3:28 AM, Maat, Herbert ter wrote: Dear GrADS-users, I have a very simple question which is probably easy to answer as well. I have temperature data in a timeseries of, let's say, 120 timesteps (30 days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to display a selection of this data, namely only temperature at 00Z and skipping the other timesteps, thus display times 1, 4, 9, and so on. Is there a simple statement to do this in GrADS or do I need to end up writing a script for this? Kind regards, Herbert ter Maat *************************************************** 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/20080703/6316cdb5/attachment.html From skimball at USOUTHAL.EDU Thu Jul 3 14:47:15 2008 From: skimball at USOUTHAL.EDU (Sytske Kimball) Date: Thu, 3 Jul 2008 13:47:15 -0500 Subject: skipping time - another question Message-ID: Now I have a question: what if you want to plot the entire timeseries, but with a different increment displayed on the x (time) axis then the one grads picks? I normally use 'set xlint', but that didn't work for real dates (like in the exmaple below). A look in the grads documentation confirmed that 'set xlint' does not work for a date/time axis. Does anyone know what to do instead? Thanks in advance, Sytske Sytske Kimball Associate Professor of Meteorology Dept. of Earth Sciences University of South Alabama Mobile, AL phone (251) 460-7031 fax (251) 460-7886 email skimball at usouthal.edu >>> Stephen R McMillan 07/03/08 12:08 PM >>> Herbert, You could also try using the "skip" command. Perhaps something like this will also work for you: set time 00z01jun2007 00z01jul2007 set missconn on (if connecting lines desired, as Mike suggested) d skip(air,4) Above should display data only at 00z timesteps (assuming 4x daily as in your sample). Stephen McMillan Mike Bosilovich Sent by: GRADSUSR at LIST.CINECA.IT 07/03/2008 11:24 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: skipping time I came up with a way to do this. It is not elegant, but works, and I would be happy to see another option. The idea is to maskout all other times except the one you want to plot. However, grads won't let you take a single time variable, and expand that into a time series so that it can mask the other times. So, here goes, for 00Z: set time 18Z1Jul2004 00Z2Jul2004 define m1=1 set time 00Z2Jul2004 06Z2Jul2004 define m2=1 set time 00Z2Jul2004 18Z2jul2004 define m3=m2*m1 *m3=1 at 00Z2Jul, undef all others modify m3 diurnal * now m3=1 for all 00Z set time 00Z1Jul2004 18Z31jul2004 display var*m3 The display should have a time series of 00Z points. If you want those points connected then prior to display, set missconn on. Hope it helps, Mike On Thu, Jul 3, 2008 at 3:28 AM, Maat, Herbert ter wrote: Dear GrADS-users, I have a very simple question which is probably easy to answer as well. I have temperature data in a timeseries of, let's say, 120 timesteps (30 days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to display a selection of this data, namely only temperature at 00Z and skipping the other timesteps, thus display times 1, 4, 9, and so on. Is there a simple statement to do this in GrADS or do I need to end up writing a script for this? Kind regards, Herbert ter Maat *************************************************** 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. *************************************************** From mike.bosilovich at GMAIL.COM Thu Jul 3 15:57:25 2008 From: mike.bosilovich at GMAIL.COM (Mike Bosilovich) Date: Thu, 3 Jul 2008 15:57:25 -0400 Subject: skipping time - another question In-Reply-To: <486CD861020000B900023DEE@gwia1.usouthal.edu> Message-ID: Documentation may well be harder than writing the code :-) I always wrote skip off for the time purpose discussed earlier because the docs say x,y.... However, sounds like you want set xlabs in the docs If you want a mean annual cycle (instead of the grads default) set xlabs JAN | FEB | MAR ... or diurnal cycle set xlabs 00Z | 03Z | 06Z .... Mike On Thu, Jul 3, 2008 at 2:47 PM, Sytske Kimball wrote: > Now I have a question: what if you want to plot the entire timeseries, but > with a > different increment displayed on the x (time) axis then the one grads > picks? I normally > use 'set xlint', but that didn't work for real dates (like in the exmaple > below). A look in > the grads documentation confirmed that 'set xlint' does not work for a > date/time axis. > Does anyone know what to do instead? > > Thanks in advance, > Sytske > > Sytske Kimball > Associate Professor of Meteorology > Dept. of Earth Sciences > University of South Alabama > Mobile, AL > phone (251) 460-7031 > fax (251) 460-7886 > email skimball at usouthal.edu > >>> Stephen R McMillan 07/03/08 12:08 PM >>> > Herbert, > You could also try using the "skip" command. Perhaps something like this > will also work for you: > > set time 00z01jun2007 00z01jul2007 > set missconn on (if connecting lines desired, as Mike suggested) > d skip(air,4) > > Above should display data only at 00z timesteps (assuming 4x daily as in > your sample). > > Stephen McMillan > > > > Mike Bosilovich > Sent by: GRADSUSR at LIST.CINECA.IT > 07/03/2008 11:24 AM > Please respond to > GRADSUSR at LIST.CINECA.IT > > > To > GRADSUSR at LIST.CINECA.IT > cc > > Subject > Re: skipping time > > > > > > > I came up with a way to do this. It is not elegant, but works, and I would > be happy to see another option. The idea is to maskout all other times > except the one you want to plot. However, grads won't let you take a > single time variable, and expand that into a time series so that it can > mask the other times. So, here goes, for 00Z: > > set time 18Z1Jul2004 00Z2Jul2004 > define m1=1 > set time 00Z2Jul2004 06Z2Jul2004 > define m2=1 > set time 00Z2Jul2004 18Z2jul2004 > define m3=m2*m1 *m3=1 at 00Z2Jul, undef all others > modify m3 diurnal * now m3=1 for all 00Z > > set time 00Z1Jul2004 18Z31jul2004 > display var*m3 > > The display should have a time series of 00Z points. If you want those > points connected then prior to display, set missconn on. > > Hope it helps, > > Mike > > On Thu, Jul 3, 2008 at 3:28 AM, Maat, Herbert ter > wrote: > Dear GrADS-users, > > I have a very simple question which is probably easy to answer as well. > > I have temperature data in a timeseries of, let's say, 120 timesteps (30 > days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to > display a selection of this data, namely only temperature at 00Z and > skipping the other timesteps, thus display times 1, 4, 9, and so on. Is > there a simple statement to do this in GrADS or do I need to end up > writing a script for this? > > Kind regards, > Herbert ter Maat > > > > > > *************************************************** > 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/20080703/bd3c64a4/attachment.html From sreeds82 at YAHOO.COM Fri Jul 4 05:35:38 2008 From: sreeds82 at YAHOO.COM (sreenivas desamsetti) Date: Fri, 4 Jul 2008 02:35:38 -0700 Subject: Finding terrain.deck properties using MMOUT_DOMAIN1 Message-ID: Dear sir, I previously created MMINPUT, LOWBDY and BDYOUT files using terrain.deck file for 2 domains. Now that terrain.deck file doesnot exist with me. Shall I found the co-ordinates, and some parameters which can be used in the terrain.deck file. If so please kindly disscuss with me. I will be thankful if I can get the same terrain.deck back. Thanking you D.Srinivas Desamsetti Srinivas Research Scholar, Department of Meteorology and Oceanography, Andhra University, Visakhapatnam - 530 003 91 9948530618 (M) sreeds82 at yahoo.com sreeds82 at rediffmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080704/82478c72/attachment.html From eduardo.agosta at GMAIL.COM Fri Jul 4 09:55:58 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Fri, 4 Jul 2008 10:55:58 -0300 Subject: different initial time data files Message-ID: Hi you all! How should I deal with data calculus when having two different initial time for the variable because they belong to two different nc files?? Thanks for some suggestion. Eduardo. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080704/1340d8ac/attachment.html From stegert at IFM.UNI-HAMBURG.DE Mon Jul 7 05:23:44 2008 From: stegert at IFM.UNI-HAMBURG.DE (Stegert) Date: Mon, 7 Jul 2008 11:23:44 +0200 Subject: Template option for files with Julian Day names In-Reply-To: <486CF3A6.9010907@atmos.umd.edu> Message-ID: Dear Renu, even though I work with GrADS for a longer time, there are still so many options and functions within GrADS I do not know (as the CHSUB I didn't consider yet, same as variables in the CTL file: "%y2"...). So, for GrADS itself I do not know a fine solution - and that's the point when I start a workaround. For your needs three ways come to my mind: 1) To edit the CTL with 20448 CHSUBs is impossible, yet not a hard effort within a batch script I guess. Within 3 loops the year, the JD and the time are incremented and the .bin filename is formed of these variables. 2) Opening one .bin file by another by editing the first and the 10th line of the CTL in a loop and opening (and closing!) it each time. 3) I personally prefer data sets including all necessary data in one file (as long not getting too complex). So, I would most obviously compile all .bin files in one file first A Fortran program would do it quite fast I guess... open(1,file='per.bin',access='direct',form='unformatted',recl=RecLen) Also with three loops of yr, jd and hr building up the filename in a variable "BinName" while incrementing "RecNumber" each time. open(2,file=BinName) DO yr=... DO mon=... DO hr=... read(2,*)per ! the right way to read depends on your file format... write(1,rec=RecNumber,iostat=ios)per close(2) ENDDO ENDDO ENDDO The CTL would then have TDEF 20448 LINEAR 00z01jan2001 3 hr and you could access it in normal GrADS way... Hope this helps. If you'd like detailed information on one of these skripts feel free to write again. Cheers, Christoph Renu R. Joseph schrieb: > Hello All, > > I need suggestions/solutions on dealing with the problem of opening > files with julian day included in the names. > > I have files that named with the following name format: m6s4_3hYYDDDHH.bin, > where YY represents the last two digits of the years (01 is 2001); DDD > is julian day; HH is the hour of the start of the measurement > > If grads allowed a julian day formats it could have been opened with > "m6s4_3h%y2%jd%h200.bin" in the dset file location as follows. > > dset ^m6s4_3h%y2%jd3%h2.bin > OPTIONS template > OPTIONS big_endian > OPTIONS yrev > UNDEF -9999 > TITLE Precipitation estimates > XDEF 1440 LINEAR 0.125 0.25 > YDEF 480 LINEAR -59.875 0.25 > ZDEF 01 LEVELS 1 > TDEF 1 LINEAR 00z01jan2001 3 hr > VARS 1 > per 1 99 (mm/hr) > ENDVARS > > Since I cannot do that, I know that one of the options is to use > DSET ^m6s4_3h%ch.bin > CHSUB 1 1 m6s4_3h0100100.bin > CHSUB 2 2 m6s4_3h0100100.bin > .................................................. > CHSUB 20448 20448 m6s4_3h0736521.bin > OPTIONS template > > but that includes so many lines. Is there a way around this, unless I > physically rename files to be consistent with grads format names? > > Thanks, > -Renu > > > -- Christoph Stegert: Ecological Modelling group Institut f?r Meereskunde - Universit?t Hamburg, Germany Zentrum f. Marine u. Atmosph?rische Wissenschaften ZMAW Bundesstr.53, 20146 Hamburg - Phone: +49-40/ 42838-7486 From skimball at USOUTHAL.EDU Mon Jul 7 10:22:13 2008 From: skimball at USOUTHAL.EDU (Sytske Kimball) Date: Mon, 7 Jul 2008 09:22:13 -0500 Subject: skipping time - another question In-Reply-To: <79aecae70807031257j57f5e051w62de5b509776028f@mail.gmail.com> Message-ID: Thanks! Seems you can't automate the label interval, you just have to give it the labels to write in. Oh well. Thanks for your help though. Sytske >>> Mike Bosilovich 7/3/2008 2:57 PM >>> Documentation may well be harder than writing the code :-) I always wrote skip off for the time purpose discussed earlier because the docs say x,y.... However, sounds like you want set xlabs in the docs If you want a mean annual cycle (instead of the grads default) set xlabs JAN | FEB | MAR ... or diurnal cycle set xlabs 00Z | 03Z | 06Z .... Mike On Thu, Jul 3, 2008 at 2:47 PM, Sytske Kimball wrote: Now I have a question: what if you want to plot the entire timeseries, but with a different increment displayed on the x (time) axis then the one grads picks? I normally use 'set xlint', but that didn't work for real dates (like in the exmaple below). A look in the grads documentation confirmed that 'set xlint' does not work for a date/time axis. Does anyone know what to do instead? Thanks in advance, Sytske Sytske Kimball Associate Professor of Meteorology Dept. of Earth Sciences University of South Alabama Mobile, AL phone (251) 460-7031 fax (251) 460-7886 email skimball at usouthal.edu >>> Stephen R McMillan 07/03/08 12:08 PM >>> Herbert, You could also try using the "skip" command. Perhaps something like this will also work for you: set time 00z01jun2007 00z01jul2007 set missconn on (if connecting lines desired, as Mike suggested) d skip(air,4) Above should display data only at 00z timesteps (assuming 4x daily as in your sample). Stephen McMillan Mike Bosilovich Sent by: GRADSUSR at LIST.CINECA.IT 07/03/2008 11:24 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: skipping time I came up with a way to do this. It is not elegant, but works, and I would be happy to see another option. The idea is to maskout all other times except the one you want to plot. However, grads won't let you take a single time variable, and expand that into a time series so that it can mask the other times. So, here goes, for 00Z: set time 18Z1Jul2004 00Z2Jul2004 define m1=1 set time 00Z2Jul2004 06Z2Jul2004 define m2=1 set time 00Z2Jul2004 18Z2jul2004 define m3=m2*m1 *m3=1 at 00Z2Jul, undef all others modify m3 diurnal * now m3=1 for all 00Z set time 00Z1Jul2004 18Z31jul2004 display var*m3 The display should have a time series of 00Z points. If you want those points connected then prior to display, set missconn on. Hope it helps, Mike On Thu, Jul 3, 2008 at 3:28 AM, Maat, Herbert ter wrote: Dear GrADS-users, I have a very simple question which is probably easy to answer as well. I have temperature data in a timeseries of, let's say, 120 timesteps (30 days times 4 points per day; 00z 06z 12z and 18z). Using GrADS I want to display a selection of this data, namely only temperature at 00Z and skipping the other timesteps, thus display times 1, 4, 9, and so on. Is there a simple statement to do this in GrADS or do I need to end up writing a script for this? Kind regards, Herbert ter Maat *************************************************** 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/20080707/f3862acb/attachment.html From brian.walsh at AMEC.COM Mon Jul 7 10:52:09 2008 From: brian.walsh at AMEC.COM (Walsh, Brian J) Date: Mon, 7 Jul 2008 10:52:09 -0400 Subject: Plotting Tropical Sub Domains Message-ID: Hello, I was looking for some help on plotting data from a moving hurricane/tropical domain. For instance, with the HWRF and GFDL hurricane models you can get a high resolution sub-domain centered on the storm. The domain of each time step varies, so when I concatenate the grib files together, GrADS uses the latitude/longitude domain from the 00 hour file for all of the time steps (or at least it keeps the latitude/longitude legends as the same numbers even with changing domains with each successive time step). Is there any easy solution to this problem or do I need to plot charts after each time step? Thanks, Brian Walsh The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. Its contents (including any attachments) may contain confidential and/or privileged information. If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080707/f4416049/attachment.html From dasilva at ALUM.MIT.EDU Mon Jul 7 23:43:54 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 7 Jul 2008 23:43:54 -0400 Subject: Template option for files with Julian Day names In-Reply-To: <486CF3A6.9010907@atmos.umd.edu> Message-ID: On Thu, Jul 3, 2008 at 11:43 AM, Renu R. Joseph wrote: > Hello All, > > I need suggestions/solutions on dealing with the problem of opening files > with julian day included in the names. > > I have files that named with the following name format: m6s4_3hYYDDDHH.bin, > where YY represents the last two digits of the years (01 is 2001); DDD is > julian day; HH is the hour of the start of the measurement > > If grads allowed a julian day formats it could have been opened with > "m6s4_3h%y2%jd%h200.bin" in the dset file location as follows. > > dset ^m6s4_3h%y2%jd3%h2.bin > OPTIONS template > OPTIONS big_endian > OPTIONS yrev > UNDEF -9999 > TITLE Precipitation estimates > XDEF 1440 LINEAR 0.125 0.25 > YDEF 480 LINEAR -59.875 0.25 > ZDEF 01 LEVELS 1 > TDEF 1 LINEAR 00z01jan2001 3 hr > VARS 1 > per 1 99 (mm/hr) > ENDVARS > > Since I cannot do that, I know that one of the options is to use > DSET ^m6s4_3h%ch.bin > CHSUB 1 1 m6s4_3h0100100.bin > CHSUB 2 2 m6s4_3h0100100.bin > .................................................. > CHSUB 20448 20448 m6s4_3h0736521.bin > OPTIONS template > > but that includes so many lines. Is there a way around this, unless I > physically rename files to be consistent with grads format names? > You basically outlined your options. The best solution is to include julian day support in grads templating, something that is badly needed (I also ran into such files before). In my case I had a good reason not to rename the files (other applications required the julian day type of name conventions), and sometimes making a copy of the files isn't an option (disk space considerations). Here are a couple of options which leave the original files alone: 1) Symlink each file to a grads friendly name in a separate directory, e.g., grads/m6s4_2007-02-05_12z.bin. You should be able to do this with a relatively simple script. 2) If you know C, implement %j3 support in GrADS and contribute a patch. In GrADS both v1.9 and v2, look at this function in gautil.c: /* Given a file name template and a dt structure, fill in to get the file name */ char *gafndt (char *fn, struct dt *dtim, struct dt *dtimi, float *vals, struct gachsub *pch1st, int t) { Good Luck, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080707/8cefab86/attachment.html From z3bao at UWATERLOO.CA Tue Jul 8 09:54:03 2008 From: z3bao at UWATERLOO.CA (Zhenhao Bao) Date: Tue, 8 Jul 2008 09:54:03 -0400 Subject: Template option for files with Julian Day names In-Reply-To: <77fcd6b20807072043y617810a7we7fb9d9b2bb669e6@mail.gmail.com> Message-ID: Zhenhao Bao (PhD) Dept. of Geography University of Waterloo 200 university avenue west Waterloo, Ontario, Canada n2l 3g1 z3bao at uwaterloo.ca Quoting Arlindo da Silva : > On Thu, Jul 3, 2008 at 11:43 AM, Renu R. Joseph > wrote: > >> Hello All, >> >> I need suggestions/solutions on dealing with the problem of opening files >> with julian day included in the names. >> >> I have files that named with the following name format: m6s4_3hYYDDDHH.bin, >> where YY represents the last two digits of the years (01 is 2001); DDD is >> julian day; HH is the hour of the start of the measurement >> >> If grads allowed a julian day formats it could have been opened with >> "m6s4_3h%y2%jd%h200.bin" in the dset file location as follows. >> >> dset ^m6s4_3h%y2%jd3%h2.bin >> OPTIONS template >> OPTIONS big_endian >> OPTIONS yrev >> UNDEF -9999 >> TITLE Precipitation estimates >> XDEF 1440 LINEAR 0.125 0.25 >> YDEF 480 LINEAR -59.875 0.25 >> ZDEF 01 LEVELS 1 >> TDEF 1 LINEAR 00z01jan2001 3 hr >> VARS 1 >> per 1 99 (mm/hr) >> ENDVARS >> >> Since I cannot do that, I know that one of the options is to use >> DSET ^m6s4_3h%ch.bin >> CHSUB 1 1 m6s4_3h0100100.bin >> CHSUB 2 2 m6s4_3h0100100.bin >> .................................................. >> CHSUB 20448 20448 m6s4_3h0736521.bin >> OPTIONS template >> >> but that includes so many lines. Is there a way around this, unless I >> physically rename files to be consistent with grads format names? >> > > > You basically outlined your options. The best solution is to include julian > day support in grads templating, something that is badly needed (I also ran > into such files before). In my case I had a good reason not to rename the > files (other applications required the julian day type of name conventions), > and sometimes making a copy of the files isn't an option (disk space > considerations). Here are a couple of options which leave the original files > alone: > > 1) Symlink each file to a grads friendly name in a separate directory, e.g., > grads/m6s4_2007-02-05_12z.bin. You should be able to do this with a > relatively simple script. > > 2) If you know C, implement %j3 support in GrADS and contribute a patch. In > GrADS both v1.9 and v2, look at this function in gautil.c: > > /* Given a file name template and a dt structure, fill in to get the file > name */ > char *gafndt (char *fn, struct dt *dtim, struct dt *dtimi, float *vals, > struct gachsub *pch1st, int t) { > > Good Luck, > > Arlindo > > > -- > Arlindo da Silva > dasilva at alum.mit.edu > From jhuddleston at HUGHES.NET Wed Jul 9 15:17:19 2008 From: jhuddleston at HUGHES.NET (John Huddleston) Date: Wed, 9 Jul 2008 13:17:19 -0600 Subject: Testers for the windows version of grads-2.0.a2 supporting GRIB2 Message-ID: Hi I've built a version of grads-2.0.a2 for windows. Initially it will require a Cygwin implementation on your PC. There is not much documentation for the new features; however, I've tested it against the older ones. If anyone is interested in testing this version of GrADS on Windows I can send you a copy of the EXE. John Huddleston Easy Computer Solutions -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080709/1397183a/attachment.html From antoyc at YAHOO.COM Wed Jul 9 23:07:22 2008 From: antoyc at YAHOO.COM (Antoy Chang) Date: Wed, 9 Jul 2008 20:07:22 -0700 Subject: Grads Library error? In-Reply-To: Message-ID: Hi everyone, Can someone assist me on this problem. I cannot plot the map?---> error on font*.dat, lores.dat and error in?local/usr/lib?appear.By the way, i'm using grads 1.8 win32e for MS windows PC. Your assistance on?my vague problem is highly appreciated. Thank you. KC -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080709/83e372a0/attachment.html From rafeh at DREAMSCAPEIT.COM Thu Jul 10 14:30:35 2008 From: rafeh at DREAMSCAPEIT.COM (Rafe Hargrove) Date: Thu, 10 Jul 2008 20:30:35 +0200 Subject: US and Mexico Precipitation Gridded Analysis Message-ID: I am trying to analyze recent weather precipitation posted by NOAA using the Grads program. NOAA posts daily files at the location below: ftp://ftp.cpc.ncep.noaa.gov/precip/wd52ws/us-mex/realtime/ I have been able to view the data for a given days without issue. I would like to identify accumulated precipitation by month since February 1, 2008. My ctl file is shown below: dset %y4%m2%d2.ll title US-Mexico Daily Precip Analysis options template big_endian sequential undef -9999. xdef 321 linear -140.00 0.25 ydef 201 linear 10.00 0.25 zdef 1 linear 1 1 tdef 29 linear 01feb2008 1dy vars 1 p 0 0 Total Precip (in/day) ENDVARS Here is what I get within Grads: Institute for Global Environment and Society All Rights Reserved Config: v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats printim Issue 'q config' command for more information. Landscape mode? (no for portrait): no GX Package Initialization: Size = 8.5 11 use c:\prf\griddeddatga-> use c:\prf\griddeddata\rafemonth c:\prf\griddeddata\rafemonth.ctl has been opened as File 1 and is the default file ga-> d p Cannot contour grid - all undefined values ga-> query file File 1 : US-Mexico Daily Precip Analysis Descriptor: c:\prf\griddeddata\rafemonth.ctl Binary: %y4%m2%d2.ll Type = Gridded Xsize = 321 Ysize = 201 Zsize = 1 Tsize = 29 Number of Variables = 1 p 0 0 Total Precip (in/day) ga-> query time Time = 00Z01FEB2008 to 00Z01FEB2008 Fri to Fri ga-> I can't seem to get the time sequence to work as it should and I'm a little unsure about the error message related to contour grids. My objective is to identify and document monthly recorded rainfall by grid (.25X.25 degrees) I would gladly pay for your assistance. Thanks, Rafe Hargrove 325-725-3036 From dasilva at ALUM.MIT.EDU Thu Jul 10 21:47:28 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Thu, 10 Jul 2008 21:47:28 -0400 Subject: convertion binary In-Reply-To: <154587.1263.qm@web31603.mail.mud.yahoo.com> Message-ID: On Thu, Jul 10, 2008 at 8:49 AM, wilfredo burgoa wrote: > Dear Arlindo, > > Please, with the option Set gxout fwrite manages to extract data in binary > (GrADS) and now would want to turn those data to a format for example to be > visualized in a spreadsheet (Excel). > Off hand the only thing I think of right now is to write a program to convert this binary file into ASCII. In this case, skip the fwrite altogether and use "set gxout print" to create ASCII directly from GrADS; look at the user guide for instructions. If you are an experienced C programmer, another possibility is look for a C library that writes in Excel binary format and write a conversion utility. This can also be done in Python using PyGrADS and the python excel tools. However, I've never done it myself (I do not use excel). Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080710/3f5c3dcb/attachment.html From fmashingia at YAHOO.COM Thu Jul 10 22:51:57 2008 From: fmashingia at YAHOO.COM (Fredrick Mashingia) Date: Thu, 10 Jul 2008 19:51:57 -0700 Subject: convertion binary Message-ID: On Thu, Jul 10, 2008 at 8:49 AM, wilfredo burgoa wrote: Dear Arlindo, ?Please,?with the option Set gxout fwrite manages to extract data in binary (GrADS) and now would want to turn those data to a format for example to be visualized in a?spreadsheet (Excel). Off hand the only thing I think of right now is to write a program to convert this binary file into ASCII. In this case, skip the fwrite altogether and use "set gxout print" to create ASCII directly from GrADS; look at the user guide for instructions. If you are an experienced C programmer, another possibility is look for a C library that writes in Excel binary format and write a conversion utility. This can also be done in Python using PyGrADS and the python excel tools. ?However, I've never done it myself (I do not use excel). ?Hello Wilfred, Arlindo is right with both options.However,the first option is the simplest and the second option is not necessary for your case. You can write FORTRAN program to?convert?your binary file into?ASCII then open the excel worksheet and load the ASCII file. Text imort wizard will appear,choose the file type that best describe your data,set the delimiters your data contains and finally set column data format. Otherwise if you use the "set gxout print" to create ASCII directly from GrADS then you can copy and paste the data into excell worksheet. In short,excel supports ASCII fomart. I hope that helps Mashingia ------------------------------------------------------ MASHINGIA, Fredrick A. J. B.Sc., M.Sc. Eng. (UDSM) Republic Of Tanzania. Phone: ?+255?713?619?561? E-mail: fmashingia at yahoo.com ------------------------------------------------------ -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080710/e7de5f4f/attachment.html From laznik at PRO-VREME.NET Fri Jul 11 03:00:44 2008 From: laznik at PRO-VREME.NET (Jost Laznik) Date: Fri, 11 Jul 2008 09:00:44 +0200 Subject: write function Message-ID: Dear users I would like to use write function, write (filename, record <, append>) What should be written in place where record and append is. Or i should left those strings unchanged? thanks jost From bhv at TROPMET.RES.IN Fri Jul 11 03:29:06 2008 From: bhv at TROPMET.RES.IN (Bakshi Vaid) Date: Fri, 11 Jul 2008 13:59:06 +0630 Subject: HI ALL In-Reply-To: <20080711070129.9D07F2002F@mx2.cineca.it> Message-ID: Dear All, I would like to compute velocity potential. Does anybody have a script to compute it?. Please send. Thank you in advance. Best Regards Hardeep V. From abbasmofidi at YAHOO.COM Fri Jul 11 03:42:10 2008 From: abbasmofidi at YAHOO.COM (abbas mofidi) Date: Fri, 11 Jul 2008 00:42:10 -0700 Subject: Output AScii file! Message-ID: Dear GrADS Users, I need to get the output of GrADS in text format. My problem is that the output AScii file is arranged in reverse order. I would be very grateful if someone can help me. Best regards, A. mofidi this is my script: 'reinit' 'sdfopen ncdata\Caspian-winter\AVHRR\sst4-path-eot.20011202.nc' 'sdfopen ncdata\Caspian-winter\AVHRR\sst4-path-eot.20011201.nc' ??? 'set dfile 1' ??? 'set lon 50.125 50.625' ??? 'set lat 38.125 39.125' ??? 'set t 1' ??? 'set z 1' ??? 'define S1 = sst' ??? 'set dfile 2' ??? 'set lon 50.125 50.625' ??? 'set lat 38.125 39.125' ??? 'set t 1' ??? 'set z 1' ??? 'define S2 = sst' ??? 'set gxout print' ??? 'set prnopts %g 3 1' ??? 'display S1-S2' ??? rc=write(DSST2Dec01.txt,result) ??? rc=close(DSST2Dec01.txt) return ? This is output .txt file:? Printing Grid -- 15 Values -- Undef = -999? ??????? ?-0.13? -0.12? -0.12 ???????? -0.14? -0.13? -0.13 ???????? -0.14? -0.13? -0.13 ???????? -0.139 -0.11? -0.09 ???????? -0.12? -0.24? -0.139 The latitude and Longitude of above output file is: Lat/Lon?? 50.125?? 50.375?? 50.625 ? 38.125?? -0.13????? -0.12???? -0.12 38.375?? -0.14????? -0.13???? -0.13 38.625?? -0.14????? -0.13???? -0.13 38.875?? -0.139???? -0.11???? -0.09 39.125?? -0.12????? -0.24???? -0.139 While, it should be like this: Lat/Lon?? 50.125??? 50.375???50.625 39.125?? -0.12?????? -0.24???? -0.139 38.875?? -0.139??? ?-0.11???? -0.09 38.625?? -0.14?????? -0.13???? -0.13 38.375?? -0.14????? ?-0.13???? -0.13 38.125?? -0.13????? ?-0.12???? -0.12 ? -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080711/0558b987/attachment.html From isilda at UEVORA.PT Fri Jul 11 05:16:10 2008 From: isilda at UEVORA.PT (Isilda da Cunha Menezes) Date: Fri, 11 Jul 2008 11:16:10 +0200 Subject: opening gradsda Message-ID: Hello! I installed the new version of Linux for grads (GrADS 2.0.a2) and I use it to open a file .nc (netcdf), but I am having trouble opening the gradsdap. I wrote, as I am accustomed to using the Linux version GrADS 1.8s11, gradsdap -lc "sdfopen air.nc" but got an error: gradsdap: error while loading shared libraries: libssl.so.4: cannot open shared object file: No such file or directory isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ gradsdap gradsdap: error while loading shared libraries: libssl.so.4: cannot open shared object file: No such file or directory isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ gradsdap gradsdap: error while loading shared libraries: libssl.so.4: cannot open shared object file: No such file or directory isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ grads and I have the same error when I try to open only the gradsdap. But the grads opens: isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ grads Grid Analysis and Display System (GrADS) Version 2.0.a2 Copyright (c) 1988-2008 by Brian Doty and the Institute for Global Environment and Society (IGES) GrADS comes with ABSOLUTELY NO WARRANTY See file COPYRIGHT for more information Config: v2.0.a2 little-endian readline printim grib2 netcdf hdf4-sds Issue 'q config' command for more information. Landscape mode? ('n' for portrait): -l GX Package Initialization: Size = 11 8.5 ga-> quit No hardcopy metafile open GX package terminated isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ I am committing a mistake to open the gradsdap or was poorly installed? Isilda From Patrick.Eriksson at FIMR.FI Fri Jul 11 09:37:53 2008 From: Patrick.Eriksson at FIMR.FI (Patrick B. Eriksson) Date: Fri, 11 Jul 2008 16:37:53 +0300 Subject: Patrick Eriksson/FIMR/FI is out of the office. Message-ID: I will be out of the office starting 11.07.2008 and will not return until 11.08.2008. Jag ?r p? semester under ovan n?mnda tid. Olen lomalla y.m. aikana. From mapyli at YORKU.CA Fri Jul 11 12:30:25 2008 From: mapyli at YORKU.CA (PY Li) Date: Fri, 11 Jul 2008 18:30:25 +0200 Subject: Performing remote analysis on Windows Message-ID: Hi, I'm a newbie, and just wondering if it's possible to perform remote analysis on Windows on the data on nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_DAILY I've tried the following operation found in GrADS-DODS Server - User's Guide (http://nomads.ncdc.noaa.gov:9091/user.html): ga-> sdfopen http://motherlode.ucar.edu:9090/dods/_expr_{C20C_A,C20C_B} {ave((rh.1-rh.2),time=1jan1960,time=1dec1960)} {0:360,-90:90,1000:1000,1nov1976:1nov1976} but always get this error message: Open Error: Can't open description file Thanks, PY From jma at COLA.IGES.ORG Fri Jul 11 13:45:04 2008 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Fri, 11 Jul 2008 13:45:04 -0400 Subject: Performing remote analysis on Windows In-Reply-To: <20080711163110.41756200D5@mx2.cineca.it> Message-ID: Oof! That documentation is outdated. Try this instead: ga-> sdfopen http://dataportal.ucar.edu:9191/dods/_expr_{/C20C/ C20C_A,/C20C/C20C_B} {ave((rh.1-rh.2),time=1jan1960,time=1dec1960)} {0:360,-90:90,1000:1000,1nov1976:1nov1976} (Note line should not wrap when you enter this command.) Jennifer On Jul 11, 2008, at 12:30 PM, PY Li wrote: > Hi, > > I'm a newbie, and just wondering if it's possible to perform remote > analysis on Windows on the data on > nomads.ncdc.noaa.gov:9091/dods/NCEP_NARR_DAILY > > I've tried the following operation found in GrADS-DODS Server - User's > Guide (http://nomads.ncdc.noaa.gov:9091/user.html): > > ga-> sdfopen http://motherlode.ucar.edu:9090/dods/_expr_ > {C20C_A,C20C_B} > {ave((rh.1-rh.2),time=1jan1960,time=1dec1960)} > {0:360,-90:90,1000:1000,1nov1976:1nov1976} > > but always get this error message: > > Open Error: Can't open description file > > Thanks, > > PY -- 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/20080711/aae05f27/attachment.html From dasilva at ALUM.MIT.EDU Sat Jul 12 00:35:45 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 12 Jul 2008 00:35:45 -0400 Subject: PyGrADS: Platform specific notes Message-ID: Dear PyGrADS Users, I just started a section on the PyGrADS wiki with platform specific notes, mostly about installation issues: http://opengrads.org/wiki/index.php?title=Python_Interface_to_GrADS#Platform_Specific_Notes If you have installed PyGrADS on platform not yet listed and would like to share some tips, please contact me directly. Thank you, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080712/183322e6/attachment.html From dasilva at ALUM.MIT.EDU Sat Jul 12 00:47:40 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 12 Jul 2008 00:47:40 -0400 Subject: Packaging GrADS for Mac OS X Message-ID: Hi, Is anybody out there who has experience packaging Mac OS X applications and building universal binaries on Darwin? I'd like prepare native, universal Mac packages for the OpenGrADS releases and I'd be very appreciative if someone who has done this before could help us out. Please drop us a note at opengrads-devel at lists.sourceforge.net if you think you can help. Thank you, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080712/29337637/attachment.html From dasilva at ALUM.MIT.EDU Sat Jul 12 01:02:06 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 12 Jul 2008 01:02:06 -0400 Subject: opening gradsda In-Reply-To: <6536a80c5131255ea3008fb8a6e6300d.squirrel@mail.uevora.pt> Message-ID: On Fri, Jul 11, 2008 at 5:16 AM, Isilda da Cunha Menezes wrote: > Hello! > I installed the new version of Linux for grads (GrADS 2.0.a2) and I use it > to open a file .nc (netcdf), but I am having trouble opening the gradsdap. > I wrote, as I am accustomed to using the Linux version GrADS 1.8s11, > gradsdap -lc "sdfopen air.nc" but got an error: > > gradsdap: error while loading shared libraries: libssl.so.4: cannot open > shared object file: No such file or directory > isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ gradsdap > gradsdap: error while loading shared libraries: libssl.so.4: cannot open > shared object file: No such file or directory > isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ gradsdap > gradsdap: error while loading shared libraries: libssl.so.4: cannot open > shared object file: No such file or directory > isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ grads > > and I have the same error when I try to open only the gradsdap. But the > grads opens: > > isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ grads > > Grid Analysis and Display System (GrADS) Version 2.0.a2 > Copyright (c) 1988-2008 by Brian Doty and the > Institute for Global Environment and Society (IGES) > GrADS comes with ABSOLUTELY NO WARRANTY > See file COPYRIGHT for more information > > Config: v2.0.a2 little-endian readline printim grib2 netcdf hdf4-sds > Issue 'q config' command for more information. > Landscape mode? ('n' for portrait): -l > GX Package Initialization: Size = 11 8.5 > ga-> quit > No hardcopy metafile open > GX package terminated > isilda at bram:~/BRAMS-prev/dprep/Gera_DP$ > > I am committing a mistake to open the gradsdap or was poorly installed? > > As the messages says, you are missing the libssl.so.4 on your system. Ask your sys admin to verify whether SSL is installed. It has been my experience that SSL is not very portable outside the US and I have omitted it from the OpenGrADS builds. Linking SSL statically is a technical solution, but unfortunately it violates US export laws since SSL has encryption technology; I have additional info at the Supplibs wiki: http://opengrads.org/wiki/index.php?title=Supplemental_Libraries_%28Supplibs%29#Why_is_gradsdap_.28formerly_gradsdods.29_built_without_SSL_support.3F Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080712/ed734d9a/attachment.html From dasilva at ALUM.MIT.EDU Sat Jul 12 01:07:41 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sat, 12 Jul 2008 01:07:41 -0400 Subject: HI ALL In-Reply-To: <20080711072325.M65291@mail.tropmet.res.in> Message-ID: On Fri, Jul 11, 2008 at 3:29 AM, Bakshi Vaid wrote: > Dear All, > > I would like to compute velocity potential. Does anybody have a script to > compute it?. > Please send. Thank you in advance. > Search the list for the keywords "fish, psi, chi, extensions". Some info here as well: http://opengrads.org/doc/udxt/fish/ Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080712/fdd103d2/attachment.html From bhv at TROPMET.RES.IN Sat Jul 12 03:26:13 2008 From: bhv at TROPMET.RES.IN (Bakshi Vaid) Date: Sat, 12 Jul 2008 13:56:13 +0630 Subject: HI ALL In-Reply-To: <77fcd6b20807112207s512371c3hc51955d989a222fb@mail.gmail.com> Message-ID: Dear Arlindo Sir, Thank you very much for kind help. I will try yours said route to compute required variable. Best Regards Bakshi Hardeep Vaid Research Associate Theoretical Studies Division Indian Institute of Tropical Meteorology Pashan, Pune, INDIA PIN:411008 Email: hardeep_vaid at yahoo.co.in Phone: +91-20-25893600 Ext.277 (office) +91-9850140881 (mobile) -- Open WebMail Project (http://openwebmail.org) ---------- Original Message ----------- From: da Silva To: GRADSUSR at LIST.CINECA.IT Sent: Sat, 12 Jul 2008 01:07:41 -0400 Subject: Re: HI ALL > On Fri, Jul 11, 2008 at 3:29 AM, Bakshi Vaid wrote: > > > Dear All, > > > > I would like to compute velocity potential. Does anybody have a script to > > compute it?. > > Please send. Thank you in advance. > > > > Search the list for the keywords "fish, psi, chi, extensions". Some info > here as well: > > http://opengrads.org/doc/udxt/fish/ > > Arlindo > > -- > Arlindo da Silva > dasilva at alum.mit.edu ------- End of Original Message ------- From ricardom at ACD.UFRJ.BR Sat Jul 12 13:18:10 2008 From: ricardom at ACD.UFRJ.BR (Ricardo Marcelo) Date: Sat, 12 Jul 2008 14:18:10 -0300 Subject: Output AScii file! In-Reply-To: <36465.24401.qm@web30404.mail.mud.yahoo.com> Message-ID: as a quick work around may be you could consider using 'sort' unix command. On Fri, Jul 11, 2008 at 4:42 AM, abbas mofidi wrote: > Dear GrADS Users, > I need to get the output of GrADS in text format. > My problem is that the output AScii file is arranged in reverse order. > I would be very grateful if someone can help me. > > Best regards, > A. mofidi > > > this is my script: > 'reinit' > 'sdfopen ncdata\Caspian-winter\AVHRR\sst4-path-eot.20011202.nc' > 'sdfopen ncdata\Caspian-winter\AVHRR\sst4-path-eot.20011201.nc' > 'set dfile 1' > 'set lon 50.125 50.625' > 'set lat 38.125 39.125' > 'set t 1' > 'set z 1' > 'define S1 = sst' > 'set dfile 2' > 'set lon 50.125 50.625' > 'set lat 38.125 39.125' > 'set t 1' > 'set z 1' > 'define S2 = sst' > 'set gxout print' > 'set prnopts %g 3 1' > 'display S1-S2' > rc=write(DSST2Dec01.txt,result) > rc=close(DSST2Dec01.txt) > return > > > > This is output .txt file: > > Printing Grid -- 15 Values -- Undef = -999 > -0.13 -0.12 -0.12 > -0.14 -0.13 -0.13 > -0.14 -0.13 -0.13 > -0.139 -0.11 -0.09 > -0.12 -0.24 -0.139 > > > The latitude and Longitude of above output file is: > > Lat/Lon 50.125 50.375 50.625 > 38.125 -0.13 -0.12 -0.12 > 38.375 -0.14 -0.13 -0.13 > 38.625 -0.14 -0.13 -0.13 > 38.875 -0.139 -0.11 -0.09 > 39.125 -0.12 -0.24 -0.139 > > > While, it should be like this: > > Lat/Lon 50.125 50.375 50.625 > 39.125 -0.12 -0.24 -0.139 > 38.875 -0.139 -0.11 -0.09 > 38.625 -0.14 -0.13 -0.13 > 38.375 -0.14 -0.13 -0.13 > 38.125 -0.13 -0.12 -0.12 > > > > -- ============================================================ Ricardo Marcelo da Silva LAMMA - Laborat?rio de Modelagem de Processos Marinhos e Atmosf?ricos Depto de Meteorologia - Instituto de Geoci?ncias Centro de Ci?ncias Matem?ticas e da Natureza Universidade Federal do Rio de Janeiro e-mail: ricardom at acd.ufrj.br phone: 55 21 2598-9470 r 24 ============================================================ -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080712/b43fc42f/attachment.html From pemck at GMX.NET Mon Jul 14 16:18:43 2008 From: pemck at GMX.NET (Paul emck) Date: Mon, 14 Jul 2008 22:18:43 +0200 Subject: Bug with MINing in script? Message-ID: Hi, This script will not return a value, hence neither write something sensible to the text file. When executing the script until *Here and then manually enter 'd min(minair,lon=290,lon=300)' into the GrADS command line it does return the desired value perfectly (-1.70758). Somebody a hint? Thanks, I would really appreciate some! (This script is simplified; it calculates zonal maximum negativ anomalies. I stripped the loops etc. for clearness) ------script begin--------- 'reinit' 'c' file='lat100N.txt' 'sdfopen air.day.sig995.ltm.nc' 'sdfopen air.sig995.1948.nc' 'set lon 290 300' 'set lat 10' 'define minair=air.2(t=1)-273.15-air.1(t=1)' *Here 'd min(minair,lon=290,lon=300)' val=subwrd(result,4) OK=write(lat100N.txt,val,append) ------script end--------- Config Grads: v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats printim From hoop at COLORADO.EDU Mon Jul 14 18:53:44 2008 From: hoop at COLORADO.EDU (Don Hooper) Date: Mon, 14 Jul 2008 16:53:44 -0600 Subject: Bug with MINing in script? Message-ID: Seems to me that the issue here is that you want word four of the second line in result. So, adding the line: line2=sublin(result,2) and changing the other assignment to: val=subwrd(line2,4) Should get you what you want. -Hoop > From owner-GRADSUSR at LIST.CINECA.IT Mon Jul 14 16:08:35 2008 > Date: Tue, 15 Jul 2008 00:07:48 +0200 > Reply-To: GRADSUSR at LIST.CINECA.IT > From: Automatic digest processor > Subject: GRADSUSR Digest - 12 Jul 2008 to 14 Jul 2008 (#2008-195) > To: Recipients of GRADSUSR digests > > There is one message totalling 40 lines in this issue. > > Topics of the day: > > 1. Bug with MINing in script? > > ---------------------------------------------------------------------- > > Date: Mon, 14 Jul 2008 22:18:43 +0200 > From: Paul emck > Subject: Bug with MINing in script? > > Hi, > This script will not return a value, hence neither write something sensible > to the text file. When executing the script until *Here and then manually > enter 'd min(minair,lon=290,lon=300)' into the GrADS command line it does > return the desired value perfectly (-1.70758). > Somebody a hint? Thanks, I would really appreciate some! > > (This script is simplified; it calculates zonal maximum negativ anomalies. > I stripped the loops etc. for clearness) > > ------script begin--------- > 'reinit' > 'c' > file='lat100N.txt' > 'sdfopen air.day.sig995.ltm.nc' > 'sdfopen air.sig995.1948.nc' > > 'set lon 290 300' > 'set lat 10' > > 'define minair=air.2(t=1)-273.15-air.1(t=1)' > > *Here > > 'd min(minair,lon=290,lon=300)' > > val=subwrd(result,4) > OK=write(lat100N.txt,val,append) > > ------script end--------- > > Config Grads: v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats > printim > > ------------------------------ > > End of GRADSUSR Digest - 12 Jul 2008 to 14 Jul 2008 (#2008-195) > *************************************************************** > From leebyerle at COMCAST.NET Mon Jul 14 22:50:24 2008 From: leebyerle at COMCAST.NET (Lee Byerle) Date: Mon, 14 Jul 2008 21:50:24 -0500 Subject: Template option for files with Julian Day names In-Reply-To: <20080708095403.8jx1zs9dcs8440c0@www.nexusmail.uwaterloo.ca> Message-ID: Zhenhao Bao- I agree a %j3 option would be great. I had a similar issue as yours once and used a script from an open forum (attached) that has worked okay to convert to Gregorian calendar. But in that case I needed the mmdd info to create a new grib record and needed naming convention changed., and no longer needed the julian in the file name. Best wishes, Lee Zhenhao Bao wrote: > Zhenhao Bao (PhD) > Dept. of Geography > University of Waterloo > 200 university avenue west > Waterloo, Ontario, Canada n2l 3g1 > z3bao at uwaterloo.ca > > > Quoting Arlindo da Silva : > >> On Thu, Jul 3, 2008 at 11:43 AM, Renu R. Joseph >> wrote: >> >>> Hello All, >>> >>> I need suggestions/solutions on dealing with the problem of opening >>> files >>> with julian day included in the names. >>> >>> I have files that named with the following name format: >>> m6s4_3hYYDDDHH.bin, >>> where YY represents the last two digits of the years (01 is 2001); >>> DDD is >>> julian day; HH is the hour of the start of the measurement >>> >>> If grads allowed a julian day formats it could have been opened with >>> "m6s4_3h%y2%jd%h200.bin" in the dset file location as follows. >>> >>> dset ^m6s4_3h%y2%jd3%h2.bin >>> OPTIONS template >>> OPTIONS big_endian >>> OPTIONS yrev >>> UNDEF -9999 >>> TITLE Precipitation estimates >>> XDEF 1440 LINEAR 0.125 0.25 >>> YDEF 480 LINEAR -59.875 0.25 >>> ZDEF 01 LEVELS 1 >>> TDEF 1 LINEAR 00z01jan2001 3 hr >>> VARS 1 >>> per 1 99 (mm/hr) >>> ENDVARS >>> >>> Since I cannot do that, I know that one of the options is to use >>> DSET ^m6s4_3h%ch.bin >>> CHSUB 1 1 m6s4_3h0100100.bin >>> CHSUB 2 2 m6s4_3h0100100.bin >>> .................................................. >>> CHSUB 20448 20448 m6s4_3h0736521.bin >>> OPTIONS template >>> >>> but that includes so many lines. Is there a way around this, unless I >>> physically rename files to be consistent with grads format names? >>> >> >> >> You basically outlined your options. The best solution is to include >> julian >> day support in grads templating, something that is badly needed (I >> also ran >> into such files before). In my case I had a good reason not to >> rename the >> files (other applications required the julian day type of name >> conventions), >> and sometimes making a copy of the files isn't an option (disk space >> considerations). Here are a couple of options which leave the >> original files >> alone: >> >> 1) Symlink each file to a grads friendly name in a separate >> directory, e.g., >> grads/m6s4_2007-02-05_12z.bin. You should be able to do this with a >> relatively simple script. >> >> 2) If you know C, implement %j3 support in GrADS and contribute a >> patch. In >> GrADS both v1.9 and v2, look at this function in gautil.c: >> >> /* Given a file name template and a dt structure, fill in to get the >> file >> name */ >> char *gafndt (char *fn, struct dt *dtim, struct dt *dtimi, float *vals, >> struct gachsub *pch1st, int t) { >> >> Good Luck, >> >> Arlindo >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> > -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: julian_to_gregorian.csh Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080714/7e1283f6/attachment.pl From teddyallen at YAHOO.COM Tue Jul 15 09:38:51 2008 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Tue, 15 Jul 2008 06:38:51 -0700 Subject: hovmoller question Message-ID: Dear GrADS community, I am a new GrADS user and am progressing quickly with the software thanks to the generosity of this group and the supporting documentation. However, I am currently stuck on one application that I just can't seem to figure out. I would like to create a shaded Hovmoller plot that describes monthly average values along the y-axis (12 tics) and longitude along the x-axis. I am using the GPCP monthly dataset and can easily generate a plot describing all of the years of the dataset along the y-axis (29 years), but can't seem how to figure out how to average the dataset to display monthly climatologies. Any suggestions on how to display monthly averages along the y-axis on a hovmoller plot? Below is my script for the hovmoller plot describing all of the years along the y-axis. I am using GrADS NC version 1.8SL11 on a windows computer. 'reinit' 'sdfopen c:\testdata\netcdf\precip.mon.mean.nc' 'set grid off' 'set grads off' 'set gxout shaded' 'set t 253 348' *time above displays GPCP years from 2000-2007 'set y 1' 'set lon 270 300' 'set clevs 0 .5 1 1.5 2 2.5 3 3.5 4 4.5 6 6.5 7 7.5 8 8.5 9' 'set grads off' 'set yaxis 253 348 4' 'set csmooth on' 'd ave(precip,lat=15,lat=17)' 'draw xlab longitude' 'draw ylab time' 'draw title Meridional Average (lat=15N-20N) of GPCP 2000-2007' 'cbarn' 'print' 'disable print' Thank you very much for your help. Sincerely, teddy allen To move quickly, go alone. To move far, go together. --- African Proverb From bernd.becker at METOFFICE.GOV.UK Tue Jul 15 09:50:41 2008 From: bernd.becker at METOFFICE.GOV.UK (Bernd Becker) Date: Tue, 15 Jul 2008 14:50:41 +0100 Subject: Ensemble gridded obs data set in netcdf Message-ID: Hello, has anyone had success on how to read the gridded obs data set in netcdf with Grads? http://eca.knmi.nl/download/ensembles/download.php 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.com - http://www.metoffice.com From davide.sacchetti at ARPAL.ORG Tue Jul 15 10:17:51 2008 From: davide.sacchetti at ARPAL.ORG (Davide Sacchetti) Date: Tue, 15 Jul 2008 16:17:51 +0200 Subject: hovmoller question In-Reply-To: <618381.82350.qm@web51511.mail.re2.yahoo.com> Message-ID: on my opinion the easiest way is to create an ad hoc dataset with the means you are interested on, eg: open files for 2000 ... 2007 'set gxout fwrite' loop for 12 months: 'd (precip.1+precip.2+...+precip.7)/7' 'disable fwrite' create ctl file for grads.fwrite open it and use your script have nice diagrams ... bye bye Davide p.s.: ... your reply-to was set to yourself instead of the list On Tue, 2008-07-15 at 06:38 -0700, Teddy Allen wrote: > Dear GrADS community, > I am a new GrADS user and am progressing quickly with the software thanks to the generosity of this group and the supporting documentation. However, I am currently stuck on one application that I just can't seem to figure out. > I would like to create a shaded Hovmoller plot that describes monthly average values along the y-axis (12 tics) and longitude along the x-axis. I am using the GPCP monthly dataset and can easily generate a plot describing all of the years of the dataset along the y-axis (29 years), but can't seem how to figure out how to average the dataset to display monthly climatologies. Any suggestions on how to display monthly averages along the y-axis on a hovmoller plot? > Below is my script for the hovmoller plot describing all of the years along the y-axis. I am using GrADS NC version 1.8SL11 on a windows computer. > > 'reinit' > 'sdfopen c:\testdata\netcdf\precip.mon.mean.nc' > 'set grid off' > 'set grads off' > 'set gxout shaded' > 'set t 253 348' > *time above displays GPCP years from 2000-2007 > 'set y 1' > 'set lon 270 300' > 'set clevs 0 .5 1 1.5 2 2.5 3 3.5 4 4.5 6 6.5 7 7.5 8 8.5 9' > 'set grads off' > 'set yaxis 253 348 4' > 'set csmooth on' > 'd ave(precip,lat=15,lat=17)' > 'draw xlab longitude' > 'draw ylab time' > 'draw title Meridional Average (lat=15N-20N) of GPCP 2000-2007' > 'cbarn' > 'print' > 'disable print' > > Thank you very much for your help. > Sincerely, > teddy allen > > To move quickly, go alone. To move far, go together. --- African Proverb > -- 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 Charles.Seman at NOAA.GOV Tue Jul 15 10:42:25 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Tue, 15 Jul 2008 10:42:25 -0400 Subject: hovmoller question In-Reply-To: <618381.82350.qm@web51511.mail.re2.yahoo.com> Message-ID: Teddy, Please find attached a script pptn_climo_Hovmuller.gs and a sample plot which I hope will help... Please let me know if you have any questions about the script... Thanks, Chuck Teddy Allen wrote: > Dear GrADS community, > I am a new GrADS user and am progressing quickly with the software thanks to the generosity of this group and the supporting documentation. However, I am currently stuck on one application that I just can't seem to figure out. > I would like to create a shaded Hovmoller plot that describes monthly average values along the y-axis (12 tics) and longitude along the x-axis. I am using the GPCP monthly dataset and can easily generate a plot describing all of the years of the dataset along the y-axis (29 years), but can't seem how to figure out how to average the dataset to display monthly climatologies. Any suggestions on how to display monthly averages along the y-axis on a hovmoller plot? > Below is my script for the hovmoller plot describing all of the years along the y-axis. I am using GrADS NC version 1.8SL11 on a windows computer. > > 'reinit' > 'sdfopen c:\testdata\netcdf\precip.mon.mean.nc' > 'set grid off' > 'set grads off' > 'set gxout shaded' > 'set t 253 348' > *time above displays GPCP years from 2000-2007 > 'set y 1' > 'set lon 270 300' > 'set clevs 0 .5 1 1.5 2 2.5 3 3.5 4 4.5 6 6.5 7 7.5 8 8.5 9' > 'set grads off' > 'set yaxis 253 348 4' > 'set csmooth on' > 'd ave(precip,lat=15,lat=17)' > 'draw xlab longitude' > 'draw ylab time' > 'draw title Meridional Average (lat=15N-20N) of GPCP 2000-2007' > 'cbarn' > 'print' > 'disable print' > > Thank you very much for your help. > Sincerely, > teddy allen > > To move quickly, go alone. To move far, go together. --- African Proverb > -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: pptn_climo_Hovmuller.gs Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080715/0290900a/attachment.pl -------------- next part -------------- A non-text attachment was scrubbed... Name: pptn_climo_Hovmuller.eps Type: application/ps Size: 24258 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080715/0290900a/attachment.bin From Charles.Seman at NOAA.GOV Tue Jul 15 11:57:57 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Tue, 15 Jul 2008 11:57:57 -0400 Subject: hovmoller question In-Reply-To: <487CB751.4050404@noaa.gov> Message-ID: Listserv, The script I sent before contains an error: for "t+13" the averaging starts at t=14 (which is actually an error! sorry about that, it should be t+12 so that the averaging starts at t=13 which is the January of the second year in the dataset for the script... I've attached a new script and plot...) Thanks, Chuck Charles Seman wrote: > Teddy, > > Please find attached a script pptn_climo_Hovmuller.gs and a sample plot > which I hope will help... Please let me know if you have any questions > about the script... > > Thanks, > Chuck > > Teddy Allen wrote: >> Dear GrADS community, >> I am a new GrADS user and am progressing quickly with the software >> thanks to the generosity of this group and the supporting >> documentation. However, I am currently stuck on one application that >> I just can't seem to figure out. >> I would like to create a shaded Hovmoller plot that describes monthly >> average values along the y-axis (12 tics) and longitude along the >> x-axis. I am using the GPCP monthly dataset and can easily generate a >> plot describing all of the years of the dataset along the y-axis (29 >> years), but can't seem how to figure out how to average the dataset >> to display monthly climatologies. Any suggestions on how to display >> monthly averages along the y-axis on a hovmoller plot? >> Below is my script for the hovmoller plot describing all of the years >> along the y-axis. I am using GrADS NC version 1.8SL11 on a windows >> computer. >> >> 'reinit' >> 'sdfopen c:\testdata\netcdf\precip.mon.mean.nc' >> 'set grid off' >> 'set grads off' >> 'set gxout shaded' >> 'set t 253 348' >> *time above displays GPCP years from 2000-2007 >> 'set y 1' >> 'set lon 270 300' >> 'set clevs 0 .5 1 1.5 2 2.5 3 3.5 4 4.5 6 6.5 7 7.5 8 8.5 9' >> 'set grads off' >> 'set yaxis 253 348 4' >> 'set csmooth on' >> 'd ave(precip,lat=15,lat=17)' >> 'draw xlab longitude' >> 'draw ylab time' >> 'draw title Meridional Average (lat=15N-20N) of GPCP 2000-2007' >> 'cbarn' >> 'print' >> 'disable print' >> >> Thank you very much for your help. >> Sincerely, >> teddy allen >> >> To move quickly, go alone. To move far, go together. --- African Proverb >> > > -- > > Please note that Charles.Seman at noaa.gov should be considered my NOAA > email address, not cjs at gfdl.noaa.gov. > > ******************************************************************** > Charles Seman Charles.Seman at noaa.gov > U.S. Department of Commerce / NOAA / OAR > Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 > 201 Forrestal Road fax: (609) 987-5063 > Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ > ******************************************************************** > > "The contents of this message are mine personally and do not necessarily > reflect any position of the Government or NOAA." -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: pptn_climo_Hovmuller.gs Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080715/38a224ea/attachment.pl -------------- next part -------------- A non-text attachment was scrubbed... Name: pptn_climo_Hovmuller.eps Type: application/ps Size: 24000 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080715/38a224ea/attachment.bin From pemck at GMX.NET Tue Jul 15 16:48:58 2008 From: pemck at GMX.NET (Paul emck) Date: Tue, 15 Jul 2008 22:48:58 +0200 Subject: Bug with MINing in script? Message-ID: Thank you so much. It did the job. Cheers, Paul From sreeds82 at YAHOO.COM Wed Jul 16 04:01:06 2008 From: sreeds82 at YAHOO.COM (sreenivas desamsetti) Date: Wed, 16 Jul 2008 01:01:06 -0700 Subject: How to work with GRIB2 in GrADS Message-ID: Dear GrADS Users, Could you please tell me how to work the GRIB2 files in GrADS directly. Means how to create ctl, idx files and please let me know the necessary libraries need to be installed Thanking you D Srinivas Desamsetti Srinivas Research Scholar, Department of Meteorology and Oceanography, Andhra University, Visakhapatnam - 530 003 91 9948530618 (M) sreeds82 at yahoo.com sreeds82 at rediffmail.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080716/59b075bb/attachment.html From Charles.Seman at NOAA.GOV Wed Jul 16 10:32:52 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Wed, 16 Jul 2008 10:32:52 -0400 Subject: new pptn_climo_Hovmuller.gs... Message-ID: Dear GrADS Users, The pptn_climo_Hovmuller.gs script I sent to the Listserv contains an error: for "t+13" the averaging starts at t=14 (which is actually an error! sorry about that, it should be t+12 so that the averaging starts at t=13 which is the January of the second year in the dataset for the script... I've attached a new script and plot...) Sorry about that, Chuck -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: pptn_climo_Hovmuller.gs Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080716/c7ba7546/attachment.pl -------------- next part -------------- A non-text attachment was scrubbed... Name: pptn_climo_Hovmuller.eps Type: application/ps Size: 24000 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080716/c7ba7546/attachment.bin From guojingtian at OUC.EDU.CN Wed Jul 16 23:41:29 2008 From: guojingtian at OUC.EDU.CN (=?gb2312?B?ufm+tMzs?=) Date: Thu, 17 Jul 2008 11:41:29 +0800 Subject: full bar represents 10 ? Message-ID: Dear Grads users, Who knows the full bar representing values when use the command "set gxout barb", it's 10 or other?? regards jingtian From arabonis at EGS.UCT.AC.ZA Thu Jul 17 12:09:29 2008 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Thu, 17 Jul 2008 18:09:29 +0200 Subject: http://www.monsoondata.org:9090/dods/ Message-ID: Hi I can't seem to connect to http://www.monsoondata.org:9090/dods/ , is the dods server down or am I experiencing a local problem? 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 jma at COLA.IGES.ORG Thu Jul 17 12:40:13 2008 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Thu, 17 Jul 2008 12:40:13 -0400 Subject: http://www.monsoondata.org:9090/dods/ In-Reply-To: <487F6EB9.8010902@egs.uct.ac.za> Message-ID: I have restarted the 9090 server -- please try again. --Jennifer On Jul 17, 2008, at 12:09 PM, Jean Pierre Arabonis wrote: > Hi > > I can't seem to connect to http://www.monsoondata.org:9090/ > dods/ , is > the dods server down or am I experiencing a local problem? > > 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 -- 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/20080717/59a8669d/attachment.html From jhuddleston at HUGHES.NET Thu Jul 17 15:53:54 2008 From: jhuddleston at HUGHES.NET (John Huddleston) Date: Thu, 17 Jul 2008 13:53:54 -0600 Subject: Update to the grib2 and w3 libraries and the Windows Cygwin build Message-ID: Hi all, The http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/ site contains grib2 and w3 library source code. The G2 library has been updated to 1.1.6. There are some local parameter G2 changes increasing the MAXPARAM=416 to MAXPARAM=454 and it includes changes to 351, 352, and 353. The W3 library has been updated to 1.6. A new grid 195 has been added and grid 196 has been changed. I've built the cnvgrib utility for Windows Cygwin. All it requires is the cygwin1.dll file from http://sources.redhat.com/cygwin. I'll rebuild the GrADS 2.0a2 binary and send it out to those who have previously asked for a copy. Jennifer, I have not heard any feedback on problems with this Cygwin build. Can you put it on you site for others to download? Thank you, John Huddleston, PhD, PE Easy Computer Solutions -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080717/82af482a/attachment.html From anelisebaur at YAHOO.COM.BR Thu Jul 17 16:04:11 2008 From: anelisebaur at YAHOO.COM.BR (Anelise Baur) Date: Thu, 17 Jul 2008 13:04:11 -0700 Subject: From grib file to binary file Message-ID: Hello grads users! ? I?m trying to write a single binary file, the first records from two different grib files (nww3.UGRD.grb and nww3.VGRD.grb, from NOAA?s site). When I open this new file in GRADS I keep getting errors, particularly for the second variable. Can someone help me? My procedures were: ? 1? I?ve created the .ctl files (nww3.UGRD.ctl and nww3.VGRD.ctl) with grib2ctl.pl; I?ve obtained the following files: ? dset ^nww3.UGRD.grb index ^nww3.UGRD.grb.idx undef 9.999E+20 title nww3.UGRD.grb * produced by grib2ctl v0.9.12.5p39c dtype grib 233 options yrev ydef 157 linear -78.000000 1 xdef 288 linear 0.000000 1.250000 tdef 61 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 1 UGRDsfc 0 33,1,1 ** surface u wind [m/s] ENDVARS ? ? dset ^nww3.vgrd.grb index ^nww3.vgrd.grb.idx undef 9.999E+20 title nww3.vgrd.grb * produced by grib2ctl v0.9.12.5p39c dtype grib 233 options yrev ydef 157 linear -78.000000 1 xdef 288 linear 0.000000 1.250000 tdef 61 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 1 VGRDsfc 0 34,1,1 ** surface v wind [m/s] ENDVARS ? 2? I?ve created the grb.idx files (nww3.UGRD.grb.idx and nww3.VGRD.grb.idx) with grib2ctl.pl and gribmap. ? 3? I?ve opened nww3.UGRD.ctl and nww3.VGRD.ctl in the GRADS. ? 4? I did ?set gxout fwrite?, then ?d UGRDsfc? and ?d VGRDsfc?. I?ve obtained the binary file grads.fwrite (with 354kb). ? 5? Then, I?ve created a .ctl file for the file grads.fwrite, called teste.ctl: ? dset ^grads.fwrite undef 9.999E+20 title teste20_06 ydef 158 linear -78.000000 1 xdef 289 linear 0.000000 1.250000 tdef 1 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 2 u 0 99 ** surface v wind [m/s] v 0 99 ** surface v wind [m/s] ENDVARS ? 6? I?ve opened teste.ctl in the GRADS, and I?ve written ?set gxout shaded?, then ?d u?, without problems. The problem appears when I display the next variable,v. GRADS doesn?t recognize v. The error is: ? ? Low Level I/O Error: Read error on data file Data file name = e:/swan_pendrive/02_07/25_06/25_06/grads.fwrite Error reading 289 bytes at location 90457 Data Request Error: Error for variable 'v' Error ocurred at column 1 DISPLAY error: Invalid expression Expression = v ? ? It indicates a problem with the binary file grads.fwrite. ? Thanks in advance, ? Anelise Pereira Baur Universidade Federal do Rio Grande do Sul Instituto de Matem?tica anelisebaur at yahoo.com.br Novos endere?os, o Yahoo! que voc? conhece. Crie um email novo com a sua cara @ymail.com ou @rocketmail.com. http://br.new.mail.yahoo.com/addresses -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080717/45ae56e9/attachment.html From ela at COLA.IGES.ORG Thu Jul 17 19:16:04 2008 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Thu, 17 Jul 2008 19:16:04 -0400 Subject: From grib file to binary file In-Reply-To: <2032096642.82151216336377716.JavaMail.root@mail.iges.org> Message-ID: Hello Anelise, The dimensions for ydef (and possibly for xdef) in your teste.ctl file are wrong. ydef should have a dimension of 157, the same as it is in the ctl files for the grib datasets. The situation for xdef is more complicated because the datasets are global in longitude. If you open such a dataset in grads, and it is the first dataset you open in that grads session, the dimension environment will automatically "wrap around" in longitude. What this means is that data at the starting longitude (0.0 degrees in this case) is duplicated internally by grads, i.e. if the longitudes of the grid points in the dataset are: 0.0 1.25 2.5 ... 357.5 358.75 (288 longitudes) then grads will duplicate the data along 0.0 degrees: 0.0 1.25 2.5 ... 357.5 358.75 0.0 (289 longitudes) so that the x-dimension of the internally stored data is actually 289 instead of 288. If you write data in 'fwrite' mode without changing the range of longitude with 'set x' or 'set lon', the resulting fwrite file will have an x-dimension of 289. You probably don't want this. To prevent the wrap-around, do the following before writing the fwrite file: set x 1 288 or, alternatively, set lon 0 358.75 The first option is easier since it's necessary to calculate the easternmost longitude (358.75 in this case). Now, your fwrite file should be slightly smaller than before (288*157*4*2 = 361728 bytes). In your teste.ctl file, the dimensions of xdef and ydef should be 288 and 157. It should work now. Best regards, Eric L. Altshuler Assistant Research Scientist Center for Ocean-Land-Atmosphere Studies 4041 Powder Mill Road, Suite 302 Calverton, MD 20705-3106 USA E-mail: ela at cola.iges.org Phone: (301) 902-1257 Fax: (301) 595-9793 ----- Original Message ----- From: "Anelise Baur" To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, July 17, 2008 4:04:11 PM GMT -05:00 US/Canada Eastern Subject: From grib file to binary file Hello grads users! ? I???m trying to write a single binary file, the first records from two different grib files (nww3.UGRD.grb and nww3.VGRD.grb, from NOAA???s site). When I open this new file in GRADS I keep getting errors, particularly for the second variable. Can someone help me? My procedures were: ? 1?? I???ve created the .ctl files (nww3.UGRD.ctl and nww3.VGRD.ctl) with grib2ctl.pl; I???ve obtained the following files: ? dset ^nww3.UGRD.grb index ^nww3.UGRD.grb.idx undef 9.999E+20 title nww3.UGRD.grb * produced by grib2ctl v0.9.12.5p39c dtype grib 233 options yrev ydef 157 linear -78.000000 1 xdef 288 linear 0.000000 1.250000 tdef 61 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 1 UGRDsfc 0 33,1,1 ** surface u wind [m/s] ENDVARS ? ? dset ^nww3.vgrd.grb index ^nww3.vgrd.grb.idx undef 9.999E+20 title nww3.vgrd.grb * produced by grib2ctl v0.9.12.5p39c dtype grib 233 options yrev ydef 157 linear -78.000000 1 xdef 288 linear 0.000000 1.250000 tdef 61 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 1 VGRDsfc 0 34,1,1 ** surface v wind [m/s] ENDVARS ? 2?? I???ve created the grb.idx files (nww3.UGRD.grb.idx and nww3.VGRD.grb.idx) with grib2ctl.pl and gribmap. ? 3?? I???ve opened nww3.UGRD.ctl and nww3.VGRD.ctl in the GRADS. ? 4?? I did ???set gxout fwrite???, then ???d UGRDsfc??? and ???d VGRDsfc???. I???ve obtained the binary file grads.fwrite (with 354kb). ? 5?? Then, I???ve created a .ctl file for the file grads.fwrite, called teste.ctl: ? dset ^grads.fwrite undef 9.999E+20 title teste20_06 ydef 158 linear -78.000000 1 xdef 289 linear 0.000000 1.250000 tdef 1 linear 12Z28apr2008 3hr zdef 1 linear 1 1 vars 2 u 0 99 ** surface v wind [m/s] v 0 99 ** surface v wind [m/s] ENDVARS ? 6?? I???ve opened teste.ctl in the GRADS, and I???ve written ???set gxout shaded???, then ???d u???, without problems. The problem appears when I display the next variable,v. GRADS doesn???t recognize v. The error is: ? ? Low Level I/O Error: Read error on data file Data file name = e:/swan_pendrive/02_07/25_06/25_06/grads.fwrite Error reading 289 bytes at location 90457 Data Request Error: Error for variable 'v' Error ocurred at column 1 DISPLAY error: Invalid expression Expression = v ? ? It indicates a problem with the binary file grads.fwrite. ? Thanks in advance, ? Anelise Pereira Baur Universidade Federal do Rio Grande do Sul Instituto de Matem??tica anelisebaur at yahoo.com.br Novos endere??os, o Yahoo! que voc?? conhece. Crie um email novo com a sua cara @ymail.com ou @rocketmail.com. http://br.new.mail.yahoo.com/addresses From axl419 at YAHOO.COM Fri Jul 18 02:17:23 2008 From: axl419 at YAHOO.COM (See Hai Ooi) Date: Thu, 17 Jul 2008 23:17:23 -0700 Subject: Converting ECMWF .pp data file into GrADS-readable data format In-Reply-To: <1038842649.82171216336564661.JavaMail.root@mail.iges.org> Message-ID: Dear all, Grateful if someone can provide me the methodology in converting the ECMWF Reanalysis .pp data file (grib2 file ?) into GrADS-readable (window version) data format (such as .nc). Thanking you in anticipation, Ooi Send instant messages to your online friends http://uk.messenger.yahoo.com From arabonis at EGS.UCT.AC.ZA Fri Jul 18 03:15:40 2008 From: arabonis at EGS.UCT.AC.ZA (Jean Pierre Arabonis) Date: Fri, 18 Jul 2008 09:15:40 +0200 Subject: http://www.monsoondata.org:9090/dods/ In-Reply-To: <5B346327-970B-4AA8-8215-1EB7579DDD54@cola.iges.org> Message-ID: Thanks - it works now - JP Jennifer Adams wrote: > I have restarted the 9090 server -- please try again. --Jennifer > > On Jul 17, 2008, at 12:09 PM, Jean Pierre Arabonis wrote: > >> Hi >> >> I can't seem to connect to http://www.monsoondata.org:9090/ >> dods/ , is >> the dods server down or am I experiencing a local problem? >> >> 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 > > -- > 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 From lincoln at CPTEC.INPE.BR Fri Jul 18 10:07:01 2008 From: lincoln at CPTEC.INPE.BR (Lincoln Muniz Alves) Date: Fri, 18 Jul 2008 11:07:01 -0300 Subject: Converting ECMWF .pp data file into GrADS-readable data format In-Reply-To: <591309.36972.qm@web31005.mail.mud.yahoo.com> Message-ID: Dear Ooi, I think that xconv it is a good tool. http://badc.nerc.ac.uk/help/software/xconv/ Lincoln Quoting See Hai Ooi : > Dear all, > > Grateful if someone can provide me the methodology in converting the > ECMWF Reanalysis .pp data file (grib2 file ?) into GrADS-readable (window > version) data format (such as .nc). > > Thanking you in anticipation, > > Ooi > > Send instant messages to your online friends http://uk.messenger.yahoo.com > -- Lincoln Muniz Alves Meteorologista - Grupo de Previs?o Clim?tica CPTEC/INPE +55 12 3186-8580 From ffalcini at QUIPO.IT Fri Jul 18 11:09:56 2008 From: ffalcini at QUIPO.IT (Falcini) Date: Fri, 18 Jul 2008 17:09:56 +0200 Subject: No subject Message-ID: Dear all, I am writing you for a Segmentation Fault break in a multiple file analysis. I have more than 30000 .nc files. Could I skip the problem by defining (probably by means of fwrite) smaller x, y, and z range?? An thus, to work with this new single file? Any suggestion and examples would be more than useful Many Thanks in advance Federico Falcini From axl419 at YAHOO.COM Fri Jul 18 18:26:56 2008 From: axl419 at YAHOO.COM (See Hai Ooi) Date: Fri, 18 Jul 2008 15:26:56 -0700 Subject: Converting ECMWF .pp data file into GrADS-readable data format In-Reply-To: <1216390021.4880a385bcecd@webmail.cptec.inpe.br> Message-ID: Dear Mr Lincoln, Thank you. By using the downloaded xconv (version 1.91), I notice that there are two files in the attached .pp data set. Both are 32-bit ieee byteswaped. However, after selecting one of them to be converted as test.nc file, errors occurred when I opened it in GrADS (window version 1.9.0-rc1). Can you guide me further, and if possible, provide me the correct step-by-step procedures ? Best regards. Ooi --- On Fri, 18/7/08, Lincoln Muniz Alves wrote: > From: Lincoln Muniz Alves > Subject: Re: Converting ECMWF .pp data file into GrADS-readable data format > To: GRADSUSR at LIST.CINECA.IT > Date: Friday, 18 July, 2008, 10:07 PM > Dear Ooi, > I think that xconv it is a good tool. > > http://badc.nerc.ac.uk/help/software/xconv/ > > Lincoln Send instant messages to your online friends http://uk.messenger.yahoo.com -------------- next part -------------- A non-text attachment was scrubbed... Name: 000000000000.15.201.000000.2000.05.01.00.00.pp Type: application/octet-stream Size: 1188672 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080718/5234ec96/attachment.obj -------------- next part -------------- A non-text attachment was scrubbed... Name: test.nc Type: application/octet-stream Size: 10092 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080718/5234ec96/attachment-0001.obj -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: test_msg.txt Url: http://gradsusr.org/pipermail/gradsusr/attachments/20080718/5234ec96/attachment.txt From aspratt2000 at YAHOO.COM Sat Jul 19 01:09:11 2008 From: aspratt2000 at YAHOO.COM (Mr Aaron Pratt) Date: Fri, 18 Jul 2008 22:09:11 -0700 Subject: Using PCGrads to examine TRMM HDF files Message-ID: Hello, I have downloaded TRMM orbit data files; I believe the file format is HDF. I wanted to use PCGrADS to open this file, but I am having problems. I know that there is a GrADSHDF executable within the PCGrADS that should be able to open the files directly, but that doesn't seem to work. I also try to use the TSDIS viewer to convert the HDF file (or at least a subsection) into a standard GrADS file. However, the .ctl file that was crreted still did not open (there seemed to be problems with the .dat file that was created). Is the Windows version of GrADS able to examine these TRMM orbit files. Also, are there any scripts available? Thanks in advance, Aaron From ranyere.nobrega at SIPAM.GOV.BR Sat Jul 19 07:12:30 2008 From: ranyere.nobrega at SIPAM.GOV.BR (=?ISO-8859-1?Q?Rany=E9re_N=F3brega?=) Date: Sat, 19 Jul 2008 07:12:30 -0400 Subject: Using PCGrads to examine TRMM HDF files In-Reply-To: <636837.84284.qm@web62408.mail.re1.yahoo.com> Message-ID: Dear Aaron, Use tthis ftp address disc2.nascom.nasa.gov/data/TRMM/ . In this address,the files can read directly in GrADS and yet own CTL. Ranyere Nobrega SIPAM - Sistema de Prote??o da Amaz?nia Mr Aaron Pratt escreveu: > Hello, > I have downloaded TRMM orbit data files; I believe the file format is HDF. I wanted to use PCGrADS to open this file, but I am having problems. I know that there is a GrADSHDF executable within the PCGrADS that should be able to open the files directly, but that doesn't seem to work. I also try to use the TSDIS viewer to convert the HDF file (or at least a subsection) into a standard GrADS file. However, the .ctl file that was crreted still did not open (there seemed to be problems with the .dat file that was created). Is the Windows version of GrADS able to examine these TRMM orbit files. Also, are there any scripts available? Thanks in advance, > Aaron > > -- __________________________________________________________ Rany?re Silva N?brega (Doutorando em meteorologia) Chefe Interino - Divis?o de Meteorologia e Climatologia Presid?ncia da Rep?blica - Casa Civil Sistema de Prote??o da Amaz?nia - Sipam Centro T?cnico e Operacional de Porto Velho Fone: (69)3217-6291 email: ranyere.nobrega at sipam.gov.br From polarlow at GMAIL.COM Sat Jul 19 12:50:38 2008 From: polarlow at GMAIL.COM (Polar Low) Date: Sat, 19 Jul 2008 13:50:38 -0300 Subject: New User Message-ID: Hello All, I`m a complete novice to grads. I have started to read the documentation but it was a little challenging, it seems more geared up for a university student or somone who already has a good understanding of it. I would like to use grads for creating meteorological products, such as nwp charts from grib data. Where would be a best place for me to learn how to do this, bearing in mind my "novice" status? I have the option to install grads on a windows or linux based pc, whichever it easiest for me. Thanks. Robert. From Charles.Seman at NOAA.GOV Sat Jul 19 14:39:50 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Sat, 19 Jul 2008 14:39:50 -0400 Subject: New User In-Reply-To: Message-ID: Robert, There is a tutorial at: http://grads.iges.org/grads/gadoc/tutorial.html Have you visited it? I hope this helps, Chuck Polar Low wrote: > Hello All, > > I`m a complete novice to grads. I have started to read the > documentation but it was a little challenging, it seems more geared up > for a university student or somone who already has a good > understanding of it. I would like to use grads for creating > meteorological products, such as nwp charts from grib data. > > Where would be a best place for me to learn how to do this, bearing in > mind my "novice" status? > > I have the option to install grads on a windows or linux based pc, > whichever it easiest for me. > > Thanks. > Robert. > -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." From aspratt2000 at YAHOO.COM Sat Jul 19 21:18:06 2008 From: aspratt2000 at YAHOO.COM (Mr Aaron Pratt) Date: Sat, 19 Jul 2008 18:18:06 -0700 Subject: Using Windows-based GrADS-HDF to examine TRMM orbit data files Message-ID: Hello, I am trying to?create plots derived from?the TRMM 2A12 orbital product (it's the hydrometeor profile data)?. It's not gridded data, but rather in orbital swaths. The files are HDF files, so that's why I thought that GrADS might be able to open them. However, I wasn't sure if the HDF files needed to be converted into a GrADS type file (like a .ctl file and the corresponding .dat file). I was told that the files could be converted through the following link: http://disc2.nascom.nasa.gov/data/TRMM/; however I don't see anywhere on the website where this is possible. Thanks in advance for any help that can be provided. Aaron -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080719/16ac04d0/attachment.html From ksrao at TROPMET.RES.IN Mon Jul 21 03:07:17 2008 From: ksrao at TROPMET.RES.IN (Seshagiri Rao) Date: Mon, 21 Jul 2008 13:37:17 +0630 Subject: Hi Message-ID: Hi, I'm wondering if anyone can help me in plotting the seasonal climatology say JJAS using the ERA40 accumlated rainfall dataset. Please share if anyone has scirpt.. Thanks in advance!!! With regards, K.S.Rao **************************************** K.Seshagiri Rao, Research Fellow, Climatology & Hydrometeorology Division, Indian Institute of Tropical Meteorology, Dr.Homibhabha Road, Pashan, Pune-411 008, India. Phone:- 020-25893600Ext.363 Mobile:- 09850957963. Alt.Email:- seshu_meteo16 at yahoo.com. **************************************** -- Open WebMail Project (http://openwebmail.org) From madhusoodanan at GMAIL.COM Mon Jul 21 03:39:38 2008 From: madhusoodanan at GMAIL.COM (Dr. Madhusoodanan M.S.) Date: Mon, 21 Jul 2008 13:09:38 +0530 Subject: Hi In-Reply-To: <20080721070358.M18870@mail.tropmet.res.in> Message-ID: Hi, pls refer in grads index : modify climatology seasonal http://www.iges.org/grads/gadoc/gadocindex.html regards, Madhu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080721/4d7b1219/attachment.html From ksrao at TROPMET.RES.IN Mon Jul 21 06:01:45 2008 From: ksrao at TROPMET.RES.IN (Seshagiri Rao) Date: Mon, 21 Jul 2008 16:31:45 +0630 Subject: Hi In-Reply-To: <726d83190807210039x4e11004crf4c191b91b5e3fb8@mail.gmail.com> Message-ID: Hi, Thanks for advice. Sorry i forgot to mention that i would like to convert total accumlated rainfall to total rainfall.. Thanks and regards, K. S. Rao **************************************** K.Seshagiri Rao, Research Fellow, Climatology & Hydrometeorology Division, Indian Institute of Tropical Meteorology, Dr.Homibhabha Road, Pashan, Pune-411 008, India. Phone:- 020-25893600Ext.363 Mobile:- 09850957963. Alt.Email:- seshu_meteo16 at yahoo.com. **************************************** -- Open WebMail Project (http://openwebmail.org) ---------- Original Message ----------- From: "Dr. Madhusoodanan M.S." To: GRADSUSR at LIST.CINECA.IT Sent: Mon, 21 Jul 2008 13:09:38 +0530 Subject: Re: Hi > Hi, > > pls refer in grads index : modify climatology seasonal > > http://www.iges.org/grads/gadoc/gadocindex.html > > regards, > Madhu ------- End of Original Message ------- From paul.earnshaw at METOFFICE.GOV.UK Mon Jul 21 06:17:18 2008 From: paul.earnshaw at METOFFICE.GOV.UK (Paul Earnshaw) Date: Mon, 21 Jul 2008 11:17:18 +0100 Subject: Converting ECMWF .pp data file into GrADS-readable data format In-Reply-To: <282632.19423.qm@web31001.mail.mud.yahoo.com> Message-ID: Ooi, It would be useful to know what the source of your data was. If this is supposed to be ECMWF Reanalysis data then can you not get it direct from ECMWF in GRIB format? The problems you are facing are in the original .pp file. Note that this is not a GRIB file, it is a file format (called PP) used by the UK Met Office. It looks like all the floating point data in the file has the wrong byte order, or is corrupted. I can only suggest that you query the source of the data and make sure what they give you is correct. Cheers, Paul On Fri, 2008-07-18 at 15:26 -0700, See Hai Ooi wrote: > Dear Mr Lincoln, > > Thank you. > > By using the downloaded xconv (version 1.91), I notice that there are two files in the attached .pp data set. Both are 32-bit ieee byteswaped. However, after selecting one of them to be converted as test.nc file, errors occurred when I opened it in GrADS (window version 1.9.0-rc1). Can you guide me further, and if possible, provide me the correct step-by-step procedures ? > > Best regards. > > Ooi > > > --- On Fri, 18/7/08, Lincoln Muniz Alves wrote: > > > From: Lincoln Muniz Alves > > Subject: Re: Converting ECMWF .pp data file into GrADS-readable data format > > To: GRADSUSR at LIST.CINECA.IT > > Date: Friday, 18 July, 2008, 10:07 PM > > Dear Ooi, > > I think that xconv it is a good tool. > > > > http://badc.nerc.ac.uk/help/software/xconv/ > > > > Lincoln > > Send instant messages to your online friends http://uk.messenger.yahoo.com -- Dr. Paul Earnshaw Global Model Development and Diagnostics Met Office FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel: +44 (0)1392 884525 Fax: +44 (0)1392 885681 email: paul.earnshaw at metoffice.gov.uk website: www.metoffice.gov.uk Met Office climate change predictions can now be viewed on Google Earth http://www.metoffice.gov.uk/research/hadleycentre/google/ From axl419 at YAHOO.COM Mon Jul 21 06:39:52 2008 From: axl419 at YAHOO.COM (See Hai Ooi) Date: Mon, 21 Jul 2008 03:39:52 -0700 Subject: Converting ECMWF .pp data file into GrADS-readable data format In-Reply-To: <1216635438.11357.11.camel@eld188.desktop.frd.metoffice.com> Message-ID: Dear Dr Paul Earnshaw, Thank you. I will follow your advice to do that. Regards. Ooi --- On Mon, 21/7/08, Paul Earnshaw wrote: > From: Paul Earnshaw > Subject: Re: Converting ECMWF .pp data file into GrADS-readable data format > To: GRADSUSR at LIST.CINECA.IT > Date: Monday, 21 July, 2008, 6:17 PM > Ooi, > > It would be useful to know what the source of your data was. If this is supposed to be ECMWF Reanalysis data then can you not get it direct from ECMWF in GRIB format? > > The problems you are facing are in the original .pp file. Note that this is not a GRIB file, it is a file format (called PP) used by the UK Met Office. > > It looks like all the floating point data in the file has the wrong byte order, or is corrupted. I can only suggest that you query the source of the data and make sure what they give you is correct. > > Cheers, > Paul > > On Fri, 2008-07-18 at 15:26 -0700, See Hai Ooi wrote: > > Dear Mr Lincoln, > > > > Thank you. > > > > By using the downloaded xconv (version 1.91), I notice that there are two files in the attached .pp data set. Both are 32-bit ieee byteswaped. However, after selecting one of them to be converted as test.nc file, errors occurred when I opened it in GrADS (window version 1.9.0-rc1). Can you guide me further, and if possible, provide me the correct step-by-step procedures ? > > > > Best regards. > > > > Ooi > > > > > > --- On Fri, 18/7/08, Lincoln Muniz Alves > wrote: > > > > > From: Lincoln Muniz Alves > > > > Subject: Re: Converting ECMWF .pp data file into > GrADS-readable data format > > > To: GRADSUSR at LIST.CINECA.IT > > > Date: Friday, 18 July, 2008, 10:07 PM > > > Dear Ooi, > > > I think that xconv it is a good tool. > > > > > > http://badc.nerc.ac.uk/help/software/xconv/ > > > > > > Lincoln > > > -- > Dr. Paul Earnshaw Global Model Development and > Diagnostics > Met Office FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel: +44 (0)1392 884525 Fax: +44 (0)1392 885681 > email: paul.earnshaw at metoffice.gov.uk website: > www.metoffice.gov.uk > > Met Office climate change predictions can now be viewed on > Google Earth > http://www.metoffice.gov.uk/research/hadleycentre/google/ Send instant messages to your online friends http://uk.messenger.yahoo.com From senya at ATMOS.UMD.EDU Mon Jul 21 07:34:21 2008 From: senya at ATMOS.UMD.EDU (Semyon Grodsky) Date: Mon, 21 Jul 2008 07:34:21 -0400 Subject: landsea mask 2x2 In-Reply-To: <48329DB7.2000001@gofferje.homelinux.org> Message-ID: Folks, if you have a landsea mask on a ---------------------- xdef 180 linear 0 2 ydef 90 linear -89 2 ---------------------- grid, please share it with me. Regards, --Senya ----------------------------------------------------------- Semyon Grodsky Computer and Space Science Building (Bldg. #224), Room 2409 Department of Atmospheric and Oceanic Science University of Maryland College Park, MD 20742 Phone: 301-405-5330 Fax: 301-314-9482 E-mail: senya at atmos.umd.edu From aspratt2000 at YAHOO.COM Mon Jul 21 08:11:00 2008 From: aspratt2000 at YAHOO.COM (Mr Aaron Pratt) Date: Mon, 21 Jul 2008 05:11:00 -0700 Subject: Using Windows-based GrADS to examine TRMM orbit data files Message-ID: Hello, ? I am trying to?create plots derived from?the TRMM 2A12 orbital product (it's the hydrometeor profile data)?. It's not gridded data, but rather in orbital swaths. The files are HDF files, so that's why I thought that GrADS might be able to open them. However, I wasn't sure if the HDF files needed to be converted into a GrADS type file (like a .ctl file and the corresponding .dat file). I was told that the files could be converted through the following link: http://disc2.nascom.nasa.gov/data/TRMM/; however I don't see anywhere on the website where this is possible. Thanks in advance for any help that can be provided. ? Aaron -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080721/fea7cb01/attachment.html From mitchell at ATMOS.WASHINGTON.EDU Mon Jul 21 12:27:02 2008 From: mitchell at ATMOS.WASHINGTON.EDU (Todd MItchell) Date: Mon, 21 Jul 2008 09:27:02 -0700 Subject: landsea mask 2x2 In-Reply-To: Message-ID: Hi Semyon, Monday 21 July 2008 I have posted various elevation and fractional-land datasets at http://jisao.washington.edu/data/elevation/ Todd Mitchell From dasilva at ALUM.MIT.EDU Mon Jul 21 12:27:21 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 21 Jul 2008 12:27:21 -0400 Subject: Converter grib2 -> grib1 for Freebsd In-Reply-To: <20080718084938.sdda7di4r48gok88@webmail.pagoria.it> Message-ID: On Fri, Jul 18, 2008 at 2:49 AM, wrote: > Arlindo, > > I need to convert grib2 -> grib1 because I want to use "clhilo" and > "regrid2" that only work with GRIB1. > I'd suggest you download the GrADS extensions for v1.9 from sf.net (FreeBSD version): http://sourceforge.net/project/showfiles.php?group_id=161773&package_id=256757 and install according to the instructions here: http://opengrads.org/wiki/index.php?title=User_Defined_Extensions (Basically you untar the file and set a couple of environment variables - I've tested this on FreeBSD). Right now these extensions only work with GrADS v1.9. The "g2" extensions allows you to read grib2 in v1.9. You can find an example of using lats4d to convert grib2 to other formats here: http://opengrads.org/doc/udxt/grib2/ My most recent lats4d scripts are attached. Also, these extensions include re() a much faster implementation of regrid2(). Documentation is here: http://opengrads.org/doc/udxt/re/ BTW, since you can read grib2 directly in GrADS v1.9, do you still need the grib2->grib1 file conversion? I'm cc'ing gradsusr as this issue is not specific of FreeBSD. Good Luck, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- A non-text attachment was scrubbed... Name: lats4d.sh Type: application/x-sh Size: 1740 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080721/d9fcf2ea/attachment.sh -------------- next part -------------- A non-text attachment was scrubbed... Name: lats4d.gs Type: application/octet-stream Size: 104733 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080721/d9fcf2ea/attachment.obj From aspratt2000 at YAHOO.COM Mon Jul 21 12:46:54 2008 From: aspratt2000 at YAHOO.COM (Mr Aaron Pratt) Date: Mon, 21 Jul 2008 09:46:54 -0700 Subject: Using PCGrads to examine TRMM HDF files Message-ID: Hello, I have downloaded TRMM orbit data files; I believe the file format is HDF. The files are th 2A12 files (this is the TMI hydrometeor profile). I wanted to use PCGrADS to open this file, but I am having problems. I know that there is a GrADSHDF executable within the PCGrADS folder that should be able to open the files directly, but that doesn't seem to work. I also try to use the TSDIS viewer to convert the HDF file (or at least a subsection) into a standard GrADS file. However, the .ctl file that was crreted still did not open (there seemed to be problems with the .dat file that was created). Is the Windows version of GrADS able to examine these TRMM orbit files. Also, are there any scripts available? Or should I try sticking with using GrADS on a Linux machine. Thanks in advance, ?Aaron -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080721/a79eb3a5/attachment.html From dasilva at ALUM.MIT.EDU Mon Jul 21 13:10:11 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 21 Jul 2008 13:10:11 -0400 Subject: Using Windows-based GrADS to examine TRMM orbit data files In-Reply-To: <81201.10219.qm@web62412.mail.re1.yahoo.com> Message-ID: On Mon, Jul 21, 2008 at 8:11 AM, Mr Aaron Pratt wrote: > Hello, > > I am trying to create plots derived from the TRMM 2A12 orbital product (it's > the hydrometeor profile data) . It's not gridded data, but rather in orbital > swaths. Well, in a way swath data is gridded, just the I/J axis do not correspond to lon/lat... The files are HDF files, so that's why I thought that GrADS might be> > able to open them. However, I wasn't sure if the HDF files needed to be > converted into a GrADS type file (like a .ctl file and the corresponding > .dat file). I was told that the files could be converted through the > following link: http://disc2.nascom.nasa.gov/data/TRMM/; however I don't see > anywhere on the website where this is possible. Thanks in advance for any > help that can be provided. > Reading HDF-EOS swatch data is a common need, but not something that can be directly done in GrADS at this point. The bottom line is, the data needs to be remapped to a regular lat/lon grid (I stress *remap* because linear interpolation is not something that should be applied to precipitation data). This can be accomplished in2 ways, 1) Implicitly by GrADS if you create a proper PDEF file. Nearest neighbor interpolation is the safest, but mass conserving schemes could be implemented as well. I believe a generic utility could be written to automatically create PDEF files for swath data, or at at least for the HDF-EOS variety. One could also write a GrADS extension along the lines of the g2() extension for doing the same thing on the fly. 2) Explicitly, by writing a file conversion program that does the remapping and writes out, say, a grads binary file. You can also convert the swath data to station data. Such utilities may already exist. For a quick inspection of the data you can always "cheat.". Swath data is after all given on a rectangular array. You can pretend the crosstrack axis is longitude and along track is latitude, and write a ctl/ddf that lies to GrADS. You will be able to look at your numbers, but the map overlay will not make any sense. As I said, quick and dirty. Brian, Jennifer: what are your plans for swath data? Do you see it as gridded or station data? Both representations can be quite useful, depending on the nature of the dataset. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From jma at COLA.IGES.ORG Mon Jul 21 13:11:36 2008 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Mon, 21 Jul 2008 13:11:36 -0400 Subject: Using Windows-based GrADS-HDF to examine TRMM orbit data files In-Reply-To: <694993.72694.qm@web62401.mail.re1.yahoo.com> Message-ID: No version of GrADS handles HDF-formatted swath data. You can try to convert it from HDF into the GrADS station-data binary format, but this would require a considerable amount of customized programming. Jennifer -- Jennifer M. Adams IGES/COLA 4041 Powder Mill Road, Suite 302 Calverton, MD 20705 jma at cola.iges.org On Jul 19, 2008, at 9:18 PM, Mr Aaron Pratt wrote: > Hello, > > I am trying to create plots derived from the TRMM 2A12 orbital > product (it's the hydrometeor profile data) . It's not gridded > data, but rather in orbital swaths. The files are HDF files, so > that's why I thought that GrADS might be able to open them. > However, I wasn't sure if the HDF files needed to be converted into > a GrADS type file (like a .ctl file and the corresponding .dat > file). I was told that the files could be converted through the > following link: http://disc2.nascom.nasa.gov/data/TRMM/; however I > don't see anywhere on the website where this is possible. Thanks in > advance for any help that can be provided. > > Aaron -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080721/2f1954f0/attachment.html From dasilva at ALUM.MIT.EDU Mon Jul 21 13:23:18 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 21 Jul 2008 13:23:18 -0400 Subject: New User In-Reply-To: Message-ID: On Sat, Jul 19, 2008 at 12:50 PM, Polar Low wrote: > Hello All, > > I`m a complete novice to grads. I have started to read the > documentation but it was a little challenging, it seems more geared up > for a university student or somone who already has a good > understanding of it. I would like to use grads for creating > meteorological products, such as nwp charts from grib data. > > Where would be a best place for me to learn how to do this, bearing in > mind my "novice" status? > > I have the option to install grads on a windows or linux based pc, > whichever it easiest for me. > I'd say *installation* is easier on Windows: http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows Now, on the Linux/Windows issue let's not start a holly war. All depends on your comfort level with them, and whether you can justify the investment in learning something new. You can also install cygwin on Windows for a very Linux like environment. Depending on the age of your hardware, Linux can have a much smaller front print and feel much faster, than say, Window XP. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From Chris.Lynnes at NASA.GOV Mon Jul 21 13:39:18 2008 From: Chris.Lynnes at NASA.GOV (Christopher Lynnes) Date: Mon, 21 Jul 2008 13:39:18 -0400 Subject: Using Windows-based GrADS to examine TRMM orbit data files In-Reply-To: <77fcd6b20807211010q7ecb8004j7f31a418d14bb07@mail.gmail.com> Message-ID: On Jul 21, 2008, at 1:10 PM, Arlindo da Silva wrote: > On Mon, Jul 21, 2008 at 8:11 AM, Mr Aaron Pratt > wrote: >> Hello, >> >> I am trying to create plots derived from the TRMM 2A12 orbital >> product (it's >> the hydrometeor profile data) . It's not gridded data, but rather >> in orbital >> swaths. > > Well, in a way swath data is gridded, just the I/J axis do not > correspond to lon/lat... > > > Brian, Jennifer: what are your plans for swath data? Do you see it as > gridded or station data? Both representations can be quite useful, > depending on the nature of the dataset. > > Arlindo Along those lines, we have been experimenting at the GES DISC with reformatting HDF-EOS swath data into NetCDF CF-1 compliant format, using CF-1's two-dimensional latitude and longitude variables, i.e., WITHOUT explicitly reprojecting the data. Tools such as IDV and McIDAS-V have been able to understand this to the extent that they can ingest the netCDF versions of the data and plot the data correctly on a map. Ferret can read and plot them, but with simple x-y axes, not on a map. Unfortunately, I have been unable to use sdfopen in GrADS to read these files ("SDF file has no discernable X coordinate."). It would be nice if... Here is an example of the on-the-fly reformatting of AIRS Level 2 swath to NetCDF: http://aurapar2u.ecs.nasa.gov/airspar1/daac-bin/OTF/HTTP_services.cgi?FILENAME=%2Fdata%2Fs4pa%2FAqua_AIRS_Level2%2FAIRX2RET.005%2F2008%2F202%2FAIRS.2008.07.20.094.L2.RetStd.v5.2.2.0.G08203122008.hdf&LABEL=AIRS.2008.07.20.094.L2.RetStd.v5.2.2.0.G08203122008.nc&SERVICE=NetCDF&VERSION=1 -- Christopher Lynnes NASA/GSFC, Code 610.2 301-614-5185 From smcmillan at PLANALYTICS.COM Mon Jul 21 15:05:21 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Mon, 21 Jul 2008 14:05:21 -0500 Subject: Using Shape Files in GrADS Message-ID: I am able to display shape file polygons in GrADS. However, I would like to be able to color-fill the component polygons according to a value, say, temperature. Is this possible? If so, please share how to do it. I can draw all the lines with the 'shp_lines....' command, and an overall/inclusive filled polygon with the 'shp_polyf...' command, specifying a single color with the 'set line...' command. However, I want to be able to control individual shape poly fills within the inclusive boundaries, such as states/provinces, of the shape file. I am currently experimenting with the "admin98" set. The associated .dbf file contains 14 columns/fields of geographical data, including what appear to be unique identifiers (US20, USA-KAN, etc.). The command line output indicates the following syntax, with 1 to 3 arguments: Usage: shp_polyf shapefile [object] How do I specify the object? Is the object the 2nd argument? What is the 3rd argument, and how do I specify it? The idea is that my script would read a text file as input, where individual shapes (such as state/province boundaries) are associated with temperature or some other variable. I am using GrADS v1.9.0-rc1 (installed from OpenGrADS), on a Windows XP platform. Any help on this would be appreciated. Stephen McMillan *************************************************** 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/20080721/2b3acf94/attachment.html From mike.bosilovich at GMAIL.COM Mon Jul 21 15:44:42 2008 From: mike.bosilovich at GMAIL.COM (Mike Bosilovich) Date: Mon, 21 Jul 2008 15:44:42 -0400 Subject: Using Windows-based GrADS to examine TRMM orbit data files In-Reply-To: <81201.10219.qm@web62412.mail.re1.yahoo.com> Message-ID: If you really want the orbital data, rather than a higher level gridded product, I'm pretty sure that you need an intermediate step to put the swath(s) in a grid that grads is friendly with. For MODIS L2 swaths, I used hdflook (on linux) to save lat-lon grids from the swaths. However it was very cumbersome, and you would not want to do that with large amounts of data. I haven't seen grads open swaths (and the other replies confirm that), but I might make use of such a utility. So, you need an HDF tool, for the PC that will open swath, and write a grid. That grid would go into PCgrads. Good Luck, Mike On Mon, Jul 21, 2008 at 8:11 AM, Mr Aaron Pratt wrote: > Hello, > > I am trying to create plots derived from the TRMM 2A12 orbital product > (it's the hydrometeor profile data) . It's not gridded data, but rather in > orbital swaths. The files are HDF files, so that's why I thought that GrADS > might be able to open them. However, I wasn't sure if the HDF files needed > to be converted into a GrADS type file (like a .ctl file and the > corresponding .dat file). I was told that the files could be converted > through the following link: http://disc2.nascom.nasa.gov/data/TRMM/; > however I don't see anywhere on the website where this is possible. Thanks > in advance for any help that can be provided. > > Aaron > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080721/f163136d/attachment.html From stergios.misios at ZMAW.DE Mon Jul 21 15:36:06 2008 From: stergios.misios at ZMAW.DE (Stergios Misios) Date: Mon, 21 Jul 2008 21:36:06 +0200 Subject: Producing a movie Message-ID: Hello all, i want to produce an animation using GRAds. I have a grads output meta file with several frames included, and i miss the last step in order to have an animation. Can i have a help on that? which program should i use to have a gif or an mpeg animation? Or my way to have a single grads metafile with all frames included in the same file is totally wrong? Thank you in advance, Stergios From jim at ATMOS.COLOSTATE.EDU Mon Jul 21 16:07:16 2008 From: jim at ATMOS.COLOSTATE.EDU (Jim Benedict) Date: Mon, 21 Jul 2008 14:07:16 -0600 Subject: Current time Message-ID: GrADS users, I can't find an answer to this anywhere: In my simple GrADS script, I define a few variables and create multiple plots (the 'set grads off' command has been used for these plots). On the final plot, though, I'd like a time stamp of the current time (that is, the time at which the plot was generated) to appear. I've tried using "set grads on" for this final plot, but with that command only the "GrADS: COLA/IGES" appears and not the timestamp. Any help would be appreciated! Thanks in advance, Jim From esertel at ENVSCI.RUTGERS.EDU Mon Jul 21 16:14:54 2008 From: esertel at ENVSCI.RUTGERS.EDU (Elif SERTEL) Date: Mon, 21 Jul 2008 16:14:54 -0400 Subject: Current time In-Reply-To: <3203BD9E-816B-49EE-A353-0E30CC0B8C8E@atmos.colostate.edu> Message-ID: You should try : set timelab on On Mon, July 21, 2008 4:07 pm, Jim Benedict wrote: > GrADS users, > > I can't find an answer to this anywhere: In my simple GrADS script, I > define a few variables and create multiple plots (the 'set grads off' > command has been used for these plots). On the final plot, though, > I'd like a time stamp of the current time (that is, the time at which > the plot was generated) to appear. I've tried using "set grads on" > for this final plot, but with that command only the "GrADS: COLA/IGES" > appears and not the timestamp. Any help would be appreciated! > > Thanks in advance, > Jim > -- Elif SERTEL, Ph. D. Visiting Fulbright Researcher Rutgers University Department of Environmental Sciences http://www.envsci.rutgers.edu/~esertel/ Elif SERTEL, Ph. D. Research Associate Istanbul Technical University Civil Engineering Faculty, Remote Sensing Division http://www.ins.itu.edu.tr/esertel/ From dd_ren at ROSSBY.METR.OU.EDU Mon Jul 21 16:35:41 2008 From: dd_ren at ROSSBY.METR.OU.EDU (Ren Diandong) Date: Mon, 21 Jul 2008 15:35:41 -0500 Subject: Producing a movie In-Reply-To: <4884E526.5080209@zmaw.de> Message-ID: Matlab read in the frames and create an avi. Diandong On Mon, 2008-07-21 at 21:36 +0200, Stergios Misios wrote: > Hello all, > i want to produce an animation using GRAds. I have a grads output meta > file with several frames included, and i miss the last step in order to > have an animation. Can i have a help on that? which program should i use > to have a gif or an mpeg animation? > Or my way to have a single grads metafile with all frames included in > the same file is totally wrong? > > Thank you in advance, > Stergios > From rhart at FSU.EDU Mon Jul 21 16:41:36 2008 From: rhart at FSU.EDU (Robert Hart) Date: Mon, 21 Jul 2008 16:41:36 -0400 Subject: Producing a movie In-Reply-To: <1216672541.27541.2.camel@romulus.metr.ou.edu> Message-ID: Stergios, If you create individual gif images from the frames in GrADS, you can quite easily create an animated gif using gifsicle: http://www.lcdf.org/gifsicle/ Also includes compressed animated gifs to keep file size down, although the cpu overhead grows. Bob ----- Original Message ----- From: Ren Diandong Date: Monday, July 21, 2008 4:36 pm Subject: Re: Producing a movie To: GRADSUSR at LIST.CINECA.IT > Matlab read in the frames and create an avi. > > Diandong > > On Mon, 2008-07-21 at 21:36 +0200, Stergios Misios wrote: > > Hello all, > > i want to produce an animation using GRAds. I have a grads output > meta> file with several frames included, and i miss the last step > in order to > > have an animation. Can i have a help on that? which program > should i use > > to have a gif or an mpeg animation? > > Or my way to have a single grads metafile with all frames > included in > > the same file is totally wrong? > > > > Thank you in advance, > > Stergios > > > From dd_ren at ROSSBY.METR.OU.EDU Mon Jul 21 18:27:54 2008 From: dd_ren at ROSSBY.METR.OU.EDU (Ren Diandong) Date: Mon, 21 Jul 2008 17:27:54 -0500 Subject: image files Message-ID: Hi, Stergios: If you already created the gifs, you can also use UNIX utility to create on single image file. For example, CONVERT may well suffice. Diandong From sibachirnordine at HOTMAIL.COM Tue Jul 22 02:08:06 2008 From: sibachirnordine at HOTMAIL.COM (sibachir nordine) Date: Tue, 22 Jul 2008 06:08:06 +0000 Subject: Current time In-Reply-To: <52304.128.6.30.197.1216671294.squirrel@webmail.envsci.rutgers.edu> Message-ID: Hi all; There is the way to display date time base of the chars . 'q time' res = subwrd(result,3) basedate = substr(res,4,12) basetime = substr(res,1,3) res = subwrd(result,6) baseday=substr(res,1,3) draw title CFO MIN 2m Temperature \ Valid To: ' baseday ', ' basedate ' ' basetime 'printim tmin2m-'basedate'-'basetime'.jpg x800 y600 white I hope it will help you. Sibachir Nordine. System Engineer. Presidency of Meteorology & Environment. Saudi Arabia. _________________________________________________________________ With Windows Live for mobile, your contacts travel with you. http://www.windowslive.com/mobile/overview.html?ocid=TXT_TAGLM_WL_mobile_072008 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080722/30759e40/attachment.html From sibachirnordine at HOTMAIL.COM Tue Jul 22 02:13:41 2008 From: sibachirnordine at HOTMAIL.COM (sibachir nordine) Date: Tue, 22 Jul 2008 06:13:41 +0000 Subject: Current time In-Reply-To: <52304.128.6.30.197.1216671294.squirrel@webmail.envsci.rutgers.edu> Message-ID: Hi all, Here is the way to display date time base of the charts. ''q time' res = subwrd(result,3) basedate = substr(res,4,12) basetime = substr(res,1,3) res = subwrd(result,6) baseday=substr(res,1,3) draw title CFO MIN 2m Temperature \ Valid To: ' baseday ', ' basedate ' ' basetime 'printim tmin2m-'basedate'-'basetime'.jpg x800 y600 white Best regards. Sibachir Nordine System Engineer. Central Forecasting Office Presidency of Meteorology & Environment Saudi Arabia _________________________________________________________________ Stay in touch when you're away with Windows Live Messenger. http://www.windowslive.com/messenger/overview.html?ocid=TXT_TAGLM_WL_messenger2_072008 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080722/02f477b2/attachment.html From Isla.Gilmour at DESHAW.COM Tue Jul 22 03:33:24 2008 From: Isla.Gilmour at DESHAW.COM (Gilmour, Isla) Date: Tue, 22 Jul 2008 03:33:24 -0400 Subject: Reading GFS GRIB2 files: issues with gribmap Message-ID: Using GrADS (v2.0.a2), wgrib2 (v1.7.7c) , g2ctl (v0.0.3h) and with perl installed, on UNIX (SunOS i386) Having downloaded grib2 GFS files of the form gfs.tccz.pgrb2.fxxx.grib2 from NCEP (as described at http://www.nco.ncep.noaa.gov/pmb/products/gfs/ ) I created a ctl file using g2ctl. This looks relatively ok (see end of email for header), other than the tdef which is listed as 1mo rather than something more reasonable like 3h or 12h (the forecast leadtime varies from 3hrly 0-180h, then 12hrly out to 384h) so that's the second question... can GrADS ctl files deal with non-uniform time steps? Leaving tdef aside, the next step is to use gribmap (the one which comes with grads-2.0.a2) , but this failed: > g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl > gribmap -i gfs.ctl Open Error: Data file type invalid --> The invalid description file record is: --> dtype grib2 The data file was not opened. Is there an alternative version of gribmap to use for GRIB2 files? Or a workaround which preserves the data in grib2 format? (Or are these the "strange" gfs files mentioned in a previous thread?) Any help would be appreciated. Isla Gilmour ---------------------------------------- gfs.ctl header dset /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 index /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2.idx undef 9.999E+20 title /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 * produced by g2ctl v0.0.3h * griddef=1:0:(360 x 181):grid_template=0: lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160 dtype grib2 ydef 181 linear -90.000000 1 xdef 360 linear 0.000000 1.000000 tdef 1 linear 00Z14jul2008 1mo * PROFILE hPa zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000 options pascals vars 139 [list of vars, which tally with wgrib2 listing] From Daniele.Gandini at ARPA.PIEMONTE.IT Tue Jul 22 03:58:36 2008 From: Daniele.Gandini at ARPA.PIEMONTE.IT (Daniele Gandini) Date: Tue, 22 Jul 2008 09:58:36 +0200 Subject: Reading GFS GRIB2 files: issues with gribmap Message-ID: I think you have to separate the original file in two ones, the first from 0 to 180hr and the second from 192 to 384 hr. You can manually change time increment in ctl file produced by g2ctl if isn't ok and then you can run gribmap. Bye Daniel ----- Original Message ----- From: "Gilmour, Isla" To: Sent: Tuesday, July 22, 2008 9:33 AM Subject: Reading GFS GRIB2 files: issues with gribmap Using GrADS (v2.0.a2), wgrib2 (v1.7.7c) , g2ctl (v0.0.3h) and with perl installed, on UNIX (SunOS i386) Having downloaded grib2 GFS files of the form gfs.tccz.pgrb2.fxxx.grib2 from NCEP (as described at http://www.nco.ncep.noaa.gov/pmb/products/gfs/ ) I created a ctl file using g2ctl. This looks relatively ok (see end of email for header), other than the tdef which is listed as 1mo rather than something more reasonable like 3h or 12h (the forecast leadtime varies from 3hrly 0-180h, then 12hrly out to 384h) so that's the second question... can GrADS ctl files deal with non-uniform time steps? Leaving tdef aside, the next step is to use gribmap (the one which comes with grads-2.0.a2) , but this failed: > g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl > gribmap -i gfs.ctl Open Error: Data file type invalid --> The invalid description file record is: --> dtype grib2 The data file was not opened. Is there an alternative version of gribmap to use for GRIB2 files? Or a workaround which preserves the data in grib2 format? (Or are these the "strange" gfs files mentioned in a previous thread?) Any help would be appreciated. Isla Gilmour ---------------------------------------- gfs.ctl header dset /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 index /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2.idx undef 9.999E+20 title /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 * produced by g2ctl v0.0.3h * griddef=1:0:(360 x 181):grid_template=0: lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160 dtype grib2 ydef 181 linear -90.000000 1 xdef 360 linear 0.000000 1.000000 tdef 1 linear 00Z14jul2008 1mo * PROFILE hPa zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000 options pascals vars 139 [list of vars, which tally with wgrib2 listing] From jma at COLA.IGES.ORG Tue Jul 22 09:58:45 2008 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Tue, 22 Jul 2008 09:58:45 -0400 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: <42E13EA8EEB5EE40B800168F8CA164AD016A6E78@mailnyc4.nyc.deshaw.com> Message-ID: On Jul 22, 2008, at 3:33 AM, Gilmour, Isla wrote: > Using GrADS (v2.0.a2), wgrib2 (v1.7.7c) , g2ctl (v0.0.3h) and with > perl > installed, on UNIX (SunOS i386) > > Having downloaded grib2 GFS files of the form > gfs.tccz.pgrb2.fxxx.grib2 > from NCEP > (as described at http://www.nco.ncep.noaa.gov/pmb/products/gfs/ ) > I created a ctl file using g2ctl. This looks relatively ok (see > end of > email for header), other than the tdef which is listed as 1mo rather > than something more reasonable like 3h or 12h (the forecast leadtime > varies from 3hrly 0-180h, then 12hrly out to 384h) so that's the > second > question... can GrADS ctl files deal with non-uniform time steps? You need two descriptors, as Daniele described. tdef 33 linear 00z22jul08 12hr tdef 61 linear 00z22jul08 3hr > > Leaving tdef aside, the next step is to use gribmap (the one which > comes > with grads-2.0.a2) , but this failed: > >> g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl >> gribmap -i gfs.ctl > Open Error: Data file type invalid > --> The invalid description file record is: > --> dtype grib2 > The data file was not opened. The 'dtype grib2' option should work with gribmap (the one that comes with 2.0.a2). Please check your $PATH or use a complete path to gribmap to make sure you're using the right one. Jennifer > Is there an alternative version of gribmap to use for GRIB2 files? > Or a > workaround which preserves the data in grib2 format? (Or are these > the > "strange" gfs files mentioned in a previous thread?) > Any help would be appreciated. > Isla Gilmour > > ---------------------------------------- > > gfs.ctl header > > dset /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 > index /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2.idx > undef 9.999E+20 > title /data/weather/gfs/20080714/00/gfs.t00z.pgrbf84.grib2 > * produced by g2ctl v0.0.3h > * griddef=1:0:(360 x 181):grid_template=0: lat-lon grid:(360 x 181) > units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to > -90.000000 > by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160 > > dtype grib2 > ydef 181 linear -90.000000 1 > xdef 360 linear 0.000000 1.000000 > tdef 1 linear 00Z14jul2008 1mo > * PROFILE hPa > zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 > 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 > 10000 > 7000 5000 3000 2000 1000 > options pascals > vars 139 > [list of vars, which tally with wgrib2 listing] -- 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/20080722/56ca0356/attachment.html From Isla.Gilmour at DESHAW.COM Tue Jul 22 10:03:03 2008 From: Isla.Gilmour at DESHAW.COM (Gilmour, Isla) Date: Tue, 22 Jul 2008 10:03:03 -0400 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: A<5A4C6D49-52BA-467F-87BA-D652BB80435A@cola.iges.org> Message-ID: Jennifer Thanks for your response. Indeed, I checked the $PATH, and tried using a complete path to gribmap, but continued to get the same error. Has anyone used gribmap from the grads v2.a.02 "package" successfully with GFS pgrbf.grib2 files? Isla ________________________________ Leaving tdef aside, the next step is to use gribmap (the one which comes with grads-2.0.a2) , but this failed: g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl gribmap -i gfs.ctl Open Error: Data file type invalid --> The invalid description file record is: --> dtype grib2 The data file was not opened. The 'dtype grib2' option should work with gribmap (the one that comes with 2.0.a2). Please check your $PATH or use a complete path to gribmap to make sure you're using the right one. Jennifer -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080722/43ac7d7f/attachment.html From klevey at CUSTOMWEATHER.COM Tue Jul 22 11:52:51 2008 From: klevey at CUSTOMWEATHER.COM (Kevin M Levey) Date: Tue, 22 Jul 2008 08:52:51 -0700 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: <42E13EA8EEB5EE40B800168F8CA164AD016A6E86@mailnyc4.nyc.deshaw.com> Message-ID: TUE 22JUL08: 0850pDT Hi All I've had success using GRADS 2.a.02 to create CTL and IDX files, using gribmap (which really is "grib2map" in version 2.a.02 but still called "gribmap"). Here is my process (I am showing the 00hour as an example - do the same for ALL hours) perl g2ctl.pl $REMOTE_DIR1/gfs.t12z.pgrbf00.grb2 $TEMP_DIR/ gfs2.t12z.pgrbf00.idx >$TEMP_DIR/gfs2_000.ctl This will create both the CTL and IDX files in *my* TEMP_DIR using data from *my* data source folder (it can be whatever you want it to be) Then the following: $GRADS2_DIR/gribmap -0 -i $TEMP_DIR/gfs2_000.ctl where GRADS2_DIR=/usr/local/bin/grads-2.0.a2/bin (*my* location of the latest version of grads) Fire up grads 2 and open gfs2_000.ctl This process works perfectly for me. Regards Kevin M Levey, MSc in Oceans and Atmospheric Sciences (University of Cape Town) Director of Meteorological Operations CustomWeather, Inc. San Francisco, California, USA "Taking the World by Storm!" http://www.customweather.com http://www.myforecast.com http://www.1stweather.com cell: 415-794-0411 work: 415-777-3566 email: klevey at customweather.com On Jul 22, 2008, at 07/22/08 - 7:03 AM, Gilmour, Isla wrote: > Jennifer > > Thanks for your response. Indeed, I checked the $PATH, and tried > using a complete path to gribmap, but continued to get the same > error. Has anyone used gribmap from the grads v2.a.02 "package" > successfully with GFS pgrbf.grib2 files? > > Isla > >> Leaving tdef aside, the next step is to use gribmap (the one which >> comes >> with grads-2.0.a2) , but this failed: >> >>> g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl >>> gribmap -i gfs.ctl >> Open Error: Data file type invalid >> --> The invalid description file record is: >> --> dtype grib2 >> The data file was not opened. > The 'dtype grib2' option should work with gribmap (the one that > comes with 2.0.a2). Please check your $PATH or use a complete path > to gribmap to make sure you're using the right one. > > Jennifer > >> -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080722/f1914191/attachment.html From jma at COLA.IGES.ORG Tue Jul 22 12:40:27 2008 From: jma at COLA.IGES.ORG (Jennifer Adams) Date: Tue, 22 Jul 2008 12:40:27 -0400 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: Message-ID: On Jul 22, 2008, at 11:52 AM, Kevin M Levey wrote: > gribmap (which really is "grib2map" in version 2.a.02 but still > called "gribmap"). This is not exactly correct. There are two separate utilities for scanning a grib file (gribscan and grib2scan). But the work of creating the index files is done by one utility, gribmap. Inside gribmap, the code for grib1 and grib2 files is completely separate, and the messages it prints to the users will say 'grib1map' or grib2map' depending on the which version of grib is being indexed. Jennifer > Here is my process (I am showing the 00hour as an example - do the > same for ALL hours) > > perl g2ctl.pl $REMOTE_DIR1/gfs.t12z.pgrbf00.grb2 $TEMP_DIR/ > gfs2.t12z.pgrbf00.idx >$TEMP_DIR/gfs2_000.ctl > > This will create both the CTL and IDX files in *my* TEMP_DIR using > data from *my* data source folder (it can be whatever you want it > to be) > > Then the following: > > $GRADS2_DIR/gribmap -0 -i $TEMP_DIR/gfs2_000.ctl > > where GRADS2_DIR=/usr/local/bin/grads-2.0.a2/bin (*my* location of > the latest version of grads) > > Fire up grads 2 and open gfs2_000.ctl > > This process works perfectly for me. > > Regards > > Kevin M Levey, MSc in Oceans and Atmospheric Sciences (University > of Cape Town) > Director of Meteorological Operations > CustomWeather, Inc. > San Francisco, California, USA > > "Taking the World by Storm!" > > http://www.customweather.com > http://www.myforecast.com > http://www.1stweather.com > > cell: 415-794-0411 > work: 415-777-3566 > email: klevey at customweather.com > > > > > On Jul 22, 2008, at 07/22/08 - 7:03 AM, Gilmour, Isla wrote: > >> Jennifer >> >> Thanks for your response. Indeed, I checked the $PATH, and tried >> using a complete path to gribmap, but continued to get the same >> error. Has anyone used gribmap from the grads v2.a.02 "package" >> successfully with GFS pgrbf.grib2 files? >> >> Isla >> >>> Leaving tdef aside, the next step is to use gribmap (the one >>> which comes >>> with grads-2.0.a2) , but this failed: >>> >>>> g2ctl -0 /data/weather/gfs/20080714/00/gfs.t00z.pgrbf > gfs.ctl >>>> gribmap -i gfs.ctl >>> Open Error: Data file type invalid >>> --> The invalid description file record is: >>> --> dtype grib2 >>> The data file was not opened. >> The 'dtype grib2' option should work with gribmap (the one that >> comes with 2.0.a2). Please check your $PATH or use a complete path >> to gribmap to make sure you're using the right one. >> >> 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/20080722/32587a92/attachment.html From smcmillan at PLANALYTICS.COM Tue Jul 22 12:54:25 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Tue, 22 Jul 2008 11:54:25 -0500 Subject: Using Shape Files in GrADS In-Reply-To: Message-ID: This is a re-send of inquiry sent yesterday. Perhaps no one in user community has actually used the shape file feature. If so, then this is intended for the person who contributed this function to v1.9.0-rc1 for Windows. All I need is an EXAMPLE or two of the correct syntax when an "object" is involved, any additional required command arguments, and any pre-requisite commands. I had done a search of the user archives regarding shape files, but found only my previous relevant questions and partial responses. Thanks-- Stephen McMillan Stephen R. McMillan/Planalytics 07/21/2008 02:05 PM To GRADSUSR at LIST.CINECA.IT cc Subject Using Shape Files in GrADS I am able to display shape file polygons in GrADS. However, I would like to be able to color-fill the component polygons according to a value, say, temperature. Is this possible? If so, please share how to do it. I can draw all the lines with the 'shp_lines....' command, and an overall/inclusive filled polygon with the 'shp_polyf...' command, specifying a single color with the 'set line...' command. However, I want to be able to control individual shape poly fills within the inclusive boundaries, such as states/provinces, of the shape file. I am currently experimenting with the "admin98" set. The associated .dbf file contains 14 columns/fields of geographical data, including what appear to be unique identifiers (US20, USA-KAN, etc.). The command line output indicates the following syntax, with 1 to 3 arguments: Usage: shp_polyf shapefile [object] How do I specify the object? Is the object the 2nd argument? What is the 3rd argument, and how do I specify it? The idea is that my script would read a text file as input, where individual shapes (such as state/province boundaries) are associated with temperature or some other variable. I am using GrADS v1.9.0-rc1 (installed from OpenGrADS), on a Windows XP platform. Any help on this would be appreciated. Stephen McMillan *************************************************** 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/20080722/4a01395e/attachment.html From hersala at GMAIL.COM Tue Jul 22 14:04:49 2008 From: hersala at GMAIL.COM (hersala) Date: Tue, 22 Jul 2008 15:04:49 -0300 Subject: a question Message-ID: Hello Grads users: Could someone help me with this: How can I show this "animation" more slowly or by steps, i.e., by clicking the mouse o pressing Enter ? Thanks in advance Hern?n Sala Instituto Ant?rtico Argentino c draw title 1000mb Heights set mproj sps set lat -90 -50 set lon 0 360 set gxout shaded set t 1 31 d hgt run cbar -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080722/f32e7588/attachment.html From ela at COLA.IGES.ORG Tue Jul 22 15:58:29 2008 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Tue, 22 Jul 2008 15:58:29 -0400 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: <807534308.94751216756676547.JavaMail.root@mail.iges.org> Message-ID: Hi Kevin, If you can use g2ctl to generate the index file, why do you run gribmap at all? Eric L. Altshuler Assistant Research Scientist Center for Ocean-Land-Atmosphere Studies 4041 Powder Mill Road, Suite 302 Calverton, MD 20705-3106 USA E-mail: ela at cola.iges.org Phone: (301) 902-1257 Fax: (301) 595-9793 ----- Original Message ----- perl g2ctl.pl $REMOTE_DIR1/gfs.t12z.pgrbf00.grb2 $TEMP_DIR/ gfs2.t12z.pgrbf00.idx >$TEMP_DIR/gfs2_000.ctl This will create both the CTL and IDX files in *my* TEMP_DIR using data from *my* data source folder (it can be whatever you want it to be) Then the following: $GRADS2_DIR/gribmap -0 -i $TEMP_DIR/gfs2_000.ctl ---------------------------- From klevey at CUSTOMWEATHER.COM Tue Jul 22 16:47:22 2008 From: klevey at CUSTOMWEATHER.COM (Kevin M Levey) Date: Tue, 22 Jul 2008 13:47:22 -0700 Subject: Reading GFS GRIB2 files: issues with gribmap In-Reply-To: <1948737129.94771216756709703.JavaMail.root@mail.iges.org> Message-ID: To create the IDX file. On Jul 22, 2008, at 07/22/08 - 12:58 PM, Eric Altshuler wrote: > Hi Kevin, > > If you can use g2ctl to generate the index file, why do you run > gribmap at all? > > > Eric L. Altshuler > Assistant Research Scientist > Center for Ocean-Land-Atmosphere Studies > 4041 Powder Mill Road, Suite 302 > Calverton, MD 20705-3106 > USA > > E-mail: ela at cola.iges.org > Phone: (301) 902-1257 > Fax: (301) 595-9793 > > ----- Original Message ----- > > perl g2ctl.pl $REMOTE_DIR1/gfs.t12z.pgrbf00.grb2 $TEMP_DIR/ > gfs2.t12z.pgrbf00.idx >$TEMP_DIR/gfs2_000.ctl > > This will create both the CTL and IDX files in *my* TEMP_DIR using > data from *my* data source folder (it can be whatever you want it to > be) > > Then the following: > > $GRADS2_DIR/gribmap -0 -i $TEMP_DIR/gfs2_000.ctl > > ---------------------------- From Charles.Seman at NOAA.GOV Tue Jul 22 19:00:05 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Tue, 22 Jul 2008 19:00:05 -0400 Subject: a question In-Reply-To: <5b40c90807221104h73987111m26aaf794d100fc0a@mail.gmail.com> Message-ID: Hersala, A couple more links for documentation on animation in GrADS: http://grads.iges.org/grads/gadoc/animation.html and http://ccr.meteor.wisc.edu/model/processing/grads/animation.gs Each of the above web sites show ways you can slow down the speed of the animation. I hope this helps, Chuck hersala wrote: > Hello Grads users: > Could someone help me with this: > How can I show this "animation" more slowly or by steps, i.e., by > clicking the mouse o pressing Enter ? > Thanks in advance > > Hern?n Sala > Instituto Ant?rtico Argentino > > > > c > draw title 1000mb Heights > set mproj sps > set lat -90 -50 > set lon 0 360 > set gxout shaded > set t 1 31 > d hgt > run cbar > -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." From dasilva at ALUM.MIT.EDU Wed Jul 23 21:12:18 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Wed, 23 Jul 2008 21:12:18 -0400 Subject: Using Shape Files in GrADS In-Reply-To: Message-ID: On Tue, Jul 22, 2008 at 12:54 PM, Stephen R McMillan wrote: > > This is a re-send of inquiry sent yesterday. Perhaps no one in user > community has actually used the shape file feature. If so, then this is > intended for the person who contributed this function to v1.9.0-rc1 for > Windows. All I need is an EXAMPLE or two of the correct syntax when an > "object" is involved, any additional required command arguments, and any > pre-requisite commands. > > I had done a search of the user archives regarding shape files, but found > only my previous relevant questions and partial responses. > I just found a silly bug that was preventing the shp_polyf feature from working. I am reconstructing my Win32 build environment and I will post a revised shape.dll as soon as I can. I believe Graziano has fixed the issue with line thickness, etc. The fixed sources have been checked in at the OpenGrADS CVS repository if you rather buid from sources. The syntax should be something like: ga-> shp_polyg admin98 8 to plot object number "8"; In this particular case, the admin98.shp file as 2,604 shapes (objects), so you can enter a number from 0 to 2,603 (it is 0-based as in C). Each shape (object) corresponds to a single small polygon, although some objects may not be plotted at all (it must be closed for it to be filled). Arlindo > Thanks-- > Stephen McMillan > > > Stephen R. McMillan/Planalytics > > 07/21/2008 02:05 PM > > To > GRADSUSR at LIST.CINECA.IT > cc > Subject > Using Shape Files in GrADS > > > > I am able to display shape file polygons in GrADS. However, I would like to > be able to color-fill the component polygons according to a value, say, > temperature. Is this possible? If so, please share how to do it. > > I can draw all the lines with the 'shp_lines....' command, and an > overall/inclusive filled polygon with the 'shp_polyf...' command, specifying > a single color with the 'set line...' command. However, I want to be able to > control individual shape poly fills within the inclusive boundaries, such as > states/provinces, of the shape file. I am currently experimenting with the > "admin98" set. The associated .dbf file contains 14 columns/fields of > geographical data, including what appear to be unique identifiers (US20, > USA-KAN, etc.). > > The command line output indicates the following syntax, with 1 to 3 > arguments: > Usage: shp_polyf shapefile [object] > > How do I specify the object? Is the object the 2nd argument? What is the 3rd > argument, and how do I specify it? > > The idea is that my script would read a text file as input, where individual > shapes (such as state/province boundaries) are associated with temperature > or some other variable. > > I am using GrADS v1.9.0-rc1 (installed from OpenGrADS), on a Windows XP > platform. > > Any help on this would be appreciated. > > Stephen McMillan > > > *************************************************** > 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. > *************************************************** > -- Arlindo da Silva dasilva at alum.mit.edu From smcmillan at PLANALYTICS.COM Thu Jul 24 12:32:58 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Thu, 24 Jul 2008 11:32:58 -0500 Subject: Using Shape Files in GrADS In-Reply-To: <77fcd6b20807231812g13455ed9jf74146f99920688e@mail.gmail.com> Message-ID: Arlindo, I'm not yet familiar with how to build from sources, but I appreciate your mentioning that option. Regarding the usage syntax, did you mean shp_polyf instead of shp_polyg? Otherwise, what is "polyg" supposed to do? Also, in their present state (on current version I'm using), both shp_polyf and shp_lines require either 1 or 3 arguments. If I try two args as in your example, I get the GrADS error "invalid number of arguments; need 1 or 3 but got 2". If this error is valid, then the user will need to know the following for proper use: Arg1 Shapefile name such as "admin98" Arg2 Presumably some object identifier, such as the "8" you mentioned. Also, how to determine the range of acceptable object identifiers. Right now, even with admin98, how do I know which sub-polygon 8, 100, 2603, etc. applies to? For now, I'll just assume they correspond to the order listed in the .dbf component of the shapefile. Arg3 What is its purpose? How do I identify the sutable range of Arg3? I'll look forward to your fix of that "silly little bug"! Thanks for your response, and keep up the great work! Stephen McMillan Arlindo da Silva Sent by: GRADSUSR at LIST.CINECA.IT 07/23/2008 08:12 PM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: Using Shape Files in GrADS On Tue, Jul 22, 2008 at 12:54 PM, Stephen R McMillan wrote: > > This is a re-send of inquiry sent yesterday. Perhaps no one in user > community has actually used the shape file feature. If so, then this is > intended for the person who contributed this function to v1.9.0-rc1 for > Windows. All I need is an EXAMPLE or two of the correct syntax when an > "object" is involved, any additional required command arguments, and any > pre-requisite commands. > > I had done a search of the user archives regarding shape files, but found > only my previous relevant questions and partial responses. > I just found a silly bug that was preventing the shp_polyf feature from working. I am reconstructing my Win32 build environment and I will post a revised shape.dll as soon as I can. I believe Graziano has fixed the issue with line thickness, etc. The fixed sources have been checked in at the OpenGrADS CVS repository if you rather buid from sources. The syntax should be something like: ga-> shp_polyg admin98 8 to plot object number "8"; In this particular case, the admin98.shp file as 2,604 shapes (objects), so you can enter a number from 0 to 2,603 (it is 0-based as in C). Each shape (object) corresponds to a single small polygon, although some objects may not be plotted at all (it must be closed for it to be filled). Arlindo > Thanks-- > Stephen McMillan > > > Stephen R. McMillan/Planalytics > > 07/21/2008 02:05 PM > > To > GRADSUSR at LIST.CINECA.IT > cc > Subject > Using Shape Files in GrADS > > > > I am able to display shape file polygons in GrADS. However, I would like to > be able to color-fill the component polygons according to a value, say, > temperature. Is this possible? If so, please share how to do it. > > I can draw all the lines with the 'shp_lines....' command, and an > overall/inclusive filled polygon with the 'shp_polyf...' command, specifying > a single color with the 'set line...' command. However, I want to be able to > control individual shape poly fills within the inclusive boundaries, such as > states/provinces, of the shape file. I am currently experimenting with the > "admin98" set. The associated .dbf file contains 14 columns/fields of > geographical data, including what appear to be unique identifiers (US20, > USA-KAN, etc.). > > The command line output indicates the following syntax, with 1 to 3 > arguments: > Usage: shp_polyf shapefile [object] > > How do I specify the object? Is the object the 2nd argument? What is the 3rd > argument, and how do I specify it? > > The idea is that my script would read a text file as input, where individual > shapes (such as state/province boundaries) are associated with temperature > or some other variable. > > I am using GrADS v1.9.0-rc1 (installed from OpenGrADS), on a Windows XP > platform. > > Any help on this would be appreciated. > > Stephen McMillan > > > *************************************************** > 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. > *************************************************** > -- Arlindo da Silva dasilva at alum.mit.edu *************************************************** 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/20080724/942287d6/attachment.html From dasilva at ALUM.MIT.EDU Fri Jul 25 00:30:55 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Fri, 25 Jul 2008 00:30:55 -0400 Subject: Using Shape Files in GrADS In-Reply-To: Message-ID: On Thu, Jul 24, 2008 at 12:32 PM, Stephen R McMillan wrote: > > Arlindo, > I'm not yet familiar with how to build from sources, but I appreciate your > mentioning that option. > > Regarding the usage syntax, did you mean shp_polyf instead of shp_polyg? Yes, sorry for the typo. > Otherwise, what is "polyg" supposed to do? Also, in their present state (on > current version I'm using), both shp_polyf and shp_lines require either 1 or > 3 arguments. That was the bug. It should require 1 or 2. If I try two args as in your example, I get the GrADS error > "invalid number of arguments; need 1 or 3 but got 2". If this error is > valid, then the user will need to know the following for proper use: > > Arg1 Shapefile name such as "admin98" > Arg2 Presumably some object identifier, such as the "8" you > mentioned. Also, how to determine the range of acceptable object > identifiers. When you plot the whole things it prints this out. It would be good to add an option to query for this. I'll see if I can get a build this weekend. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From smcmillan at PLANALYTICS.COM Fri Jul 25 10:28:47 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Fri, 25 Jul 2008 09:28:47 -0500 Subject: Using Shape Files in GrADS In-Reply-To: <77fcd6b20807242130k51ba9d98i214bfa9fe521570a@mail.gmail.com> Message-ID: Thanks for the update--I look forward to the build. Stephen Arlindo da Silva Sent by: GRADSUSR at LIST.CINECA.IT 07/24/2008 11:30 PM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Re: Using Shape Files in GrADS On Thu, Jul 24, 2008 at 12:32 PM, Stephen R McMillan wrote: > > Arlindo, > I'm not yet familiar with how to build from sources, but I appreciate your > mentioning that option. > > Regarding the usage syntax, did you mean shp_polyf instead of shp_polyg? Yes, sorry for the typo. > Otherwise, what is "polyg" supposed to do? Also, in their present state (on > current version I'm using), both shp_polyf and shp_lines require either 1 or > 3 arguments. That was the bug. It should require 1 or 2. If I try two args as in your example, I get the GrADS error > "invalid number of arguments; need 1 or 3 but got 2". If this error is > valid, then the user will need to know the following for proper use: > > Arg1 Shapefile name such as "admin98" > Arg2 Presumably some object identifier, such as the "8" you > mentioned. Also, how to determine the range of acceptable object > identifiers. When you plot the whole things it prints this out. It would be good to add an option to query for this. I'll see if I can get a build this weekend. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu *************************************************** 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/20080725/b856e8dc/attachment.html From madan at EES.HOKUDAI.AC.JP Sat Jul 26 21:41:57 2008 From: madan at EES.HOKUDAI.AC.JP (madan sigdel) Date: Sat, 26 Jul 2008 18:41:57 -0700 Subject: about output data Message-ID: Dear Users I have calculated the out put result using fwrite and made separate binary file and read the output data by control file. But the results seem some strange. Is the way of doing correct? My file looks like this ********************* Open shum.mon.mean.ctl Open uwnd.mon.mean.ctl Open vwnd.mon.mean.ctl Set lev 1000 Set lat -90 90 Set lon 0 360 Set t 277 672 Define qu=u.2*sh Define qv=v.3*sh Define flx1=vint(lev,qu,500) Define flx2=vint(lev,qv,500) Set gxout fwrite Set fwrite data.bin d flx1 d flx2 disable fwrite ******************** Thank you .............................. Madan -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080726/1fe53d4b/attachment.html From bruyerec at UCAR.EDU Sat Jul 26 05:38:58 2008 From: bruyerec at UCAR.EDU (Cindy Bruyere) Date: Sat, 26 Jul 2008 03:38:58 -0600 Subject: about output data In-Reply-To: <000001c8ef89$f8089850$9ddb3285@hokkaido1s2x9u> Message-ID: I will be out of the office till August 4, 2008. From dasilva at ALUM.MIT.EDU Sun Jul 27 22:52:57 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Sun, 27 Jul 2008 22:52:57 -0400 Subject: help In-Reply-To: <720636.24344.qm@web31607.mail.mud.yahoo.com> Message-ID: On Mon, Jul 21, 2008 at 10:47 PM, wilfredo burgoa wrote: > Dear Arlindo, > > Help me, please. I have two data files with extension *.nc and would want > to carry out the difference between the archives soon to be read with > GrADS???. For example: uwnd20C.nc-uwndA2.nc. > I am not quite sure what you mean by "soon to be read with GrADS"... If each one of your files are GrADS readable, then you can compute the difference trivially with lats4d. From the OS command line: % lats4d.sh -nc -i uwnd20C.nc -j uwndA2.nc -func '@.1- at .2' -o uwnd_diff -v will create file uwnd_diff.nc with the desired difference. If the individual files are not GrADS readable you will need to address that first, either with a DDF (easier) or with a CTL (more flexible). Consult the User's Guide for details. BTW, please post questions like this to gradsusr and *not* directly to me. I am much more likely to answer a GrADS related question on gradsusr then on my personal inbox. Besides, you will have the benefit of many other knowledgeable users out there, and the answer may benefit other users as well. Thank you, Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From davide.sacchetti at ARPAL.ORG Mon Jul 28 02:19:27 2008 From: davide.sacchetti at ARPAL.ORG (Davide Sacchetti) Date: Mon, 28 Jul 2008 08:19:27 +0200 Subject: about output data In-Reply-To: <000001c8ef89$f8089850$9ddb3285@hokkaido1s2x9u> Message-ID: at a first look you define QU, QV only for level 1000 ... you can't integrate then from 1000 to 500. Use directly: Define flx1=vint(lev,u.2*sh.1,500) or set lev 1000 500 before definig QU, QV Bye bye Davide On Sat, 2008-07-26 at 18:41 -0700, madan sigdel wrote: > > > Dear Users > > > > I have calculated the out put result using fwrite and made separate > binary file and read the output data by control file. > > > > But the results seem some strange. Is the way of doing correct? > > > > My file looks like this > > ********************* > > Open shum.mon.mean.ctl > > Open uwnd.mon.mean.ctl > > Open vwnd.mon.mean.ctl > > > > Set lev 1000 > > Set lat -90 90 > > Set lon 0 360 > > > > Set t 277 672 > > > > Define qu=u.2*sh > > Define qv=v.3*sh > > > > Define flx1=vint(lev,qu,500) > > Define flx2=vint(lev,qv,500) > > > > Set gxout fwrite > > Set fwrite data.bin > > d flx1 > > d flx2 > > disable fwrite > > > > ******************** > > Thank you > > .............................. > > Madan > > -- 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 pablo.santos at NOAA.GOV Mon Jul 28 09:04:20 2008 From: pablo.santos at NOAA.GOV (Pablo Santos) Date: Mon, 28 Jul 2008 15:04:20 +0200 Subject: Grads and Grib2 g2ctl.pl Question Message-ID: Hello: I have a grb2 file (http://www.srh.noaa.gov/mfl/research/cwarf_test.grb2) that I create a ctl file for using g2ctl.pl. The ctl file looks like this: [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ cat cwarf_test.ctl dset ^cwarf_test.grb2 index ^cwarf_test.grb2.idx undef 9.999E+20 title cwarf_test.grb2 * produced by g2ctl v0.0.3h * griddef=1.1:4:(571 x 591):grid_template=30: Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 LoV 315.000000 LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP -90.000000 LonSP 0.000000 North Pole (571 dtype grib2 xdef 571 linear 0 0.1 ydef 591 linear 0 0.1 tdef 17 linear 00Z08jun2008 3hr zdef 1 linear 1 1 vars 3 DIRPWsfc 0,1,0 10,0,10 ** surface none Primary Wave Direction [deg] HTSGWsfc 0,1,0 10,0,3 ** surface none Significant Height of Combined Wind Waves and Swell [m] PERPWsfc 0,1,0 10,0,11 ** surface none Primary Wave Mean Period [s] ENDVARS But this seems to be a preprojected grib2 I am trying to display which means I need the pdef line added to the ctl file. Is there a way to make the g2ctl.pl utility do that? If not, could someone help me figure out for this particular grid what the pdef line is supposed to look like? Below is the syntax of what the pdef line I am supposed to add to the ctl file is supposed to look like. I am trying to figure out what the line is supposed to look like based on the output of the wgrib2 command below: [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ wgrib2 -grid cwarf_test.grb2 1.1:4:grid_template=30: Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 LoV 315.000000 LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP -90.000000 LonSP 0.000000 North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 Does anyone know? In any case, if someone can help it would be appreciated. Thanks. Pablo PDEF isize jsize LCCR latref lonref iref jref Struelat Ntruelat slon dx dy PDEF isize jsize LCC latref lonref iref jref Struelat Ntruelat slon dx dy Example: PDEF 103 69 lccr 30 -88 51.5 34.5 20 40 -88 90000 90000 Args: isize The size of the native grid in the x direction jsize The size of the native grid in the y direction latref reference latitude lonref reference longitude (in degrees, E is positive, W is negative) iref i of ref point jref j of ref point Struelat S true lat Ntruelat N true lat slon standard longitude dx grid X increment in meters dy grid Y increment in meters From Oliver.Fuhrer at METEOSWISS.CH Mon Jul 28 10:56:37 2008 From: Oliver.Fuhrer at METEOSWISS.CH (Oliver Fuhrer) Date: Mon, 28 Jul 2008 16:56:37 +0200 Subject: Compare diurnal cycles In-Reply-To: A<20080728130506.57AF92041F@mx2.cineca.it> Message-ID: Hi all! I have two variables in Grads containing the following diurnal cycles... flux1 (defined at 00Z, 03Z, 06Z, 09Z, 12Z, 15Z, 18Z, 21Z) flux2 (defined at 00Z, 01Z, 02Z, ..., 23Z) f1 at 00Z is the mean flux from 00Z to 03Z. f2 at 00Z is the mena flux from 00Z to 01Z. Now when I use Grads to plot the two on the same line chart, they are phase-shifted in time, since the flux1 at 00Z should be plotted at 01:30Z and the flux2 at 00Z should be plotted at 00:30Z, but they are not. How can I get Grads to time shift the two variables? Thanks for you help! Oliver > -----Original Message----- > From: GRADSUSR at LIST.CINECA.IT > [mailto:GRADSUSR at LIST.CINECA.IT] On Behalf Of Pablo Santos > Sent: Montag, 28. Juli 2008 15:04 > To: GRADSUSR at LIST.CINECA.IT > Subject: Grads and Grib2 g2ctl.pl Question > > Hello: > > I have a grb2 file > (http://www.srh.noaa.gov/mfl/research/cwarf_test.grb2) > that I create a ctl file for using g2ctl.pl. > > The ctl file looks like this: > > [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ cat > cwarf_test.ctl > dset ^cwarf_test.grb2 > index ^cwarf_test.grb2.idx > undef 9.999E+20 > title cwarf_test.grb2 > * produced by g2ctl v0.0.3h > * griddef=1.1:4:(571 x 591):grid_template=30: Lambert > Conformal: (571 x 591) > input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 > LoV 315.000000 > LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP > -90.000000 LonSP > 0.000000 North Pole (571 > dtype grib2 > xdef 571 linear 0 0.1 > ydef 591 linear 0 0.1 > tdef 17 linear 00Z08jun2008 3hr > zdef 1 linear 1 1 > vars 3 > DIRPWsfc 0,1,0 10,0,10 ** surface none Primary Wave Direction [deg] > HTSGWsfc 0,1,0 10,0,3 ** surface none Significant Height > of Combined Wind > Waves and Swell [m] > PERPWsfc 0,1,0 10,0,11 ** surface none Primary Wave Mean Period [s] > ENDVARS > > But this seems to be a preprojected grib2 I am trying to > display which means > I need the pdef line added to the ctl file. Is there a way to make the > g2ctl.pl utility do that? If not, could someone help me > figure out for this > particular grid what the pdef line is supposed to look like? > > Below is the syntax of what the pdef line I am supposed to > add to the ctl > file is supposed to look like. I am trying to figure out what > the line is > supposed to look like based on the output of the wgrib2 command below: > > [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ wgrib2 -grid > cwarf_test.grb2 > 1.1:4:grid_template=30: > Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 > Lat1 28.201020 Lon1 279.309660 LoV 315.000000 > LatD 28.500000 Latin1 30.000000 Latin2 30.000000 > LatSP -90.000000 LonSP 0.000000 > North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 > > Does anyone know? In any case, if someone can help it would > be appreciated. > Thanks. > > Pablo > > PDEF isize jsize LCCR latref lonref iref jref Struelat > Ntruelat slon dx dy > PDEF isize jsize LCC latref lonref iref jref Struelat > Ntruelat slon dx dy > Example: PDEF 103 69 lccr 30 -88 51.5 34.5 20 40 -88 90000 90000 > Args: > isize The size of the native grid in the x direction > jsize The size of the native grid in the y direction > latref reference latitude > lonref reference longitude (in degrees, E is positive, W is negative) > iref i of ref point > jref j of ref point > Struelat S true lat > Ntruelat N true lat > slon standard longitude > dx grid X increment in meters > dy grid Y increment in meters > From hersala at GMAIL.COM Mon Jul 28 11:09:01 2008 From: hersala at GMAIL.COM (hersala) Date: Mon, 28 Jul 2008 12:09:01 -0300 Subject: a consult In-Reply-To: <5b40c90807280739g42e3d132t310d7631c4908258@mail.gmail.com> Message-ID: Hello Grads users: Could you help me with this doubt: is it possible that some Grads versions do not recognize the instruction '!sleep 1' in a .gs file ? Because the same .gs runs properly in one computer but not it does not in another, in the last one appears sleep: not known. Thank you ! Hern?n Sala Instituto Ant?rtico Argentino -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080728/db16f725/attachment.html From ela at COLA.IGES.ORG Mon Jul 28 12:14:26 2008 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Mon, 28 Jul 2008 12:14:26 -0400 Subject: a consult In-Reply-To: <125951949.106951217261561671.JavaMail.root@mail.iges.org> Message-ID: Hernan, The sleep command is a unix command. Using ! in a grads script tells grads to have the shell interpret whatever follows the ! symbol. So your problem is not related to grads. Check your PATH environment variable and see if /bin and /usr/bin are in your PATH. Most common unix commands are found in /bin or /usr/bin. Best regards, Eric L. Altshuler Assistant Research Scientist Center for Ocean-Land-Atmosphere Studies 4041 Powder Mill Road, Suite 302 Calverton, MD 20705-3106 USA E-mail: ela at cola.iges.org Phone: (301) 902-1257 Fax: (301) 595-9793 ----- Original Message ----- From: "hersala" To: GRADSUSR at LIST.CINECA.IT Sent: Monday, July 28, 2008 11:09:01 AM GMT -05:00 US/Canada Eastern Subject: a consult Hello Grads users: Could you help me with this doubt: is it possible that some Grads versions do not recognize the instruction '!sleep 1' in a .gs file ? Because the same .gs runs properly in one computer but not it does not in another, in the last one appears sleep: not known. Thank you ! Hern?n Sala Instituto Ant?rtico Argentino From ecarter at WEATHEREXTREME.COM Mon Jul 28 14:51:27 2008 From: ecarter at WEATHEREXTREME.COM (Elizabeth Austin) Date: Mon, 28 Jul 2008 11:51:27 -0700 Subject: viewing NOAA bin files - question Message-ID: Hello Grads Users, I ordered some wind/wave data from the NOAA site: http://hurricane.ncdc.noaa.gov/pls/plhas/HAS.FileAppSelect?datasetname=9957GRA The files are .bin files and I am not sure how to view them in grads. grib2ctl.pl does not work. Do I need to use degrib or ? I am attaching one of the files. Also, on this same site, some of the data were .pcx files (I know what these are but am unable to view them in ANY program, e.g., Quark, Photoshop, Paint, picture viewer, illustrator, etc.). I am also attaching one of these just on the off-chance that someone knows how to view these files. Thanks very much. Elizabeth Austin ____________________________ Elizabeth J. Austin, Ph.D., C.C.M. President WeatherExtreme Ltd. Lake Tahoe, Nevada USA Phone: 775-636-8553 Fax: 775-636-8430 E-fax: 775-562-9333 -------------- next part -------------- A non-text attachment was scrubbed... Name: ds.waveh.bin Type: application/macbinary Size: 24314 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080728/44588533/attachment.bin -------------- next part -------------- A non-text attachment was scrubbed... Name: 2007081802_SRRS0136.PCX Type: application/octet-stream Size: 283771 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080728/44588533/attachment.obj From Wesley.Ebisuzaki at NOAA.GOV Mon Jul 28 14:51:08 2008 From: Wesley.Ebisuzaki at NOAA.GOV (Wesley Ebisuzaki) Date: Mon, 28 Jul 2008 14:51:08 -0400 Subject: Grads and Grib2 g2ctl.pl Question In-Reply-To: <20080728130506.57AF92041F@mx2.cineca.it> Message-ID: Pablo, I couldn't display the file because the "dx" and "dy" parameters were not given. 1.1:4:grid_template=30: Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 LoV 315.000000 LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP -90.000000 LonSP 0.000000 North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 This is common problem with NCEP grib2 files. The dx and dy values are actually encoded but a flag indicates that the dx/dy values are missing. You can change the flag value by wgrib2 -set_flag_table_3.3 48 in.grb -grib out.grb Once I did this, I could display the fields. You may need to upgrade wgrib2 to get the set_flag_table_3.3 option. Wesley Ebisuzaki Pablo Santos wrote: > Hello: > > I have a grb2 file (http://www.srh.noaa.gov/mfl/research/cwarf_test.grb2) > that I create a ctl file for using g2ctl.pl. > > The ctl file looks like this: > > [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ cat cwarf_test.ctl > dset ^cwarf_test.grb2 > index ^cwarf_test.grb2.idx > undef 9.999E+20 > title cwarf_test.grb2 > * produced by g2ctl v0.0.3h > * griddef=1.1:4:(571 x 591):grid_template=30: Lambert Conformal: (571 x 591) > input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 LoV 315.000000 > LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP -90.000000 LonSP > 0.000000 North Pole (571 > dtype grib2 > xdef 571 linear 0 0.1 > ydef 591 linear 0 0.1 > tdef 17 linear 00Z08jun2008 3hr > zdef 1 linear 1 1 > vars 3 > DIRPWsfc 0,1,0 10,0,10 ** surface none Primary Wave Direction [deg] > HTSGWsfc 0,1,0 10,0,3 ** surface none Significant Height of Combined Wind > Waves and Swell [m] > PERPWsfc 0,1,0 10,0,11 ** surface none Primary Wave Mean Period [s] > ENDVARS > > But this seems to be a preprojected grib2 I am trying to display which means > I need the pdef line added to the ctl file. Is there a way to make the > g2ctl.pl utility do that? If not, could someone help me figure out for this > particular grid what the pdef line is supposed to look like? > > Below is the syntax of what the pdef line I am supposed to add to the ctl > file is supposed to look like. I am trying to figure out what the line is > supposed to look like based on the output of the wgrib2 command below: > > [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ wgrib2 -grid > cwarf_test.grb2 > 1.1:4:grid_template=30: > Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 > Lat1 28.201020 Lon1 279.309660 LoV 315.000000 > LatD 28.500000 Latin1 30.000000 Latin2 30.000000 > LatSP -90.000000 LonSP 0.000000 > North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 > > Does anyone know? In any case, if someone can help it would be appreciated. > Thanks. > > Pablo > > PDEF isize jsize LCCR latref lonref iref jref Struelat Ntruelat slon dx dy > PDEF isize jsize LCC latref lonref iref jref Struelat Ntruelat slon dx dy > Example: PDEF 103 69 lccr 30 -88 51.5 34.5 20 40 -88 90000 90000 > Args: > isize The size of the native grid in the x direction > jsize The size of the native grid in the y direction > latref reference latitude > lonref reference longitude (in degrees, E is positive, W is negative) > iref i of ref point > jref j of ref point > Struelat S true lat > Ntruelat N true lat > slon standard longitude > dx grid X increment in meters > dy grid Y increment in meters > From ameliagraziano at TISCALINET.IT Mon Jul 28 15:47:19 2008 From: ameliagraziano at TISCALINET.IT (amelia graziano) Date: Mon, 28 Jul 2008 21:47:19 +0200 Subject: lats4d2.sh - convert GRIB2 ==> GRIB1 In-Reply-To: <488E151C.7050305@noaa.gov> Message-ID: Hello, I'm trying to convert GRIB2 ==> GRIB1 using lats4d2.sh. Unfortunately I cannot find the correct script to do it. I've a GRIB2 file gfs.t00z.pgrb2f03 and I want to convert to GRIB1 I would apprecciate some help Thanks in advance Amelia From pablo.santos at NOAA.GOV Mon Jul 28 17:40:28 2008 From: pablo.santos at NOAA.GOV (Pablo Santos) Date: Mon, 28 Jul 2008 17:40:28 -0400 Subject: Grads and Grib2 g2ctl.pl Question In-Reply-To: <488E151C.7050305@noaa.gov> Message-ID: Wesley: Got it. Creating the grb file with the option below and running it again through g2ctl did the trick creating the pdef line in the ctl file. I was able to do it with the version of grads (2) and wgrib2 I have in linux. However, I also often work with grads in windows/cygwin. I downloaded the latest version of wgrib2 for windows/cygwin I am using from ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/Windows_XP but it does not have the flag you are referring to below. Also, does anybody know if there is a compiled version of grads2 with support for grib2 for windows in the work? Thanks, Pablo Wesley Ebisuzaki wrote: > Pablo, > > I couldn't display the file because the "dx" and "dy" parameters were > not given. > > 1.1:4:grid_template=30: > Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 > Lat1 28.201020 Lon1 279.309660 LoV 315.000000 > LatD 28.500000 Latin1 30.000000 Latin2 30.000000 > LatSP -90.000000 LonSP 0.000000 > North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 > > This is common problem with NCEP grib2 files. The dx and dy values are > actually > encoded but a flag indicates that the dx/dy values are missing. You can > change the > flag value by > > wgrib2 -set_flag_table_3.3 48 in.grb -grib out.grb > > Once I did this, I could display the fields. You may need to upgrade > wgrib2 to get the > set_flag_table_3.3 option. > > Wesley Ebisuzaki > > > > Pablo Santos wrote: >> Hello: >> >> I have a grb2 file >> (http://www.srh.noaa.gov/mfl/research/cwarf_test.grb2) >> that I create a ctl file for using g2ctl.pl. >> >> The ctl file looks like this: >> >> [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ cat cwarf_test.ctl >> dset ^cwarf_test.grb2 >> index ^cwarf_test.grb2.idx >> undef 9.999E+20 >> title cwarf_test.grb2 >> * produced by g2ctl v0.0.3h >> * griddef=1.1:4:(571 x 591):grid_template=30: Lambert Conformal: (571 >> x 591) >> input WE:SN output WE:SN res 0 Lat1 28.201020 Lon1 279.309660 LoV >> 315.000000 >> LatD 28.500000 Latin1 30.000000 Latin2 30.000000 LatSP -90.000000 LonSP >> 0.000000 North Pole (571 >> dtype grib2 >> xdef 571 linear 0 0.1 >> ydef 591 linear 0 0.1 >> tdef 17 linear 00Z08jun2008 3hr >> zdef 1 linear 1 1 >> vars 3 >> DIRPWsfc 0,1,0 10,0,10 ** surface none Primary Wave Direction [deg] >> HTSGWsfc 0,1,0 10,0,3 ** surface none Significant Height of >> Combined Wind >> Waves and Swell [m] >> PERPWsfc 0,1,0 10,0,11 ** surface none Primary Wave Mean Period [s] >> ENDVARS >> >> But this seems to be a preprojected grib2 I am trying to display >> which means >> I need the pdef line added to the ctl file. Is there a way to make the >> g2ctl.pl utility do that? If not, could someone help me figure out >> for this >> particular grid what the pdef line is supposed to look like? >> >> Below is the syntax of what the pdef line I am supposed to add to the >> ctl >> file is supposed to look like. I am trying to figure out what the >> line is >> supposed to look like based on the output of the wgrib2 command below: >> >> [nawips at MFL-LW-WRFPOST ~/FIT_CSTAR_Test_data_files]$ wgrib2 -grid >> cwarf_test.grb2 >> 1.1:4:grid_template=30: >> Lambert Conformal: (571 x 591) input WE:SN output WE:SN res 0 >> Lat1 28.201020 Lon1 279.309660 LoV 315.000000 >> LatD 28.500000 Latin1 30.000000 Latin2 30.000000 >> LatSP -90.000000 LonSP 0.000000 >> North Pole (571 x 591) Dx 0.000000 m Dy 0.000000 m mode 0 >> >> Does anyone know? In any case, if someone can help it would be >> appreciated. >> Thanks. >> >> Pablo >> >> PDEF isize jsize LCCR latref lonref iref jref Struelat Ntruelat slon >> dx dy >> PDEF isize jsize LCC latref lonref iref jref Struelat Ntruelat slon >> dx dy >> Example: PDEF 103 69 lccr 30 -88 51.5 34.5 20 40 -88 90000 90000 >> Args: >> isize The size of the native grid in the x direction >> jsize The size of the native grid in the y direction >> latref reference latitude >> lonref reference longitude (in degrees, E is positive, W is negative) >> iref i of ref point >> jref j of ref point >> Struelat S true lat >> Ntruelat N true lat >> slon standard longitude >> dx grid X increment in meters >> dy grid Y increment in meters >> From dasilva at ALUM.MIT.EDU Mon Jul 28 20:13:43 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 28 Jul 2008 20:13:43 -0400 Subject: Grads and Grib2 g2ctl.pl Question In-Reply-To: <488E3CCC.2040601@noaa.gov> Message-ID: On Mon, Jul 28, 2008 at 5:40 PM, Pablo Santos wrote: > Wesley: > > Got it. Creating the grb file with the option below and running it again > through g2ctl did the trick creating the pdef line in the ctl file. I > was able to do it with the version of grads (2) and wgrib2 I have in linux. > > However, I also often work with grads in windows/cygwin. I downloaded > the latest version of wgrib2 for windows/cygwin I am using from > > ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2/Windows_XP > > but it does not have the flag you are referring to below. > > Also, does anybody know if there is a compiled version of grads2 with > support for grib2 for windows in the work? > The next release of Win32 GrADS due sometime this summer will very likely include grads2/gradsdap binaries. In the meantime, search the list as someone has contributed a grads 2.0.a2 binaries. Also, Win32 GrADS v1.9.0-rc1 comes with an extension that allows you to read grib2. At this point you are better off using it with perl from cygwin. Info on Win32 GrADS is here: http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows The g2 extension is documented here: http://opengrads.org/doc/udxt/grib2/ and comes bundled with Win32 GrADS. Arlindo -- Arlindo da Silva dasilva at alum.mit.edu From dasilva at ALUM.MIT.EDU Mon Jul 28 20:18:11 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Mon, 28 Jul 2008 20:18:11 -0400 Subject: lats4d2.sh - convert GRIB2 ==> GRIB1 In-Reply-To: <200807281947.m6SJlQ8Z017458@smtpi1.ngi.it> Message-ID: On Mon, Jul 28, 2008 at 3:47 PM, amelia graziano wrote: > Hello, > > I'm trying to convert GRIB2 ==> GRIB1 using lats4d2.sh. Unfortunately > I cannot find the correct script to do it. > > I've a GRIB2 file gfs.t00z.pgrb2f03 and I want to convert to GRIB1 > > I would apprecciate some help > 1) Search this list for my recent posting with las4d.sh/lats4d.gs attachment 2) Make sure you have the g2() extension installed: http://opengrads.org/doc/udxt/grib2/ http://opengrads.org/wiki/index.php?title=User_Defined_Extensions 3) You must have grads 1.9.0-rc1available from http://sourceforge.net/project/showfiles.php?group_id=161773 Arlindo > Thanks in advance > > Amelia > -- Arlindo da Silva dasilva at alum.mit.edu From anwibuma at YAHOO.COM Mon Jul 28 21:45:40 2008 From: anwibuma at YAHOO.COM (wilfredo burgoa) Date: Mon, 28 Jul 2008 18:45:40 -0700 Subject: ttest In-Reply-To: <77fcd6b20807281718i466101e7qf55b16c916453d5e@mail.gmail.com> Message-ID: Dear all users, ? Please. Statistical significance scripts calculate the t-test statistic for a two-sample difference of means test and shade the grid cells that are significant at a chosen level. I have two data file with extensions *.nc. When run scripts in GrADS the error message is?undefined grid contours. My firts conclusion is that them point of grid of the archives is not identical, in such sense, That I can do????????????????? Yours Faithfully, Andreas ? ____________________________________________________________________________________ Yahoo! MTV Blog & Rock >?Cu?ntanos tu historia, inspira una canci?n y g?nate un viaje a los Premios MTV! Participa aqu? http://mtvla.yahoo.com/ -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080728/dffd3b6c/attachment.html From anwibuma at YAHOO.COM Mon Jul 28 21:57:33 2008 From: anwibuma at YAHOO.COM (wilfredo burgoa) Date: Mon, 28 Jul 2008 18:57:33 -0700 Subject: ttest Message-ID: Dear all users, ? Please. Statistical significance scripts calculate the t-test statistic for a two-sample difference of means test and shade the grid cells that are significant at a chosen level. I have two data file with extensions *.nc. When run scripts in GrADS the error message is?undefined grid contours. My firts conclusion is that them point of grid of the archives is not identical, in such sense, That I can do????????????????? Yours Faithfully, Andreas ? ____________________________________________________________________________________ Yahoo! MTV Blog & Rock >?Cu?ntanos tu historia, inspira una canci?n y g?nate un viaje a los Premios MTV! Participa aqu? http://mtvla.yahoo.com/ -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080728/fe5edaa5/attachment.html From mike.bosilovich at GMAIL.COM Mon Jul 28 22:10:26 2008 From: mike.bosilovich at GMAIL.COM (Mike Bosilovich) Date: Mon, 28 Jul 2008 22:10:26 -0400 Subject: Compare diurnal cycles In-Reply-To: <19E843D905DC46438D74D64CAF17340905C7FE32@zueex01.meteoswiss.ch> Message-ID: Hi Oliver, It should be possible, but you didn't say what format the data is in. If it is flat binary with a ctl file, make a new ctl with your center point time instead of the report time (in other words, in the TDEF line, change 00Z to 01:30Z, and in the other 00Z to 00:30Z). It all should work, unless you have the hour templated in the file name. If that is the case, or if you have grib or a self defined format, it may not be possible to center the average directly. I haven't tried changing the time in an DDF file, so I'm not sure that that would work. One way would be to write a binary with lats4d (lats4d -i infile -o outfile -format stream), then create the ctl with the center time point. The lats4d option may not work for single point time series in a self defined format. It typically hasn't liked -lat and -lon being points. I suggest creating a new ctl so that in the future, after you haven't used the data in a while, it will be apparent in your scripts which uses the centered average, versus the reported time. Good Luck, Mike On Mon, Jul 28, 2008 at 10:56 AM, Oliver Fuhrer wrote: > > Hi all! > > I have two variables in Grads containing the following diurnal cycles... > > flux1 (defined at 00Z, 03Z, 06Z, 09Z, 12Z, 15Z, 18Z, 21Z) > flux2 (defined at 00Z, 01Z, 02Z, ..., 23Z) > > f1 at 00Z is the mean flux from 00Z to 03Z. > f2 at 00Z is the mena flux from 00Z to 01Z. > > Now when I use Grads to plot the two on the same line chart, they are > phase-shifted in time, since the flux1 at 00Z should be plotted at > 01:30Z and the flux2 at 00Z should be plotted at 00:30Z, but they are > not. > > How can I get Grads to time shift the two variables? > > Thanks for you help! > Oliver > > From zaheer_a_babar at HOTMAIL.COM Tue Jul 29 07:15:35 2008 From: zaheer_a_babar at HOTMAIL.COM (Zaheer Babar) Date: Tue, 29 Jul 2008 11:15:35 +0000 Subject: Problem of Grads on Windows Vista Message-ID: Hi! I have newly installed vista, but Im unable to run Grads on it. Vista is installed but when I run grads dods a window pops up saying "windows is trying to solve the problem" and then it fails. Any suggestions Zaheer _________________________________________________________________ Keep your kids safer online with Windows Live Family Safety. http://www.windowslive.com/family_safety/overview.html?ocid=TXT_TAGLM_WL_family_safety_072008 -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080729/90ea6ea1/attachment.html From zaheer_a_babar at HOTMAIL.COM Tue Jul 29 07:15:51 2008 From: zaheer_a_babar at HOTMAIL.COM (Zaheer Babar) Date: Tue, 29 Jul 2008 11:15:51 +0000 Subject: Problem of Grads on Windows Vista Message-ID: Hi! I have newly installed vista, but Im unable to run Grads on it. Vista is installed but when I run grads dods a window pops up saying "windows is trying to solve the problem" and then it fails. Any suggestions Zaheer _________________________________________________________________ Time for vacation? WIN what you need- enter now! http://www.gowindowslive.com/summergiveaway/?ocid=tag_jlyhm -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080729/50662c6e/attachment.html From qdai at MATH.UALBERTA.CA Tue Jul 29 14:26:36 2008 From: qdai at MATH.UALBERTA.CA (Qingfang Dai) Date: Tue, 29 Jul 2008 12:26:36 -0600 Subject: How to plot the contour shaded with a color bar? Message-ID: Dear all users, I want to plot EOF data into contour shaded with a color bar. Now my GS file is as follows: 'reinit' 'open ../ABMapsForReport/EOF_Jan.ctl' 'set grads off' 'set mproj nps' 'set grid off' 'set frame off' 'set gxout shaded' 'set clevs -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 1' 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6' 'set cstyle 1' 'set csmooth off' 'set t 1 1' 'set mpdset ../province' 'd prep' 'set font 1' 'set string 1 tc 3' 'set strsiz 0.15 0.15' 'draw title EOF1 1961-2000 January' 'cbarn' 'printim ../EOF1_Jan.gif gif white' But I don't know why the color bar I got without data. I am wondering how to plot a coutour shaded with a color bar which display each color denotes each data interval. Thanks. Yours Faithfully, Fang From mmacleod at SCOTIAWEATHER.COM Tue Jul 29 14:38:47 2008 From: mmacleod at SCOTIAWEATHER.COM (mmacleod) Date: Tue, 29 Jul 2008 18:38:47 +0000 Subject: How to plot the contour shaded with a color bar? In-Reply-To: <20080729181645.M31456@math.ualberta.ca> Message-ID: Hi Fang If you can send a copy of the GIF file generated so I can see what is produced I might be able to help. We do this type of output all the time. I am attaching a copy of an output (it is in PNG but the same thing appears in GIF if I specify that format) to see if this is the type of thing you are looking for. Regards, Mac Qingfang Dai wrote: > Dear all users, > > I want to plot EOF data into contour shaded with a color bar. Now my GS file > is as follows: > > 'reinit' > 'open ../ABMapsForReport/EOF_Jan.ctl' > 'set grads off' > 'set mproj nps' > 'set grid off' > 'set frame off' > 'set gxout shaded' > 'set clevs -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 1' > 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6' > 'set cstyle 1' > 'set csmooth off' > 'set t 1 1' > 'set mpdset ../province' > 'd prep' > 'set font 1' > 'set string 1 tc 3' > 'set strsiz 0.15 0.15' > 'draw title EOF1 1961-2000 January' > 'cbarn' > 'printim ../EOF1_Jan.gif gif white' > > But I don't know why the color bar I got without data. I am wondering how to > plot a coutour shaded with a color bar which display each color denotes each > data interval. Thanks. > > > Yours Faithfully, > > Fang > -- M.A. (Mac) MacLeod President and General Manager Scotia Weather Services Inc 192 Wyse Road, Suite 8, Dartmouth, N.S. B3A 1M9 Tele: 902-468-3866 Fax: 902-461-1768 E-mail: mmacleod at scotiaweather.com Visit us: www.scotiaweather.com -------------- next part -------------- A non-text attachment was scrubbed... Name: GEMRgn_Ern_NA_12HrTotal-PcpnH48.png Type: image/png Size: 44992 bytes Desc: not available Url : http://gradsusr.org/pipermail/gradsusr/attachments/20080729/a7c42c87/attachment.png From eduardo.agosta at GMAIL.COM Tue Jul 29 14:38:35 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Tue, 29 Jul 2008 15:38:35 -0300 Subject: How to plot the contour shaded with a color bar? In-Reply-To: <20080729181645.M31456@math.ualberta.ca> Message-ID: Fang, Verify whether you data are whithin the interval range you proposed (usually if your data interval is out of the settings, one-coloured shade is drawn).....or try to put the cbar after the display command...or use it as run c:\...path...\cbarn.gs I hope it helps. Eduardo. 2008/7/29 Qingfang Dai > Dear all users, > > I want to plot EOF data into contour shaded with a color bar. Now my GS > file > is as follows: > > 'reinit' > 'open ../ABMapsForReport/EOF_Jan.ctl' > 'set grads off' > 'set mproj nps' > 'set grid off' > 'set frame off' > 'set gxout shaded' > 'set clevs -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 1' > 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6' > 'set cstyle 1' > 'set csmooth off' > 'set t 1 1' > 'set mpdset ../province' > 'd prep' > 'set font 1' > 'set string 1 tc 3' > 'set strsiz 0.15 0.15' > 'draw title EOF1 1961-2000 January' > 'cbarn' > 'printim ../EOF1_Jan.gif gif white' > > But I don't know why the color bar I got without data. I am wondering how > to > plot a coutour shaded with a color bar which display each color denotes > each > data interval. Thanks. > > > Yours Faithfully, > > Fang > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080729/77d23821/attachment.html From eduardo.agosta at GMAIL.COM Tue Jul 29 14:45:12 2008 From: eduardo.agosta at GMAIL.COM (Eduardo Agosta Scarel) Date: Tue, 29 Jul 2008 15:45:12 -0300 Subject: ttest In-Reply-To: <352360.42790.qm@web31604.mail.mud.yahoo.com> Message-ID: Andreas, If the gridded resolution between the data files are different, few than less is left to do. Try to specify your request showing the script you are using...perhaps the problem is other thing. Good luck, Eduardo. 2008/7/28 wilfredo burgoa > Dear all users, > > > > Please. Statistical significance scripts calculate the *t-test statistic*for a > *two-sample difference of means test *and shade the grid cells that are > significant at a chosen level. I have two data file with extensions *.nc. > When run scripts in GrADS the error message is undefined grid contours. My > firts conclusion is that them point of grid of the archives is not > identical, in such sense, That I can do????????????????? > > Yours Faithfully, > > Andreas > > > > ------------------------------ > > Yahoo! MTV Blog & Rock > ?Cu?ntanos tu historia, inspira una canci?n y g?nate un viaje a los Premios > MTV! > Participa aqu? http://mtvla.yahoo.com/ > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080729/550a0ece/attachment.html From mmacleod at SCOTIAWEATHER.COM Tue Jul 29 15:14:41 2008 From: mmacleod at SCOTIAWEATHER.COM (mmacleod) Date: Tue, 29 Jul 2008 19:14:41 +0000 Subject: How to plot the contour shaded with a color bar? In-Reply-To: <20080729184455.M80802@math.ualberta.ca> Message-ID: Hi Fang If you want the contour values as I had in the PNG file I sent then after 'd prep' and before 'set font 1' add 'set gxout contour' 'set clevs -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 1' 'set ccols 1 1 1 1 1 1 1 1 1 1 1 1' 'set cstyle 2 ' 'set cthick 4 ' 'd prep' This will draw contours in black as the selected level then number the contour in black as well. If you want to use the a different colour then change the 1 (black to the number of the colour you want) Note: if you use the same colour for the contours as the shading then you will not see the contours nor the value of the contours displayed.) Hope this helps. Mac Qingfang Dai wrote: > Hi Mac, > > The attachment is the GIF file I got. Thanks. > > Fang > > > On Tue, 29 Jul 2008 18:38:47 +0000, mmacleod wrote > >> Hi Fang >> >> If you can send a copy of the GIF file generated so I can see what is >> produced I might be able to help. We do this type of output all the >> time. I am attaching a copy of an output (it is in PNG but the same >> thing appears in GIF if I specify that format) to see if this is the >> type of thing you are looking for. >> >> Regards, >> >> Mac >> >> Qingfang Dai wrote: >> >>> Dear all users, >>> >>> I want to plot EOF data into contour shaded with a color bar. Now my GS >>> > file > >>> is as follows: >>> >>> 'reinit' >>> 'open ../ABMapsForReport/EOF_Jan.ctl' >>> 'set grads off' >>> 'set mproj nps' >>> 'set grid off' >>> 'set frame off' >>> 'set gxout shaded' >>> 'set clevs -20 -18 -16 -14 -12 -10 -8 -6 -4 -2 0 1' >>> 'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6' >>> 'set cstyle 1' >>> 'set csmooth off' >>> 'set t 1 1' >>> 'set mpdset ../province' >>> 'd prep' >>> 'set font 1' >>> 'set string 1 tc 3' >>> 'set strsiz 0.15 0.15' >>> 'draw title EOF1 1961-2000 January' >>> 'cbarn' >>> 'printim ../EOF1_Jan.gif gif white' >>> >>> But I don't know why the color bar I got without data. I am wondering >>> > how to > >>> plot a coutour shaded with a color bar which display each color denotes >>> > each > >>> data interval. Thanks. >>> >>> >>> Yours Faithfully, >>> >>> Fang >>> >>> >> -- >> M.A. (Mac) MacLeod >> President and General Manager >> Scotia Weather Services Inc >> 192 Wyse Road, Suite 8, >> Dartmouth, N.S. B3A 1M9 >> >> Tele: 902-468-3866 >> Fax: 902-461-1768 >> E-mail: mmacleod at scotiaweather.com >> Visit us: www.scotiaweather.com >> > > > -- > Department of Math and Stat Sciences (http://www.math.ualberta.ca) > > > > ------------------------------------------------------------------------ > -- M.A. (Mac) MacLeod President and General Manager Scotia Weather Services Inc 192 Wyse Road, Suite 8, Dartmouth, N.S. B3A 1M9 Tele: 902-468-3866 Fax: 902-461-1768 E-mail: mmacleod at scotiaweather.com Visit us: www.scotiaweather.com From dasilva at ALUM.MIT.EDU Tue Jul 29 15:58:15 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Tue, 29 Jul 2008 15:58:15 -0400 Subject: Problem of Grads on Windows Vista In-Reply-To: Message-ID: On Tue, Jul 29, 2008 at 7:15 AM, Zaheer Babar wrote: > Hi! > > I have newly installed vista, but Im unable to run Grads on it. Vista is > installed but when I run grads dods a window pops up saying "windows is > trying to solve the problem" and then it fails. > Any suggestions > Which version of GrADS are you using? Try v1.9.0-rc1: http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows -- Arlindo da Silva dasilva at alum.mit.edu From brian.walsh at AMEC.COM Wed Jul 30 10:25:37 2008 From: brian.walsh at AMEC.COM (Walsh, Brian J) Date: Wed, 30 Jul 2008 10:25:37 -0400 Subject: Topographical Maps Message-ID: Hello, I'm looking to plot topographical maps (resolution of 1 km or less) using GrADS. Can anyone direct me on where to obtain this data (specifically in gridded format) for North America? Thanks, Brian The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. Its contents (including any attachments) may contain confidential and/or privileged information. If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080730/a1bf8264/attachment.html From teddyallen at YAHOO.COM Wed Jul 30 10:50:40 2008 From: teddyallen at YAHOO.COM (Teddy Allen) Date: Wed, 30 Jul 2008 07:50:40 -0700 Subject: Topographical Maps In-Reply-To: Message-ID: Brian, Maybe this may help you: http://jisao.washington.edu/data/elevation/ teddy To move quickly, go alone. To move far, go together. --- African Proverb --- On Wed, 7/30/08, Walsh, Brian J wrote: > From: Walsh, Brian J > Subject: Topographical Maps > To: GRADSUSR at LIST.CINECA.IT > Date: Wednesday, July 30, 2008, 7:25 AM > Hello, > > I'm looking to plot topographical maps (resolution of 1 > km or less) > using GrADS. Can anyone direct me on where to obtain this > data > (specifically in gridded format) for North America? > > Thanks, > > Brian > > > The information contained in this e-mail is intended only > for the individual or entity to whom it is addressed. > Its contents (including any attachments) may contain > confidential and/or privileged information. > If you are not an intended recipient you must not use, > disclose, disseminate, copy or print its contents. > If you receive this e-mail in error, please notify the > sender by reply e-mail and delete and destroy the message. From smcmillan at PLANALYTICS.COM Wed Jul 30 11:11:58 2008 From: smcmillan at PLANALYTICS.COM (Stephen R McMillan) Date: Wed, 30 Jul 2008 10:11:58 -0500 Subject: Topographical Maps In-Reply-To: Message-ID: Brian, You might also try the following: http://monsoondata.org:9090/dods/topo Stephen "Walsh, Brian J" Sent by: GRADSUSR at LIST.CINECA.IT 07/30/2008 09:25 AM Please respond to GRADSUSR at LIST.CINECA.IT To GRADSUSR at LIST.CINECA.IT cc Subject Topographical Maps Hello, I'm looking to plot topographical maps (resolution of 1 km or less) using GrADS. Can anyone direct me on where to obtain this data (specifically in gridded format) for North America? Thanks, Brian *************************************************** 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/20080730/7acc1cba/attachment.html From john.bacon at WORSTEAD.CO.UK Wed Jul 30 11:15:24 2008 From: john.bacon at WORSTEAD.CO.UK (John_Bacon) Date: Wed, 30 Jul 2008 16:15:24 +0100 Subject: File close error Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080730/3b25cff8/attachment.html From michael.l.sestak at NAVY.MIL Wed Jul 30 11:32:22 2008 From: michael.l.sestak at NAVY.MIL (Sestak, Michael L CIV 63134) Date: Wed, 30 Jul 2008 08:32:22 -0700 Subject: File close error In-Reply-To: <4890858C.8080802@worstead.co.uk> Message-ID: An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080730/4fa534ba/attachment.html From rafeh at DREAMSCAPEIT.COM Wed Jul 30 11:46:18 2008 From: rafeh at DREAMSCAPEIT.COM (Rafe Hargrove) Date: Wed, 30 Jul 2008 10:46:18 -0500 Subject: UNSUBSCRIBE In-Reply-To: <000901c8b653$9129b6b0$91b1a8c0@arpa.veneto.it> Message-ID: UNSUBSCRIBE This email and any attached files are confidential and intended solely for the intended recipient(s). If you are not the named recipient you should not read, distribute, copy or alter this email. Any views or opinions expressed in this email are those of the author and do not represent those of the company. Warning: Although precautions have been taken to make sure no viruses are present in this email, the company cannot accept responsibility for any loss or damage that arise from the use of this email or attachments. -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080730/f3f9a727/attachment.html From Charles.Seman at NOAA.GOV Wed Jul 30 11:40:17 2008 From: Charles.Seman at NOAA.GOV (Charles Seman) Date: Wed, 30 Jul 2008 11:40:17 -0400 Subject: File close error In-Reply-To: <4890858C.8080802@worstead.co.uk> Message-ID: John, I've made the mistake of trying to close a file # other than the last one opened... so the files must be closed in reverse numerical order... for example, if files 1, 2, 3 were opened, then they must be closed as 3, 2, 1... (in a script: 'close 3' ; 'close 2' ; 'close 1'). I hope this helps, Chuck John_Bacon wrote: > Hi > Having opened several files and used the /set dfile/ # command to make > them active, where # is the file descriptor, am having difficulty in > closing files again them? However I try to use close file# I get an > invalid file number error. > Can anyone point out what I'm doing wrong. > Cheers > John Bacon -- Please note that Charles.Seman at noaa.gov should be considered my NOAA email address, not cjs at gfdl.noaa.gov. ******************************************************************** Charles Seman Charles.Seman at noaa.gov U.S. Department of Commerce / NOAA / OAR Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547 201 Forrestal Road fax: (609) 987-5063 Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/ ******************************************************************** "The contents of this message are mine personally and do not necessarily reflect any position of the Government or NOAA." From dasilva at ALUM.MIT.EDU Wed Jul 30 20:52:23 2008 From: dasilva at ALUM.MIT.EDU (Arlindo da Silva) Date: Wed, 30 Jul 2008 20:52:23 -0400 Subject: Creating Shaded Plots from Text File In-Reply-To: Message-ID: On Wed, Jul 30, 2008 at 9:47 AM, Walsh, Brian J wrote: > Do you have an example of how to use a text file with this type of data and > plot it in GrADS. I'm having some trouble understanding the plotting of > station data in GrADS. There is no ready made utility for reading station data in ASCII, you have to write a program to convert it to the proper binary format. Actually, B.-J. Tsuang may have something... -- Arlindo da Silva dasilva at alum.mit.edu From isilda at UEVORA.PT Thu Jul 31 06:27:14 2008 From: isilda at UEVORA.PT (Isilda da Cunha Menezes) Date: Thu, 31 Jul 2008 12:27:14 +0200 Subject: problems read binary grads file with a fortran file Message-ID: Hello! I have a problem when I try to read with a fortran program the binary data generated by the grads that extracts with the fwrite. I extracts the data as follows: 'open prec-temp-mocamb_g3.ctl' 'set gxout fwrite' 'set fwrite -le mocambique.gra' 'set lon -8.967 -7.077' 'set lat 36.999 38.511' ************precipitacao****** t=2 while(t<746) 'set t 't 'set z 1' 'd (totpcp-totpcp(t-1))+(acccon-acccon(t-1))' ************temperature******** 'set z 1' 'd tempc' * Fim do loop no tempo t=t+1 endwhile *'disable fwrite' 'quit' And I read the data created in a fortran program as follows: program ler real a(42,42),x(42,42,2,24,31) integer var,dia,hora,lat,long,soma,j,id,imi,imf,hour open(11,file='mocambique.gra',status='old',form='unformatted') open(12,file='dados-mocamb.dat',status='new',form='formatted') do hour=1,744 do var=1,2 read(11) ((a(long,lat),long=1,42),lat=1,42) enddo enddo write(*,*)'estou aqui' But get an error of core dump. I tried reading the matrix of different ways, but always get the same error. This file .gra is read directly with a fortran program or will have to use another method? If not needed, someone could indicate me where I have the error? Best regards, Isilda From ela at COLA.IGES.ORG Thu Jul 31 15:21:56 2008 From: ela at COLA.IGES.ORG (Eric Altshuler) Date: Thu, 31 Jul 2008 15:21:56 -0400 Subject: problems read binary grads file with a fortran file In-Reply-To: <27015055.120931217531831613.JavaMail.root@mail.iges.org> Message-ID: Hi Isilda, The problem may be that your fortran program opens unit 11 for sequential access but grads writes the fwrite file as flat binary (no control words) which would have to be read in the fortran program by direct access. Opening unit 11 for direct access is one way to solve the problem, but it is easier to tell grads to write the fwrite file sequentially by adding the -sq flag to 'set fwrite': 'set fwrite -sq -le mocambique.gra' Also, you have specified the dimensions of the array 'a' as (42,42). Be sure that the X-Y dimensions in the grads dimension environment are also (42,42) when the fwrite file is written. It's not possible to tell from your 'set lon' and 'set lat' commands how many grid points are in the lat-lon region you have specified. Best regards, Eric L. Altshuler Assistant Research Scientist Center for Ocean-Land-Atmosphere Studies 4041 Powder Mill Road, Suite 302 Calverton, MD 20705-3106 USA E-mail: ela at cola.iges.org Phone: (301) 902-1257 Fax: (301) 595-9793 ----- Original Message ----- From: "Isilda da Cunha Menezes" To: GRADSUSR at LIST.CINECA.IT Sent: Thursday, July 31, 2008 6:27:14 AM GMT -05:00 US/Canada Eastern Subject: problems read binary grads file with a fortran file Hello! I have a problem when I try to read with a fortran program the binary data generated by the grads that extracts with the fwrite. I extracts the data as follows: 'open prec-temp-mocamb_g3.ctl' 'set gxout fwrite' 'set fwrite -le mocambique.gra' 'set lon -8.967 -7.077' 'set lat 36.999 38.511' ************precipitacao****** t=2 while(t<746) 'set t 't 'set z 1' 'd (totpcp-totpcp(t-1))+(acccon-acccon(t-1))' ************temperature******** 'set z 1' 'd tempc' * Fim do loop no tempo t=t+1 endwhile *'disable fwrite' 'quit' And I read the data created in a fortran program as follows: program ler real a(42,42),x(42,42,2,24,31) integer var,dia,hora,lat,long,soma,j,id,imi,imf,hour open(11,file='mocambique.gra',status='old',form='unformatted') open(12,file='dados-mocamb.dat',status='new',form='formatted') do hour=1,744 do var=1,2 read(11) ((a(long,lat),long=1,42),lat=1,42) enddo enddo write(*,*)'estou aqui' But get an error of core dump. I tried reading the matrix of different ways, but always get the same error. This file .gra is read directly with a fortran program or will have to use another method? If not needed, someone could indicate me where I have the error? Best regards, Isilda