[Dart-dev] DART/branches Revision: 11974

[dart at ucar.edu]

nancy at ucar.edu
2017-10-11 17:57:08 -0600 (Wed, 11 Oct 2017)
116
this program became the basis for test_quad_irreg_interp and
then the fully regular quad case was split from that.




Deleted: DART/branches/quad_interp/models/utilities/test_quad_interp.f90
===================================================================
--- DART/branches/quad_interp/models/utilities/test_quad_interp.f90	2017-10-11 23:56:21 UTC (rev 11973)
+++ DART/branches/quad_interp/models/utilities/test_quad_interp.f90	2017-10-11 23:57:08 UTC (rev 11974)
@@ -1,239 +0,0 @@
-! DART software - Copyright 2004 - 2013 UCAR. This open source software is
-! provided by UCAR, "as is", without charge, subject to all terms of use at
-! http://www.image.ucar.edu/DAReS/DART/DART_download
-!
-! DART $Id$
-
-program test_quad_interp
-
-! intended to show how the state structure and quad code can be used
-! together.  start with a simple regional grid and work out from there.
-
-! Modules that are absolutely required for use are listed
-use        types_mod, only : r8, i8, MISSING_R8, deg2rad, rad2deg
-use    utilities_mod, only : error_handler, open_file, close_file
-use   random_seq_mod, only : init_random_seq, random_seq_type, &
-                             random_uniform, random_gaussian
-
-use quad_utils_mod, only : quad_interp_handle, init_quad_interp, finalize_quad_interp, set_quad_coords,       &
-                           quad_lon_lat_locate, quad_lon_lat_evaluate, GRID_QUAD_FULLY_REGULAR,               &
-                           GRID_QUAD_IRREG_SPACED_REGULAR, GRID_QUAD_FULLY_IRREGULAR, GRID_QUAD_UNKNOWN_TYPE, &
-                           QUAD_LOCATED_UNKNOWN, QUAD_LOCATED_CELL_CENTERS, QUAD_LOCATED_LON_EDGES,           &
-                           QUAD_LOCATED_LAT_EDGES, QUAD_LOCATED_CELL_CORNERS
-
-
-implicit none
-
-
-type(quad_interp_handle) :: h
-
-! data grid size
-integer, parameter :: nx = 9
-integer, parameter :: ny = 5
-
-! locations of data grid corners
-real(r8) :: data_lons(nx, ny) = MISSING_R8
-real(r8) :: data_lats(nx, ny) = MISSING_R8
-
-! extents of the data grid (these mimic a regional model's grid)
-real(r8) :: start_lon = 100.0_r8
-real(r8) :: end_lon   = 150.5_r8
-real(r8) :: start_lat = -11.4_r8
-real(r8) :: end_lat   =  34.1_r8
-
-! angle to rotate data grid in degrees
-! positive is counterclockwise; will rotate
-! around lower left grid point (start lon/lat).
-!real(r8) :: angle = 45.0_r8
-!real(r8) :: angle = 30.0_r8
-!real(r8) :: angle =  90.0_r8
-!real(r8) :: angle = -30.0_r8
-!real(r8) :: angle = -10.0_r8
- real(r8) :: angle = 0.0_r8
-
-! deform grid by this fraction of the deltas
-real(r8) :: lon_def = 0.3_r8
-real(r8) :: lat_def = 0.3_r8
-
-! data values on the grid
-real(r8) :: grid_data(nx, ny) = MISSING_R8
-
-
-! sampling grid size
-integer, parameter :: nrx = 210
-integer, parameter :: nry = 150
-
-! locations of sampling grid
-real(r8) :: reg_lons(nrx) = MISSING_R8
-real(r8) :: reg_lats(nrx) = MISSING_R8
-
-! extents of the sampling grid
-real(r8) :: reg_start_lon = 110.0_r8
-real(r8) :: reg_end_lon   = 140.0_r8
-real(r8) :: reg_start_lat = -20.0_r8
-real(r8) :: reg_end_lat   =  30.0_r8
-
-! where interpolated values are stored on reg grid
-real(r8) :: interp_data(nrx, nry) = MISSING_R8
-
-
-type(random_seq_type) :: ran
-
-integer  :: i, j
-real(r8) :: lon_del, lat_del, reg_lon_del, reg_lat_del
-integer  :: lon_bot, lat_bot, lon_top, lat_top
-integer  :: istatus
-real(r8) :: invals(4), outval
-integer  :: iunit_orig, iunit_interp
-
-call init_random_seq(ran)
-
-iunit_orig = open_file('original_data.txt', action='write')
-iunit_interp = open_file('interp_data.txt', action='write')
-
-lon_del = (end_lon - start_lon) / (nx-1)
-lat_del = (end_lat - start_lat) / (ny-1)
-