Changeset 226 for trunk/SRC/ToBeReviewed/GRILLE/decoupeterre.pro
- Timestamp:
- 03/16/07 10:22:26 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/GRILLE/decoupeterre.pro
r163 r226 48 48 ; Used to pass your keywords 49 49 ; 50 ; @keyword TRI 50 ; @keyword TRI 51 51 ; This keyword serve to obtain, thanks to grille, the triangulation which 52 52 ; refer to the grid but only on the part of the zoom. This of triangulation … … 56 56 ; 57 57 ; @keyword WDEPTH 58 ; To specify that the field is at W depth instead of T 58 ; To specify that the field is at W depth instead of T 59 59 ; depth (automatically activated if vargrid eq 'W') 60 60 ; … … 184 184 if vargrid EQ 'V' OR vargrid EQ 'F' then $ 185 185 glam = glamv[firstx:lastx, firsty:lasty] 186 ENDIF ELSE BEGIN ; to drawsectionbottom 186 ENDIF ELSE BEGIN ; to drawsectionbottom 187 187 if vargrid EQ 'V' OR vargrid EQ 'F' OR finite(glamu[0]) EQ 0 then $ 188 188 glam = glamf[0 > (firstx-1):lastx, firsty:lasty] $ … … 194 194 if vargrid EQ 'U' OR vargrid EQ 'F' then $ 195 195 gphi = gphiu[firstx:lastx, firsty:lasty] 196 ENDIF ELSE BEGIN ; to drawsectionbottom 196 ENDIF ELSE BEGIN ; to drawsectionbottom 197 197 if vargrid EQ 'U' OR vargrid EQ 'F' OR finite(gphiv[0]) EQ 0 then $ 198 198 gphi = gphif[firstx:lastx, 0 > (firsty-1):lasty] $ … … 211 211 firstz = firstzt 212 212 lastz = lastzt 213 ENDELSE 213 ENDELSE 214 214 nz = lastz-firstz+1 215 215 ;------------------------------------------------------------ … … 221 221 profond = firstz NE 0 222 222 END 223 ; for the ver ical section, we have to choose the right mask according223 ; for the vertical section, we have to choose the right mask according 224 224 ; to the grid point and to the direction of the section 225 225 'xz':BEGIN … … 267 267 END 268 268 keyword_set(key_partialstep) and type EQ 'yz' $ 269 AND nx EQ 1 AND keyword_set(realsection):BEGIN 269 AND nx EQ 1 AND keyword_set(realsection):BEGIN 270 270 bottom = total(mask, 3) 271 271 good = where(bottom NE 0 AND bottom NE nz+1) … … 290 290 tri = triangule(mask, coinmonte = coinmonte, coindescend = coindescend, _extra = ex) 291 291 indicezoom = (lindgen(jpi, jpj))[firstx:lastx, firsty:lasty] 292 ENDIF ELSE BEGIN 292 ENDIF ELSE BEGIN 293 293 ; Otherwise, we recuperate the part of triangulation that interest us and we number them well!! 294 if nx EQ jpi AND ny EQ jpj then tri = triangles_list ELSE BEGIN 294 if nx EQ jpi AND ny EQ jpj then tri = triangles_list ELSE BEGIN 295 295 msk = bytarr(jpi, jpj) 296 296 msk[firstx:lastx, firsty:lasty] = 1 … … 306 306 ENDELSE 307 307 ;------------------------------------------------------------------- 308 if keyword_set(key_performance) THEN print, 'temps decoupeterre', systime(1)-tempsun 308 if keyword_set(key_performance) THEN print, 'temps decoupeterre', systime(1)-tempsun 309 309 ;------------------------------------------------------------ 310 310 return
Note: See TracChangeset
for help on using the changeset viewer.