[2] | 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+premiery*jpi) |
---|
| 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 | |
---|