Message-ID: Steve, Thanks so much - I didn't even notice this. I thought I was going crazy! -Joe On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan wrote: > Hello Joe, > > Your time query is introducing a space into the "day" variable. PRINTIM > doesn't like spaces in the filename string. Try this instead: > day=substr(result, 11, 9) > > It should work with above change. > > Stephen McMillan > On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: > >> Greetings! >> >> I'm just getting started with GrADS and I'm still a beginner, but I'm >> having a problem that I can't seem to resolve from reading the documentation >> and searching online, including the archives of this mailing list. >> >> Problem: When using a variable in the filename to save the output, I'm >> getting an error. >> >> Here's the relevant code: >> >> 'printim out/RTMA_TMP_'day'_'hr'.png' >> >> where I am defining the variable day and hour using the following: >> 'q time' >> hr=substr(result, 8, 2) >> day=substr(result, 11, 10) >> >> The result of which, when I print as a string or in the title, comes out >> to be 23 for the hour and 16AUG2011. >> >> Expect result: I should have a file output called >> RTMA_TMP_16AUG2011_23.png >> >> What actually is happening: The file output is just called >> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >> error: PRINTIM error: Invalid option; ignored >> >> I've been able to sometimes get this to work, like if I put a few >> characters (not including an underscore or period) in the output filename - >> but I can't figure it out. I'd like to use either underscores, dashes or >> periods to separate the variables in the output file name, and I've seen it >> done in example scripts I've found on the web. >> >> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >> full grads config info: >> ga-> q config >> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds >> hdf5 opendap-grids,stn athena geotiff shapefile >> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >> Copyright (c) 1988-2010 by Brian Doty and the >> Institute for Global Environment and Society (IGES) >> This program is distributed WITHOUT ANY WARRANTY >> See file COPYRIGHT for more information. >> >> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >> >> This version of GrADS has been configured with the following options: >> o Built on a LITTLE ENDIAN machine >> o Athena Widget GUI ENABLED >> o Command line editing ENABLED >> http://tiswww.case.edu/php/chet/readline/rltop.html >> o printim command for image output ENABLED >> http://www.zlib.net >> http://www.libpng.org/pub/png/libpng.html >> http://www.libgd.org/Main_Page >> o GRIB2 interface ENABLED >> http://www.ijg.org >> http://www.ece.uvic.ca/~mdadams/jasper >> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >> g2clib-1.1.8 >> o NetCDF interface ENABLED >> http://www.unidata.ucar.edu/software/netcdf >> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >> o OPeNDAP gridded data interface ENABLED >> o OPeNDAP station data interface ENABLED >> http://iges.org/grads/gadoc/supplibs.html >> libgadap 2.0.oga.1 >> o HDF4 and HDF5 interfaces ENABLED >> http://hdfgroup.org >> HDF 4.2r3 >> HDF5 1.8.4 >> o GeoTIFF and KML/TIFF output ENABLED >> http://www.libtiff.org >> http://geotiff.osgeo.org >> o KML contour output ENABLED >> o Shapefile interface ENABLED >> http://shapelib.maptools.org >> >> For additional information please consult http://iges.org/grads >> >> Not sure if its needed, but here is also the code I am using: >> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >> things as this has been my sandbox script to learn new things.) >> >> Thanks in advance for any and all assistance! >> >> -Joe Moore >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > *************************************************** > 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. > *************************************************** > > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/d853dfa2/attachment.html From dasilva at alum.mit.edu Sun Sep 4 13:22:11 2011 From: dasilva at alum.mit.edu (Arlindo da Silva) Date: Sun, 4 Sep 2011 13:22:11 -0400 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: On Sun, Sep 4, 2011 at 12:46 PM, Joe Moore wrote: > Steve, > > Thanks so much - I didn't even notice this. I thought I was going crazy! > > BTW, since you are using an opengrads build, trying replacing "printim" with "gxyat" in your script for better quality images. Unless you have a performance critical application, the fact that gxyat takes a little longer may not matter. Arlindo -Joe > > > On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan < > smcmillan at planalytics.com> wrote: > >> Hello Joe, >> >> Your time query is introducing a space into the "day" variable. PRINTIM >> doesn't like spaces in the filename string. Try this instead: >> day=substr(result, 11, 9) >> >> It should work with above change. >> >> Stephen McMillan >> On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: >> >>> Greetings! >>> >>> I'm just getting started with GrADS and I'm still a beginner, but I'm >>> having a problem that I can't seem to resolve from reading the documentation >>> and searching online, including the archives of this mailing list. >>> >>> Problem: When using a variable in the filename to save the output, I'm >>> getting an error. >>> >>> Here's the relevant code: >>> >>> 'printim out/RTMA_TMP_'day'_'hr'.png' >>> >>> where I am defining the variable day and hour using the following: >>> 'q time' >>> hr=substr(result, 8, 2) >>> day=substr(result, 11, 10) >>> >>> The result of which, when I print as a string or in the title, comes out >>> to be 23 for the hour and 16AUG2011. >>> >>> Expect result: I should have a file output called >>> RTMA_TMP_16AUG2011_23.png >>> >>> What actually is happening: The file output is just called >>> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >>> error: PRINTIM error: Invalid option; ignored >>> >>> I've been able to sometimes get this to work, like if I put a few >>> characters (not including an underscore or period) in the output filename - >>> but I can't figure it out. I'd like to use either underscores, dashes or >>> periods to separate the variables in the output file name, and I've seen it >>> done in example scripts I've found on the web. >>> >>> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >>> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >>> full grads config info: >>> ga-> q config >>> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf >>> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile >>> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >>> Copyright (c) 1988-2010 by Brian Doty and the >>> Institute for Global Environment and Society (IGES) >>> This program is distributed WITHOUT ANY WARRANTY >>> See file COPYRIGHT for more information. >>> >>> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >>> >>> This version of GrADS has been configured with the following options: >>> o Built on a LITTLE ENDIAN machine >>> o Athena Widget GUI ENABLED >>> o Command line editing ENABLED >>> http://tiswww.case.edu/php/chet/readline/rltop.html >>> o printim command for image output ENABLED >>> http://www.zlib.net >>> http://www.libpng.org/pub/png/libpng.html >>> http://www.libgd.org/Main_Page >>> o GRIB2 interface ENABLED >>> http://www.ijg.org >>> http://www.ece.uvic.ca/~mdadams/jasper >>> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >>> g2clib-1.1.8 >>> o NetCDF interface ENABLED >>> http://www.unidata.ucar.edu/software/netcdf >>> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >>> o OPeNDAP gridded data interface ENABLED >>> o OPeNDAP station data interface ENABLED >>> http://iges.org/grads/gadoc/supplibs.html >>> libgadap 2.0.oga.1 >>> o HDF4 and HDF5 interfaces ENABLED >>> http://hdfgroup.org >>> HDF 4.2r3 >>> HDF5 1.8.4 >>> o GeoTIFF and KML/TIFF output ENABLED >>> http://www.libtiff.org >>> http://geotiff.osgeo.org >>> o KML contour output ENABLED >>> o Shapefile interface ENABLED >>> http://shapelib.maptools.org >>> >>> For additional information please consult http://iges.org/grads >>> >>> Not sure if its needed, but here is also the code I am using: >>> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >>> things as this has been my sandbox script to learn new things.) >>> >>> Thanks in advance for any and all assistance! >>> >>> -Joe Moore >>> >>> _______________________________________________ >>> gradsusr mailing list >>> gradsusr at gradsusr.org >>> http://gradsusr.org/mailman/listinfo/gradsusr >>> >>> >> *************************************************** >> 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. >> *************************************************** >> >> >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/0d8974ee/attachment-0001.html From instantweathermaps at gmail.com Sun Sep 4 13:30:00 2011 From: instantweathermaps at gmail.com (Joey Woodson) Date: Sun, 4 Sep 2011 12:30:00 -0500 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/b63340dd/attachment.html From hallak at model.iag.usp.br Sat Sep 10 07:25:09 2011 From: hallak at model.iag.usp.br (Ricardo Hallak) Date: Sat, 10 Sep 2011 09:25:09 -0200 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: <20110910111216.M81399@model.iag.usp.br> Dear Bernd, I would use the fwrite command: For each computation, you can fwrite the new defined fields in a single binary file, if they have the same horizontal dimensions. With an appropriate ctl that describe correctly the new binary file as the defined fields were at the same time, you can do calculations. Ricardo On Sat, 10 Sep 2011 09:30:06 +0100, Becker, Bernd wrote > Thanks Arlindo, > ? > That would be nice for single numbers. > ? > It as one notch more elaborate: > ? > ? 'set dfile 2' > ? 'set t 2? 33' > ? 'define obs???? = 'var'.2' > > ?'set dfile 3' > ? 'set t 2? 33' > ?'define fx???? = 'var'.3' > > file 2 and 3 have different time frames. > ? > d fx-obs > ? > fails. > And it would be a real bummer to have to calculate the times from one file relative > to the other file. > ? > Compositing in this way would be really hard. > ? > Enjoy the weekend! > Bernd. > ? > -- > Bernd Becker?? Climate Impacts Product Development > Met Office Hadley Centre??? FitzRoy Road?Exeter?? Devon EX1 3PB? United Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > ? > > ----------------------------------------------------------------------- From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva > Sent: 09 September 2011 17:40 > To: GrADS Users Forum > Subject: Re: [gradsusr] define variable independent of time > > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: > > > Simple question, > complicated answer: > > I would like to composite stuff from several ?files. > > Each file has its own time and date info in its own .ctl file. > > How can I add elements from different such files? > > How can I break away from the strict treatment of the time dimension in > grads? > > open file1.ctl > define a=ave(b,t=1,t=34) > close 1 > open file2.ctl > define c= ave(b,t=3,t=9) > close 2 > etc. > > d a+c > Cannot contour grid - all undefined values > > > > Try: > > d a(t=1)+c(t=1) > > -- > Arlindo da Silva > dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/61de2c96/attachment.html From rafanoharana at gmail.com Sat Sep 10 11:57:36 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Sat, 10 Sep 2011 08:57:36 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear John, Thank you for your answer. I do have ogr2ogr and a shapefile of the study area. Thank you anyway for the file. In your command, you only deal with .shp files, and ogr2ogr cannot be used within GrADS. The problem is that I have a .nc file which contains a variable. I want to clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can deal with .nc format (I also tried but not working). Thanks, Serge On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John < Huddleston at cira.colostate.edu> wrote: > Serge**** > > ** ** > > Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to clip > your data.**** > > ** ** > > You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux > use ?apt-get install gdal-bin?.**** > > ** ** > > Use the following command:**** > > ** ** > > ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp**** > > ** ** > > where Study_Area.shp is a shapefile in the attached zip file and > Your_input_data.shp is the shapefile you want to clip.**** > > ** ** > > John Huddleston, PhD**** > > Cooperative Institute for Research in the Atmosphere**** > > ** ** > > ** ** > > ** ** > > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Rafanoharana Serge Claudio > *Sent:* Friday, September 09, 2011 12:59 PM > *To:* GrADS Users Forum > *Subject:* [gradsusr] Crop country boundary or masking an area - still not > working**** > > ** ** > > Dear all, > > I've already posted this message before but I can't get it worked. > If any of you have any idea to crop the study area, please let me know. One > said that it is possible to grab a shapefile of the top part and fill it > with black (I guess he meant re-setting the value of the variable covered by > this shapefile to 0). If you have any experience or any tips on how to do > that, thanks in advance for your help. > > Kind regards, > > Serge**** > > ---------- Forwarded message ---------- > From: *Rafanoharana Serge Claudio* > Date: Fri, Aug 26, 2011 at 7:53 AM > Subject: Crop country boundary > To: GrADS Users Forum > > > Dear all, > > I would like to ask if there is a way to crop the data to a specific area. > In GrADS, if we use 'set mpdset hires', we have the country boundary. > However for my case I only need the study area (please refer to the attached > file). > > I found that http://www.diva-gis.org/datadown provides Grid data. Can we > use the data from this Website to crop it or...? > > Thank you, > > -- > Serge **** > > > > > -- > Serge**** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/2a873fa2/attachment-0001.html From jdduda at iastate.edu Sat Sep 10 13:40:56 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Sat, 10 Sep 2011 12:40:56 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Thanks Arlindo, > > That would be nice for single numbers. > > It as one notch more elaborate: > > 'set dfile 2' > 'set t 2 33' > 'define obs = 'var'.2' > 'set dfile 3' > 'set t 2 33' > 'define fx = 'var'.3' > file 2 and 3 have different time frames. > > d fx-obs > > fails. > And it would be a real bummer to have to calculate the times from one file > relative > to the other file. > > Compositing in this way would be really hard. > > Enjoy the weekend! > Bernd. > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva > *Sent:* 09 September 2011 17:40 > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> >> Simple question, >> complicated answer: >> >> I would like to composite stuff from several files. >> >> Each file has its own time and date info in its own .ctl file. >> >> How can I add elements from different such files? >> >> How can I break away from the strict treatment of the time dimension in >> grads? >> >> open file1.ctl >> define a=ave(b,t=1,t=34) >> close 1 >> open file2.ctl >> define c= ave(b,t=3,t=9) >> close 2 >> etc. >> >> d a+c >> Cannot contour grid - all undefined values >> >> > Try: > > d a(t=1)+c(t=1) > > > -- > Arlindo da Silva > dasilva at alum.mit.edu > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/acae5974/attachment.html From yaqiang.wang at gmail.com Sat Sep 10 22:57:20 2011 From: yaqiang.wang at gmail.com (Yaqiang Wang) Date: Sun, 11 Sep 2011 10:57:20 +0800 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: You can use maskout function to crop your grid data. Or use basemap script to 'crop' your figure. On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio wrote: > Dear John, > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > area. Thank you anyway for the file. > In your command, you only deal with .shp files, and ogr2ogr cannot be used > within GrADS. > The problem is that I have a .nc file which contains a variable. I want to > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > deal with .nc format (I also tried but not working). > > Thanks, > > Serge > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > wrote: >> >> Serge >> >> >> >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to >> clip your data. >> >> >> >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux >> use ?apt-get install gdal-bin?. >> >> >> >> Use the following command: >> >> >> >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp >> >> >> >> where Study_Area.shp is a shapefile in the attached zip file and >> Your_input_data.shp is the shapefile you want to clip. >> >> >> >> John Huddleston, PhD >> >> Cooperative Institute for Research in the Atmosphere >> >> >> >> >> >> >> >> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] >> On Behalf Of Rafanoharana Serge Claudio >> Sent: Friday, September 09, 2011 12:59 PM >> To: GrADS Users Forum >> Subject: [gradsusr] Crop country boundary or masking an area - still not >> working >> >> >> >> Dear all, >> >> I've already posted this message before but I can't get it worked. >> If any of you have any idea to crop the study area, please let me know. >> One said that it is possible to grab a shapefile of the top part? and fill >> it with black (I guess he meant re-setting the value of the variable covered >> by this shapefile to 0). If you have any experience or any tips on how to do >> that, thanks in advance for your help. >> >> Kind regards, >> >> Serge >> >> ---------- Forwarded message ---------- >> From: Rafanoharana Serge Claudio >> Date: Fri, Aug 26, 2011 at 7:53 AM >> Subject: Crop country boundary >> To: GrADS Users Forum >> >> >> Dear all, >> >> I would like to ask if there is a way to crop the data to a specific area. >> In GrADS, if we use 'set mpdset hires', we have the country boundary. >> However for my case I only need the study area (please refer to the attached >> file). >> >> I found that http://www.diva-gis.org/datadown provides Grid data. Can we >> use the data from this Website to crop it or...? >> >> Thank you, >> >> -- >> Serge >> >> >> -- >> Serge >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> > > > > -- > Serge > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- ************************************************* Dr. Yaqiang Wang Chinese Academy of Meteorological Sciences (CAMS) 46, Zhong-Guan-Cun South Avenue Beijing, 100081 China yaqiang.wang at gmail.com ************************************************** From rafanoharana at gmail.com Mon Sep 12 03:27:09 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Mon, 12 Sep 2011 00:27:09 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang wrote: > You can use maskout function to crop your grid data. Or use basemap > script to 'crop' your figure. > > On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > > Dear John, > > > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > > area. Thank you anyway for the file. > > In your command, you only deal with .shp files, and ogr2ogr cannot be > used > > within GrADS. > > The problem is that I have a .nc file which contains a variable. I want > to > > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > > deal with .nc format (I also tried but not working). > > > > Thanks, > > > > Serge > > > > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: > >> > >> Serge > >> > >> > >> > >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to > >> clip your data. > >> > >> > >> > >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on > Linux > >> use ?apt-get install gdal-bin?. > >> > >> > >> > >> Use the following command: > >> > >> > >> > >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp > >> > >> > >> > >> where Study_Area.shp is a shapefile in the attached zip file and > >> Your_input_data.shp is the shapefile you want to clip. > >> > >> > >> > >> John Huddleston, PhD > >> > >> Cooperative Institute for Research in the Atmosphere > >> > >> > >> > >> > >> > >> > >> > >> From: gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] > >> On Behalf Of Rafanoharana Serge Claudio > >> Sent: Friday, September 09, 2011 12:59 PM > >> To: GrADS Users Forum > >> Subject: [gradsusr] Crop country boundary or masking an area - still not > >> working > >> > >> > >> > >> Dear all, > >> > >> I've already posted this message before but I can't get it worked. > >> If any of you have any idea to crop the study area, please let me know. > >> One said that it is possible to grab a shapefile of the top part and > fill > >> it with black (I guess he meant re-setting the value of the variable > covered > >> by this shapefile to 0). If you have any experience or any tips on how > to do > >> that, thanks in advance for your help. > >> > >> Kind regards, > >> > >> Serge > >> > >> ---------- Forwarded message ---------- > >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM > >> Subject: Crop country boundary > >> To: GrADS Users Forum > >> > >> > >> Dear all, > >> > >> I would like to ask if there is a way to crop the data to a specific > area. > >> In GrADS, if we use 'set mpdset hires', we have the country boundary. > >> However for my case I only need the study area (please refer to the > attached > >> file). > >> > >> I found that http://www.diva-gis.org/datadown provides Grid data. Can > we > >> use the data from this Website to crop it or...? > >> > >> Thank you, > >> > >> -- > >> Serge > >> > >> > >> -- > >> Serge > >> > >> _______________________________________________ > >> gradsusr mailing list > >> gradsusr at gradsusr.org > >> http://gradsusr.org/mailman/listinfo/gradsusr > >> > > > > > > > > -- > > Serge > > > > _______________________________________________ > > gradsusr mailing list > > gradsusr at gradsusr.org > > http://gradsusr.org/mailman/listinfo/gradsusr > > > > > > > > -- > ************************************************* > Dr. Yaqiang Wang > Chinese Academy of Meteorological Sciences (CAMS) > 46, Zhong-Guan-Cun South Avenue > Beijing, 100081 > China > > yaqiang.wang at gmail.com > ************************************************** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/545035f7/attachment.html From bernd.becker at metoffice.gov.uk Mon Sep 12 04:03:26 2011 From: bernd.becker at metoffice.gov.uk (Becker, Bernd) Date: Mon, 12 Sep 2011 09:03:26 +0100 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Jeff, That is is the really hard part. Each item fx and obs describes a chunk of a time series. One would have to manually calculate the start and end dates with respect to the reference dfile. If only (and I think Mike explained it a long time ago, and I cannot find the email) one could set file 1 set t 2 34 define whatever set dfile 2 set t 6 39 define theother and happily combine whatever with theother. Think about multimodel ensemble hindcast analysis, trying to composit observed phenomena with modelled phenomena that simply don't happen simultaneously. All the best, Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Jeffrey Duda Sent: 10 September 2011 18:41 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd wrote: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/44da6546/attachment-0001.html From Huddleston at cira.colostate.edu Mon Sep 12 08:00:55 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Mon, 12 Sep 2011 05:00:55 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

, Message-ID: Serge See http://www.ncl.ucar.edu/Document/Tools/ncl_convert2nc.shtml and convert your study area shapefile to a masking netCDF file. Then write a small C/FORTRAN application to read your input netCDF file and the masking netCDF and write out a new netCDF file where the area outside your study area is the missing "undef" value. John Huddleston ________________________________ From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of Rafanoharana Serge Claudio [rafanoharana at gmail.com] Sent: Monday, September 12, 2011 1:27 AM To: GrADS Users Forum Subject: Re: [gradsusr] Crop country boundary or masking an area - still not working Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang > wrote: You can use maskout function to crop your grid data. Or use basemap script to 'crop' your figure. On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > Dear John, > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > area. Thank you anyway for the file. > In your command, you only deal with .shp files, and ogr2ogr cannot be used > within GrADS. > The problem is that I have a .nc file which contains a variable. I want to > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > deal with .nc format (I also tried but not working). > > Thanks, > > Serge > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: >> >> Serge >> >> >> >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to >> clip your data. >> >> >> >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux >> use ?apt-get install gdal-bin?. >> >> >> >> Use the following command: >> >> >> >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp >> >> >> >> where Study_Area.shp is a shapefile in the attached zip file and >> Your_input_data.shp is the shapefile you want to clip. >> >> >> >> John Huddleston, PhD >> >> Cooperative Institute for Research in the Atmosphere >> >> >> >> >> >> >> >> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] >> On Behalf Of Rafanoharana Serge Claudio >> Sent: Friday, September 09, 2011 12:59 PM >> To: GrADS Users Forum >> Subject: [gradsusr] Crop country boundary or masking an area - still not >> working >> >> >> >> Dear all, >> >> I've already posted this message before but I can't get it worked. >> If any of you have any idea to crop the study area, please let me know. >> One said that it is possible to grab a shapefile of the top part and fill >> it with black (I guess he meant re-setting the value of the variable covered >> by this shapefile to 0). If you have any experience or any tips on how to do >> that, thanks in advance for your help. >> >> Kind regards, >> >> Serge >> >> ---------- Forwarded message ---------- >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM >> Subject: Crop country boundary >> To: GrADS Users Forum > >> >> >> Dear all, >> >> I would like to ask if there is a way to crop the data to a specific area. >> In GrADS, if we use 'set mpdset hires', we have the country boundary. >> However for my case I only need the study area (please refer to the attached >> file). >> >> I found that http://www.diva-gis.org/datadown provides Grid data. Can we >> use the data from this Website to crop it or...? >> >> Thank you, >> >> -- >> Serge >> >> >> -- >> Serge >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> > > > > -- > Serge > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- ************************************************* Dr. Yaqiang Wang Chinese Academy of Meteorological Sciences (CAMS) 46, Zhong-Guan-Cun South Avenue Beijing, 100081 China yaqiang.wang at gmail.com ************************************************** _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/52babf1b/attachment.html From jdduda at iastate.edu Mon Sep 12 12:16:59 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Mon, 12 Sep 2011 11:16:59 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Bernd, To help me understand, could you provide a simple example? Jeff On Mon, Sep 12, 2011 at 3:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/59590758/attachment-0001.html From rafanoharana at gmail.com Mon Sep 12 12:29:21 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Mon, 12 Sep 2011 09:29:21 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: On Mon, Sep 12, 2011 at 10:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > Well, have you thought about creating a ctl with each member corresponding to a different file? In grads ensemble implementation you can specify different beginning/ending times for each member, with TDEF specifying the time interval that covers all the members. I typically do this to handle twice daily forecasts: each forecast is an ensemble member, each starting at the initial condition and going out 5 days. Here is how it looks on our GDS: http://portal.nccs.nasa.gov/cgi-lats4d/opendap.cgi?&path=/GEOS-5/fp/0.25_deg/fcast/inst1_2d_hwl_Nx.all.info Of course, this assumes that the files have the same variables and grids, just different time span. Arlindo > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/5ce7c317/attachment-0001.html From jma at cola.iges.org Mon Sep 12 13:57:19 2011 From: jma at cola.iges.org (Jennifer Adams) Date: Mon, 12 Sep 2011 13:57:19 -0400 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: I work with a Debian oon Linux. 2011/9/14 Huddleston, John > CL > > What is your operating system? > > John Huddleston > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On > Behalf Of Christelle LACROIX [cl at meteoconsult.fr] > *Sent:* Tuesday, September 13, 2011 4:42 AM > *To:* GrADS Users Forum > *Subject:* [gradsusr] problem with installation of grads > > Hi, > > > I have a problem with grads, when i run a script i have this error: > "/opt/grads-2.0.a9/bin/grads: error while loading shared libraries: > libX11.so.6: cannot open shared object file: No such file or directory" > > whereas libX11.so.6 exists, when i do a locate libX11 i have this result: > /usr/lib/libX11.so.6 > /usr/lib/libX11.so.6.3.0 > > I don't understand this problem, if you have an idea, it'll be great!! > > Thanks > > CL > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110914/cca92e02/attachment.html From smcmillan at planalytics.com Wed Sep 14 09:10:22 2011 From: smcmillan at planalytics.com (Stephen McMillan) Date: Wed, 14 Sep 2011 09:10:22 -0400 Subject: [gradsusr] Problem accessing NCEP (NOMADS) GrADS-DODS server (from a web server) In-Reply-To: References: Message-ID: Joe, I realize this is 12 hours later, but I just now successfully opened the same dataset. Here's the command output: Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3 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.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff Issue 'q config' command for more information. Loading User Defined Extensions table ... ok. GX Package Initialization: Size = 11 8.5 wrote: > All, > > I'm having a problem using the GrADS Data Server (GrADS-DODS) from NCEP > (NOMADS). I am using the sdfopen command to access the dataset (RTMA 2.5km). > When I run the script locally on my Ubuntu machine, everything works great! > However, when I transfered the script to my web host I get an error. > > Here's the line from my script: > 'sdfopen > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p5'dapdate'/rtma2p5_anl_'currhr'z > ' > (where dapdate=YYYYMMDD and currhr=HH - yes, I've double checked to make > sure the data exists.) > > Here's the output: > ... > Running in Batch mode > No hardcopy metafile open > All files closed; all defined objects released; > All GrADS attributes have been reinitialized > Error: nc_open failed to open file > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p520110913/rtma2p5_anl_23z > NetCDF: I/O failure > gadsdf: Couldn't ingest SDF metadata. > ... > > I've tried searching the web and this (excellent) mailing list assistance > but I can't seem to figure out what's going wrong. I suspect it may have > something to do with my web host. (I use Bluehost- I know they're a somewhat > subpar host, but I'm just a student experimenting with different things!) > > Here's what uname -a gives on my web server: > Linux box566.bluehost.com 2.6.32-46.1.BHsmp #1 SMP Tue Sep 6 12:18:02 MDT > 2011 x86_64 x86_64 x86_64 GNU/Linux > > If there's other tests I can run to probe the issue, let me know. Also, is > this the best way to use the GrADS-DODS server? Should I just download the > grib2 locally instead? > > I was planning to download the grib2 locally, but when I discovered the > GrADS-DODS server I found out that this might not be necessary. I plan to > run this script hourly, so please advise as to which method would be best. > > Thanks a lot! > > -Joe Moore > > The script in question: > http://pastebin.com/MDQ9vk3y > > My GrADS config on my web server: > Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds > hdf5 opendap-grids,stn athena geotiff shapefile > Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 > Copyright (c) 1988-2010 by Brian Doty and the > Institute for Global Environment and Society (IGES) > This program is distributed WITHOUT ANY WARRANTY > See file COPYRIGHT for more information. > > Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu > > This version of GrADS has been configured with the following options: > o Built on a LITTLE ENDIAN machine > o Athena Widget GUI ENABLED > o Command line editing ENABLED > http://tiswww.case.edu/php/chet/readline/rltop.html > o printim command for image output ENABLED > http://www.zlib.net > http://www.libpng.org/pub/png/libpng.html > http://www.libgd.org/Main_Page > o GRIB2 interface ENABLED > http://www.ijg.org > http://www.ece.uvic.ca/~mdadams/jasper > http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 > g2clib-1.1.8 > o NetCDF interface ENABLED > http://www.unidata.ucar.edu/software/netcdf > netcdf 4.1.1 of Jul 10 2010 23:36:25 $ > o OPeNDAP gridded data interface ENABLED > o OPeNDAP station data interface ENABLED > http://iges.org/grads/gadoc/supplibs.html > libgadap 2.0.oga.1 > o HDF4 and HDF5 interfaces ENABLED > http://hdfgroup.org > HDF 4.2r3 > HDF5 1.8.4 > o GeoTIFF and KML/TIFF output ENABLED > http://www.libtiff.org > http://geotiff.osgeo.org > o KML contour output ENABLED > o Shapefile interface ENABLED > http://shapelib.maptools.org > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > *************************************************** 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/20110914/5c694fbd/attachment.html From Huddleston at cira.colostate.edu Wed Sep 14 09:37:46 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Wed, 14 Sep 2011 06:37:46 -0700 Subject: [gradsusr] problem with installation of grads In-Reply-To: References:

Message-ID: CL Try a fresh install of the OpenGrADS. If a re-install of the OpenGrads does not help, then try http://vista.cira.colostate.edu/nco/GrADS/grads-2.0.rc4-i686-linux-ubuntu32.tar.gz or http://vista.cira.colostate.edu/nco/GrADS/grads-2.0.rc4-x86_64-linux-ubuntu64.tar.gz The GrADS executables have been statically linked so there is no dependency on any shared object libraries. John From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Christelle LACROIX Sent: Wednesday, September 14, 2011 6:34 AM To: GrADS Users Forum Subject: Re: [gradsusr] problem with installation of grads I work with a Debian oon Linux. 2011/9/14 Huddleston, John > CL What is your operating system? John Huddleston ________________________________ From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of Christelle LACROIX [cl at meteoconsult.fr] Sent: Tuesday, September 13, 2011 4:42 AM To: GrADS Users Forum Subject: [gradsusr] problem with installation of grads Hi, I have a problem with grads, when i run a script i have this error: "/opt/grads-2.0.a9/bin/grads: error while loading shared libraries: libX11.so.6: cannot open shared object file: No such file or directory" whereas libX11.so.6 exists, when i do a locate libX11 i have this result: /usr/lib/libX11.so.6 /usr/lib/libX11.so.6.3.0 I don't understand this problem, if you have an idea, it'll be great!! Thanks CL _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org

Message-ID: Steve, Thanks so much - I didn't even notice this. I thought I was going crazy! -Joe On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan wrote: > Hello Joe, > > Your time query is introducing a space into the "day" variable. PRINTIM > doesn't like spaces in the filename string. Try this instead: > day=substr(result, 11, 9) > > It should work with above change. > > Stephen McMillan > On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: > >> Greetings! >> >> I'm just getting started with GrADS and I'm still a beginner, but I'm >> having a problem that I can't seem to resolve from reading the documentation >> and searching online, including the archives of this mailing list. >> >> Problem: When using a variable in the filename to save the output, I'm >> getting an error. >> >> Here's the relevant code: >> >> 'printim out/RTMA_TMP_'day'_'hr'.png' >> >> where I am defining the variable day and hour using the following: >> 'q time' >> hr=substr(result, 8, 2) >> day=substr(result, 11, 10) >> >> The result of which, when I print as a string or in the title, comes out >> to be 23 for the hour and 16AUG2011. >> >> Expect result: I should have a file output called >> RTMA_TMP_16AUG2011_23.png >> >> What actually is happening: The file output is just called >> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >> error: PRINTIM error: Invalid option; ignored >> >> I've been able to sometimes get this to work, like if I put a few >> characters (not including an underscore or period) in the output filename - >> but I can't figure it out. I'd like to use either underscores, dashes or >> periods to separate the variables in the output file name, and I've seen it >> done in example scripts I've found on the web. >> >> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >> full grads config info: >> ga-> q config >> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds >> hdf5 opendap-grids,stn athena geotiff shapefile >> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >> Copyright (c) 1988-2010 by Brian Doty and the >> Institute for Global Environment and Society (IGES) >> This program is distributed WITHOUT ANY WARRANTY >> See file COPYRIGHT for more information. >> >> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >> >> This version of GrADS has been configured with the following options: >> o Built on a LITTLE ENDIAN machine >> o Athena Widget GUI ENABLED >> o Command line editing ENABLED >> http://tiswww.case.edu/php/chet/readline/rltop.html >> o printim command for image output ENABLED >> http://www.zlib.net >> http://www.libpng.org/pub/png/libpng.html >> http://www.libgd.org/Main_Page >> o GRIB2 interface ENABLED >> http://www.ijg.org >> http://www.ece.uvic.ca/~mdadams/jasper >> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >> g2clib-1.1.8 >> o NetCDF interface ENABLED >> http://www.unidata.ucar.edu/software/netcdf >> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >> o OPeNDAP gridded data interface ENABLED >> o OPeNDAP station data interface ENABLED >> http://iges.org/grads/gadoc/supplibs.html >> libgadap 2.0.oga.1 >> o HDF4 and HDF5 interfaces ENABLED >> http://hdfgroup.org >> HDF 4.2r3 >> HDF5 1.8.4 >> o GeoTIFF and KML/TIFF output ENABLED >> http://www.libtiff.org >> http://geotiff.osgeo.org >> o KML contour output ENABLED >> o Shapefile interface ENABLED >> http://shapelib.maptools.org >> >> For additional information please consult http://iges.org/grads >> >> Not sure if its needed, but here is also the code I am using: >> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >> things as this has been my sandbox script to learn new things.) >> >> Thanks in advance for any and all assistance! >> >> -Joe Moore >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > *************************************************** > 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. > *************************************************** > > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/d853dfa2/attachment-0002.html From dasilva at alum.mit.edu Sun Sep 4 13:22:11 2011 From: dasilva at alum.mit.edu (Arlindo da Silva) Date: Sun, 4 Sep 2011 13:22:11 -0400 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: On Sun, Sep 4, 2011 at 12:46 PM, Joe Moore wrote: > Steve, > > Thanks so much - I didn't even notice this. I thought I was going crazy! > > BTW, since you are using an opengrads build, trying replacing "printim" with "gxyat" in your script for better quality images. Unless you have a performance critical application, the fact that gxyat takes a little longer may not matter. Arlindo -Joe > > > On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan < > smcmillan at planalytics.com> wrote: > >> Hello Joe, >> >> Your time query is introducing a space into the "day" variable. PRINTIM >> doesn't like spaces in the filename string. Try this instead: >> day=substr(result, 11, 9) >> >> It should work with above change. >> >> Stephen McMillan >> On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: >> >>> Greetings! >>> >>> I'm just getting started with GrADS and I'm still a beginner, but I'm >>> having a problem that I can't seem to resolve from reading the documentation >>> and searching online, including the archives of this mailing list. >>> >>> Problem: When using a variable in the filename to save the output, I'm >>> getting an error. >>> >>> Here's the relevant code: >>> >>> 'printim out/RTMA_TMP_'day'_'hr'.png' >>> >>> where I am defining the variable day and hour using the following: >>> 'q time' >>> hr=substr(result, 8, 2) >>> day=substr(result, 11, 10) >>> >>> The result of which, when I print as a string or in the title, comes out >>> to be 23 for the hour and 16AUG2011. >>> >>> Expect result: I should have a file output called >>> RTMA_TMP_16AUG2011_23.png >>> >>> What actually is happening: The file output is just called >>> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >>> error: PRINTIM error: Invalid option; ignored >>> >>> I've been able to sometimes get this to work, like if I put a few >>> characters (not including an underscore or period) in the output filename - >>> but I can't figure it out. I'd like to use either underscores, dashes or >>> periods to separate the variables in the output file name, and I've seen it >>> done in example scripts I've found on the web. >>> >>> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >>> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >>> full grads config info: >>> ga-> q config >>> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf >>> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile >>> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >>> Copyright (c) 1988-2010 by Brian Doty and the >>> Institute for Global Environment and Society (IGES) >>> This program is distributed WITHOUT ANY WARRANTY >>> See file COPYRIGHT for more information. >>> >>> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >>> >>> This version of GrADS has been configured with the following options: >>> o Built on a LITTLE ENDIAN machine >>> o Athena Widget GUI ENABLED >>> o Command line editing ENABLED >>> http://tiswww.case.edu/php/chet/readline/rltop.html >>> o printim command for image output ENABLED >>> http://www.zlib.net >>> http://www.libpng.org/pub/png/libpng.html >>> http://www.libgd.org/Main_Page >>> o GRIB2 interface ENABLED >>> http://www.ijg.org >>> http://www.ece.uvic.ca/~mdadams/jasper >>> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >>> g2clib-1.1.8 >>> o NetCDF interface ENABLED >>> http://www.unidata.ucar.edu/software/netcdf >>> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >>> o OPeNDAP gridded data interface ENABLED >>> o OPeNDAP station data interface ENABLED >>> http://iges.org/grads/gadoc/supplibs.html >>> libgadap 2.0.oga.1 >>> o HDF4 and HDF5 interfaces ENABLED >>> http://hdfgroup.org >>> HDF 4.2r3 >>> HDF5 1.8.4 >>> o GeoTIFF and KML/TIFF output ENABLED >>> http://www.libtiff.org >>> http://geotiff.osgeo.org >>> o KML contour output ENABLED >>> o Shapefile interface ENABLED >>> http://shapelib.maptools.org >>> >>> For additional information please consult http://iges.org/grads >>> >>> Not sure if its needed, but here is also the code I am using: >>> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >>> things as this has been my sandbox script to learn new things.) >>> >>> Thanks in advance for any and all assistance! >>> >>> -Joe Moore >>> >>> _______________________________________________ >>> gradsusr mailing list >>> gradsusr at gradsusr.org >>> http://gradsusr.org/mailman/listinfo/gradsusr >>> >>> >> *************************************************** >> 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. >> *************************************************** >> >> >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/0d8974ee/attachment-0002.html From instantweathermaps at gmail.com Sun Sep 4 13:30:00 2011 From: instantweathermaps at gmail.com (Joey Woodson) Date: Sun, 4 Sep 2011 12:30:00 -0500 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/b63340dd/attachment-0002.html From hallak at model.iag.usp.br Sat Sep 10 07:25:09 2011 From: hallak at model.iag.usp.br (Ricardo Hallak) Date: Sat, 10 Sep 2011 09:25:09 -0200 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: <20110910111216.M81399@model.iag.usp.br> Dear Bernd, I would use the fwrite command: For each computation, you can fwrite the new defined fields in a single binary file, if they have the same horizontal dimensions. With an appropriate ctl that describe correctly the new binary file as the defined fields were at the same time, you can do calculations. Ricardo On Sat, 10 Sep 2011 09:30:06 +0100, Becker, Bernd wrote > Thanks Arlindo, > ? > That would be nice for single numbers. > ? > It as one notch more elaborate: > ? > ? 'set dfile 2' > ? 'set t 2? 33' > ? 'define obs???? = 'var'.2' > > ?'set dfile 3' > ? 'set t 2? 33' > ?'define fx???? = 'var'.3' > > file 2 and 3 have different time frames. > ? > d fx-obs > ? > fails. > And it would be a real bummer to have to calculate the times from one file relative > to the other file. > ? > Compositing in this way would be really hard. > ? > Enjoy the weekend! > Bernd. > ? > -- > Bernd Becker?? Climate Impacts Product Development > Met Office Hadley Centre??? FitzRoy Road?Exeter?? Devon EX1 3PB? United Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > ? > > ----------------------------------------------------------------------- From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva > Sent: 09 September 2011 17:40 > To: GrADS Users Forum > Subject: Re: [gradsusr] define variable independent of time > > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: > > > Simple question, > complicated answer: > > I would like to composite stuff from several ?files. > > Each file has its own time and date info in its own .ctl file. > > How can I add elements from different such files? > > How can I break away from the strict treatment of the time dimension in > grads? > > open file1.ctl > define a=ave(b,t=1,t=34) > close 1 > open file2.ctl > define c= ave(b,t=3,t=9) > close 2 > etc. > > d a+c > Cannot contour grid - all undefined values > > > > Try: > > d a(t=1)+c(t=1) > > -- > Arlindo da Silva > dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/61de2c96/attachment-0002.html From rafanoharana at gmail.com Sat Sep 10 11:57:36 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Sat, 10 Sep 2011 08:57:36 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear John, Thank you for your answer. I do have ogr2ogr and a shapefile of the study area. Thank you anyway for the file. In your command, you only deal with .shp files, and ogr2ogr cannot be used within GrADS. The problem is that I have a .nc file which contains a variable. I want to clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can deal with .nc format (I also tried but not working). Thanks, Serge On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John < Huddleston at cira.colostate.edu> wrote: > Serge**** > > ** ** > > Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to clip > your data.**** > > ** ** > > You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux > use ?apt-get install gdal-bin?.**** > > ** ** > > Use the following command:**** > > ** ** > > ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp**** > > ** ** > > where Study_Area.shp is a shapefile in the attached zip file and > Your_input_data.shp is the shapefile you want to clip.**** > > ** ** > > John Huddleston, PhD**** > > Cooperative Institute for Research in the Atmosphere**** > > ** ** > > ** ** > > ** ** > > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Rafanoharana Serge Claudio > *Sent:* Friday, September 09, 2011 12:59 PM > *To:* GrADS Users Forum > *Subject:* [gradsusr] Crop country boundary or masking an area - still not > working**** > > ** ** > > Dear all, > > I've already posted this message before but I can't get it worked. > If any of you have any idea to crop the study area, please let me know. One > said that it is possible to grab a shapefile of the top part and fill it > with black (I guess he meant re-setting the value of the variable covered by > this shapefile to 0). If you have any experience or any tips on how to do > that, thanks in advance for your help. > > Kind regards, > > Serge**** > > ---------- Forwarded message ---------- > From: *Rafanoharana Serge Claudio* > Date: Fri, Aug 26, 2011 at 7:53 AM > Subject: Crop country boundary > To: GrADS Users Forum > > > Dear all, > > I would like to ask if there is a way to crop the data to a specific area. > In GrADS, if we use 'set mpdset hires', we have the country boundary. > However for my case I only need the study area (please refer to the attached > file). > > I found that http://www.diva-gis.org/datadown provides Grid data. Can we > use the data from this Website to crop it or...? > > Thank you, > > -- > Serge **** > > > > > -- > Serge**** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/2a873fa2/attachment-0002.html From jdduda at iastate.edu Sat Sep 10 13:40:56 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Sat, 10 Sep 2011 12:40:56 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Thanks Arlindo, > > That would be nice for single numbers. > > It as one notch more elaborate: > > 'set dfile 2' > 'set t 2 33' > 'define obs = 'var'.2' > 'set dfile 3' > 'set t 2 33' > 'define fx = 'var'.3' > file 2 and 3 have different time frames. > > d fx-obs > > fails. > And it would be a real bummer to have to calculate the times from one file > relative > to the other file. > > Compositing in this way would be really hard. > > Enjoy the weekend! > Bernd. > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva > *Sent:* 09 September 2011 17:40 > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> >> Simple question, >> complicated answer: >> >> I would like to composite stuff from several files. >> >> Each file has its own time and date info in its own .ctl file. >> >> How can I add elements from different such files? >> >> How can I break away from the strict treatment of the time dimension in >> grads? >> >> open file1.ctl >> define a=ave(b,t=1,t=34) >> close 1 >> open file2.ctl >> define c= ave(b,t=3,t=9) >> close 2 >> etc. >> >> d a+c >> Cannot contour grid - all undefined values >> >> > Try: > > d a(t=1)+c(t=1) > > > -- > Arlindo da Silva > dasilva at alum.mit.edu > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/acae5974/attachment-0002.html From yaqiang.wang at gmail.com Sat Sep 10 22:57:20 2011 From: yaqiang.wang at gmail.com (Yaqiang Wang) Date: Sun, 11 Sep 2011 10:57:20 +0800 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang wrote: > You can use maskout function to crop your grid data. Or use basemap > script to 'crop' your figure. > > On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > > Dear John, > > > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > > area. Thank you anyway for the file. > > In your command, you only deal with .shp files, and ogr2ogr cannot be > used > > within GrADS. > > The problem is that I have a .nc file which contains a variable. I want > to > > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > > deal with .nc format (I also tried but not working). > > > > Thanks, > > > > Serge > > > > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: > >> > >> Serge > >> > >> > >> > >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to > >> clip your data. > >> > >> > >> > >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on > Linux > >> use ?apt-get install gdal-bin?. > >> > >> > >> > >> Use the following command: > >> > >> > >> > >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp > >> > >> > >> > >> where Study_Area.shp is a shapefile in the attached zip file and > >> Your_input_data.shp is the shapefile you want to clip. > >> > >> > >> > >> John Huddleston, PhD > >> > >> Cooperative Institute for Research in the Atmosphere > >> > >> > >> > >> > >> > >> > >> > >> From: gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] > >> On Behalf Of Rafanoharana Serge Claudio > >> Sent: Friday, September 09, 2011 12:59 PM > >> To: GrADS Users Forum > >> Subject: [gradsusr] Crop country boundary or masking an area - still not > >> working > >> > >> > >> > >> Dear all, > >> > >> I've already posted this message before but I can't get it worked. > >> If any of you have any idea to crop the study area, please let me know. > >> One said that it is possible to grab a shapefile of the top part and > fill > >> it with black (I guess he meant re-setting the value of the variable > covered > >> by this shapefile to 0). If you have any experience or any tips on how > to do > >> that, thanks in advance for your help. > >> > >> Kind regards, > >> > >> Serge > >> > >> ---------- Forwarded message ---------- > >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM > >> Subject: Crop country boundary > >> To: GrADS Users Forum > >> > >> > >> Dear all, > >> > >> I would like to ask if there is a way to crop the data to a specific > area. > >> In GrADS, if we use 'set mpdset hires', we have the country boundary. > >> However for my case I only need the study area (please refer to the > attached > >> file). > >> > >> I found that http://www.diva-gis.org/datadown provides Grid data. Can > we > >> use the data from this Website to crop it or...? > >> > >> Thank you, > >> > >> -- > >> Serge > >> > >> > >> -- > >> Serge > >> > >> _______________________________________________ > >> gradsusr mailing list > >> gradsusr at gradsusr.org > >> http://gradsusr.org/mailman/listinfo/gradsusr > >> > > > > > > > > -- > > Serge > > > > _______________________________________________ > > gradsusr mailing list > > gradsusr at gradsusr.org > > http://gradsusr.org/mailman/listinfo/gradsusr > > > > > > > > -- > ************************************************* > Dr. Yaqiang Wang > Chinese Academy of Meteorological Sciences (CAMS) > 46, Zhong-Guan-Cun South Avenue > Beijing, 100081 > China > > yaqiang.wang at gmail.com > ************************************************** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/545035f7/attachment-0002.html From bernd.becker at metoffice.gov.uk Mon Sep 12 04:03:26 2011 From: bernd.becker at metoffice.gov.uk (Becker, Bernd) Date: Mon, 12 Sep 2011 09:03:26 +0100 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Jeff, That is is the really hard part. Each item fx and obs describes a chunk of a time series. One would have to manually calculate the start and end dates with respect to the reference dfile. If only (and I think Mike explained it a long time ago, and I cannot find the email) one could set file 1 set t 2 34 define whatever set dfile 2 set t 6 39 define theother and happily combine whatever with theother. Think about multimodel ensemble hindcast analysis, trying to composit observed phenomena with modelled phenomena that simply don't happen simultaneously. All the best, Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Jeffrey Duda Sent: 10 September 2011 18:41 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd wrote: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/44da6546/attachment-0002.html From Huddleston at cira.colostate.edu Mon Sep 12 08:00:55 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Mon, 12 Sep 2011 05:00:55 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

, Message-ID: Serge See http://www.ncl.ucar.edu/Document/Tools/ncl_convert2nc.shtml and convert your study area shapefile to a masking netCDF file. Then write a small C/FORTRAN application to read your input netCDF file and the masking netCDF and write out a new netCDF file where the area outside your study area is the missing "undef" value. John Huddleston ________________________________ From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of Rafanoharana Serge Claudio [rafanoharana at gmail.com] Sent: Monday, September 12, 2011 1:27 AM To: GrADS Users Forum Subject: Re: [gradsusr] Crop country boundary or masking an area - still not working Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang > wrote: You can use maskout function to crop your grid data. Or use basemap script to 'crop' your figure. On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > Dear John, > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > area. Thank you anyway for the file. > In your command, you only deal with .shp files, and ogr2ogr cannot be used > within GrADS. > The problem is that I have a .nc file which contains a variable. I want to > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > deal with .nc format (I also tried but not working). > > Thanks, > > Serge > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: >> >> Serge >> >> >> >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to >> clip your data. >> >> >> >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux >> use ?apt-get install gdal-bin?. >> >> >> >> Use the following command: >> >> >> >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp >> >> >> >> where Study_Area.shp is a shapefile in the attached zip file and >> Your_input_data.shp is the shapefile you want to clip. >> >> >> >> John Huddleston, PhD >> >> Cooperative Institute for Research in the Atmosphere >> >> >> >> >> >> >> >> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] >> On Behalf Of Rafanoharana Serge Claudio >> Sent: Friday, September 09, 2011 12:59 PM >> To: GrADS Users Forum >> Subject: [gradsusr] Crop country boundary or masking an area - still not >> working >> >> >> >> Dear all, >> >> I've already posted this message before but I can't get it worked. >> If any of you have any idea to crop the study area, please let me know. >> One said that it is possible to grab a shapefile of the top part and fill >> it with black (I guess he meant re-setting the value of the variable covered >> by this shapefile to 0). If you have any experience or any tips on how to do >> that, thanks in advance for your help. >> >> Kind regards, >> >> Serge >> >> ---------- Forwarded message ---------- >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM >> Subject: Crop country boundary >> To: GrADS Users Forum > >> >> >> Dear all, >> >> I would like to ask if there is a way to crop the data to a specific area. >> In GrADS, if we use 'set mpdset hires', we have the country boundary. >> However for my case I only need the study area (please refer to the attached >> file). >> >> I found that http://www.diva-gis.org/datadown provides Grid data. Can we >> use the data from this Website to crop it or...? >> >> Thank you, >> >> -- >> Serge >> >> >> -- >> Serge >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> > > > > -- > Serge > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- ************************************************* Dr. Yaqiang Wang Chinese Academy of Meteorological Sciences (CAMS) 46, Zhong-Guan-Cun South Avenue Beijing, 100081 China yaqiang.wang at gmail.com ************************************************** _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/52babf1b/attachment-0002.html From jdduda at iastate.edu Mon Sep 12 12:16:59 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Mon, 12 Sep 2011 11:16:59 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Bernd, To help me understand, could you provide a simple example? Jeff On Mon, Sep 12, 2011 at 3:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/59590758/attachment-0002.html From rafanoharana at gmail.com Mon Sep 12 12:29:21 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Mon, 12 Sep 2011 09:29:21 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: On Mon, Sep 12, 2011 at 10:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > Well, have you thought about creating a ctl with each member corresponding to a different file? In grads ensemble implementation you can specify different beginning/ending times for each member, with TDEF specifying the time interval that covers all the members. I typically do this to handle twice daily forecasts: each forecast is an ensemble member, each starting at the initial condition and going out 5 days. Here is how it looks on our GDS: http://portal.nccs.nasa.gov/cgi-lats4d/opendap.cgi?&path=/GEOS-5/fp/0.25_deg/fcast/inst1_2d_hwl_Nx.all.info Of course, this assumes that the files have the same variables and grids, just different time span. Arlindo > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/5ce7c317/attachment-0002.html From jma at cola.iges.org Mon Sep 12 13:57:19 2011 From: jma at cola.iges.org (Jennifer Adams) Date: Mon, 12 Sep 2011 13:57:19 -0400 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: I work with a Debian oon Linux. 2011/9/14 Huddleston, John > CL > > What is your operating system? > > John Huddleston > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On > Behalf Of Christelle LACROIX [cl at meteoconsult.fr] > *Sent:* Tuesday, September 13, 2011 4:42 AM > *To:* GrADS Users Forum > *Subject:* [gradsusr] problem with installation of grads > > Hi, > > > I have a problem with grads, when i run a script i have this error: > "/opt/grads-2.0.a9/bin/grads: error while loading shared libraries: > libX11.so.6: cannot open shared object file: No such file or directory" > > whereas libX11.so.6 exists, when i do a locate libX11 i have this result: > /usr/lib/libX11.so.6 > /usr/lib/libX11.so.6.3.0 > > I don't understand this problem, if you have an idea, it'll be great!! > > Thanks > > CL > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110914/cca92e02/attachment-0002.html From smcmillan at planalytics.com Wed Sep 14 09:10:22 2011 From: smcmillan at planalytics.com (Stephen McMillan) Date: Wed, 14 Sep 2011 09:10:22 -0400 Subject: [gradsusr] Problem accessing NCEP (NOMADS) GrADS-DODS server (from a web server) In-Reply-To: References: Message-ID: Joe, I realize this is 12 hours later, but I just now successfully opened the same dataset. Here's the command output: Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3 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.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff Issue 'q config' command for more information. Loading User Defined Extensions table ... ok. GX Package Initialization: Size = 11 8.5 wrote: > All, > > I'm having a problem using the GrADS Data Server (GrADS-DODS) from NCEP > (NOMADS). I am using the sdfopen command to access the dataset (RTMA 2.5km). > When I run the script locally on my Ubuntu machine, everything works great! > However, when I transfered the script to my web host I get an error. > > Here's the line from my script: > 'sdfopen > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p5'dapdate'/rtma2p5_anl_'currhr'z > ' > (where dapdate=YYYYMMDD and currhr=HH - yes, I've double checked to make > sure the data exists.) > > Here's the output: > ... > Running in Batch mode > No hardcopy metafile open > All files closed; all defined objects released; > All GrADS attributes have been reinitialized > Error: nc_open failed to open file > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p520110913/rtma2p5_anl_23z > NetCDF: I/O failure > gadsdf: Couldn't ingest SDF metadata. > ... > > I've tried searching the web and this (excellent) mailing list assistance > but I can't seem to figure out what's going wrong. I suspect it may have > something to do with my web host. (I use Bluehost- I know they're a somewhat > subpar host, but I'm just a student experimenting with different things!) > > Here's what uname -a gives on my web server: > Linux box566.bluehost.com 2.6.32-46.1.BHsmp #1 SMP Tue Sep 6 12:18:02 MDT > 2011 x86_64 x86_64 x86_64 GNU/Linux > > If there's other tests I can run to probe the issue, let me know. Also, is > this the best way to use the GrADS-DODS server? Should I just download the > grib2 locally instead? > > I was planning to download the grib2 locally, but when I discovered the > GrADS-DODS server I found out that this might not be necessary. I plan to > run this script hourly, so please advise as to which method would be best. > > Thanks a lot! > > -Joe Moore > > The script in question: > http://pastebin.com/MDQ9vk3y > > My GrADS config on my web server: > Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds > hdf5 opendap-grids,stn athena geotiff shapefile > Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 > Copyright (c) 1988-2010 by Brian Doty and the > Institute for Global Environment and Society (IGES) > This program is distributed WITHOUT ANY WARRANTY > See file COPYRIGHT for more information. > > Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu > > This version of GrADS has been configured with the following options: > o Built on a LITTLE ENDIAN machine > o Athena Widget GUI ENABLED > o Command line editing ENABLED > http://tiswww.case.edu/php/chet/readline/rltop.html > o printim command for image output ENABLED > http://www.zlib.net > http://www.libpng.org/pub/png/libpng.html > http://www.libgd.org/Main_Page > o GRIB2 interface ENABLED > http://www.ijg.org > http://www.ece.uvic.ca/~mdadams/jasper > http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 > g2clib-1.1.8 > o NetCDF interface ENABLED > http://www.unidata.ucar.edu/software/netcdf > netcdf 4.1.1 of Jul 10 2010 23:36:25 $ > o OPeNDAP gridded data interface ENABLED > o OPeNDAP station data interface ENABLED > http://iges.org/grads/gadoc/supplibs.html > libgadap 2.0.oga.1 > o HDF4 and HDF5 interfaces ENABLED > http://hdfgroup.org > HDF 4.2r3 > HDF5 1.8.4 > o GeoTIFF and KML/TIFF output ENABLED > http://www.libtiff.org > http://geotiff.osgeo.org > o KML contour output ENABLED > o Shapefile interface ENABLED > http://shapelib.maptools.org > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > *************************************************** 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/20110914/5c694fbd/attachment-0002.html From Huddleston at cira.colostate.edu Wed Sep 14 09:37:46 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Wed, 14 Sep 2011 06:37:46 -0700 Subject: [gradsusr] problem with installation of grads In-Reply-To: References:

Message-ID: Steve, Thanks so much - I didn't even notice this. I thought I was going crazy! -Joe On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan wrote: > Hello Joe, > > Your time query is introducing a space into the "day" variable. PRINTIM > doesn't like spaces in the filename string. Try this instead: > day=substr(result, 11, 9) > > It should work with above change. > > Stephen McMillan > On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: > >> Greetings! >> >> I'm just getting started with GrADS and I'm still a beginner, but I'm >> having a problem that I can't seem to resolve from reading the documentation >> and searching online, including the archives of this mailing list. >> >> Problem: When using a variable in the filename to save the output, I'm >> getting an error. >> >> Here's the relevant code: >> >> 'printim out/RTMA_TMP_'day'_'hr'.png' >> >> where I am defining the variable day and hour using the following: >> 'q time' >> hr=substr(result, 8, 2) >> day=substr(result, 11, 10) >> >> The result of which, when I print as a string or in the title, comes out >> to be 23 for the hour and 16AUG2011. >> >> Expect result: I should have a file output called >> RTMA_TMP_16AUG2011_23.png >> >> What actually is happening: The file output is just called >> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >> error: PRINTIM error: Invalid option; ignored >> >> I've been able to sometimes get this to work, like if I put a few >> characters (not including an underscore or period) in the output filename - >> but I can't figure it out. I'd like to use either underscores, dashes or >> periods to separate the variables in the output file name, and I've seen it >> done in example scripts I've found on the web. >> >> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >> full grads config info: >> ga-> q config >> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds >> hdf5 opendap-grids,stn athena geotiff shapefile >> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >> Copyright (c) 1988-2010 by Brian Doty and the >> Institute for Global Environment and Society (IGES) >> This program is distributed WITHOUT ANY WARRANTY >> See file COPYRIGHT for more information. >> >> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >> >> This version of GrADS has been configured with the following options: >> o Built on a LITTLE ENDIAN machine >> o Athena Widget GUI ENABLED >> o Command line editing ENABLED >> http://tiswww.case.edu/php/chet/readline/rltop.html >> o printim command for image output ENABLED >> http://www.zlib.net >> http://www.libpng.org/pub/png/libpng.html >> http://www.libgd.org/Main_Page >> o GRIB2 interface ENABLED >> http://www.ijg.org >> http://www.ece.uvic.ca/~mdadams/jasper >> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >> g2clib-1.1.8 >> o NetCDF interface ENABLED >> http://www.unidata.ucar.edu/software/netcdf >> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >> o OPeNDAP gridded data interface ENABLED >> o OPeNDAP station data interface ENABLED >> http://iges.org/grads/gadoc/supplibs.html >> libgadap 2.0.oga.1 >> o HDF4 and HDF5 interfaces ENABLED >> http://hdfgroup.org >> HDF 4.2r3 >> HDF5 1.8.4 >> o GeoTIFF and KML/TIFF output ENABLED >> http://www.libtiff.org >> http://geotiff.osgeo.org >> o KML contour output ENABLED >> o Shapefile interface ENABLED >> http://shapelib.maptools.org >> >> For additional information please consult http://iges.org/grads >> >> Not sure if its needed, but here is also the code I am using: >> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >> things as this has been my sandbox script to learn new things.) >> >> Thanks in advance for any and all assistance! >> >> -Joe Moore >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > *************************************************** > 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. > *************************************************** > > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/d853dfa2/attachment-0003.html From dasilva at alum.mit.edu Sun Sep 4 13:22:11 2011 From: dasilva at alum.mit.edu (Arlindo da Silva) Date: Sun, 4 Sep 2011 13:22:11 -0400 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: On Sun, Sep 4, 2011 at 12:46 PM, Joe Moore wrote: > Steve, > > Thanks so much - I didn't even notice this. I thought I was going crazy! > > BTW, since you are using an opengrads build, trying replacing "printim" with "gxyat" in your script for better quality images. Unless you have a performance critical application, the fact that gxyat takes a little longer may not matter. Arlindo -Joe > > > On Sun, Sep 4, 2011 at 11:51 AM, Stephen McMillan < > smcmillan at planalytics.com> wrote: > >> Hello Joe, >> >> Your time query is introducing a space into the "day" variable. PRINTIM >> doesn't like spaces in the filename string. Try this instead: >> day=substr(result, 11, 9) >> >> It should work with above change. >> >> Stephen McMillan >> On Sat, Sep 3, 2011 at 4:53 PM, Joe Moore wrote: >> >>> Greetings! >>> >>> I'm just getting started with GrADS and I'm still a beginner, but I'm >>> having a problem that I can't seem to resolve from reading the documentation >>> and searching online, including the archives of this mailing list. >>> >>> Problem: When using a variable in the filename to save the output, I'm >>> getting an error. >>> >>> Here's the relevant code: >>> >>> 'printim out/RTMA_TMP_'day'_'hr'.png' >>> >>> where I am defining the variable day and hour using the following: >>> 'q time' >>> hr=substr(result, 8, 2) >>> day=substr(result, 11, 10) >>> >>> The result of which, when I print as a string or in the title, comes out >>> to be 23 for the hour and 16AUG2011. >>> >>> Expect result: I should have a file output called >>> RTMA_TMP_16AUG2011_23.png >>> >>> What actually is happening: The file output is just called >>> RTMA_TMP_23AUG2011 (that's it- no .png) and when I run I get the following >>> error: PRINTIM error: Invalid option; ignored >>> >>> I've been able to sometimes get this to work, like if I put a few >>> characters (not including an underscore or period) in the output filename - >>> but I can't figure it out. I'd like to use either underscores, dashes or >>> periods to separate the variables in the output file name, and I've seen it >>> done in example scripts I've found on the web. >>> >>> Configuration: I am running OpenGrADS 2.0.a9.oga.1 on Ubuntu 11.04. The >>> data is coming from the RTMA grib2 from the NOAA NOMADS server. Here is my >>> full grads config info: >>> ga-> q config >>> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf >>> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile >>> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 >>> Copyright (c) 1988-2010 by Brian Doty and the >>> Institute for Global Environment and Society (IGES) >>> This program is distributed WITHOUT ANY WARRANTY >>> See file COPYRIGHT for more information. >>> >>> Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu >>> >>> This version of GrADS has been configured with the following options: >>> o Built on a LITTLE ENDIAN machine >>> o Athena Widget GUI ENABLED >>> o Command line editing ENABLED >>> http://tiswww.case.edu/php/chet/readline/rltop.html >>> o printim command for image output ENABLED >>> http://www.zlib.net >>> http://www.libpng.org/pub/png/libpng.html >>> http://www.libgd.org/Main_Page >>> o GRIB2 interface ENABLED >>> http://www.ijg.org >>> http://www.ece.uvic.ca/~mdadams/jasper >>> http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 >>> g2clib-1.1.8 >>> o NetCDF interface ENABLED >>> http://www.unidata.ucar.edu/software/netcdf >>> netcdf 4.1.1 of Jul 10 2010 23:36:25 $ >>> o OPeNDAP gridded data interface ENABLED >>> o OPeNDAP station data interface ENABLED >>> http://iges.org/grads/gadoc/supplibs.html >>> libgadap 2.0.oga.1 >>> o HDF4 and HDF5 interfaces ENABLED >>> http://hdfgroup.org >>> HDF 4.2r3 >>> HDF5 1.8.4 >>> o GeoTIFF and KML/TIFF output ENABLED >>> http://www.libtiff.org >>> http://geotiff.osgeo.org >>> o KML contour output ENABLED >>> o Shapefile interface ENABLED >>> http://shapelib.maptools.org >>> >>> For additional information please consult http://iges.org/grads >>> >>> Not sure if its needed, but here is also the code I am using: >>> http://pastebin.com/hHuXR6ae (There's a few comments and irrelevant >>> things as this has been my sandbox script to learn new things.) >>> >>> Thanks in advance for any and all assistance! >>> >>> -Joe Moore >>> >>> _______________________________________________ >>> gradsusr mailing list >>> gradsusr at gradsusr.org >>> http://gradsusr.org/mailman/listinfo/gradsusr >>> >>> >> *************************************************** >> 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. >> *************************************************** >> >> >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110904/0d8974ee/attachment-0003.html From instantweathermaps at gmail.com Sun Sep 4 13:30:00 2011 From: instantweathermaps at gmail.com (Joey Woodson) Date: Sun, 4 Sep 2011 12:30:00 -0500 Subject: [gradsusr] printim error due to variable in filename In-Reply-To: References:

Message-ID: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/b63340dd/attachment-0003.html From hallak at model.iag.usp.br Sat Sep 10 07:25:09 2011 From: hallak at model.iag.usp.br (Ricardo Hallak) Date: Sat, 10 Sep 2011 09:25:09 -0200 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: <20110910111216.M81399@model.iag.usp.br> Dear Bernd, I would use the fwrite command: For each computation, you can fwrite the new defined fields in a single binary file, if they have the same horizontal dimensions. With an appropriate ctl that describe correctly the new binary file as the defined fields were at the same time, you can do calculations. Ricardo On Sat, 10 Sep 2011 09:30:06 +0100, Becker, Bernd wrote > Thanks Arlindo, > ? > That would be nice for single numbers. > ? > It as one notch more elaborate: > ? > ? 'set dfile 2' > ? 'set t 2? 33' > ? 'define obs???? = 'var'.2' > > ?'set dfile 3' > ? 'set t 2? 33' > ?'define fx???? = 'var'.3' > > file 2 and 3 have different time frames. > ? > d fx-obs > ? > fails. > And it would be a real bummer to have to calculate the times from one file relative > to the other file. > ? > Compositing in this way would be really hard. > ? > Enjoy the weekend! > Bernd. > ? > -- > Bernd Becker?? Climate Impacts Product Development > Met Office Hadley Centre??? FitzRoy Road?Exeter?? Devon EX1 3PB? United Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > ? > > ----------------------------------------------------------------------- From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva > Sent: 09 September 2011 17:40 > To: GrADS Users Forum > Subject: Re: [gradsusr] define variable independent of time > > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: > > > Simple question, > complicated answer: > > I would like to composite stuff from several ?files. > > Each file has its own time and date info in its own .ctl file. > > How can I add elements from different such files? > > How can I break away from the strict treatment of the time dimension in > grads? > > open file1.ctl > define a=ave(b,t=1,t=34) > close 1 > open file2.ctl > define c= ave(b,t=3,t=9) > close 2 > etc. > > d a+c > Cannot contour grid - all undefined values > > > > Try: > > d a(t=1)+c(t=1) > > -- > Arlindo da Silva > dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/61de2c96/attachment-0003.html From rafanoharana at gmail.com Sat Sep 10 11:57:36 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Sat, 10 Sep 2011 08:57:36 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear John, Thank you for your answer. I do have ogr2ogr and a shapefile of the study area. Thank you anyway for the file. In your command, you only deal with .shp files, and ogr2ogr cannot be used within GrADS. The problem is that I have a .nc file which contains a variable. I want to clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can deal with .nc format (I also tried but not working). Thanks, Serge On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John < Huddleston at cira.colostate.edu> wrote: > Serge**** > > ** ** > > Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to clip > your data.**** > > ** ** > > You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux > use ?apt-get install gdal-bin?.**** > > ** ** > > Use the following command:**** > > ** ** > > ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp**** > > ** ** > > where Study_Area.shp is a shapefile in the attached zip file and > Your_input_data.shp is the shapefile you want to clip.**** > > ** ** > > John Huddleston, PhD**** > > Cooperative Institute for Research in the Atmosphere**** > > ** ** > > ** ** > > ** ** > > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Rafanoharana Serge Claudio > *Sent:* Friday, September 09, 2011 12:59 PM > *To:* GrADS Users Forum > *Subject:* [gradsusr] Crop country boundary or masking an area - still not > working**** > > ** ** > > Dear all, > > I've already posted this message before but I can't get it worked. > If any of you have any idea to crop the study area, please let me know. One > said that it is possible to grab a shapefile of the top part and fill it > with black (I guess he meant re-setting the value of the variable covered by > this shapefile to 0). If you have any experience or any tips on how to do > that, thanks in advance for your help. > > Kind regards, > > Serge**** > > ---------- Forwarded message ---------- > From: *Rafanoharana Serge Claudio* > Date: Fri, Aug 26, 2011 at 7:53 AM > Subject: Crop country boundary > To: GrADS Users Forum > > > Dear all, > > I would like to ask if there is a way to crop the data to a specific area. > In GrADS, if we use 'set mpdset hires', we have the country boundary. > However for my case I only need the study area (please refer to the attached > file). > > I found that http://www.diva-gis.org/datadown provides Grid data. Can we > use the data from this Website to crop it or...? > > Thank you, > > -- > Serge **** > > > > > -- > Serge**** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/2a873fa2/attachment-0003.html From jdduda at iastate.edu Sat Sep 10 13:40:56 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Sat, 10 Sep 2011 12:40:56 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Thanks Arlindo, > > That would be nice for single numbers. > > It as one notch more elaborate: > > 'set dfile 2' > 'set t 2 33' > 'define obs = 'var'.2' > 'set dfile 3' > 'set t 2 33' > 'define fx = 'var'.3' > file 2 and 3 have different time frames. > > d fx-obs > > fails. > And it would be a real bummer to have to calculate the times from one file > relative > to the other file. > > Compositing in this way would be really hard. > > Enjoy the weekend! > Bernd. > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva > *Sent:* 09 September 2011 17:40 > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> >> Simple question, >> complicated answer: >> >> I would like to composite stuff from several files. >> >> Each file has its own time and date info in its own .ctl file. >> >> How can I add elements from different such files? >> >> How can I break away from the strict treatment of the time dimension in >> grads? >> >> open file1.ctl >> define a=ave(b,t=1,t=34) >> close 1 >> open file2.ctl >> define c= ave(b,t=3,t=9) >> close 2 >> etc. >> >> d a+c >> Cannot contour grid - all undefined values >> >> > Try: > > d a(t=1)+c(t=1) > > > -- > Arlindo da Silva > dasilva at alum.mit.edu > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110910/acae5974/attachment-0003.html From yaqiang.wang at gmail.com Sat Sep 10 22:57:20 2011 From: yaqiang.wang at gmail.com (Yaqiang Wang) Date: Sun, 11 Sep 2011 10:57:20 +0800 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang wrote: > You can use maskout function to crop your grid data. Or use basemap > script to 'crop' your figure. > > On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > > Dear John, > > > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > > area. Thank you anyway for the file. > > In your command, you only deal with .shp files, and ogr2ogr cannot be > used > > within GrADS. > > The problem is that I have a .nc file which contains a variable. I want > to > > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > > deal with .nc format (I also tried but not working). > > > > Thanks, > > > > Serge > > > > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: > >> > >> Serge > >> > >> > >> > >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to > >> clip your data. > >> > >> > >> > >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on > Linux > >> use ?apt-get install gdal-bin?. > >> > >> > >> > >> Use the following command: > >> > >> > >> > >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp > >> > >> > >> > >> where Study_Area.shp is a shapefile in the attached zip file and > >> Your_input_data.shp is the shapefile you want to clip. > >> > >> > >> > >> John Huddleston, PhD > >> > >> Cooperative Institute for Research in the Atmosphere > >> > >> > >> > >> > >> > >> > >> > >> From: gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] > >> On Behalf Of Rafanoharana Serge Claudio > >> Sent: Friday, September 09, 2011 12:59 PM > >> To: GrADS Users Forum > >> Subject: [gradsusr] Crop country boundary or masking an area - still not > >> working > >> > >> > >> > >> Dear all, > >> > >> I've already posted this message before but I can't get it worked. > >> If any of you have any idea to crop the study area, please let me know. > >> One said that it is possible to grab a shapefile of the top part and > fill > >> it with black (I guess he meant re-setting the value of the variable > covered > >> by this shapefile to 0). If you have any experience or any tips on how > to do > >> that, thanks in advance for your help. > >> > >> Kind regards, > >> > >> Serge > >> > >> ---------- Forwarded message ---------- > >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM > >> Subject: Crop country boundary > >> To: GrADS Users Forum > >> > >> > >> Dear all, > >> > >> I would like to ask if there is a way to crop the data to a specific > area. > >> In GrADS, if we use 'set mpdset hires', we have the country boundary. > >> However for my case I only need the study area (please refer to the > attached > >> file). > >> > >> I found that http://www.diva-gis.org/datadown provides Grid data. Can > we > >> use the data from this Website to crop it or...? > >> > >> Thank you, > >> > >> -- > >> Serge > >> > >> > >> -- > >> Serge > >> > >> _______________________________________________ > >> gradsusr mailing list > >> gradsusr at gradsusr.org > >> http://gradsusr.org/mailman/listinfo/gradsusr > >> > > > > > > > > -- > > Serge > > > > _______________________________________________ > > gradsusr mailing list > > gradsusr at gradsusr.org > > http://gradsusr.org/mailman/listinfo/gradsusr > > > > > > > > -- > ************************************************* > Dr. Yaqiang Wang > Chinese Academy of Meteorological Sciences (CAMS) > 46, Zhong-Guan-Cun South Avenue > Beijing, 100081 > China > > yaqiang.wang at gmail.com > ************************************************** > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/545035f7/attachment-0003.html From bernd.becker at metoffice.gov.uk Mon Sep 12 04:03:26 2011 From: bernd.becker at metoffice.gov.uk (Becker, Bernd) Date: Mon, 12 Sep 2011 09:03:26 +0100 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Jeff, That is is the really hard part. Each item fx and obs describes a chunk of a time series. One would have to manually calculate the start and end dates with respect to the reference dfile. If only (and I think Mike explained it a long time ago, and I cannot find the email) one could set file 1 set t 2 34 define whatever set dfile 2 set t 6 39 define theother and happily combine whatever with theother. Think about multimodel ensemble hindcast analysis, trying to composit observed phenomena with modelled phenomena that simply don't happen simultaneously. All the best, Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Jeffrey Duda Sent: 10 September 2011 18:41 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time You can reference the time of a variable when you use it in a statement. For example, you could do 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' This should help you do what you want to do. Jeff Duda On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd wrote: Thanks Arlindo, That would be nice for single numbers. It as one notch more elaborate: 'set dfile 2' 'set t 2 33' 'define obs = 'var'.2' 'set dfile 3' 'set t 2 33' 'define fx = 'var'.3' file 2 and 3 have different time frames. d fx-obs fails. And it would be a real bummer to have to calculate the times from one file relative to the other file. Compositing in this way would be really hard. Enjoy the weekend! Bernd. -- Bernd Becker Climate Impacts Product Development Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United Kingdom Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 E-mail:bernd.becker at metoffice.gov.uk - http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/ bernd-becker http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game ________________________________ From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Arlindo da Silva Sent: 09 September 2011 17:40 To: GrADS Users Forum Subject: Re: [gradsusr] define variable independent of time On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd wrote: Simple question, complicated answer: I would like to composite stuff from several files. Each file has its own time and date info in its own .ctl file. How can I add elements from different such files? How can I break away from the strict treatment of the time dimension in grads? open file1.ctl define a=ave(b,t=1,t=34) close 1 open file2.ctl define c= ave(b,t=3,t=9) close 2 etc. d a+c Cannot contour grid - all undefined values Try: d a(t=1)+c(t=1) -- Arlindo da Silva dasilva at alum.mit.edu _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/44da6546/attachment-0003.html From Huddleston at cira.colostate.edu Mon Sep 12 08:00:55 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Mon, 12 Sep 2011 05:00:55 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

, Message-ID: Serge See http://www.ncl.ucar.edu/Document/Tools/ncl_convert2nc.shtml and convert your study area shapefile to a masking netCDF file. Then write a small C/FORTRAN application to read your input netCDF file and the masking netCDF and write out a new netCDF file where the area outside your study area is the missing "undef" value. John Huddleston ________________________________ From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On Behalf Of Rafanoharana Serge Claudio [rafanoharana at gmail.com] Sent: Monday, September 12, 2011 1:27 AM To: GrADS Users Forum Subject: Re: [gradsusr] Crop country boundary or masking an area - still not working Dear Yaqiang, Thank you for your email. I will work on it and let you know ASAP. Best Regards, Serge On Sat, Sep 10, 2011 at 7:57 PM, Yaqiang Wang > wrote: You can use maskout function to crop your grid data. Or use basemap script to 'crop' your figure. On Sat, Sep 10, 2011 at 11:57 PM, Rafanoharana Serge Claudio > wrote: > Dear John, > > Thank you for your answer. I do have ogr2ogr and a shapefile of the study > area. Thank you anyway for the file. > In your command, you only deal with .shp files, and ogr2ogr cannot be used > within GrADS. > The problem is that I have a .nc file which contains a variable. I want to > clip this .nc file (but not a .shp file). I am not sure that ogr2ogr can > deal with .nc format (I also tried but not working). > > Thanks, > > Serge > > > On Fri, Sep 9, 2011 at 1:09 PM, Huddleston, John > > wrote: >> >> Serge >> >> >> >> Attached is your study area shapefile (Study_Area.shp). Use ogr2ogr to >> clip your data. >> >> >> >> You can get ogr2ogr from http://www.gdal.org/ogr2ogr.html or if on Linux >> use ?apt-get install gdal-bin?. >> >> >> >> Use the following command: >> >> >> >> ogr2ogr -clipsrc Study_Area.shp output.shp Your_input_data.shp >> >> >> >> where Study_Area.shp is a shapefile in the attached zip file and >> Your_input_data.shp is the shapefile you want to clip. >> >> >> >> John Huddleston, PhD >> >> Cooperative Institute for Research in the Atmosphere >> >> >> >> >> >> >> >> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] >> On Behalf Of Rafanoharana Serge Claudio >> Sent: Friday, September 09, 2011 12:59 PM >> To: GrADS Users Forum >> Subject: [gradsusr] Crop country boundary or masking an area - still not >> working >> >> >> >> Dear all, >> >> I've already posted this message before but I can't get it worked. >> If any of you have any idea to crop the study area, please let me know. >> One said that it is possible to grab a shapefile of the top part and fill >> it with black (I guess he meant re-setting the value of the variable covered >> by this shapefile to 0). If you have any experience or any tips on how to do >> that, thanks in advance for your help. >> >> Kind regards, >> >> Serge >> >> ---------- Forwarded message ---------- >> From: Rafanoharana Serge Claudio > >> Date: Fri, Aug 26, 2011 at 7:53 AM >> Subject: Crop country boundary >> To: GrADS Users Forum > >> >> >> Dear all, >> >> I would like to ask if there is a way to crop the data to a specific area. >> In GrADS, if we use 'set mpdset hires', we have the country boundary. >> However for my case I only need the study area (please refer to the attached >> file). >> >> I found that http://www.diva-gis.org/datadown provides Grid data. Can we >> use the data from this Website to crop it or...? >> >> Thank you, >> >> -- >> Serge >> >> >> -- >> Serge >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> > > > > -- > Serge > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- ************************************************* Dr. Yaqiang Wang Chinese Academy of Meteorological Sciences (CAMS) 46, Zhong-Guan-Cun South Avenue Beijing, 100081 China yaqiang.wang at gmail.com ************************************************** _______________________________________________ gradsusr mailing list gradsusr at gradsusr.org http://gradsusr.org/mailman/listinfo/gradsusr -- Serge -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/52babf1b/attachment-0003.html From jdduda at iastate.edu Mon Sep 12 12:16:59 2011 From: jdduda at iastate.edu (Jeffrey Duda) Date: Mon, 12 Sep 2011 11:16:59 -0500 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: Bernd, To help me understand, could you provide a simple example? Jeff On Mon, Sep 12, 2011 at 3:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Jeff Duda Iowa State University Meteorology Graduate Student www.meteor.iastate.edu/~jdduda -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/59590758/attachment-0003.html From rafanoharana at gmail.com Mon Sep 12 12:29:21 2011 From: rafanoharana at gmail.com (Rafanoharana Serge Claudio) Date: Mon, 12 Sep 2011 09:29:21 -0700 Subject: [gradsusr] Crop country boundary or masking an area - still not working In-Reply-To: References:

Message-ID: On Mon, Sep 12, 2011 at 10:03 AM, Becker, Bernd < bernd.becker at metoffice.gov.uk> wrote: > ** > Jeff, > > That is is the really hard part. > Each item fx and obs describes a chunk of a time series. > One would have to manually calculate the start and end dates > with respect to the reference dfile. > > If only (and I think Mike explained it a long time ago, and I cannot find > the > email) one could > set file 1 > set t 2 34 > define whatever > set dfile 2 > set t 6 39 > define theother > > and happily combine whatever with theother. > > Think about multimodel ensemble hindcast analysis, trying to composit > observed phenomena with modelled phenomena that simply > don't happen simultaneously. > > Well, have you thought about creating a ctl with each member corresponding to a different file? In grads ensemble implementation you can specify different beginning/ending times for each member, with TDEF specifying the time interval that covers all the members. I typically do this to handle twice daily forecasts: each forecast is an ensemble member, each starting at the initial condition and going out 5 days. Here is how it looks on our GDS: http://portal.nccs.nasa.gov/cgi-lats4d/opendap.cgi?&path=/GEOS-5/fp/0.25_deg/fcast/inst1_2d_hwl_Nx.all.info Of course, this assumes that the files have the same variables and grids, just different time span. Arlindo > All the best, > Bernd. > > > > -- > Bernd Becker Climate Impacts Product Development > Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United > Kingdom > Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 > E-mail:bernd.becker at metoffice.gov.uk - > http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker > http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game > > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [mailto: > gradsusr-bounces at gradsusr.org] *On Behalf Of *Jeffrey Duda > *Sent:* 10 September 2011 18:41 > > *To:* GrADS Users Forum > *Subject:* Re: [gradsusr] define variable independent of time > > You can reference the time of a variable when you use it in a statement. > For example, you could do > > 'd var.2(time=12Z01JAN1996) - var.3(time=18Z13AUG2010)' > > This should help you do what you want to do. > > Jeff Duda > > On Sat, Sep 10, 2011 at 3:30 AM, Becker, Bernd < > bernd.becker at metoffice.gov.uk> wrote: > >> ** >> Thanks Arlindo, >> >> That would be nice for single numbers. >> >> It as one notch more elaborate: >> >> 'set dfile 2' >> 'set t 2 33' >> 'define obs = 'var'.2' >> 'set dfile 3' >> 'set t 2 33' >> 'define fx = 'var'.3' >> file 2 and 3 have different time frames. >> >> d fx-obs >> >> fails. >> And it would be a real bummer to have to calculate the times from one file >> relative >> to the other file. >> >> Compositing in this way would be really hard. >> >> Enjoy the weekend! >> Bernd. >> >> >> -- >> Bernd Becker Climate Impacts Product Development >> Met Office Hadley Centre FitzRoy Road Exeter Devon EX1 3PB United >> Kingdom >> Tel.: +44 (0) 1392 884511 Fax: +44 (0)870 900 5050 >> E-mail:bernd.becker at metoffice.gov.uk - >> http://www.metoffice.gov.uk/research/our-scientists/seasonal-to-decadal/bernd-becker >> http://www.metoffice.gov.uk/news/releases/archive/2011/weather-game >> >> >> ------------------------------ >> *From:* gradsusr-bounces at gradsusr.org [mailto: >> gradsusr-bounces at gradsusr.org] *On Behalf Of *Arlindo da Silva >> *Sent:* 09 September 2011 17:40 >> *To:* GrADS Users Forum >> *Subject:* Re: [gradsusr] define variable independent of time >> >> On Fri, Sep 9, 2011 at 12:05 PM, Becker, Bernd < >> bernd.becker at metoffice.gov.uk> wrote: >> >>> >>> Simple question, >>> complicated answer: >>> >>> I would like to composite stuff from several files. >>> >>> Each file has its own time and date info in its own .ctl file. >>> >>> How can I add elements from different such files? >>> >>> How can I break away from the strict treatment of the time dimension in >>> grads? >>> >>> open file1.ctl >>> define a=ave(b,t=1,t=34) >>> close 1 >>> open file2.ctl >>> define c= ave(b,t=3,t=9) >>> close 2 >>> etc. >>> >>> d a+c >>> Cannot contour grid - all undefined values >>> >>> >> Try: >> >> d a(t=1)+c(t=1) >> >> >> -- >> Arlindo da Silva >> dasilva at alum.mit.edu >> >> _______________________________________________ >> gradsusr mailing list >> gradsusr at gradsusr.org >> http://gradsusr.org/mailman/listinfo/gradsusr >> >> > > > -- > Jeff Duda > Iowa State University > Meteorology Graduate Student > www.meteor.iastate.edu/~jdduda > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -- Arlindo da Silva dasilva at alum.mit.edu -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110912/5ce7c317/attachment-0003.html From jma at cola.iges.org Mon Sep 12 13:57:19 2011 From: jma at cola.iges.org (Jennifer Adams) Date: Mon, 12 Sep 2011 13:57:19 -0400 Subject: [gradsusr] define variable independent of time In-Reply-To: References:

Message-ID: I work with a Debian oon Linux. 2011/9/14 Huddleston, John > CL > > What is your operating system? > > John Huddleston > > ------------------------------ > *From:* gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] On > Behalf Of Christelle LACROIX [cl at meteoconsult.fr] > *Sent:* Tuesday, September 13, 2011 4:42 AM > *To:* GrADS Users Forum > *Subject:* [gradsusr] problem with installation of grads > > Hi, > > > I have a problem with grads, when i run a script i have this error: > "/opt/grads-2.0.a9/bin/grads: error while loading shared libraries: > libX11.so.6: cannot open shared object file: No such file or directory" > > whereas libX11.so.6 exists, when i do a locate libX11 i have this result: > /usr/lib/libX11.so.6 > /usr/lib/libX11.so.6.3.0 > > I don't understand this problem, if you have an idea, it'll be great!! > > Thanks > > CL > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > -------------- next part -------------- An HTML attachment was scrubbed... URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110914/cca92e02/attachment-0003.html From smcmillan at planalytics.com Wed Sep 14 09:10:22 2011 From: smcmillan at planalytics.com (Stephen McMillan) Date: Wed, 14 Sep 2011 09:10:22 -0400 Subject: [gradsusr] Problem accessing NCEP (NOMADS) GrADS-DODS server (from a web server) In-Reply-To: References: Message-ID: Joe, I realize this is 12 hours later, but I just now successfully opened the same dataset. Here's the command output: Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3 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.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff Issue 'q config' command for more information. Loading User Defined Extensions table ... ok. GX Package Initialization: Size = 11 8.5 wrote: > All, > > I'm having a problem using the GrADS Data Server (GrADS-DODS) from NCEP > (NOMADS). I am using the sdfopen command to access the dataset (RTMA 2.5km). > When I run the script locally on my Ubuntu machine, everything works great! > However, when I transfered the script to my web host I get an error. > > Here's the line from my script: > 'sdfopen > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p5'dapdate'/rtma2p5_anl_'currhr'z > ' > (where dapdate=YYYYMMDD and currhr=HH - yes, I've double checked to make > sure the data exists.) > > Here's the output: > ... > Running in Batch mode > No hardcopy metafile open > All files closed; all defined objects released; > All GrADS attributes have been reinitialized > Error: nc_open failed to open file > http://nomads.ncep.noaa.gov:9090/dods/rtma2p5/rtma2p520110913/rtma2p5_anl_23z > NetCDF: I/O failure > gadsdf: Couldn't ingest SDF metadata. > ... > > I've tried searching the web and this (excellent) mailing list assistance > but I can't seem to figure out what's going wrong. I suspect it may have > something to do with my web host. (I use Bluehost- I know they're a somewhat > subpar host, but I'm just a student experimenting with different things!) > > Here's what uname -a gives on my web server: > Linux box566.bluehost.com 2.6.32-46.1.BHsmp #1 SMP Tue Sep 6 12:18:02 MDT > 2011 x86_64 x86_64 x86_64 GNU/Linux > > If there's other tests I can run to probe the issue, let me know. Also, is > this the best way to use the GrADS-DODS server? Should I just download the > grib2 locally instead? > > I was planning to download the grib2 locally, but when I discovered the > GrADS-DODS server I found out that this might not be necessary. I plan to > run this script hourly, so please advise as to which method would be best. > > Thanks a lot! > > -Joe Moore > > The script in question: > http://pastebin.com/MDQ9vk3y > > My GrADS config on my web server: > Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds > hdf5 opendap-grids,stn athena geotiff shapefile > Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1 > Copyright (c) 1988-2010 by Brian Doty and the > Institute for Global Environment and Society (IGES) > This program is distributed WITHOUT ANY WARRANTY > See file COPYRIGHT for more information. > > Built Mon Oct 4 14:16:26 EDT 2010 for i686-pc-linux-gnu > > This version of GrADS has been configured with the following options: > o Built on a LITTLE ENDIAN machine > o Athena Widget GUI ENABLED > o Command line editing ENABLED > http://tiswww.case.edu/php/chet/readline/rltop.html > o printim command for image output ENABLED > http://www.zlib.net > http://www.libpng.org/pub/png/libpng.html > http://www.libgd.org/Main_Page > o GRIB2 interface ENABLED > http://www.ijg.org > http://www.ece.uvic.ca/~mdadams/jasper > http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 > g2clib-1.1.8 > o NetCDF interface ENABLED > http://www.unidata.ucar.edu/software/netcdf > netcdf 4.1.1 of Jul 10 2010 23:36:25 $ > o OPeNDAP gridded data interface ENABLED > o OPeNDAP station data interface ENABLED > http://iges.org/grads/gadoc/supplibs.html > libgadap 2.0.oga.1 > o HDF4 and HDF5 interfaces ENABLED > http://hdfgroup.org > HDF 4.2r3 > HDF5 1.8.4 > o GeoTIFF and KML/TIFF output ENABLED > http://www.libtiff.org > http://geotiff.osgeo.org > o KML contour output ENABLED > o Shapefile interface ENABLED > http://shapelib.maptools.org > > > _______________________________________________ > gradsusr mailing list > gradsusr at gradsusr.org > http://gradsusr.org/mailman/listinfo/gradsusr > > *************************************************** 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/20110914/5c694fbd/attachment-0003.html From Huddleston at cira.colostate.edu Wed Sep 14 09:37:46 2011 From: Huddleston at cira.colostate.edu (Huddleston, John) Date: Wed, 14 Sep 2011 06:37:46 -0700 Subject: [gradsusr] problem with installation of grads In-Reply-To: References: