source: trunk/draw.pro

Last change on this file was 48, checked in by pinsard, 10 years ago

fix thanks to coding rules

File size: 1.4 KB
Line 
1PRO draw, field, glam, gphi, zpi, zpj, nbcol, dessin, nbdessin, title = title, zmin, zmax, LAND = land
2@common
3jpi = zpi
4jpj = zpj
5glamt = glam
6IF size(glamt, /n_dimensions) EQ 1 THEN glamt = glamt#replicate(1, jpj)
7glamf = glam
8gphit = gphi
9IF size(gphit, /n_dimensions) EQ 1 THEN gphit = replicate(1, jpi)#gphit
10gphif = gphi
11gdept = [0]
12gdepw = [0]
13tmask = replicate(1, jpi, jpj)
14
15e1t = replicate(1, jpi, jpj)
16e2t = replicate(1, jpi, jpj)
17e3t = [1]
18jpk = 1
19;
20varname = title
21domdef, grille = ['T', 'F']
22key_periodique = 1
23;
24ind = where(field GE 1.e20)
25IF n_elements(ind) GT 1 THEN field(ind) = 0.
26;
27IF size(glam, /n_dimensions) EQ 1 THEN triangles = triangule() ELSE triangulate, glamt, gphit, triangles
28;
29IF keyword_set(land) THEN BEGIN
30   if dessin EQ 1 then plt, field, zmin, zmax, nite = 0, /carte, petit = [nbcol, nbdessin, dessin], /remplit, div = 4, /cont_nofill, retain = 2, /nocontour, /land $
31   ELSE plt, field, zmin, zmax, nite = 0, /carte, petit = [nbcol, nbdessin, dessin], $
32    /noerase, /remplit, div = 4, /cont_nofill, /nocontour, /land
33ENDIF ELSE BEGIN
34if dessin EQ 1 then plt, field, zmin, zmax, nite = 0, /carte, petit = [nbcol, nbdessin, dessin], /remplit, div = 4, /cont_nofill, retain = 2, /nocontour, /port $
35ELSE plt, field, zmin, zmax, nite = 0, /carte, petit = [nbcol, nbdessin, dessin], $
36 /noerase, /remplit, div = 4, /cont_nofill, /nocontour, /port
37ENDELSE
38return
39end
Note: See TracBrowser for help on using the repository browser.