Changeset 48 for trunk/divfred.pro
- Timestamp:
- 03/16/14 20:38:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/divfred.pro
r41 r48 25 25 ; - Requires SAXO tools 26 26 ; - les matrices u et v peuvent de 2 a 4 dimensions. 27 ; attention pour distinguer les diff érentes configurations de u et v27 ; attention pour distinguer les différentes configurations de u et v 28 28 ; (xy, xyz, xyt, xyzt), on regarde la variable du common 29 29 ; -time qui contient le calendrier en jour julien d''IDL auquel 30 30 ; se rapportent u et v ansi que la variable 31 ; -jpt qui est le nombre de pas de temps à considérer dans time.32 ; les tableaux u et v sont d écoupés sur le même domaine33 ; g éographique. A cause du décalage des grilles T, U, V et F il est34 ; possible que ces 2 tableaux n''aient pas la m ême taille et se35 ; reportent à des indices différents. Si tel est le cas les tableaux36 ; sont red écoupés sur les indices qu'ils ont en commun et le dommaine37 ; est red éfini pour qu'il colle àces indices communs.38 ; pour éviter ces redécoupes utiliser le mot clé/memeindice dans31 ; -jpt qui est le nombre de pas de temps à considérer dans time. 32 ; les tableaux u et v sont découpés sur le même domaine 33 ; géographique. A cause du décalage des grilles T, U, V et F il est 34 ; possible que ces 2 tableaux n''aient pas la même taille et se 35 ; reportent à des indices différents. Si tel est le cas les tableaux 36 ; sont redécoupés sur les indices qu'ils ont en commun et le dommaine 37 ; est redéfini pour qu'il colle à ces indices communs. 38 ; pour éviter ces redécoupes utiliser le mot clé /memeindice dans 39 39 ; <pro>domdef</pro> 40 40 ; 41 ; les points sur le bord du dessin sont mis à!values.f_nan41 ; les points sur le bord du dessin sont mis à !values.f_nan 42 42 ; 43 43 ; @todo 44 ; ++ pas fini de comprendre, tester (compatibilit ésaxo), adapter, commenter44 ; ++ pas fini de comprendre, tester (compatibilité saxo), adapter, commenter 45 45 ; nettoyer 46 46 ; ++ a comparer et merger avec SAXO_DIR/ToBeReviewed/CALCULS/div.pro … … 111 111 v = v[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, *] 112 112 END 113 ELSE :BEGIN113 ELSE :BEGIN 114 114 zdiv = -1 115 115 GOTO, sortie … … 142 142 ; mise a !values.f_nan de la bordure 143 143 ;------------------------------------------------------------ 144 if 144 if NOT keyword_set(key_periodic) OR nx NE jpi then begin 145 145 zdiv(0, *, *) = !values.f_nan 146 146 zdiv(nx-1, *, *) = !values.f_nan … … 161 161 varunits = '1e6*s-1' 162 162 domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 163 if keyword_set(direc) then 163 if keyword_set(direc) then zdiv = moyenne(zdiv,direc,/nan, boite = boite) 164 164 END 165 165 ;---------------------------------------------------------------------------- … … 189 189 v = v[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, *] 190 190 END 191 ELSE :return, -1191 ELSE : return, -1 192 192 endcase 193 193 ;------------------------------------------------------------ … … 216 216 ; mise a !values.f_nan de la bordure 217 217 ;------------------------------------------------------------ 218 if 218 if NOT keyword_set(key_periodic) OR nx NE jpi then begin 219 219 zdiv(0, *, *) = !values.f_nan 220 220 zdiv(nx-1, *, *) = !values.f_nan … … 232 232 varunits = '1e6*s-1' 233 233 domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 234 if keyword_set(direc) then 234 if keyword_set(direc) then zdiv = grossemoyenne(zdiv,direc,/nan, boite = boite) 235 235 END 236 236 ;---------------------------------------------------------------------------- … … 247 247 ;---------------------------------------------------------------------------- 248 248 ;---------------------------------------------------------------------------- 249 ELSE :BEGIN ;xy249 ELSE : BEGIN ;xy 250 250 indice3d = lindgen(jpi, jpj, jpk) 251 251 indice3d = indice3d[indicex[0]:indicex[0]+nx-1,indicey[0]:indicey[0]+ny-1, firstzt] … … 273 273 v = v[indice2d] 274 274 END 275 ELSE :return, -1275 ELSE : return, -1 276 276 endcase 277 277 ;------------------------------------------------------------ … … 289 289 ; mise a !values.f_nan de la bordure 290 290 ;------------------------------------------------------------ 291 if 291 if NOT keyword_set(key_periodic) OR nx NE jpi then begin 292 292 zdiv(0, *) = !values.f_nan 293 293 zdiv(nx-1, *) = !values.f_nan … … 308 308 varunits = '1e6*s-1' 309 309 domdef, (glamt[indice2d])[0, 0], (glamu[indice2d])[nx-1, 0],(gphit[indice2d])[0, 0], (gphiv[indice2d])[0, ny-1], vert1, vert2, grille = ['t'] 310 if keyword_set(direc) then 310 if keyword_set(direc) then zdiv = moyenne(zdiv,direc,/nan, boite = boite) 311 311 312 312 END
Note: See TracChangeset
for help on using the changeset viewer.