[2] | 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 | ; |
---|
[48] | 27 | IF size(glam, /n_dimensions) EQ 1 THEN triangles = triangule() ELSE triangulate, glamt, gphit, triangles |
---|
[2] | 28 | ; |
---|
[48] | 29 | IF keyword_set(land) THEN BEGIN |
---|
[2] | 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 |
---|
[48] | 33 | ENDIF ELSE BEGIN |
---|
[2] | 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 |
---|
[48] | 37 | ENDELSE |
---|
[2] | 38 | return |
---|
| 39 | end |
---|