[Dart-dev] [3305] DART/trunk/ncep_obs/work: Minor fixes to the
scripting ( mostly in comments), and an
nancy at subversion.ucar.edu
nancy at subversion.ucar.edu
Wed Apr 9 16:26:55 MDT 2008
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080409/c730372f/attachment.html
-------------- next part --------------
Modified: DART/trunk/ncep_obs/work/create_obs_seq.csh
===================================================================
--- DART/trunk/ncep_obs/work/create_obs_seq.csh 2008-04-09 22:20:02 UTC (rev 3304)
+++ DART/trunk/ncep_obs/work/create_obs_seq.csh 2008-04-09 22:26:55 UTC (rev 3305)
@@ -1,4 +1,4 @@
-#!/bin/csh
+#!/bin/csh
#
# Data Assimilation Research Testbed -- DART
# Copyright 2004-2007, Data Assimilation Research Section
@@ -23,20 +23,19 @@
# you might need other options set for LSF on your system;
# e.g. -P project_charge_code, a different queue name, etc.
-#BSUB -o create_obs_seq.out
-#BSUB -e create_obs_seq.err
-#BSUB -J create_obs_seq
+set locallog = obs_seq.log
+#BSUB -J obsseq
+#BSUB -W 2:00
+
+#BSUB -o obs_seq.out
+#BSUB -e obs_seq.err
#BSUB -q share
-#BSUB -W 6:00
+#BSUB -P NNNNNNNN
#BSUB -n 1
-if ($?LS_SUBCWD) then
- cd $LS_SUBCWD
-endif
+touch $locallog
-touch create_obs_seq.out
-
set STRING = "1,$ s#,##g"
grep year input.nml >! ensstring.$$
@@ -46,6 +45,7 @@
grep month input.nml >! ensstring.$$
set ensstring = `sed -e "$STRING" ensstring.$$`
set mo = $ensstring[3]
+if ($mo < 10) set mo = 0$mo
grep day input.nml >! ensstring.$$
set ensstring = `sed -e "$STRING" ensstring.$$`
@@ -68,9 +68,11 @@
set ensstring = `sed -e "$STRING" ensstring.$$`
set daily_file = $ensstring[3]
-echo "year, mo, day1, dayn, binary = " $year $mo $day1 $dayn $binary >> create_obs_seq.out
-echo ObsBase = $ObsBase >> create_obs_seq.out
+rm -f ensstring.$$
+echo "year, mo, day1, dayn, binary = " $year $mo $day1 $dayn $binary >> $locallog
+echo ObsBase = $ObsBase >> $locallog
+
./create_real_obs
set n = $day1
@@ -82,22 +84,27 @@
if ($daily_file == .true.) then
set obs_seq = ( ${base} )
else
- # originally this script created 2 files per day, 12 and 24. changes
- # in the create_real_obs executable have it creating 4 per day now.
- # if you change it back, change the next line to have only 2 files.
+ # originally, if not doing daily files, the default was 2 files per day.
+ # now the default matches the ncep data file blocking; 4 files per day,
+ # 6 hours per file. if you want to go back to 2 files per day, you will
+ # have to modify the following line to have only 12 and 24.
set obs_seq = ( ${base}06 ${base}12 ${base}18 ${base}24 )
endif
foreach obs ($obs_seq)
if ($binary == .false.) then
if (-e $obs) then
- echo "fixing $obs pole locations" >> create_obs_seq.out
- sed -e 's/ 1.57079632679490/ 1.57079632679488/' $obs >! fixed_pole
+ # on at least one platform (ibm power5), this value is read in
+ # and rounded to be 1 bit in the least significant digit larger
+ # than 90.000 in degrees, which causes errors in the location module.
+ # drop these values just slightly so they read in as 90/-90 exactly.
+ echo " fixing $obs pole locations" >> $locallog
+ sed -e 's/ 1.57079632679490/ 1.57079632679488/' \
+ -e 's/-1.57079632679490/-1.57079632679488/' $obs >! fixed_pole${n}
else
exit
endif
- mv fixed_pole $obs
+ mv fixed_pole${n} $obs
endif
- mv $obs $Obs_base:h &
end
@ n++
end
Modified: DART/trunk/ncep_obs/work/input.nml
===================================================================
--- DART/trunk/ncep_obs/work/input.nml 2008-04-09 22:20:02 UTC (rev 3304)
+++ DART/trunk/ncep_obs/work/input.nml 2008-04-09 22:26:55 UTC (rev 3305)
@@ -1,55 +1,45 @@
-
&ncepobs_nml
- year = 2006,
- month = 1,
- day = 1,
- tot_days = 1,
- max_num = 1000000,
- ObsBase = '../prep_bufr/work/temp_obs.',
- select_obs = 1,
- ADPUPA = .true.,
- AIRCAR = .true.,
- AIRCFT = .true.,
- SATEMP = .false.,
- SFCSHP = .false.,
- ADPSFC = .false.,
- SATWND = .true.,
- obs_U = .true.,
- obs_V = .true.,
- obs_T = .true.,
+ year = 2007,
+ month = 1,
+ day = 1,
+ tot_days = 31,
+ max_num = 800000,
+ select_obs = 0,
+ ObsBase = '../prepbufr/data/temp_obs.',
+ ADPUPA = .false.,
+ AIRCAR = .false.,
+ AIRCFT = .false.,
+ SATEMP = .false.,
+ SFCSHP = .false.,
+ ADPSFC = .false.,
+ SATWND = .false.,
+ obs_U = .false.,
+ obs_V = .false.,
+ obs_T = .false.,
obs_PS = .false.,
- obs_QV = .true.,
- daily_file = .false.,
- lat1 = -89.0,
- lat2 = 89.0,
- lon1 = 0.0,
- lat2 = 360.0 /
+ obs_QV = .false.,
+ daily_file = .true.,
+ lon1 = 0.0,
+ lon2 = 360.0,
+ lat1 = -90.0,
+ lat2 = 90.0 /
&obs_sequence_nml
write_binary_obs_sequence = .false. /
-&obs_kind_nml
- /
-
&assim_model_nml
write_binary_restart_files = .true. /
-&model_nml
- /
-
-&location_nml
- /
-
&utilities_nml
TERMLEVEL = 1,
logfilename = 'dart_log.out' /
-&preprocess_nml
- input_obs_kind_mod_file = '../../obs_kind/DEFAULT_obs_kind_mod.F90',
- output_obs_kind_mod_file = '../../obs_kind/obs_kind_mod.f90',
- input_obs_def_mod_file = '../../obs_def/DEFAULT_obs_def_mod.F90',
- output_obs_def_mod_file = '../../obs_def/obs_def_mod.f90',
- input_files = '../../obs_def/obs_def_reanalysis_bufr_mod.f90',
- '../../obs_def/obs_def_gps_mod.f90',
- '../../obs_def/obs_def_eval_mod.f90' /
+&model_nml
+ /
+&location_nml
+ /
+
+&obs_kind_nml
+ /
+
More information about the Dart-dev
mailing list