- Timestamp:
- 03/16/07 10:22:26 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro
r202 r226 5 5 ; to any grid using the bilinear method 6 6 ; 7 ; @categories 7 ; @categories 8 8 ; Interpolation 9 9 ; … … 53 53 ; June 2006: Sebastien Masson (smasson\@lodyc.jussieu.fr) 54 54 ; 55 ; @version $Id$ 55 ; @version 56 ; $Id$ 56 57 ; 57 58 ;- … … 65 66 ; 66 67 jpio = (size(olonin, /dimensions))[0] 67 jpjo = (size(olonin, /dimensions))[1] 68 jpjo = (size(olonin, /dimensions))[1] 68 69 ; mask check 69 70 IF n_elements(omsk) EQ 1 AND omsk[0] EQ -1 THEN omsk = replicate(1b, jpio, jpjo) … … 73 74 stop 74 75 ENDIF 75 IF n_elements(amsk) NE jpia*jpja THEN BEGIN 76 IF n_elements(amsk) NE jpia*jpja THEN BEGIN 76 77 print, 'output grid mask do not have the good size' 77 78 stop … … 181 182 ind = good[ind] 182 183 ; now, we morph the quadrilateral ocean cell into the reference square (0 -> 1) 183 ; in addition we get the co rrdinates of the atmospheric point in this new morphed square184 ; in addition we get the coordinates of the atmospheric point in this new morphed square 184 185 IF onsphere THEN BEGIN 185 186 ; Warning! quadrilateral2square use anticlockwise quadrilateral definition … … 213 214 xy[1] EQ 0 AND (tmpmsk[0]+tmpmsk[1]) EQ 0 : foraddr[n] = -1 214 215 xy[1] EQ 1 AND (tmpmsk[2]+tmpmsk[3]) EQ 0 : foraddr[n] = -1 215 ELSE: BEGIN 216 ELSE: BEGIN 216 217 ; we keep its address 217 218 foraddr[n] = ind
Note: See TracChangeset
for help on using the changeset viewer.