[Dart-dev] DART/branches Revision: 13208

[dart at ucar.edu]

nancy at ucar.edu
2019-07-31 08:48:41 -0600 (Wed, 31 Jul 2019)
816
start of making gitm work with manhattan.  NOT FINISHED!!

created a new executable to convert a single gitm .bin
file to a netcdf file (instead of a lanai-style restart file).
it's unclear to me if what i have is a single block of a larger
dataset (since it has ghost zones) or if this is the full dataset:w

modified all the path_names files to compile with manhattan.

started work on the model_mod - made the interfaces work.
this will NOT work for filter because get_close_obs() will
need fixing, and there isn't currently a get_close_state().
also may need to write the 2 new vertical convert routines.
for now i'm just passing through the default versions which
satisfy the compiler but do nothing.

also made the quickbuild match what we currently use in
other models (doesn't build mpi targets twice, etc)




Modified: DART/branches/rma_trunk/models/gitm/dart_to_gitm.f90
===================================================================
--- DART/branches/rma_trunk/models/gitm/dart_to_gitm.f90	2019-07-25 23:24:59 UTC (rev 13207)
+++ DART/branches/rma_trunk/models/gitm/dart_to_gitm.f90	2019-07-31 14:48:41 UTC (rev 13208)
@@ -9,15 +9,10 @@
 !----------------------------------------------------------------------
 ! purpose: interface between DART and the GITM model
 !
-! method: Read DART state vector and overwrite values in a gitm restart file.
-!         If the DART state vector has an 'advance_to_time' present, a
-!         file called DART_GITM_time_control.txt is created with
-!         information appropriate to advance gitm to the requested time.
+! method: Read DART state netcdf files and overwrite values in a gitm restart file.
 !
-!         The dart_to_gitm_nml namelist setting for advance_time_present
-!         determines whether or not the input file has an 'advance_to_time'.
-!         Typically, only temporary files like 'assim_model_state_ic' have
-!         an 'advance_to_time'.
+! NOT FINISHED!  HARDLY STARTED.
+!
 !----------------------------------------------------------------------
 
 use        types_mod, only : r8
@@ -27,10 +22,8 @@
                              logfileunit, open_file, close_file
 
 use        model_mod, only : static_init_model, get_model_size, &
-                             get_gitm_restart_dirname, statevector_to_restart_file
+                             get_gitm_restart_dirname !! , statevector_to_restart_file
 
-use  assim_model_mod, only : open_restart_read, aread_state_restart, close_restart
-
 use time_manager_mod, only : time_type, print_time, print_date, operator(-), &
                              get_time, get_date
 
@@ -46,11 +39,9 @@
 ! namelist parameters with default values.
 !-----------------------------------------------------------------------
 
-character (len = 128) :: dart_to_gitm_input_file = 'dart_restart'
-logical               :: advance_time_present    = .false.
+character (len = 128) :: dart_to_gitm_input_file = 'filter_restart.nc'
 
-namelist /dart_to_gitm_nml/ dart_to_gitm_input_file, &
-                            advance_time_present
+namelist /dart_to_gitm_nml/ dart_to_gitm_input_file
 
 !----------------------------------------------------------------------
 ! global storage
@@ -58,7 +49,6 @@
 
 integer               :: iunit, io, x_size
 type(time_type)       :: model_time, adv_to_time
-real(r8), allocatable :: statevector(:)
 character(len=256)    :: gitm_restart_dirname  = 'none'
 
 !======================================================================
@@ -86,20 +76,20 @@
 write(*,*) 'to gitm restart files in directory ', "'"//trim(gitm_restart_dirname)//"'"
 
 x_size = get_model_size()
-allocate(statevector(x_size))
+!allocate(statevector(x_size))
 
 !----------------------------------------------------------------------
 ! Reads the valid time, the state, and the target time.
 !----------------------------------------------------------------------
 
-iunit = open_restart_read(dart_to_gitm_input_file)
+!iunit = open_restart_read(dart_to_gitm_input_file)
 
-if ( advance_time_present ) then
-   call aread_state_restart(model_time, statevector, iunit, adv_to_time)
-else
-   call aread_state_restart(model_time, statevector, iunit)
-endif
-call close_restart(iunit)
+!if ( advance_time_present ) then
+!   call aread_state_restart(model_time, statevector, iunit, adv_to_time)
+!else
+!   call aread_state_restart(model_time, statevector, iunit)
+!endif
+!call close_restart(iunit)
 
 !----------------------------------------------------------------------
 ! update the current gitm state vector
@@ -107,11 +97,11 @@
 ! time_manager_nml: stop_option, stop_count increments
 !----------------------------------------------------------------------
 
-call statevector_to_restart_file(statevector, gitm_restart_dirname, model_time)
+!call statevector_to_restart_file(statevector, gitm_restart_dirname, model_time)
 
-if ( advance_time_present ) then
-   call write_gitm_time_control(model_time, adv_to_time)
-endif
+!if ( advance_time_present ) then
+!   call write_gitm_time_control(model_time, adv_to_time)
+!endif
 
 !----------------------------------------------------------------------