[Dart-dev] [3298] DART/trunk/ncep_obs/prep_bufr/work:
Updated versions which include the right namelist
nancy at subversion.ucar.edu
nancy at subversion.ucar.edu
Wed Apr 9 15:39:11 MDT 2008
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20080409/e20d8cd0/attachment-0001.html
-------------- next part --------------
Modified: DART/trunk/ncep_obs/prep_bufr/work/input.nml
===================================================================
--- DART/trunk/ncep_obs/prep_bufr/work/input.nml 2008-04-09 21:09:52 UTC (rev 3297)
+++ DART/trunk/ncep_obs/prep_bufr/work/input.nml 2008-04-09 21:39:10 UTC (rev 3298)
@@ -1,41 +1,27 @@
&prep_bufr_nml
obs_window = 24.0,
obs_window_cw = 24.0,
-/
+ qctype_use = 1, 2, 3, 4 /
+! the above are suggested settings. the full set of
+! namelist items (and defaults) are here:
-! obs_window = 0.80,
-! obs_window = 1.50,
+! otype_use (= all) ! report types to use
+! obs_window = 0.8, ! observation time window (hours)
+! obs_window_cw = 1.5, ! cloud wind observation time window (hours)
+! qctype_use = 1,2,3,4,9,15 ! NCEP QC values to accept
+! land_temp_error = 2.5, ! assumed error for surface temp. observations (K)
+! land_wind_error = 3.5, ! assumed error for surface wind observations (m/s)
+! land_moist_error = 0.2 ! assumed error for surface moist. observations (%)
+
+! other examples of use:
+! obs_window = 2.00,
+! obs_window = 1.00,
+! otype_use = 120.0, 180.0, 181.0, 182.0,
! otype_use = 120.0, 130.0, 131.0, 132.0, 133.0, 180.0,
! 181.0, 182.0, 220.0, 221.0, 230.0, 231.0,
! 232.0, 233.0, 242.0, 243.0, 245.0, 246.0,
! 252.0, 253.0, 255.0, 280.0, 281.0, 282.0,
-! qctype_use = 1, 2, 3, 4, 9, 15 /
+! qctype_use = 1, 2, 3, 4,
+!
-&ncepobs_nml
- year = 2006,
- month = 08,
- day = 1,
- tot_days = 31,
- max_num = 1000000,
- ObsBase = 'temp_obs.',
- select_obs = 0,
- ADPUPA = .true.,
- AIRCAR = .true.,
- AIRCFT = .true.,
- SATEMP = .true.,
- SFCSHP = .true.,
- ADPSFC = .true.,
- SATWND = .true.,
- obs_U = .true.,
- obs_V = .true.,
- obs_T = .true.,
- obs_PS = .true.,
- obs_QV = .true.,
- daily_file = .true.,
- lat1 = -89.99,
- lat2 = 89.99,
- lon1 = 0.0,
- lat2 = 360.0 /
-
-
Modified: DART/trunk/ncep_obs/prep_bufr/work/prepbufr.csh
===================================================================
--- DART/trunk/ncep_obs/prep_bufr/work/prepbufr.csh 2008-04-09 21:09:52 UTC (rev 3297)
+++ DART/trunk/ncep_obs/prep_bufr/work/prepbufr.csh 2008-04-09 21:39:10 UTC (rev 3298)
@@ -1,4 +1,4 @@
-#!/bin/csh
+#!/bin/csh
#
# Data Assimilation Research Testbed -- DART
# Copyright 2004-2007, Data Assimilation Research Section
@@ -26,14 +26,11 @@
#BSUB -e prepbufr.err
#BSUB -J prepbufr
#BSUB -q share
-#BSUB -W 6:00
+#BSUB -W 2:00
+#BSUB -P NNNNNNNN
#BSUB -n 1
-if ($?LS_SUBCWD) then
- cd $LS_SUBCWD
-endif
-
#
#--------------------------------------------------------------
#
@@ -48,22 +45,23 @@
# Convert from big-endian BUFR files to little-endian for Intel chip systems.
# ('yes' or whatever)
-set daily = no
-set convert = yes
-set year = 2006
-set month = 1
+set daily = yes
+set convert = no
+set year = 1988
+set month = 12
set beginday = 1
#
-# end day (up to and including the last day of the month. Leap year Februaries are OK.
-# Remember that the prepqm###### file for hour 0 of the first day of the next
-# month is necessary for endday = last day of a month.)
+# end day (up to and including the last day of the month.
+# Leap year Februaries are OK.
+# Remember that the prepqm###### file for hour 0 of the first day of the next
+# month is necessary for endday = last day of a month.)
#
-set endday = 1
+set endday = 31
# Location of BUFR files (named prepqmYYYYMMDDHH)
# are assumed to be in subdirectories named YYYYMM of the path listed here.
# Those subdirectory names will be constructed below.
-set BUFR_dir = ../data/
+set BUFR_dir = 198812
set get_year = $year
# END USER SET PARAMETERS
@@ -75,12 +73,8 @@
if (($year % 100) == 0) @ days_in_mo[2] = $days_in_mo[2] - 1
if (($year % 400) == 0) @ days_in_mo[2] = $days_in_mo[2] + 1
-if ( $?LS_SUBCWD ) then
- cd $LS_SUBCWD
-endif
+rm -f prepqm.out *.err *.out
-rm prepqm.out temp_obs *.err *.out
-
# Loop over days
set day = $beginday
@@ -89,11 +83,11 @@
echo '-------------------------------------- '
# clear any old intermediate (text) BUFR file
- rm temp_obs
+ rm -f temp_obs
- # convert 1 "day"s worth (data from '6z to 6z of the next day) of BUFR files
- # into a single intermediate file if 'daily' is set to true; 4 6-hour files
- # otherwise.
+ # convert 1 "day"s worth (data from '3z to 3z of the next day) of BUFR files
+ # into a single intermediate file if 'daily' is set to true; 4 6-hour files
+ # otherwise.
set h = 0
set next_day = not
while ($h < 30)
@@ -114,9 +108,11 @@
set dd = 1
@ mm++
if ($mm > 12) then
+ if ($mm > 12 && $hh == 0) then
+ @ get_year ++
+ endif
# next year
set mm = 1
- @ get_year = $get_year++
@ yy = $get_year % 100
endif
endif
@@ -128,7 +124,8 @@
if ($dd < 10) set dd = 0$dd
if ($hh < 10) set hh = 0$hh
- # link(big endian) or make(little endian) input file 'prepqm' for prepbufr.x
+ # link(big endian) or make(little endian) input file 'prepqm'
+ # for prepbufr.x
set BUFR_loc = ${BUFR_dir}/${get_year}${mm}
if (! -e ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh}) then
echo "MISSING FILE ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh} and aborting"
@@ -143,8 +140,8 @@
mv prepqm.littleendian prepqm.in
rm prepqm.bigendian
else
- echo "copying prepqm.in to ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh}"
- cp ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh} prepqm.in
+ echo "linking prepqm.in to ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh}"
+ ln -f ${BUFR_loc}/prepqm${yy}${mm}${dd}${hh} prepqm.in
endif
if ($h == 30) then
More information about the Dart-dev
mailing list