Changeset 12243 for utils/tools/NESTING
- Timestamp:
- 2019-12-13T12:55:39+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
utils/tools/NESTING/src/agrif_interpolation.f90
r10381 r12243 198 198 CHARACTER(*) :: typevar 199 199 200 INTEGER :: nxf,nyf, zx,zy201 INTEGER :: ji,jj,jif,jjf,jic,jjc,j decx,jdecy200 INTEGER :: nxf,nyf,nxc,nyc,zx,zy 201 INTEGER :: ji,jj,jif,jjf,jic,jjc,jic1,jjc1,jdecx,jdecy 202 202 REAL*8 :: Ax, Bx, Ay, By 203 203 204 204 nxf = SIZE(tabout,1) 205 205 nyf = SIZE(tabout,2) 206 207 nxc = SIZE(tabin,1) 208 nyc = SIZE(tabin,2) 206 209 207 210 SELECT CASE(typevar) … … 248 251 Ay = 1. - By 249 252 250 tabout(ji,jj) = ( Bx * tabin(jic+1,jjc ) + Ax * tabin(jic,jjc ) ) * Ay + & 251 & ( Bx * tabin(jic+1,jjc+1) + Ax * tabin(jic,jjc+1) ) * By 253 jic1 = MIN( nxc, jic+1 ) ! avoid out of bounds for tabin below 254 jjc1 = MIN( nyc, jjc+1 ) ! -- 255 256 tabout(ji,jj) = ( Bx * tabin(jic1,jjc ) + Ax * tabin(jic,jjc ) ) * Ay + & 257 & ( Bx * tabin(jic1,jjc1) + Ax * tabin(jic,jjc1) ) * By 252 258 END DO 253 259 END DO
Note: See TracChangeset
for help on using the changeset viewer.