Commit 02b02cff authored by Carl Schreck's avatar Carl Schreck

Automated Nightly Commit - Thu May 9 00:01:20 EDT 2019

parent a9dc7def
......@@ -7,12 +7,12 @@
QUEUE=allnodes
echo $QUEUE
PY_SCRIPT=subset_interim
PY_SCRIPT=extract_interim
PY_DIR=`pwd`
LOG_DIR=$PY_DIR/log
mkdir -p $LOG_DIR
#for YEAR in {1979..2018}; do
for YEAR in {2000..2018}; do
JOB_NAME=$PY_SCRIPT"_"$YEAR
echo $JOB_NAME `date`
LOG_FILE=$LOG_DIR/$JOB_NAME.log
......@@ -30,5 +30,5 @@ mkdir -p $LOG_DIR
python -u $PY_DIR/$PY_SCRIPT.py $YEAR
#done
done
......@@ -28,7 +28,8 @@ for tar_file in tar_files:
tar.extractall(basedir+year)
cjs.tstamp('Open')
grb_files = glob.glob(basedir + year + '/*' + year + '*schreck368330')
grb_files = glob.glob(basedir + year + '/ei.oper.an.pl.regn128uv.' + year
+ '*schreck368330')
grb_files.sort()
ds = xr.open_mfdataset(grb_files, engine='cfgrib', concat_dim='time')
......
......@@ -7,7 +7,7 @@ fi
echo $QUEUE
NCL_SCRIPT=filter_waves
VARNAME="chi200"
VARNAME="u850"
FILTNAME="k09"
NCL_DIR=`pwd`
LOG_DIR=$NCL_DIR/log
......
......@@ -29,11 +29,14 @@ begin
ymin = -20
ymax = 20
timeUnits = "days since 1800-01-01 00:00:00"
tmin = cd_inv_calendar( 2017, 01, 01, 00, 0, 0, timeUnits, 0 )
tmax = cd_inv_calendar( 2018, 12, 31, 00, 0, 0, timeUnits, 0 )
if( .not.isvar("filtName") ) then
filtName = "wide"
end if
obsPerDay = 1
obsPerDay = 4
if( .not.isvar("varName") ) then
varName = "u850"
end if
......@@ -50,19 +53,19 @@ begin
pathOut = basePath + "waves/" + varName + ".anom.waves.nc"
calcHigh = False
calcLow = False
calcLow = True
calcMrg = False
calcMjo = False
calcMjo = True
calcKelvin = True
calcEr = False
calcEr = True
calcMtd = False
calcTd = False
calcWig = False
makeNewFile = False
makeNewFile = True
; Open the input files
fin = addfile( pathIn, "r" )
time = fin->time
time = fin->time({tmin:tmax})
lat = fin->lat({ymin:ymax})
lon = fin->lon
......@@ -84,7 +87,7 @@ begin
print_clock( "Reading the input data... " )
data = fin->$varName$(:,{ymin:ymax},:)
data = fin->$varName$({tmin:tmax},{ymin:ymax},:)
data&time@beginning_date = cd_string( data&time(0), "" )
data&time@ending_date = cd_string( data&time(dimsizes(data&time)-1), "" )
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment