Changeset 3182 for palm/trunk/UTIL/inifor/tests/test-interpolation.f90
- Timestamp:
- Jul 27, 2018 1:36:03 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/UTIL/inifor/tests/test-interpolation.f90
r2718 r3182 21 21 ! Current revisions: 22 22 ! ----------------- 23 ! 23 ! Updated test for new grid_definition 24 ! 24 25 ! 25 26 ! Former revisions: … … 76 77 xmin = -5.0 * TO_RADIANS, xmax = 5.5 * TO_RADIANS, & 77 78 ymin = -5.0 * TO_RADIANS, ymax = 6.5 * TO_RADIANS, & 78 zmin = 0.0, zmax = 10.0, &79 79 x0 = 0.0, y0 = 0.0, z0 = 0.0, & 80 80 nx = nlon-1, ny = nlat-1, nz = nlev-1) … … 86 86 res = assert_equal( (/cosmo_grid%lat(0), cosmo_grid % lon(0), & 87 87 cosmo_grid%lat(2), cosmo_grid % lon(2), & 88 cosmo_grid%dx*TO_DEGREES, cosmo_grid%dy*TO_DEGREES/),& 88 (cosmo_grid%lon(1) - cosmo_grid%lon(0))*TO_DEGREES, & 89 (cosmo_grid%lat(1) - cosmo_grid%lat(0))*TO_DEGREES/),& 89 90 (/-5.0 * TO_RADIANS, -5.0 * TO_RADIANS, & 90 91 6.5 * TO_RADIANS, 5.5 * TO_RADIANS, & … … 97 98 xmin = 0.0, xmax = 1.0, & 98 99 ymin = 0.0, ymax = 1.0, & 99 zmin = 0.0, zmax = 1.0, &100 100 x0 = 0.0, y0 = 0.0, z0 = 0.0, & 101 101 nx = 1, ny = 1, nz = 1) … … 127 127 ! Act 128 128 CALL find_horizontal_neighbours(cosmo_grid % lat, cosmo_grid % lon, & 129 cosmo_grid % dxi, cosmo_grid % dyi, palm_grid % clat, palm_grid % clon,&130 palm_grid % ii, palm_grid % jj)129 palm_grid % clat, palm_grid % clon, & 130 palm_grid % ii, palm_grid % jj) 131 131 132 132 ! Assert … … 178 178 ! Act 179 179 CALL find_horizontal_neighbours(cosmo_grid % lat, cosmo_grid % lon, & 180 cosmo_grid % dxi, cosmo_grid % dyi, palm_grid % clat, palm_grid % clon,&181 palm_grid % ii, palm_grid % jj)180 palm_grid % clat, palm_grid % clon, & 181 palm_grid % ii, palm_grid % jj) 182 182 183 183 CALL compute_horizontal_interp_weights(cosmo_grid % lat, cosmo_grid % lon, & 184 cosmo_grid % dxi, cosmo_grid % dyi, palm_grid % clat, & 185 palm_grid % clon, palm_grid % ii, palm_grid % jj, palm_grid % w_horiz) 184 palm_grid % clat, palm_grid % clon, & 185 palm_grid % ii, palm_grid % jj, & 186 palm_grid % w_horiz) 186 187 187 188 ! Assert
Note: See TracChangeset
for help on using the changeset viewer.