Changeset 226 for trunk/SRC/ToBeReviewed/CALCULS/remplit.pro
- Timestamp:
- 03/16/07 10:22:26 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/CALCULS/remplit.pro
r209 r226 22 22 ; 23 23 ; @keyword FILLXDIR 24 ; 24 ; 25 25 ; 26 26 ; @keyword FILLYDIR … … 51 51 ; $Id$ 52 52 ; 53 ;; 54 ;; Extrapole zinout[jpi,jpj] sur les continents en utilisant les 4 55 ;; plus proches voisins masques oceaniquement et construit un nouveau 56 ; masque 57 ;; contenant l'ancien masque oceanique PLUSles points extrapoles. 58 ;; Reitere le processus niter fois. 59 ;; C'est pas clair, essayez ! 60 ;; 61 ;; 53 ;; 54 ;; Extrapole zinout[jpi,jpj] sur les continents en utilisant les 4 55 ;; plus proches voisins masques oceaniquement et construit un nouveau masque 56 ;; contenant l'ancien masque oceanique PLUS les points extrapoles. 57 ;; Reitere le processus niter fois. 58 ;; C'est pas clair, essayez ! 59 ;; 60 ;; 62 61 ; 63 62 ; /Nan: to fill the point which have the value 64 ; !values.f_nan. W hitout this keyword, these point are not filling63 ; !values.f_nan. Without this keyword, these point are not filling 65 64 ; and stays at !values.f_nan. 66 65 ; … … 139 138 ;--------------------------------------------------------------- 140 139 ;--------------------------------------------------------------- 141 ; iteration 142 ;--------------------------------------------------------------- 143 ;--------------------------------------------------------------- 144 FOR n = 1, niter DO BEGIN 140 ; iteration 141 ;--------------------------------------------------------------- 142 ;--------------------------------------------------------------- 143 FOR n = 1, niter DO BEGIN 145 144 ; on trouve les points coast 146 145 tempdeux = systime(1) ; pour key_performance =2 … … 244 243 END 245 244 endcase 246 ; 245 ; 247 246 z[coast] = temporary(zcoast)/ temporary(weight) 248 247 ; we update the the boundary conditions of z … … 258 257 IF testvar(var = key_performance) EQ 2 THEN $ 259 258 print, 'temps remplit: une iteration ', systime(1)-tempdeux 260 ENDFOR 259 ENDFOR 261 260 fini: 262 261 ; … … 280 279 key_gridtype = oldkey_gridtype 281 280 ;--------------------------------------------------------------- 282 if keyword_set(key_performance) THEN print, 'temps remplit', systime(1)-tempsun 281 if keyword_set(key_performance) THEN print, 'temps remplit', systime(1)-tempsun 283 282 return, z 284 END 283 END 285 284
Note: See TracChangeset
for help on using the changeset viewer.