Changeset 309
- Timestamp:
- 11/12/07 17:12:55 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro
r307 r309 314 314 b = reform(forweight[*, 1], 1, nawater) 315 315 forweight = -1 ; free memory 316 newaweig = [(1 -a)*(1-b), (1-b)*a, a*b, (1-a)*b]316 newaweig = [(1.d - a)*(1.d - b), (1.d - b)*a, a*b, (1.d - a)*b] 317 317 a = -1 & b = -1 ; free memory 318 318 ; mask the weight to suppress the corner located on land … … 322 322 ; weights normalization 323 323 totalweig = total(newaweig, 1, /double) 324 IF min(totalweig, max = ma) LE 0.d then stop325 IF ma GT 1.dthen stop324 ;; IF min(totalweig, max = ma) LE 0.d then stop 325 ;; IF ma- 1.d GT 1.e-6 then stop 326 326 newaweig = newaweig/(replicate(1.d, 4)#totalweig) 327 327 ; weights
Note: See TracChangeset
for help on using the changeset viewer.