[Dart-dev] [6851] DART/trunk/models/pe2lyr/work: add the build files to make obs_seq_to_netcdf for this model.
nancy at ucar.edu
nancy at ucar.edu
Fri Mar 14 11:11:55 MDT 2014
Revision: 6851
Author: nancy
Date: 2014-03-14 11:11:53 -0600 (Fri, 14 Mar 2014)
Log Message:
-----------
add the build files to make obs_seq_to_netcdf for this model.
update the input.nml to add the namelists needed for obs_seq_to_netcdf,
and make some of the defaults better for actually running the model
(e.g. outlier_threshold, cutoff). also note where our tools need
minor tweaks for running with a model that isn't using a calendar.
Modified Paths:
--------------
DART/trunk/models/pe2lyr/work/input.nml
Added Paths:
-----------
DART/trunk/models/pe2lyr/work/mkmf_obs_seq_to_netcdf
DART/trunk/models/pe2lyr/work/path_names_obs_seq_to_netcdf
-------------- next part --------------
Modified: DART/trunk/models/pe2lyr/work/input.nml
===================================================================
--- DART/trunk/models/pe2lyr/work/input.nml 2014-03-13 17:20:22 UTC (rev 6850)
+++ DART/trunk/models/pe2lyr/work/input.nml 2014-03-14 17:11:53 UTC (rev 6851)
@@ -42,7 +42,7 @@
output_interval = 1,
num_groups = 1,
input_qc_threshold = 3.0,
- outlier_threshold = -1.0,
+ outlier_threshold = 4.0,
output_forward_op_errors = .true.,
output_timestamps = .false.,
output_inflation = .true.,
@@ -79,7 +79,7 @@
&assim_tools_nml
filter_kind = 1,
- cutoff = 100000.0,
+ cutoff = 0.02,
sort_obs_inc = .false.,
spread_restoration = .false.,
sampling_error_correction = .false.,
@@ -119,6 +119,7 @@
approximate_distance = .false.,
nlon = 71,
nlat = 36,
+ print_box_level = 0,
output_box_info = .false. /
&preprocess_nml
@@ -136,8 +137,7 @@
write_nml = 'file' /
&obs_sequence_tool_nml
- num_input_files = 2,
- filename_seq = 'obs_seq.one', 'obs_seq.two',
+ filename_seq = 'obs_seq.out',
filename_out = 'obs_seq.processed',
first_obs_days = -1,
first_obs_seconds = -1,
@@ -194,6 +194,10 @@
gregorian_cal = .false.
/
+# again - this model doesn't use a calendar. it starts at
+# day 0, seconds 0. that's equiv to year 1601 in gregorian,
+# so use that for now. but it should allow a no_calendar option.
+#
# The times in the namelist for the obs_diag program are vectors
# that follow the following sequence:
# year month day hour minute second
@@ -202,23 +206,29 @@
#
# Acceptable latitudes range from [-90, 90]
# Acceptable longitudes range from [ 0, Inf]
+#
+# currently this fails anyway because the observations are called
+# U_WIND_COMPONENT and V_WIND_COMPONENT and in trying to match up
+# corresponding wind observations into a single vector, the obs_diag
+# program expects a platform name at the start of the obs type and
+# fails without one.
&obs_diag_nml
obs_sequence_name = 'obs_seq.final',
obs_sequence_list = '',
- first_bin_center = 2003, 1, 1, 0, 0, 0 ,
- last_bin_center = 2003, 1, 7, 0, 0, 0 ,
- bin_separation = 0, 0, 0, 6, 0, 0 ,
- bin_width = 0, 0, 0, 6, 0, 0 ,
- time_to_skip = 0, 0, 1, 0, 0, 0 ,
+ first_bin_center = 1601, 1, 1, 0, 0, 0 ,
+ last_bin_center = 9999, 1, 1, 0, 0, 0 ,
+ bin_separation = 0, 0, 0, 0, 0, 0 ,
+ bin_width = 0, 0, 1, 0, 0, 0 ,
+ time_to_skip = 0, 0, 0, 0, 0, 0 ,
max_num_bins = 1000,
trusted_obs = 'null',
- Nregions = 4,
- lonlim1 = 0.0, 0.0, 0.0, 235.0,
- lonlim2 = 360.0, 360.0, 360.0, 295.0,
- latlim1 = 20.0, -80.0, -20.0, 25.0,
- latlim2 = 80.0, -20.0, 20.0, 55.0,
- reg_names = 'Northern Hemisphere', 'Southern Hemisphere', 'Tropics', 'North America',
+ Nregions = 3,
+ lonlim1 = 0.0, 0.0, 0.0,
+ lonlim2 = 360.0, 360.0, 360.0,
+ latlim1 = 20.0, -80.0, -20.0,
+ latlim2 = 80.0, -20.0, 20.0,
+ reg_names = 'Northern Hemisphere', 'Southern Hemisphere', 'Tropics',
print_mismatched_locs = .false.,
create_rank_histogram = .true.,
outliers_in_histogram = .true.,
@@ -226,3 +236,42 @@
verbose = .true.
/
+# this is what the namelist should probably look like, because
+# this model doesn't really have a calendar. it just increments
+# days and seconds. but the schedule module currently only
+# understands models with calendars, so lie to it about using
+# a gregorian calendar. it works just fine - day 0 is year 1601.
+#&schedule_nml
+# calendar = 'no_calendar',
+# first_bin_start = 0, 0,
+# first_bin_end = 9999, 0,
+# last_bin_end = 9999, 0,
+# bin_interval_days = 1000000,
+# bin_interval_seconds = 0,
+# max_num_bins = 1000,
+# print_table = .true.
+# /
+
+&schedule_nml
+ calendar = 'Gregorian',
+ first_bin_start = 1601, 1, 1, 0, 0, 0,
+ first_bin_end = 2999, 1, 1, 0, 0, 0,
+ last_bin_end = 2999, 1, 1, 0, 0, 0,
+ bin_interval_days = 1000000,
+ bin_interval_seconds = 0,
+ max_num_bins = 1000,
+ print_table = .true.
+ /
+
+
+&obs_seq_to_netcdf_nml
+ obs_sequence_name = 'obs_seq.final',
+ obs_sequence_list = '',
+ append_to_netcdf = .false.,
+ lonlim1 = 0.0,
+ lonlim2 = 360.0,
+ latlim1 = -90.0,
+ latlim2 = 90.0,
+ verbose = .false.
+ /
+
Added: DART/trunk/models/pe2lyr/work/mkmf_obs_seq_to_netcdf
===================================================================
--- DART/trunk/models/pe2lyr/work/mkmf_obs_seq_to_netcdf (rev 0)
+++ DART/trunk/models/pe2lyr/work/mkmf_obs_seq_to_netcdf 2014-03-14 17:11:53 UTC (rev 6851)
@@ -0,0 +1,18 @@
+#!/bin/csh
+#
+# 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$
+
+../../../mkmf/mkmf -p obs_seq_to_netcdf -t ../../../mkmf/mkmf.template \
+ -a "../../.." path_names_obs_seq_to_netcdf
+
+exit $status
+
+# <next few lines under version control, do not edit>
+# $URL$
+# $Revision$
+# $Date$
+
Property changes on: DART/trunk/models/pe2lyr/work/mkmf_obs_seq_to_netcdf
___________________________________________________________________
Added: svn:executable
+ *
Added: svn:mime-type
+ text/plain
Added: svn:keywords
+ Date Rev Author HeadURL Id
Added: svn:eol-style
+ native
Added: DART/trunk/models/pe2lyr/work/path_names_obs_seq_to_netcdf
===================================================================
--- DART/trunk/models/pe2lyr/work/path_names_obs_seq_to_netcdf (rev 0)
+++ DART/trunk/models/pe2lyr/work/path_names_obs_seq_to_netcdf 2014-03-14 17:11:53 UTC (rev 6851)
@@ -0,0 +1,15 @@
+assim_model/assim_model_mod.f90
+common/types_mod.f90
+location/threed_sphere/location_mod.f90
+models/pe2lyr/model_mod.f90
+models/pe2lyr/pe2lyr_mod.f90
+models/pe2lyr/spharmt_mod.f90
+mpi_utilities/null_mpi_utilities_mod.f90
+obs_def/obs_def_mod.f90
+obs_kind/obs_kind_mod.f90
+obs_sequence/obs_seq_to_netcdf.f90
+obs_sequence/obs_sequence_mod.f90
+random_seq/random_seq_mod.f90
+time_manager/schedule_mod.f90
+time_manager/time_manager_mod.f90
+utilities/utilities_mod.f90
More information about the Dart-dev
mailing list