Changeset 114 for trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro
- Timestamp:
- 06/19/06 16:14:56 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/PLOTS/VECTEUR/ajoutvect.pro
r41 r114 54 54 ;------------------------------------------------------------ 55 55 pro ajoutvect,vecteur, vectlegende, UNVECTSUR=unvectsur,VECTMIN=vectmin, VECTMAX=vectmax, _EXTRA = ex 56 ; 57 compile_opt idl2, strictarrsubs 58 ; 56 59 @common 57 60 tempsun = systime(1) ; pour key_performance … … 128 131 ; construction de u et v aux pts T 129 132 ;----------------------------------------------------------- 130 a=u (0,*)133 a=u[0,*] 131 134 u=(u+shift(u,1,0))/2. 132 if NOT keyword_set(key_periodic) OR nx NE jpi then u (0,*)=a133 a=v (*,0)135 if NOT keyword_set(key_periodic) OR nx NE jpi then u[0,*]=a 136 a=v[*,0] 134 137 v=(v+shift(v,0,1))/2. 135 if NOT keyword_set(key_periodic) OR nx NE jpi then v (*,0)=a138 if NOT keyword_set(key_periodic) OR nx NE jpi then v[*,0]=a 136 139 ;---------------------------------------------------------------------------- 137 140 ; attribution du mask et des tableau de longitude et latitude … … 191 194 if interpolle then t2 = msku*shift(msku,1,0)*mskv*shift(mskv,0,1) $ 192 195 ELSE t2 = tmask[firstxt:lastxt,firstyt:lastyt,firstzt] 193 if NOT keyword_set(key_periodic) OR nx NE jpi then t2 (0, *)=0.194 t2 (*,0)=0.196 if NOT keyword_set(key_periodic) OR nx NE jpi then t2[0, *]=0. 197 t2[*,0]=0. 195 198 terre=where(t2 eq 0) 196 199 if terre[0] ne -1 then begin 197 u (terre)=1e5198 v (terre)=1e5200 u[terre]=1e5 201 v[terre]=1e5 199 202 ENDIF 200 203 ;-----------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.