Context Navigation

grille.pro @ 2

Last change on this file since 2 was 2, checked in by opalod, 22 years ago
Initial revision
Property svn:eol-style set to `native` Property svn:executable set to ``* Property svn:keywords set to `Author Date Id Revision`
File size: 10.1 KB

Line
1	;------------------------------------------------------------
2	;------------------------------------------------------------
3	;------------------------------------------------------------
4	;+
5	; NAME:grille
6	;
7	; PURPOSE: choisit la grille qui doit etre utilisee pour faire le graphe en
8	; fonction de vargrid et renvoie les parametres correspondants calcules ds
9	; domdef.pro et reduit au domaine definit par domdef (contrairement a
10	; grandegrille.pro)
11	;
12	; CATEGORY:
13	;
14	; CALLING SEQUENCE:
15	; grille,mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,dernierx,derniery,dernierz,e1,e2,e3
16	;
17	; INPUTS:rien. ATTENTION les choix de la grille se fait a partir de la
18	; valeur de la variable globale vargrid, qui peut etre egale a 'T',
19	; 'U', 'V', 'W' ou 'F'.
20	;
21	; KEYWORD PARAMETERS:
22	; TRI si ce mot clef sert a obtenir grace a grille la
23	; triangulation qui se rapporte a la grille mais uniquement
24	; sur la partie du zoom. ce tableau de triangulation reduit
25	; est passe ds la variable que l''on a egalee a tri.par ex:
26	; grille,...,tri=triangulation_reduite. ne mot clef est
27	; utilise dans plt.pro
28	;
29	; /FORPLT: ds plt on veut que sur les points terres, glam et
30	; gphi soit egale a glamt et gphit quelle que soit la grille.
31	;
32	; /NOTRI: utile seulement qd TRI est active. dans ce cas
33	; grille retourne -1 ds la variable tri meme si la variable du
34	; common triangles est definie et differente de -1
35	;
36	; OUTPUTS:mask,glam,gphi,gdep,nx,ny,nz,premierx,premiery,premierz,
37	; dernierx,derniery,dernierz,e1,e2,e3
38	;
39	; pour leur definition cf domdef et la gestion des sous
40	; domaines sur le web
41	;
42	; Rq: ces outputs sont optionnels, si je veux recuperer que
43	; mask, glam et gphi il suffit de taper grille, mask, glam, gphi
44	;
45	; COMMON BLOCKS:
46	; common.pro congridseb.pro
47	;
48	; SIDE EFFECTS: utilise la variable globale vargird
49	;
50	; RESTRICTIONS: vargrid doit etre 'T', 'W', 'U', 'V' ou 'F'
51	;
52	; EXAMPLE:
53	;
54	; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
55	; 12/2/1999
56	; 10/11/1999 /forplt
57	;-
58	;------------------------------------------------------------
59	;------------------------------------------------------------
60	;------------------------------------------------------------
61	pro grille, mask, glam, gphi, gdep, nx, ny,nz,premierx,premiery,premierz,dernierx, derniery, dernierz, e1,e2,e3, TRI = tri, NOTRI = notri, TOUT = tout, FORPLT = forplt, _EXTRA = ex
62	@common
63	tempsun = systime(1) ; pour key_performance
64	;------------------------------------------------------------
65	if keyword_set(tout) then begin
66	oldboite = [lon1, lon2, lat1, lat2, prof1, prof2]
67	domdef, grille = vargrid, _EXTRA = ex
68	endif
69	tempdeux = systime(1) ; pour key_performance =2
70	CASE 1 OF
71	;------------------------------------------------------------
72	; grille T
73	;------------------------------------------------------------
74	vargrid eq 'OPAPTDHT' or vargrid eq 'OPAPT3DT' $
75	or vargrid eq 'T': begin
76	;scalaires
77	nx=nxt
78	ny=nyt
79	nz=nzt
80	premierx = premierxt
81	premiery = premieryt
82	premierz = premierzt
83	dernierx = dernierxt
84	derniery = dernieryt
85	dernierz = dernierzt
86	;vecteurs 2d
87	glam=glamt[premierx:dernierx, premiery:derniery]
88	gphi=gphit[premierx:dernierx, premiery:derniery]
89	e1 =e1t[premierx:dernierx, premiery:derniery]
90	e2 =e2t[premierx:dernierx, premiery:derniery]
91	;vecteurs 3d
92	mask=tmask[premierx:dernierx, premiery:derniery, premierz:dernierz]
93	end
94	;------------------------------------------------------------
95	; grille W
96	;------------------------------------------------------------
97	vargrid eq 'OPAPT3DW' or vargrid eq 'W': begin
98	;scalaires
99	nx=nxt
100	ny=nyt
101	nz=nzw
102	premierx = premierxt
103	premiery = premieryt
104	premierz = premierzw
105	dernierx = dernierxt
106	derniery = dernieryt
107	dernierz = dernierzw
108	;vecteurs 2d
109	terre = where(tmask[, , 0] EQ 0)
110	glam=glamt[premierx:dernierx, premiery:derniery]
111	gphi=gphit[premierx:dernierx, premiery:derniery]
112	e1 =e1t[premierx:dernierx, premiery:derniery]
113	e2 =e2t[premierx:dernierx, premiery:derniery]
114	;vecteurs 3d
115	mask=tmask[premierx:dernierx, premiery:derniery, premierz:dernierz]
116	end
117	;------------------------------------------------------------
118	; grille U
119	;------------------------------------------------------------
120	vargrid eq 'OPAPTDHU' or vargrid eq 'OPAPT3DU' $
121	or vargrid eq 'U': begin
122	;scalaires
123	nx=nxu
124	ny=nyu
125	nz=nzt
126	premierx = premierxu
127	premiery = premieryu
128	premierz = premierzt
129	dernierx = dernierxu
130	derniery = dernieryu
131	dernierz = dernierzt
132	;vecteurs 2d
133	glam=glamu[premierx:dernierx, premiery:derniery]
134	gphi=gphiu[premierx:dernierx, premiery:derniery]
135	if keyword_set(forplt) then BEGIN
136	mask = tmask[premierx:dernierx, premiery:derniery, niveau-1]
137	terre = where(mask EQ 0)
138	if terre[0] NE -1 then begin
139	glam[terre] = (glamt[premierx:dernierx, premiery:derniery])[terre]
140	gphi[terre] = (gphit[premierx:dernierx, premiery:derniery])[terre]
141	endif
142	ENDIF
143	e1 =e1u[premierx:dernierx, premiery:derniery]
144	e2 =e2u[premierx:dernierx, premiery:derniery]
145	;vecteurs 3d
146	mask=(umask())[premierx:dernierx, premiery:derniery, premierz:dernierz]
147	end
148	;------------------------------------------------------------
149	; grille V
150	;------------------------------------------------------------
151	vargrid eq 'OPAPTDHV' or vargrid eq 'OPAPT3DV' $
152	or vargrid eq 'V': begin
153	;scalaires
154	nx=nxv
155	ny=nyv
156	nz=nzt
157	premierx = premierxv
158	premiery = premieryv
159	premierz = premierzt
160	dernierx = dernierxv
161	derniery = dernieryv
162	dernierz = dernierzt
163	;vecteurs 2d
164	glam=glamv[premierx:dernierx, premiery:derniery]
165	gphi=gphiv[premierx:dernierx, premiery:derniery]
166	if keyword_set(forplt) then BEGIN
167	mask = tmask[premierx:dernierx, premiery:derniery, niveau-1]
168	terre = where(mask EQ 0)
169	if terre[0] NE -1 then begin
170	glam[terre] = (glamt[premierx:dernierx, premiery:derniery])[terre]
171	gphi[terre] = (gphit[premierx:dernierx, premiery:derniery])[terre]
172	endif
173	ENDIF
174	e1 =e1v[premierx:dernierx, premiery:derniery]
175	e2 =e2v[premierx:dernierx, premiery:derniery]
176	;vecteurs 3d
177	mask=(vmask())[premierx:dernierx, premiery:derniery, premierz:dernierz]
178	end
179	;------------------------------------------------------------
180	; grille F
181	;------------------------------------------------------------
182	vargrid eq 'OPAPTDHF' or vargrid eq 'OPAPT3DF' $
183	or vargrid eq 'F': begin
184	;scalaires
185	nx=nxf
186	ny=nyf
187	nz=nzt
188	premierx = premierxf
189	premiery = premieryf
190	premierz = premierzt
191	dernierx = dernierxf
192	derniery = dernieryf
193	dernierz = dernierzt
194	;vecteurs 2d
195	glam=glamf[premierx:dernierx, premiery:derniery]
196	gphi=gphif[premierx:dernierx, premiery:derniery]
197	if keyword_set(forplt) then BEGIN
198	mask = tmask[premierx:dernierx, premiery:derniery, niveau-1]
199	terre = where(mask EQ 0)
200	if terre[0] NE -1 then begin
201	glam[terre] = (glamt[premierx:dernierx, premiery:derniery])[terre]
202	gphi[terre] = (gphit[premierx:dernierx, premiery:derniery])[terre]
203	endif
204	ENDIF
205	e1 =e1f[premierx:dernierx, premiery:derniery]
206	e2 =e2f[premierx:dernierx, premiery:derniery]
207	;vecteurs 3d
208	mask=(fmask())[premierx:dernierx, premiery:derniery, premierz:dernierz]
209	END
210	;------------------------------------------------------------
211	ELSE:BEGIN
212	ras = report('Vauvaise definition de la variable vargrid: '+vargrid+'ceete variable doit etre T, U, V, W ou F')
213	stop
214	END
215	ENDCASE
216	IF testvar(var = key_performance) EQ 2 THEN $
217	print, 'temps grille: attribution des scalaires, vecteurs et tableaux ', systime(1)-tempdeux
218	;------------------------------------------------------------
219	; Variables se rapportant a la dimension verticale
220	;------------------------------------------------------------
221	tempdeux = systime(1) ; pour key_performance =2
222	if vargrid eq 'OPAPT3DW' or vargrid eq 'W' then begin
223	gdep = gdepw[premierz:dernierz]
224	e3=e3w[premierz:dernierz]
225	endif else begin
226	gdep = gdept[premierz:dernierz]
227	e3=e3t[premierz:dernierz]
228	endelse
229	IF testvar(var = key_performance) EQ 2 THEN $
230	print, 'temps grille: Variables se rapportant a la dimension verticale ', systime(1)-tempdeux
231	;------------------------------------------------------------
232	; vecteur triangulation Qd TRI est active
233	;------------------------------------------------------------
234	if arg_present(TRI) then $
235	if triangles[0] EQ -1 OR keyword_set(notri) then tri = -1 ELSE BEGIN
236	tempdeux = systime(1) ; pour key_performance =2
237	msk = bytarr(jpi, jpj)
238	msk[premierx:dernierx,premiery:derniery] = 1
239	ind = where( msk[triangles[0, ]]msk[triangles[1, ]]msk[triangles[2, *]] EQ 1 )
240	tri =triangles[, ind]-(premierx+premieryjpi)
241	y = tri/jpi
242	x = tri-y*jpi
243	tri = x+y*nx
244	IF testvar(var = key_performance) EQ 2 THEN $
245	print, 'temps grille: decoupage de la triangulation ', systime(1)-tempdeux
246	ENDELSE
247	;------------------------------------------------------------------
248	; pour s'assurer qu'il n'y a pas de dimension degenerees (=1)
249	;-------------------------------------------------------------------
250	; mask=reform(mask, /over)
251	; glam=reform(glam, /over)
252	; gphi=reform(gphi, /over)
253	; gdep=reform(gdep, /over)
254	; e1=reform(e1, /over)
255	; e2=reform(e2, /over)
256	; e3=reform(e3, /over)
257
258	if keyword_set(tout) then domdef, oldboite, grille = vargrid
259	if keyword_set(key_performance) THEN print, 'temps grille', systime(1)-tempsun
260
261	return
262
263	end
264
265
266
267
268
269
270

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

Download in other formats:

Original Format