[gradsusr] Grads earth's radius

Jennifer Adams jma at cola.iges.org
Tue Dec 13 09:38:52 EST 2011


The values (there are more than one) of the Earth radius are hard-coded into GrADS in the PDEF routines that were provided by NCEP to deal with their pre-projected data. Changing them would break the display of just about every regional product that comes out of NCEP. If you require a different value, then what is needed is a new PDEF keyword and a new routine to deal with it in gaddes.c. You can take a look at the source code and submit a patch if you want, but it might be simpler to do the calculations yourself and create your own PDEF file. Or, you could change the values in your own copy of the GrADS source and then create a local build with the values that suit the data you want to display. 
--Jennifer


On Dec 12, 2011, at 3:06 PM, Jofre Janué wrote:

>  
> I have some problems with the borders of countries and the representation of the WRF model, because don't match( you can see it with 3km resolution). I think theWRF works with a earth's radius of 6370 km and Grads is 6371.2 km.
> Is it possible to change the Grads radius of the earth to match borders with WRF?
>  
>  
>  
>  
>  
>  
>  
>  
> 
>  
> > From: gradsusr-request at gradsusr.org
> > Subject: gradsusr Digest, Vol 22, Issue 11
> > To: gradsusr at gradsusr.org
> > Date: Sun, 11 Dec 2011 18:14:57 -0500
> > 
> > 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: OpenGrads Installation (Arlindo da Silva)
> > 2. Re: OpenGrads Installation (Andrew Revering)
> > 
> > 
> > ----------------------------------------------------------------------
> > 
> > Message: 1
> > Date: Sun, 11 Dec 2011 13:47:27 -0800
> > From: Arlindo da Silva <dasilva at alum.mit.edu>
> > Subject: Re: [gradsusr] OpenGrads Installation
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CALz7xZd3=aUQw0-jEYOZJ5BnBAAGL+H+56ebCwUPS49ve=uMXA at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> > 
> > On Sun, Dec 11, 2011 at 12:44 PM, Andrew Revering <andy at f5data.com> wrote:
> > 
> > > Arlindo,
> > >
> > > I'm sending this to you privately
> > >
> > 
> > Not so, you replied to the whole list...
> > 
> > 
> > 
> > > because I dont want to accidentally send some server-sensitive information
> > > to the entire list accidentally...
> > >
> > >>
> > >> Simply type "env" at the command line and you will get a list of all
> > >> environment variables defined.
> > >>
> > >
> > > root at server [~]# env
> > > HOSTNAME=*********************
> > > /env
> > >
> > ...
> > 
> > Good, no grads related environment variable is being set, so we can rule
> > out any conflict there.
> > 
> > 
> > >
> > >> Make sure you move the Contents/ directory as a whole, keeping directory
> > >> tree intact.
> > >>
> > >
> > > So for example, should it be:
> > > /opt/opengrads/Linux.... etc.
> > > or
> > > /opt/opengrads/Contents/Linux.... etc.
> > >
> > > I currently have it as the first scenario.
> > >
> > >
> > 
> > Either is fine, provided the whole Contets/ directory (and all
> > sub-directories are moved over). In your case, make sure /opt/opengrads is
> > ahead in your path (and not another directory further down under Linux/).
> > 
> > 
> > 
> > >
> > > I don't think I'm able to do this because the server is a web server and
> > > I'm connecting through SSH.
> > >
> > 
> > Sometimes you need the -X option to enable X windows through ssh, for
> > example
> > 
> > % ssh -X myserver.domain.com
> > 
> > 
> > 
> > > So when I run opengrads I get "error in GXSTRT: unable to connect to X
> > > server"
> > >
> > > So I have to run everything with -lbc 'run ' scripts. Is there a way to do
> > > q udx and send that output to a text file?
> > >
> > >
> > You can run grads interatively without having X Windows connecivity. Just
> > start opengrads in batch mode (-b option) :
> > 
> > % opengrads -bl
> > 
> > What I am looking for is whether you see a line like this:
> > 
> > *Loading User Defined Extensions table
> > </Applications/OpenGrADS/Darwin/Versions/2.0.0.oga.1/i386/gex/udxt> ... ok.*
> > 
> > and whether
> > 
> > ga-> q udx
> > 
> > gives you a full list of extensions.
> > 
> > Arlindo
> > 
> > 
> > -- 
> > Arlindo da Silva
> > dasilva at alum.mit.edu
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111211/438e4d22/attachment-0001.html 
> > 
> > ------------------------------
> > 
> > Message: 2
> > Date: Sun, 11 Dec 2011 17:07:41 -0600
> > From: Andrew Revering <andy at f5data.com>
> > Subject: Re: [gradsusr] OpenGrads Installation
> > To: GrADS Users Forum <gradsusr at gradsusr.org>
> > Message-ID:
> > <CAHXC6uvYwyniLHKc7efqOSxA2+k1pzgn60KucTSLixAcoiDnHw at mail.gmail.com>
> > Content-Type: text/plain; charset="iso-8859-1"
> > 
> > >
> > > Not so, you replied to the whole list...
> > >>
> > >
> > Yah... I'm a damn fool.
> > 
> > 
> > > ga-> q udx
> > >
> > > gives you a full list of extensions.
> > >
> > 
> > ga-> q udx
> > 
> > Command Environment Variables Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > printenv Expand environment variables c_xenv@^env.gex
> > runenv Expand env vars and run command c_env@^env.gex
> > @ Expand env vars and run command c_env@^env.gex
> > getenv Get value of environment variable c_getenv@^env.gex
> > setenv Set value of environment variable c_setenv@^env.gex
> > 
> > Command GxYAT: Image Generation Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > gxyat Save images in PNG/SVG/PDF/PS c_gxyat@^gxyat.gex
> > set_rgba Set color red/green/blue/alpha c_rgba@^gxyat.gex
> > 
> > Command Simple Extension Example Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > hello Hello, World! sample command c_hello@^libhello.gex
> > 
> > Command Interprocess Communication Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > ipc_verb IPC verbose toggle c_Verb@^libipc.gex
> > ipc_open Open stream for save/load c_Open@^libipc.gex
> > ipc_close Close stream c_Close@^libipc.gex
> > ipc_save Save expression to stream c_Save@^libipc.gex
> > ipc_define Define variable (obsolete) c_Define@^libipc.gex
> > ipc_error Print IPC error message c_Error@^libipc.gex
> > 
> > Command LATS: File Reformatting and Subsetting Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > set_lats Set LATS parameters c_lats_set@
> > ^liblats.gex
> > query_lats Query LATS parameters c_lats_query@
> > ^liblats.gex
> > lats_grid Define LATS grid c_lats_grid@
> > ^liblats.gex
> > lats_data Write data to LATS file c_lats_data@
> > ^liblats.gex
> > reinit_lats Re-init LATS parameters c_lats_reinit@
> > ^liblats.gex
> > 
> > Command Mike Fiorino's Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > mfhilo Find max/min or H/L in 2D field c_mfhilo@^libmf.gex
> > tcprop Properties radinf from central lon/lat c_tcprop@^libmf.gex
> > tcxy2rt cartesian -> cylindrical coords results only scriptvar RESULT
> > c_tcxy2rt@^libmf.gex
> > grhist histogram of expr c_grhist@^libmf.gex
> > 
> > Command Satellite Orbit Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > set_orb Set ORB parameters c_set@^orb.gex
> > orb_track Satellite ground tracks c_track@^orb.gex
> > 
> > Command Shapefile Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > shp_lines Draw lines from shapefile c_lines@^shape.gex
> > shp_polyf Draw polygons from shapefile c_polyf@^shape.gex
> > 
> > Function Ben-Jei Tsuang's Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > lt Local time f_bjt@^libbjt.gex
> > jd Julian day f_bjt@^libbjt.gex
> > cosz Cosine solar zenith angle f_bjt@^libbjt.gex
> > dayratio Daylight ratio f_bjt@^libbjt.gex
> > if Conditional function f_bjt@^libbjt.gex
> > maxv Maximum value f_bjt@^libbjt.gex
> > minv Minimum value f_bjt@^libbjt.gex
> > which Label gridpoints f_bjt@^libbjt.gex
> > ftest F-test f_bjt@^libbjt.gex
> > ttest T-test f_bjt@^libbjt.gex
> > tfit Point linear regression f_bjt@^libbjt.gex
> > fit Global linear regression f_bjt@^libbjt.gex
> > tcorr2 Time correlation f_bjt@^libbjt.gex
> > tregr2 Point linear regression f_bjt@^libbjt.gex
> > tmave2 Time averaging w/masking f_bjt@^libbjt.gex
> > madvu Calculates -d(u*EXPR)/dx f_bjt@^libbjt.gex
> > madvv Calculates -d(V*EXPR)/dy f_bjt@^libbjt.gex
> > madvw Calculates -d(W*EXPR)/dp f_bjt@^libbjt.gex
> > muadv Zonal advection f_bjt@^libbjt.gex
> > mvadv Meridional advection f_bjt@^libbjt.gex
> > mwadv Vertical advection f_bjt@^libbjt.gex
> > satvap Saturated vapor pressure f_bjt@^libbjt.gex
> > dew Dew point temperature f_bjt@^libbjt.gex
> > lw Thermal infrared fluxes f_bjt@^libbjt.gex
> > lw2 Thermal infrared fluxes v2 f_bjt@^libbjt.gex
> > pinterp Pressure interpolation f_bjt@^libbjt.gex
> > zinterp Height interpolation f_bjt@^libbjt.gex
> > line Draws a line f_bjt@^libbjt.gex
> > vint2 Mass-weighted vertical integral f_bjt@^libbjt.gex
> > 
> > Function Streamfunction & Velocity Potential Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > fish Poisson solver f_fish@^fish.gex
> > fish_psi Compute streamfunction f_psichi@^fish.gex
> > fish_chi Compute velocity potential f_psichi@^fish.gex
> > fish_vor Compute vorticity f_psichi@^fish.gex
> > fish_div Compute divergence f_psichi@^fish.gex
> > 
> > Function Saakeskus Wx Function Collection Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > dewpt Dew-point temperature f_gsudf@^gsudf.gex
> > ept Equivalent potential temperature f_gsudf@^gsudf.gex
> > epi Potential Instability f_gsudf@^gsudf.gex
> > plcl Pressure at LCL f_gsudf@^gsudf.gex
> > shear Wind Shear f_gsudf@^gsudf.gex
> > srh1km Storm-relative Helicity 0-1km' f_gsudf@^gsudf.gex
> > srh3km Storm-relative Helicity 0-3km' f_gsudf@^gsudf.gex
> > tlcl Temperature at LCL f_gsudf@^gsudf.gex
> > tsindex Chance of T-storm Initiation f_gsudf@^gsudf.gex
> > ttindex Total Totals Index f_gsudf@^gsudf.gex
> > wchill Wind Chill Index (New) f_gsudf@^gsudf.gex
> > 
> > Function Simple Extension Example Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > hello Hello, World! sample function f_hello@^libhello.gex
> > 
> > Function Interprocess Communication Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > ipc_save Save expression to stream f_Save@^libipc.gex
> > ipc_load Load variable from file f_Load@^libipc.gex
> > 
> > Function Mike Fiorino's Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > smth2d Shuman smoother/de-smoother f_smth2d@^libmf.gex
> > xy2rt cartesian -> cylindrical coord f_xy2rt@^libmf.gex
> > uv2trw Find radial/tangential velocity f_uv2trw@^libmf.gex
> > re2 General interpolator (regrid2) f_re2@^libmf.gex
> > esmrf Saturation vapor pressure (old MRF) f_esmrf@^libmf.gex
> > linreg Linear regression: mx + b f_linreg@^libmf.gex
> > 
> > Function Satellite Orbit Extensions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > orb_mask Satellite orbit masking f_mask@^orb.gex
> > 
> > Function Regridding Functions Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > re General interpolator ffre@^re.gex
> > re_ General interpolator (verbose) ffre_dbg@^re.gex
> > regrid2 regrid2-like wrapper for re() f_gsudf@^gsudf.gex
> > geos Interpolation to GEOS-5 Grids f_gsudf@^gsudf.gex
> > reimg Interpolation for Image Generation' f_gsudf@^gsudf.gex
> > 
> > Function Spherical Harmonic Analysis Function at Library
> > ----------- ---------------------------------------
> > --------------------------
> > sh_filt Spherical harmonic filter f_shfilt@^shfilt.gex
> > sh_power Spherical harmonic spectra f_shpowr@^shfilt.gex
> > 
> > ga->
> > -------------- next part --------------
> > An HTML attachment was scrubbed...
> > URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111211/03144cd9/attachment.html 
> > 
> > ------------------------------
> > 
> > _______________________________________________
> > gradsusr mailing list
> > gradsusr at gradsusr.org
> > http://gradsusr.org/mailman/listinfo/gradsusr
> > 
> > 
> > End of gradsusr Digest, Vol 22, Issue 11
> > ****************************************
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111213/c4d558bf/attachment-0003.html 


More information about the gradsusr mailing list