Commit 1513408a authored by Carl Schreck's avatar Carl Schreck

Automated Nightly Commit - Wed Jan 30 00:00:22 EST 2019

parent 3201e46c
......@@ -19,29 +19,29 @@ load "$CJS_NCL_LIB/print_clock.ncl"
begin
print_clock( "Here we go!" )
basePath = "/home/carl/data/state_of_climate/"
basePath = "~/data/ncep1/total/"
print_clock( "Reading uLow" )
uLowVar = "uwnd850"
uLowPath = basePath + uLowVar + ".total.notc500.nc"
uLowPath = basePath + uLowVar + ".total.nc"
uLowFile = addfile( uLowPath, "r" )
uLow = uLowFile->$uLowVar$
print_clock( "Reading vLow" )
vLowVar = "vwnd850"
vLowPath = basePath + vLowVar + ".total.notc500.nc"
vLowPath = basePath + vLowVar + ".total.nc"
vLowFile = addfile( vLowPath, "r" )
vLow = vLowFile->$vLowVar$
print_clock( "Reading uHigh" )
uHighVar = "uwnd200"
uHighPath = basePath + uHighVar + ".total.notc500.nc"
uHighPath = basePath + uHighVar + ".total.nc"
uHighFile = addfile( uHighPath, "r" )
uHigh = uHighFile->$uHighVar$
print_clock( "Reading vHigh" )
vHighVar = "vwnd200"
vHighPath = basePath + vHighVar + ".total.notc500.nc"
vHighPath = basePath + vHighVar + ".total.nc"
vHighFile = addfile( vHighPath, "r" )
vHigh = vHighFile->$vHighVar$
......@@ -65,9 +65,9 @@ begin
print_clock( "Writing" )
setfileoption( "nc", "Format", "LargeFile" )
magPath = basePath + "shear.total.notc500.nc"
uShearPath = basePath + "uShear.total.notc500.nc"
vShearPath = basePath + "vShear.total.notc500.nc"
magPath = basePath + "shear.total.nc"
uShearPath = basePath + "uShear.total.nc"
vShearPath = basePath + "vShear.total.nc"
system( "rm " + magPath )
outFile = addfile(magPath,"c")
......
......@@ -20,18 +20,19 @@ load "$CJS_NCL_LIB/print_clock.ncl"
begin
print_clock( "Here we go! " )
setfileoption("nc","Format","NetCDF4")
setfileoption("nc","CompressionLevel",1)
; setfileoption("nc","Format","NetCDF4")
; setfileoption("nc","CompressionLevel",1)
setfileoption("nc","Format","LargeFile")
; These are some parameters that could be useful to have up top
if( .not.isvar("varName") ) then
varName = "t850"
varName = "skt"
end if
standardizing = False
; basePath = "/home/carl/data/cfsr/"
basePath = "/home/carl/data/nclimgrid/beta/merged/"
if( False ) then
basePath = "~/data/ncep1/"
; basePath = "/home/carl/data/nclimgrid/beta/merged/"
if( True ) then
pathIn = basePath + "total/" + varName + ".total.nc"
pathClim = basePath + "clim/" + varName + ".clim.1981.2010.nc"
pathAnom = basePath + "anom/" + varName + ".anom.nc"
......
......@@ -14,7 +14,7 @@ mkdir -p $LOG_DIR
EMPTY=\"\"
#for VAR_NAME in uShear vShear shear pmsl g500
for VAR_NAME in prcp tmin tmax tavg
for VAR_NAME in uwnd925 uShear vShear shear
do
# if [[ $VAR_NAME = "vShear" || $VAR_NAME = "uShear" || $VAR_NAME = "shear" ]]; then
ALL_LEVELS=( $EMPTY )
......
......@@ -21,23 +21,24 @@ begin
print_clock( "Here we go! " )
; These are some parameters that could be useful to have up top
varName = "skt"
varName = "sst"
newFile = True
doCalcAnom = True
; basePath = "~/data/cfsr/monthly/" + varName
; basePath = "~/data/ersst/v5/sst.mnmean"
basePath = "~/data/ersst/v5/sst.mnmean"
; basePath = "~/data/oisst/oisst.month"
; inPath = basePath + ".total.nc"
; climPath = basePath + ".clim.1981.2010.nc"
; anomPath = basePath + ".anom.nc"
; stdPath = basePath + ".std.nc"
basePath = "~/data/ncep1/monthly/"
inPath = basePath + varName + ".mon" + ".mean.nc"
climPath = basePath + varName + ".mon" + ".clim.nc"
anomPath = basePath + varName + ".mon" + ".anom.nc"
stdPath = basePath + varName + ".mon" + ".std.nc"
inPath = basePath + ".nc"
climPath = basePath + ".clim.1981.2010.nc"
anomPath = basePath + ".anom.nc"
stdPath = basePath + ".std.nc"
; basePath = "~/data/ncep1/monthly/"
; inPath = basePath + varName + ".mon" + ".mean.nc"
; climPath = basePath + varName + ".mon" + ".clim.nc"
; anomPath = basePath + varName + ".mon" + ".anom.nc"
; stdPath = basePath + varName + ".mon" + ".std.nc"
climStart = 1981001
climEnd = 2010366
......
......@@ -25,10 +25,10 @@ begin
; These are some parameters that could be useful to have up top
timeUnits = "days since 1800-01-01 00:00:00"
if( .not.isvar("inVar") ) then
inVar = (/ "uwnd" /)
inVar = (/ "skt" /)
end if
if( .not.isvar("levels") ) then
levels = (/ 850 /)
levels = (/ -1 /)
end if
if( .not.isvar("version") ) then
version = "ncep1"
......
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