Commit a9dc7def authored by Carl Schreck's avatar Carl Schreck

Automated Nightly Commit - Wed May 8 00:01:19 EDT 2019

parent 590c2629
......@@ -7,12 +7,12 @@
QUEUE=allnodes
echo $QUEUE
PY_SCRIPT=extract_interim
PY_SCRIPT=subset_interim
PY_DIR=`pwd`
LOG_DIR=$PY_DIR/log
mkdir -p $LOG_DIR
for YEAR in {1979..2018}; do
#for YEAR in {1979..2018}; do
JOB_NAME=$PY_SCRIPT"_"$YEAR
echo $JOB_NAME `date`
LOG_FILE=$LOG_DIR/$JOB_NAME.log
......@@ -30,5 +30,5 @@ for YEAR in {1979..2018}; do
python -u $PY_DIR/$PY_SCRIPT.py $YEAR
done
#done
......@@ -36,7 +36,7 @@ cjs.tstamp('Read/Write')
ds.time.encoding['units'] = 'days since 1800-01-01'
ds.time.encoding['dtype'] = 'float64'
#ds.encoding['_FillValue'] = -9999.
ds.to_netcdf(basedir + 'u850.' + year + '.nc' )
ds.to_netcdf(basedir + 'u850.' + year + '.nc')
cjs.tstamp('Closing')
ds.close()
......
......@@ -16,16 +16,19 @@ var_name = 'u850'
basedir = '/home/carl/data/era-interim/'
cjs.tstamp('Open')
in_files = glob.glob(basedir + '/' + var_name + '.201*.nc')
in_files = glob.glob(basedir + 'incoming/' + var_name + '.*.nc')
in_files.sort()
ds = xr.open_mfdataset(in_files, concat_dim='time')
cjs.tstamp('Read/Write')
ds.encoding['_FillValue'] = -9999.
ds.to_netcdf(basedir + 'total/' + var_name + '.nc' )
cjs.tstamp('Closing')
ds.close()
ds.time.encoding['units'] = 'days since 1800-01-01'
ds.time.encoding['dtype'] = 'float64'
ds = ds.rename({'u':'u850', 'latitude':'lat', 'longitude':'lon'})
ds.sel(lat=slice(15,-15)).to_netcdf(basedir + 'total/' + var_name
+ '.nc')
#cjs.tstamp('Closing')
#ds.close()
cjs.tstamp('Thank you, come again.')
......@@ -26,28 +26,28 @@ begin
; These are some parameters that could be useful to have up top
if( .not.isvar("varName") ) then
varName = "chi200"
varName = "u850"
end if
standardizing = False
basePath = "~/data/cfsr/"
basePath = "~/data/era-interim/"
; basePath = "~/data/olr/current/"
; basePath = "/home/carl/data/nclimgrid/beta/merged/"
if( True ) then
pathIn = basePath + "total/" + varName + ".total.nc"
pathClim = basePath + "clim/" + varName + ".clim.1981.2010.nc"
pathIn = basePath + "total/" + varName + ".nc"
pathClim = basePath + "clim/" + varName + ".clim.1988.2017.nc"
pathAnom = basePath + "anom/" + varName + ".anom.nc"
pathStd = basePath + "std/" + varName + ".std.nc"
else
pathIn = basePath + varName + ".total.nc"
pathClim = basePath + varName + ".clim.1981.2010.nc"
pathClim = basePath + varName + ".clim.1988.2017.nc"
pathAnom = basePath + varName + ".anom.nc"
pathStd = basePath + varName + ".std.nc"
end if
climStart = 1981001
climEnd = 2010366
climStart = 1988001
climEnd = 2017366
; Open the input files
fin = addfile( pathIn, "r" )
......
......@@ -14,7 +14,8 @@ mkdir -p $LOG_DIR
EMPTY=\"\"
#for VAR_NAME in uShear vShear shear pmsl g500
for VAR_NAME in uwnd925 uShear vShear shear
#for VAR_NAME in uwnd925 uShear vShear shear
for VAR_NAME in u850
do
# if [[ $VAR_NAME = "vShear" || $VAR_NAME = "uShear" || $VAR_NAME = "shear" ]]; then
ALL_LEVELS=( $EMPTY )
......
......@@ -35,30 +35,30 @@ begin
end if
obsPerDay = 1
if( .not.isvar("varName") ) then
varName = "rain"
varName = "u850"
end if
print( varName )
; basePath = "~/data/cfsr/"
basePath = "~/data/nasa_ymc/fcst_verif/"
basePath = "~/data/era-interim/"
; basePath = "~/data/nasa_ymc/fcst_verif/"
; basePath = "~/data/olr/current/"
; pathIn = basePath + ".std.nc"
; pathOut = basePath + ".waves.std.nc"
pathIn = basePath + varName + ".anom.nc"
pathOut = basePath + varName + "." + filtName + ".nc"
; pathIn = basePath + "anom/" + varName + ".anom.nc"
; pathOut = basePath + "waves/" + varName + ".anom.waves.nc"
; pathIn = basePath + varName + ".anom.nc"
; pathOut = basePath + varName + "." + filtName + ".nc"
pathIn = basePath + "anom/" + varName + ".anom.nc"
pathOut = basePath + "waves/" + varName + ".anom.waves.nc"
calcHigh = True
calcLow = True
calcHigh = False
calcLow = False
calcMrg = False
calcMjo = True
calcMjo = False
calcKelvin = True
calcEr = True
calcEr = False
calcMtd = False
calcTd = False
calcWig = False
makeNewFile = True
makeNewFile = False
; Open the input files
fin = addfile( pathIn, "r" )
......
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