plot hovmuller

Categories:

PNG hoevmuller AMSU OLR

Params:
dataset : in, required, type='string'

dataset to plot

canal : in, optional, type='string'

canal of AMSU

only used if AMSU dataset

iyear : in, required, type='integer'

year

Keywords:
DATASETTYPE : default='for AMSU ++, for OLR ++'

dataset type name. For each dataset, there are several file type (year vs total, global vs regional). ++

AXE : default='latitude'

geographical axe to be use for hovmuller can be either 'latitude' or 'longitude'

OVERWRITE : default=false

to overwrite the output image file if exists.

Useless if NOSAVE is set.

NOSAVE : default=false

to NOT save the output image in a file

might be usefull if several plots are needed on the same figure.

file if exist might be overwritten if OVERWRITE is set to 1.

_EXTRA :

Used to declare that this routine accepts inherited keywords

see example with SMALL keyword

Restrictions:
  • Requires SAXO
Returns:
Uses:

pltt

cm_4cal

forout.pro

Pre:

varamma_profile.sh

for AMSU dataset be sure to have cx.anyyyy.nc in the directory defined in ${VARAMMA_ID}/

for OLR dataset be sure to have olr.day.mean_yyyy_ng.nc in the directory defined in ${VARAMMA_ID}/

Post:

varamma_profile.sh

Todo:

++ parametrisation date début date fin

++ exploitation du champ info des fichiers OLR

++ gros souci avec les longitude OLR !!!!

++ split read and plot

++ saveimage or openps (for image to provide to publisher)

++ filename si small idiot car peu contenir plusieurs figures donc pas canal+year...

Examples:

To plot AMSU a4 hovmuller latitude plot for 2006 and produce a PNG file:

IDL> hovmuller, 'AMSU', 'a4', 2006L, AXE='latitude'

To plot AMSU a4 and a5 hovmuller latitude plots for 2006 on the same figure:

IDL> hovmuller, 'AMSU', 'a4', 2006L, AXE='latitude', /nosave,small=[2,1,1]
IDL> hovmuller, 'AMSU', 'a5', 2006L, AXE='latitude', /noerase,/nocolorbar,small=[2,1,2]

++ les deux dessins n'occupent pas la même taille.

To plot OLR hovmuller latitude plot for 2006 and produce a PNG file:

IDL> hovmuller, 'OLR', 'n.a.', 2006L, AXE='latitude'

To plot OLR hovmuller longitude plot for 2006 and produce a PNG file:

IDL> hovmuller, 'OLR', 'n.a', 2006L, AXE='longitude'

To plot AMSU a4 hovmuller latitude plot climatology and produce a PNG file:

IDL> hovmuller, 'AMSU', 'a4', 0000L, DATASETTYPE='varamma_t3',AXE='latitude'
History:
hovmuller latitude
canal='a7'
mois=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec']
mo_lon=[31,28,31,30,31,30,31,31,30,31,30,31]

imo=4 &im1=10  ; no du mois
file='stagiaire/'+canal+'.an2006.nc'
initncdf, file, xaxisname = 'xlon', yaxisname = 'ylat'
domdef,-10,5,-30,45  ;domdef,0,5,-30,45
j1=total(mo_lon(0:imo-2)) & J2=total(mo_lon(0:im1-1))
data=read_ncdf('moyenne_tb', j1,j2,/timestep,timevar = 'jours', file = file)
time=julday(imo,1,2006)+lindgen(jpt)  ; redefinition axe temps
pltt,window=0, data, 'yt',title='latitude - time '+canal;,min=250,max=295
saveimage, 'sorties/'+canal+'20060410-hov-10-5.png',/png ;capture d'ecran
Version:$Id: hovmuller.pro 2 2008-12-18 16:56:52Z pinsard $