Commit 6816a49a authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Fri Jul 15 00:00:12 EDT 2016

parent 2cb3603c
......@@ -6,7 +6,7 @@ else
fi
echo $QUEUE
NCL_SCRIPT=wave_meanvar
NCL_SCRIPT=wk99_diff
NCL_DIR=`pwd`
LOG_DIR=$NCL_DIR/log
mkdir -p $LOG_DIR
......
......@@ -95,6 +95,9 @@ begin
res@cnMissingValFillPattern = "SolidFill"
res@cnMissingValFillColor = "gray"
res@tiXAxisString = "HIRS"
res@tiYAxisString = "AVHRR"
; ...set the contour levels
if( shdLevels(0).eq.-999 ) then
res@cnMaxLevelCount = 8
......
......@@ -18,56 +18,76 @@ begin ; main
; These are some parameters that could be useful to have up top
plotType = "png"
plotName = "figures/wk99_diff"
plotDpi = 200
diffNotTrend = True
nTests = 1000
if( diffNotTrend ) then
plotName = "figures/wk99_diff"
aPath = "/home/carl/data/olr/compare/hirs/olr.wk99.nc"
bPath = "/home/carl/data/olr/compare/avhrr/olr.wk99.nc"
else
plotName = "figures/wk99_trend"
aPath = "/home/carl/data/olr/compare/hirs/olr.wk99.1980s.nc"
bPath = "/home/carl/data/olr/compare/hirs/olr.wk99.2000s.nc"
end if
fontHeightF = 0.02
nTests = 1000
pThresh = 0.00
calcDims = (/ 0, 1 /)
print_clock( "Reading hirs" )
inPath = "/home/carl/data/olr/compare/hirs/olr.wk99.nc"
inFile = addfile( inPath, "r" )
hirsSymLog = inFile->symLogPower(freq|:,{k|-20:20})
hirsAntiLog = inFile->antiLogPower(freq|:,{k|-20:20})
hirsSymNorm = inFile->symNormPower(freq|:,{k|-20:20})
hirsAntiNorm = inFile->antiNormPower(freq|:,{k|-20:20})
print_clock( "Reading A" )
aFile = addfile( aPath, "r" )
aSymLog = aFile->symLogPower(freq|:,{k|-20:20})
aAntiLog = aFile->antiLogPower(freq|:,{k|-20:20})
aSymNorm = aFile->symNormPower(freq|:,{k|-20:20})
aAntiNorm = aFile->antiNormPower(freq|:,{k|-20:20})
hirsSymAll = inFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
hirsAntiAll = inFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
hirsBack = conform( hirsSymAll, inFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
aSymAll = aFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
aAntiAll = aFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
aBack = conform( aSymAll, aFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
print_clock( "Reading avhrr" )
inPath = "/home/carl/data/olr/compare/avhrr/olr.wk99.nc"
inFile = addfile( inPath, "r" )
avhrrSymLog = inFile->symLogPower(freq|:,{k|-20:20})
avhrrAntiLog = inFile->antiLogPower(freq|:,{k|-20:20})
avhrrSymNorm = inFile->symNormPower(freq|:,{k|-20:20})
avhrrAntiNorm = inFile->antiNormPower(freq|:,{k|-20:20})
print_clock( "Reading B" )
bFile = addfile( bPath, "r" )
bSymLog = bFile->symLogPower(freq|:,{k|-20:20})
bAntiLog = bFile->antiLogPower(freq|:,{k|-20:20})
bSymNorm = bFile->symNormPower(freq|:,{k|-20:20})
bAntiNorm = bFile->antiNormPower(freq|:,{k|-20:20})
avhrrSymAll = inFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
avhrrAntiAll = inFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
avhrrBack = conform( avhrrSymAll, inFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
bSymAll = bFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
bAntiAll = bFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
bBack = conform( bSymAll, bFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
diffSym = bootstrap_diff( hirsSymAll, avhrrSymAll, \
diffSym = bootstrap_diff( aSymAll, bSymAll, \
nTests, pThresh, calcDims )
printMinMax( diffSym, True )
diffAnti = bootstrap_diff( hirsAntiAll, avhrrAntiAll, \
diffAnti = bootstrap_diff( aAntiAll, bAntiAll, \
nTests, pThresh, calcDims )
printMinMax( diffAnti, True )
diffSymLog = hirsSymLog - avhrrSymLog
aSymNormAll = aSymAll / aBack
copy_VarMeta( aSymAll, aSymNormAll )
bSymNormAll = bSymAll / bBack
copy_VarMeta( bSymAll, bSymNormAll )
diffSymNorm = bootstrap_diff( aSymAll, bSymAll, \
nTests, pThresh, calcDims )
printMinMax( diffSymNorm, True )
aAntiNormAll = aAntiAll / aBack
copy_VarMeta( aAntiAll, aAntiNormAll )
bAntiNormAll = bSymAll / bBack
copy_VarMeta( bAntiAll, bAntiNormAll )
diffAntiNorm = bootstrap_diff( aAntiAll, bAntiAll, \
nTests, pThresh, calcDims )
printMinMax( diffAntiNorm, True )
diffSymLog = aSymLog - bSymLog
copy_VarMeta( diffSym, diffSymLog )
diffAntiLog = hirsAntiLog - avhrrAntiLog
diffAntiLog = aAntiLog - bAntiLog
copy_VarMeta( diffAnti, diffAntiLog )
diffSymNorm = hirsSymNorm - avhrrSymNorm
copy_VarMeta( diffSym, diffSymNorm )
diffAntiNorm = hirsAntiNorm - avhrrAntiNorm
copy_VarMeta( diffAnti, diffAntiNorm )
print_clock( "Opening the workstation" )
; ...allows png or gif to work
......@@ -109,30 +129,30 @@ begin ; main
printVarSummary(diffAntiLog)
plotAntiLog = kf_draw( wks, diffAntiLog(freq|:,k|:), \
"None", False, True, res )
cjs_add_contours( wks, plotAntiLog, hirsAntiLog(freq|:,k|:), cntRes )
cjs_add_contours( wks, plotAntiLog, aAntiLog(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotAntiLog, diffAntiLog(freq|:,k|:), pvalRes )
res@gsnLeftString = "Log Symmetric"
print_clock( res@gsnLeftString )
plotSymLog = kf_draw( wks, diffSymLog(freq|:,k|:), \
"None", True, False, res )
cjs_add_contours( wks, plotSymLog, hirsSymLog(freq|:,k|:), cntRes )
cjs_add_contours( wks, plotSymLog, aSymLog(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotSymLog, diffSymLog(freq|:,k|:), pvalRes )
res@cnLevels = ispan( -5, 5, 1 ) * 0.05
res@cnLevels = ispan( -5, 5, 1 ) * 0.02
cntRes@cnLevels = cspan( 1.0, 2.0, 0.2 )
res@gsnLeftString = "Normalized Anti-symmetric"
print_clock( res@gsnLeftString )
plotAntiNorm = kf_draw( wks, diffAntiNorm(freq|:,k|:), \
"None", False, True, res )
cjs_add_contours( wks, plotAntiNorm, hirsAntiNorm(freq|:,k|:), cntRes )
cjs_add_contours( wks, plotAntiNorm, aAntiNorm(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotAntiNorm, diffAntiNorm(freq|:,k|:), pvalRes )
res@gsnLeftString = "Normalized Symmetric"
print_clock( res@gsnLeftString )
plotSymNorm = kf_draw( wks, diffSymNorm(freq|:,k|:), \
"None", True, False, res )
cjs_add_contours( wks, plotSymNorm, hirsSymNorm(freq|:,k|:), cntRes )
cjs_add_contours( wks, plotSymNorm, aSymNorm(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotSymNorm, diffSymNorm(freq|:,k|:), pvalRes )
......
Supports Markdown
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