How to make composite plots

Charles Seman Charles.Seman at NOAA.GOV
Fri Sep 25 18:59:17 EDT 2009


Sushant,

Please see http://grads.iges.org/grads/gadoc/gradcomdmodify.html which
shows how to make a climatology, and then if desired, to modify it for
calculating anomalies...

I hope this helps,
Chuck

sushant puranik wrote:
> Thank you for hint.
> I tried
> d (ave(var,t=1,t=360,12)+ave(var,t=2,t=360,12)+ave(var,t=3,t=360,12))/3
> after that i checked the result it gives average of three days only
> not for three months.
>
> Since i has ECMWF data at 6 hour interval in the monthly form (each
> file contain 1 month data).
> So i tried as per your suggestion but it gives ave of first three days
> of first file only and not for the three years file.
>
> Also i am facing the following problem
>
> [Sushant at localhost ~]$ gradsnc
>
> Grid Analysis and Display System (GrADS) Version 1.8SL11
> Copyright (c) 1988-2001 by Brian Doty
> Center for Ocean-Land-Atmosphere Studies
> Institute for Global Environment and Society
> All Rights Reserved
>
> Config: v1.8SL11 32-bit little-endian readline sdf/xdf netcdf lats printim
>
> Issue 'q config' command for more information.
>
> Landscape mode? (no for portrait):
> GX Package Initialization: Size = 11 8.5
> ga-> c
> *ga-> sdfopen /home/Sushant/Desktop/ECMWFdata/May2002data/u.nc
> <http://u.nc> *
> Scanning self-describing file:
> /home/Sushant/Desktop/ECMWFdata/May2002data/u.nc <http://u.nc>
> Found 1 data variables in SDF file.
> Found displayable variable u with 29 levels in SDF file.
> SDF file /home/Sushant/Desktop/ECMWFdata/May2002data/u.nc
> <http://u.nc> is open as file 1
> LON set to 0 360
> LAT set to -90 90
> LEV set to 50 50
> Time values set: 2002:5:1:0 2002:5:1:0
> *ga-> sdfopen /home/Sushant/Desktop/ECMWFdata/May2001data/u.nc
> <http://u.nc> *
> Scanning self-describing file:
> /home/Sushant/Desktop/ECMWFdata/May2001data/u.nc <http://u.nc>
> Found 1 data variables in SDF file.
> Found displayable variable u with 29 levels in SDF file.
> SDF file /home/Sushant/Desktop/ECMWFdata/May2001data/u.nc
> <http://u.nc> is open as file 2
> ga-> set lev 1000
> LEV set to 1000 1000
> ga-> d u
> Contouring: -20 to 20 interval 5
> ga-> c
> ga-> d u.2
> *Data Request Warning: Request beyond file limits
> Entire grid contents are set to missing data
> Dimension ranges are: X = 1 240 Y = 1 121 Z = 29 29 T = 1461 124
> Warning issued for variable = u.2
> Cannot contour grid - all undefined values *
> ga-> c
>
> Now can you please solve this problem.
>
>
>
> Thank you again
>
>
> Sushant
>
>
> On Tue, Sep 22, 2009 at 5:50 PM, Yan Libin <zjuyanlibin at 163.com
> <mailto:zjuyanlibin at 163.com>> wrote:
>
>     d ave(var,t=1,t=360,12)
>     d
>     (ave(var,t=1,t=360,12)+ave(var,t=2,t=360,12)+ave(var,t=3,t=360,12))/3
>     2009-09-22
>     ------------------------------------------------------------------------
>     Yan Libin
>     ------------------------------------------------------------------------
>     *发件人:* sushant puranik
>     *发送时间:* 2009-09-22 20:18:04
>     *收件人:* GRADSUSR
>     *抄送:*
>     *主题:* How to make composite plots
>     Hello sir
>     I wish to make composite plot for 20 year period for each month
>     (monthly mean composite) and also want to make composite for
>     season (seasonal composite)(Dec-March).
>     Can you please give me help in this regard.
>
>
>     Thanking you
>
>     Sushant
>
>

--

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
 Charles Seman                                Charles.Seman at noaa.gov
 U.S. Department of Commerce / NOAA / OAR
 Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
 201 Forrestal Road                              fax: (609) 987-5063
 Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not necessarily
reflect any position of the Government or NOAA."



More information about the gradsusr mailing list