Context Navigation

zonalmean.py @ 1014

Last change on this file since 1014 was 1002, checked in by aslmd, 12 years ago
UTIL PYTHON planetoplot_v2. possibility to load a field in a restricted interval with e.g. -t 0,1 without performing computations (self.compute=nothing). use for histograms, see histo.py. added smooth2diter. bug correction: ensure numpy array, vector, minmax equal.
Property svn:executable set to ``*
File size: 1.1 KB

Rev	Line
[910]	1	#! /usr/bin/env python
	2	from ppclass import pp
	3
	4	## AVERAGED PROFILE
	5	temp = pp()
	6	temp.file = "/home/aymeric/Big_Data/DATAPLOT/diagfired.nc"
	7	temp.var = "temp"
[1002]	8	temp.x = "-180,175"
[910]	9	temp.y = "-90,90"
	10	temp.t = "0,1"
	11	temp.get()
	12	temp.defineplot()
	13	temp.p[0].title = "This is an averaged temperature profile"
	14	temp.makeplot()
	15
	16	## ZONAL MEAN
	17	u = pp()
	18	u.file = "/home/aymeric/Big_Data/DATAPLOT/diagfired.nc"
	19	u.var = "u"
[1002]	20	u.x = "-180,175"
[910]	21	u.y = None
	22	u.t = "0.5"
[923]	23	u.filename = "zonalmean"
[910]	24	u.get()
	25	u.defineplot()
	26	u.p[0].div = 30.
	27	u.p[0].colorb = "RdBu_r"
	28	u.p[0].title = "This is a zonal mean"
	29	u.makeplot()
	30
	31	## ZONAL MINIMUM
	32	u.compute = "min"
[923]	33	u.filename = "zonalmin"
[910]	34	u.get()
	35	u.defineplot()
	36	u.p[0].div = 30.
	37	u.p[0].colorb = "cool"
	38	u.p[0].title = "This is minimum over zonal axis"
	39	u.makeplot()
	40
	41	## ZONAL MAXIMUM
	42	u.compute = "max"
[923]	43	u.filename = "myplot"
[910]	44	u.get()
	45	u.defineplot()
	46	u.p[0].div = 30.
	47	u.p[0].colorb = "hot"
	48	u.p[0].title = "This is maximum over zonal axis"
	49	u.makeplot()
	50
	51	## MAP OF MAXIMUM OVER TIME
	52	u.compute = "max"
	53	u.x = None
	54	u.y = None
	55	u.t = "0,1"
	56	u.z = 20000.
	57	u.get()
	58	u.defineplot()
	59	u.p[0].div = 30.
	60	u.p[0].title = "This is maximum over time"
	61	u.makeplot()

Note: See TracBrowser for help on using the repository browser.