[Dart-dev] DART/branches Revision: 11021

[dart at ucar.edu]

nancy at ucar.edu
2017-02-07 12:55:59 -0700 (Tue, 07 Feb 2017)
141
make this program compile, but i have questions about
the interpolate routines which don't seem to actually
return any interpolated values?




Modified: DART/branches/rma_single_file/models/ROMS/test_grid.f90
===================================================================
--- DART/branches/rma_single_file/models/ROMS/test_grid.f90	2017-02-07 18:46:44 UTC (rev 11020)
+++ DART/branches/rma_single_file/models/ROMS/test_grid.f90	2017-02-07 19:55:59 UTC (rev 11021)
@@ -29,8 +29,7 @@
 
 use       assim_model_mod, only : static_init_assim_model
 
-use  state_space_diag_mod, only : aoutput_diagnostics, init_diag_output,        &
-                                  finalize_diag_output, netcdf_file_type
+use  state_space_diag_mod, only : netcdf_file_type
 
 use      time_manager_mod, only : time_type, set_calendar_type, GREGORIAN,      &
                                   set_time, print_time, print_date, operator(-)
@@ -159,8 +158,8 @@
 call init_ensemble_manager(ens_handle, num_ens, model_size)
 
 ! Reading netcdf restart file:
-input_file_info = io_filenames_init(num_ens, .false., root_name='input')
-output_file_info = io_filenames_init(num_ens, .false., root_name='output')
+call io_filenames_init(input_file_info, num_ens, .false., .false., root_name='input')
+call io_filenames_init(output_file_info, num_ens, .false.,.false.,  root_name='output')
 
 
 !----------------------------------------------------------------------

Modified: DART/branches/rma_single_file/models/ROMS/test_roms_interpolate.f90
===================================================================
--- DART/branches/rma_single_file/models/ROMS/test_roms_interpolate.f90	2017-02-07 18:46:44 UTC (rev 11020)
+++ DART/branches/rma_single_file/models/ROMS/test_roms_interpolate.f90	2017-02-07 19:55:59 UTC (rev 11021)
@@ -10,7 +10,7 @@
 ! interpolation test routines for threed sphere locations.
 !-------------------------------------------------------------------------------
 
-use             types_mod, only : r8, i8, missing_r8, metadatalength
+use             types_mod, only : r8, i8, MISSING_R8, metadatalength
 
 use         utilities_mod, only : register_module, error_handler, E_MSG, E_ERR, &
                                   initialize_utilities, finalize_utilities,     &
@@ -308,7 +308,7 @@
                                   jval,             &
                                   kval,             &
                                   mykindindex,      &
-                                  location,         &
+                                  value,            &
                                   ios_out) result(num_passed)
 
 type(ensemble_type)   , intent(inout) :: ens_handle
@@ -317,7 +317,7 @@
 real(r8)              , intent(in)    :: jval
 real(r8)              , intent(in)    :: kval
 integer               , intent(in)    :: mykindindex
-real(r8)              , intent(out)   :: location
+real(r8)              , intent(out)   :: value
 integer               , intent(out)   :: ios_out
 integer :: num_passed
 
@@ -331,6 +331,9 @@
 ios_out   = get_location_from_ijk(ival, jval, kval, mykindindex, loc)
 my_loc    = get_location(loc) 
 
+! where is the code that actually does the interpolate??  nsc.
+value = MISSING_R8
+
 if (ios_out == 0 ) then
    if (do_output() ) &
       write(*,*) 'test_interpolate_single SUCCESS lat,lon,lev', my_loc(1),my_loc(2),my_loc(3)