[Dart-dev] DART/branches Revision: 12561

[dart at ucar.edu]

nancy at ucar.edu
2018-04-26 10:55:53 -0600 (Thu, 26 Apr 2018)
224
version from the rma_fill_inf branch.  this works for
the multiple file option but not for single file yet.
(so this is fine for cam).

reads a separate namelist because we can't currently
share the filter namelist easily.




Modified: DART/branches/recam/assimilation_code/programs/fill_inflation_restart/fill_inflation_restart.f90
===================================================================
--- DART/branches/recam/assimilation_code/programs/fill_inflation_restart/fill_inflation_restart.f90	2018-04-26 16:52:37 UTC (rev 12560)
+++ DART/branches/recam/assimilation_code/programs/fill_inflation_restart/fill_inflation_restart.f90	2018-04-26 16:55:53 UTC (rev 12561)
@@ -2,43 +2,442 @@
 ! by UCAR, "as is", without charge, subject to all terms of use at
 ! http://www.image.ucar.edu/DAReS/DART/DART_download
 !
-! $Id$ 
+! $Id$
 
-! THIS PROGRAM HAS BEEN OBSOLETED by converting our restart files to
-! NetCDF format.  you can assign a value with 'ncap2', one of the NCO
-! utility program.  for now i'm going to leave the mkmf and path_names
-! files around and the program will just print out a helpful message.
-! in a few months/years we should remove this program entirely.
-!
-! here's an example of using ncap2 to set the T,U and V inf values:
-!  ncap2 -s 'T=1.0;U=1.0;V=1.0' wrfinput_d01 prior_inf.nc
-!  ncap2 -s 'T=0.6;U=0.6;V=0.6' wrfinput_d01 prior_sd.nc 
-!
-! this works as long as you have at least version 4.4.2 of the NCO utils.
-! some earlier versions change the full 3d arrays into a single scalar.
-! if you see this, get a more recent version of the nco tools.
-!
-! nsc 27jul2016
+!----------------------------------------------------------------------
+!> purpose: generate initial inflation files so an experiment starting up
+!> can always start from a restart file without having to alter the namelist
+!> between cycles 1 and 2.
+!>
+!> an alternative to running this program is to use the nco utilities thus:
+!>
+!> Here is an example using version 4.4.2 or later of the NCO tools:
+!>   ncap2 -s "T=1.0;U=1.0;V=1.0" wrfinput_d01 prior_inflation_mean.nc
+!>   ncap2 -s "T=0.6;U=0.6;V=0.6" wrfinput_d01 prior_inflation_sd.nc'
+!>
+!----------------------------------------------------------------------
 
 program fill_inflation_restart
 
-write(*,*) ''
-write(*,*) 'This program is OBSOLETE since DART inflation files are now in NetCDF format.' 
-write(*,*) 'To fill an initial inflation file use one of the standard NCO utilities like "ncap2" with'
-write(*,*) 'a copy of a model restart file to set the initial inflation mean, and a second file'
-write(*,*) 'for the initial inflation standard deviation.  Inflation mean and sd values now'
-write(*,*) 'are formatted in the files exactly like restart values, e.g. arranged by variable'
-write(*,*) 'type like T, U, V, etc.'
-write(*,*) ''
-write(*,*) 'Here is an example using version 4.4.2 or later of the NCO tools:'
-write(*,*) '  ncap2 -s "T=1.0;U=1.0;V=1.0" wrfinput_d01 prior_inf.nc'
-write(*,*) '  ncap2 -s "T=0.6;U=0.6;V=0.6" wrfinput_d01 prior_sd.nc'
-write(*,*) ''
+use             types_mod, only : r8, i8, missing_r8, metadatalength
 
+use         utilities_mod, only : register_module, error_handler, E_MSG, E_ERR, &
+                                  find_namelist_in_file, check_namelist_read,   &
+                                  initialize_utilities, finalize_utilities
+
+use       assim_model_mod, only : static_init_assim_model
+
+use      time_manager_mod, only : time_type, set_time, print_time, print_date, operator(-), &
+                                  get_calendar_type, NO_CALENDAR
+
+use  ensemble_manager_mod, only : init_ensemble_manager, ensemble_type, get_copy, &
+                                  map_pe_to_task, map_task_to_pe
+use   state_vector_io_mod, only : state_vector_io_init, read_state, write_state
+
+use   state_structure_mod, only : state_structure_info, get_model_variable_indices, &
+                                  get_num_domains, get_num_variables, get_io_num_dims, &
+                                  get_variable_name, get_dim_name, get_dim_length
+
+use      io_filenames_mod, only : io_filenames_init, file_info_type,       &
+                                  stage_metadata_type, get_stage_metadata, &
+                                  get_restart_filename,                    &
+                                  set_file_metadata, file_info_dump,       &
+                                  set_io_copy_flag, READ_COPY, WRITE_COPY
+
+use direct_netcdf_mod,     only : finalize_single_file_io, write_augmented_state
+
+use netcdf_utilities_mod,  only : nc_open_file_readwrite, nc_add_attribute_to_variable, &
+                                  nc_begin_define_mode, nc_end_define_mode,             &
+                                  nc_define_real_variable, nc_define_dimension,         &
+                                  nc_close_file, nc_put_variable, nc_variable_exists
+
+use             model_mod, only : static_init_model, get_model_size,       &
+                                  get_state_meta_data, model_interpolate
+
+use     mpi_utilities_mod, only : initialize_mpi_utilities, finalize_mpi_utilities
+
+use netcdf
+
+implicit none
+
+! version controlled file description for error handling, do not edit
+character(len=256), parameter :: source   = &
+   "$URL$"
+character(len=32 ), parameter :: revision = "$Revision$"
+character(len=128), parameter :: revdate  = "$Date$"
+
+