[gradsusr] gradsusr Digest, Vol 43, Issue 18
Dieni Fitriani
dieniakman at gmail.com
Thu Sep 5 06:35:03 EDT 2013
i am facing problem using grads. i tried to display the result but
everytime i display it in a different time it gives me a different picture
too. can somebody help me solve this problem. my grads was set offline.
dieni
2013/9/5 <gradsusr-request at gradsusr.org>
> Send gradsusr mailing list submissions to
> gradsusr at gradsusr.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://gradsusr.org/mailman/listinfo/gradsusr
> or, via email, send a message with subject or body 'help' to
> gradsusr-request at gradsusr.org
>
> You can reach the person managing the list at
> gradsusr-owner at gradsusr.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of gradsusr digest..."
>
>
> Today's Topics:
>
> 1. Re: lterp+sdfwrite adding extra longitude values
> (Charles Seman - NOAA Federal)
> 2. Re: Decadal average plot (Hosmay Lopez)
> 3. Problems reading a netcdf file (Matt Bunkers - WFO UNR)
> 4. Re: how to merge file data using grads
> (Muhammad Yunus Ahmad Mazuki)
> 5. Re: Problems reading a netcdf file (Jeff Duda)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 04 Sep 2013 18:02:48 -0400
> From: Charles Seman - NOAA Federal <charles.seman at noaa.gov>
> Subject: Re: [gradsusr] lterp+sdfwrite adding extra longitude values
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID: <5227AE08.8000305 at noaa.gov>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Andrew,
>
> from http://grads.iges.org/grads/gadoc/gradfunclterp.html:
> ---
> "If the domain of source is larger than the domain of dest, the returned
> result will have an expanded grid to cover the requested domain."
> ---
>
> compare the longitude limits for the fine and coarse grids... it could
> be that if the specified fine grid limits are outside the coarse
> longitude grid limits, the interpolation puts data in two extra outer
> longitudes in the destination grid to satisfy the grid environment
> defined by "set X 1 [number of longitude points in fine.nc]"... so try
> "set dfile 2 ; set X 1 [number of longitude points in coarse.nc]" after
> "sdfopen coarse.nc"... this would set the global longitude limits to
> those of coarse.nc, which would then be within the limits of fine.nc
>
> Hope this helps,
> Chuck
>
> On 09/04/2013 04:42 PM, Andrew Friedman wrote:
> > Hi GrADS users,
> >
> > I am using lterp to regrid a file to a coarser resolution, and saving
> the variable to netcdf using sdfwrite.
> >
> > My problem is that the output file has 2 additional longitudes compared
> to the input grid destination description file.
> >
> > Similar to the example in the lterp documentation in the grads index:
> >
> > sdfopen fine.nc
> > set X 1 [number of longitude points in fine.nc]
> > sdfopen coarse.nc
> > output = lterp(fine,coarse.2,aave)
> > set sdfwrite output.nc
> > sdfwrite output
> >
> > I've had this issue before when there is just one input file, and
> explicitly defining longitude with the step 'set X 1 [number of longitude
> points]' works to solve it. In this case however, I'm not sure how to limit
> the longitude dimension of the output.
> > Do you have ideas how I can fix this?
> >
> > Thanks,
> > Andrew
> >
> >
> > ------
> > Andrew Friedman
> > Graduate Student, UC Berkeley Geography
> > andfried at berkeley.edu
> >
> >
> >
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
>
> --
>
> Please note that Charles.Seman at noaa.gov should be considered my NOAA
> email address, not cjs at gfdl.noaa.gov.
>
> ********************************************************************
> Charles Seman Charles.Seman at noaa.gov
> U.S. Department of Commerce / NOAA / OAR
> Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
> 201 Forrestal Road fax: (609) 987-5063
> Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
> ********************************************************************
>
> "The contents of this message are mine personally and do not reflect any
> official or unofficial position of the United States Federal Government,
> the United States Department of Commerce, or NOAA."
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 4 Sep 2013 22:57:18 +0000
> From: Hosmay Lopez <hlopez at rsmas.miami.edu>
> Subject: Re: [gradsusr] Decadal average plot
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <
> 8A902EB84795C94C8726F11DC8C574A81D456F9E at EXCHANGE1.amazon.rsmas.miami.edu>
>
> Content-Type: text/plain; charset="iso-8859-1"
>
> OK, I see,
>
> You can make use of while loops and output the averages, that is one
> record for each decaded totaling 20 as follows. Your averages will be
> output by the fwrite command, you will need a control file to read it.
>
> ***************************************
> 'set gxout fwrite'
> 'set fwrite decade_data.dat'
> is = 1
> while(is < = 2400)
>
> 'set t 'is
> 'd tloop(ave(temperature, t+0, t+119))'
> 'c'
>
> is = is + 120
> endwhile
> 'disable fwrite'
> 'set gxout contour'
> *****************************************
> Hosmay Lopez
> Rosenstiel School of Marine and Atmospheric Science
> University of Miami
> 4600 Rickenbacker Causeway Miami, FL 33149
> Office: MSC337
> Phone: (305) 298-8006
> ________________________________
> From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on
> behalf of Mubashar Ahmad [mubashardogar at gmail.com]
> Sent: Wednesday, September 04, 2013 5:18 PM
> To: GrADS Users Forum
> Subject: Re: [gradsusr] Decadal average plot
>
>
> Hi lopez,
>
> Thanks for your input but I don't want runing mean average based on 10
> years, instead I want decadal average such that after getting averaged
> value for first 10 years (average of 1 to 120 months) it should give value
> for next 10 years (average of 121 to 240th time) and so on.
>
> Thanks,
>
> Mubashar
>
> On Sep 4, 2013 9:45 PM, "Hosmay Lopez" <hlopez at rsmas.miami.edu<mailto:
> hlopez at rsmas.miami.edu>> wrote:
> Dear Mubashar,
>
> You can use the tloop function in conjunction with the ave function as
> follows:
>
> set t 61 2340
> d tloop(ave(temperature,t-60,t+60))
>
> This will create an average for decadal timescales.
> Note, the time start at month 61 = beginning of year 6 and ends at month
> 2340 = 5 years prior to year 200 of your data.
>
> At the end you should have a monthly time series of decadal averaged
> quantity.
>
> Best,
>
> Hosmay Lopez
> Rosenstiel School of Marine and Atmospheric Science
> University of Miami
> 4600 Rickenbacker Causeway Miami, FL 33149
> Office: MSC337
> Phone: (305) 298-8006
>
> ________________________________________
> From: gradsusr-bounces at gradsusr.org<mailto:gradsusr-bounces at gradsusr.org>
> [gradsusr-bounces at gradsusr.org<mailto:gradsusr-bounces at gradsusr.org>] on
> behalf of Mubashar Ahmad [mubashardogar at gmail.com<mailto:
> mubashardogar at gmail.com>]
> Sent: Wednesday, September 04, 2013 2:12 PM
> To: gradsusr
> Subject: [gradsusr] Decadal average plot
>
> Dear Grads Users,
>
> I have monthly time series data (let say surface temperature) and
> wanted to plot it on decadal average basis. I mean first value of plot
> will be average of first 120 months and the next value of my plot will
> be average of next 120 months and so on. In total for 200 years I
> should get 20 value in plot. How can I do that? Y axis will be decadal
> average temperature values and x axis can be labelled 1, 2, 3 up 20
> instead of months or years.
>
> Thanks,
>
> Mubashar
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
> http://gradsusr.org/mailman/listinfo/gradsusr
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
> http://gradsusr.org/mailman/listinfo/gradsusr
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20130904/f8a95e3e/attachment-0001.html
>
> ------------------------------
>
> Message: 3
> Date: Wed, 4 Sep 2013 17:53:42 -0600
> From: Matt Bunkers - WFO UNR <matthew.bunkers at noaa.gov>
> Subject: [gradsusr] Problems reading a netcdf file
> To: gradsusr at gradsusr.org
> Message-ID:
> <CAB_myoNDpT_8az_2Z=nW0p0dwudrtK9e9treWy=
> EJ5g2Q8RRBQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hi,
>
> I have a netcdf file that cannot be opened with sdfopen because it doesn't
> follow the COORDS convention. After creating a descriptor file with x, y,
> and t definitions I can open the file with xdfopen and display data.
> However, when I use the "set lat" and "set lon" commands in this scenario
> the program bombs out. I therefore decided to create the full descriptor
> file (with dtype netcdf, pdef set accordingly, etc.) so I could plot my
> area of interest with the "set lat" and "set lon" commands. That all seems
> to work fine; however, my data values now are all zero. This is perplexing
> because I could see non-zero data when opening the file with the xdfopen
> command, but not when I opened the netcdf file with the open command and
> full descriptor file (I'm trying to plot forecast precipitation, by the
> way). If needed I can provide a copy of my descriptor file. Thanks for
> any advice you can offer.
>
> -- Matt Bunkers, NWS Rapid City, SD
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20130904/9c2503be/attachment-0001.html
>
> ------------------------------
>
> Message: 4
> Date: Thu, 5 Sep 2013 10:25:17 +0800
> From: Muhammad Yunus Ahmad Mazuki <ukm.yunus at gmail.com>
> Subject: Re: [gradsusr] how to merge file data using grads
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <
> CAKJBA6ZJbEeWB-iJC-bLDDHtjkJSFULt+A3GmxcsWwt7D+XDvA at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> You can either use read and write back script in grads. Or you can work
> directly on the file with 'cdo'. That if you want to rewrite the files into
> one file. You can also use template option in a controller file to read all
> 10 file as one file.
>
>
> On Wed, Sep 4, 2013 at 3:57 PM, Lyndon Mark Olaguera <
> olagueralyndonmark429 at gmail.com> wrote:
>
> > what type of data file is this? an output from WRF?
> > Can you give some details..
> >
> > Lyndon
> >
> > *Lyndon Mark P. Olaguera*
> > Assistant Instructor
> > Ateneo de Manila University
> > 426 6001
> > loc 5695
> >
> >
> > On Wed, Sep 4, 2013 at 2:48 PM, Indah Lutfiyanti <
> > indahlutfiyanti387 at yahoo.com> wrote:
> >
> >> Dear Senior,
> >>
> >> please help me, I have problem how to merge data which is separate 10
> >> file to be 1 file using grads?
> >>
> >> thank you
> >>
> >> Lutfi
> >>
> >>
> >>
> >> _______________________________________________
> >> 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
> >
> >
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20130905/c6c7ea27/attachment-0001.html
>
> ------------------------------
>
> Message: 5
> Date: Wed, 4 Sep 2013 21:45:53 -0500
> From: Jeff Duda <jeffduda319 at gmail.com>
> Subject: Re: [gradsusr] Problems reading a netcdf file
> To: GrADS Users Forum <gradsusr at gradsusr.org>
> Message-ID:
> <
> CAAig09B51hgYE+YH-qJFqyuAGmJKfB3rHotkcuTYwP2GHis0-g at mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Matt,
> I would have to guess there is some sort of error in your control file. I
> frequently view netCDF files in Grads and have had similar issues in the
> past.
>
> Jeff Duda
>
>
> On Wed, Sep 4, 2013 at 6:53 PM, Matt Bunkers - WFO UNR <
> matthew.bunkers at noaa.gov> wrote:
>
> > Hi,
> >
> > I have a netcdf file that cannot be opened with sdfopen because it
> doesn't
> > follow the COORDS convention. After creating a descriptor file with x,
> y,
> > and t definitions I can open the file with xdfopen and display data.
> > However, when I use the "set lat" and "set lon" commands in this
> scenario
> > the program bombs out. I therefore decided to create the full descriptor
> > file (with dtype netcdf, pdef set accordingly, etc.) so I could plot my
> > area of interest with the "set lat" and "set lon" commands. That all
> seems
> > to work fine; however, my data values now are all zero. This is
> perplexing
> > because I could see non-zero data when opening the file with the xdfopen
> > command, but not when I opened the netcdf file with the open command and
> > full descriptor file (I'm trying to plot forecast precipitation, by the
> > way). If needed I can provide a copy of my descriptor file. Thanks for
> > any advice you can offer.
> >
> > -- Matt Bunkers, NWS Rapid City, SD
> >
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> >
> >
>
>
> --
> Jeff Duda
> Graduate research assistant
> University of Oklahoma School of Meteorology
> Center for Analysis and Prediction of Storms
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://gradsusr.org/pipermail/gradsusr/attachments/20130904/228d2488/attachment.html
>
> ------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> End of gradsusr Digest, Vol 43, Issue 18
> ****************************************
>
--
*Dieni Fitriani Buchary Akman
**Laboratory of Meteorology
Major in Applied Meteorology
Department of Geophysics and Meteorology
Faculty of Mathematics and Natural Science
Bogor Agricultural University
*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130905/82459a16/attachment-0003.html
More information about the gradsusr
mailing list