Changeset 48 for trunk/extrap.pro
- Timestamp:
- 03/16/14 20:38:39 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/extrap.pro
r41 r48 2 2 ; NAME: extrap.pro (based on remplit.pro) 3 3 ; 4 ; PURPOSE: Extrapolates ocean data values on land points 4 ; PURPOSE: Extrapolates ocean data values on land points 5 5 ; 6 6 ; CATEGORY: Subroutine … … 8 8 ; CALLING SEQUENCE: extrap, zdata, zmask, it, val 9 9 ; 10 ; INPUTS: 10 ; INPUTS: 11 11 ; zdata : field to extrapolate 12 12 ; zmask : field's mask 13 13 ; it : iteration 14 ; val : field value on masked points 14 ; val : field value on masked points 15 15 ; 16 16 ; KEYWORD PARAMETERS: None 17 17 ; 18 ; OUTPUTS: 19 ; zdata : extrapolated field 18 ; OUTPUTS: 19 ; zdata : extrapolated field 20 20 ; zmask : mask after extrapolation 21 21 ; 22 22 ; COMMON BLOCKS: 23 ; common_interp.pro23 ; common_interp.pro 24 24 ; 25 25 ; … … 38 38 ; IF keyword_set(key_performance) THEN print, systime(1)-temp 39 39 ; 40 ; 2. Extrapolating when needed 40 ; 2. Extrapolating when needed 41 41 ; ============================ 42 42 ; … … 71 71 ; IF keyword_set(key_performance) THEN print, systime(1)-temp 72 72 ; IF keyword_set(key_performance) THEN temp = systime(1) 73 ; IF n_elements(ocean2) NE jpiatm*(jpjatm+4) THEN BEGIN 73 ; IF n_elements(ocean2) NE jpiatm*(jpjatm+4) THEN BEGIN 74 74 ; land = where(zmask EQ 0.) 75 75 ; zremplit(land) = val 76 ; ENDIF 76 ; ENDIF 77 77 ; IF keyword_set(key_performance) THEN print, systime(1)-temp 78 78 ; … … 95 95 terre = where(mmmask EQ 0) 96 96 if terre[0] EQ -1 then GOTO, fini 97 ; les points du bord du cadre doivent maintenan etre dans la terre97 ; les points du bord du cadre doivent maintenant etre dans la terre 98 98 mmmask[0, *] = 0 99 99 mmmask[jpiatm+1, *] = 0 … … 119 119 +1./sqrt(2)*(z[(jpiatm+2)+1+cote]+z[(jpiatm+2)-1+cote]+z[-(jpiatm+2)+1+cote]+z[-(jpiatm+2)-1+cote]) 120 120 poids = voisins[cote] 121 121 122 122 z[cote] = zcote/poids 123 123 ;---------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.