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