Commit 2cb3603c authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Wed Jul 13 00:00:20 EDT 2016

parent c9814985
......@@ -36,8 +36,8 @@ begin
minTime = cd_inv_calendar( 1979, 01, 01, 00, 0, 0, timeUnits, 0 )
maxTime = cd_inv_calendar( 2013, 12, 31, 18, 0, 0, timeUnits, 0 )
hirsBin = ispan( 210, 350, 5 )
avhrrBin = ispan( 210, 350, 5 )
hirsBin = ispan( 210, 310, 2 )
avhrrBin = ispan( 210, 310, 2 )
print_clock( "Reading the data" )
......@@ -59,6 +59,10 @@ begin
copy_VarAtts( hirsData, hirs1d )
copy_VarAtts( avhrrData, avhrr1d )
hirsGood = ( hirs1d.ge.min(hirsBin) ).and.( hirs1d.le.max(hirsBin) )
avhrrGood = ( avhrr1d.ge.min(avhrrBin) ).and.( avhrr1d.le.max(avhrrBin) )
goodInd = ind( hirsGood.and.avhrrGood )
; Customize base plot
res = True
......@@ -112,7 +116,7 @@ begin
; res@pmLabelBarOrthogonalPosF = fontHeightF * 5
; res@pmLabelBarHeightF = fontHeightF * 2
res@plotPercent = True
res@plotPercent = False
res@setZeroToMissing = True
print_clock( "Drawing the plot" )
......@@ -128,8 +132,9 @@ begin
wks = gsn_open_wks( plotTypeLocal, plotName )
gsn_merge_colormaps( wks, "posneg_carl", "default" )
plot = cjs_draw_hist2d( wks, hirs1d, avhrr1d, hirsBin, avhrrBin, res )
cjs_add_regline( wks, plot, hirs1d, avhrr1d, True )
plot = cjs_draw_hist2d( wks, hirs1d(goodInd), avhrr1d(goodInd), \
hirsBin, avhrrBin, res )
cjs_add_regline( wks, plot, hirs1d(goodInd), avhrr1d(goodInd), True )
draw( plot )
frame( wks )
......
......@@ -20,7 +20,7 @@ function draw_meanvar( \
begin ; draw_meanvar
; These are some parameters that could be useful to have up top
nTests = 1000
pThresh = 0.95
pThresh = 0.0
calcDims = 0
......@@ -102,8 +102,11 @@ begin ; draw_meanvar
copy_VarAtts( i_cntRes, cntRes )
end if
pvalRes = True
retVal = cjs_draw_shaded_map( io_wks, diffData, res )
cjs_add_contours( io_wks, retVal, plotData, cntRes )
cjs_add_signif( io_wks, retVal, diffData, pvalRes )
printMinMax( plotData, True )
printMinMax( diffData, True )
......
......@@ -23,8 +23,8 @@ begin ; main
fontHeightF = 0.02
nTests = 10
pThresh = 0.95
nTests = 1000
pThresh = 0.00
calcDims = (/ 0, 1 /)
print_clock( "Reading hirs" )
......@@ -101,6 +101,7 @@ begin ; main
cntRes@cnLineLabelFontHeightF = 0.01
pvalRes = True
; pvalRes@cnLinesOn = False
res@cnLevels = ispan( -5, 5, 1 ) * 0.02
res@gsnLeftString = "Log Anti-symmetric"
......@@ -152,7 +153,7 @@ begin ; main
+ "+repage -density " + plotDpi + " " \\
+ plotName + ".eps " + plotName + "." + plotType )
if( .not.isStrSubset( plotType, "e" ) ) then
system( "rm -f " + plotName + ".eps" )
; system( "rm -f " + plotName + ".eps" )
end if
end if
......
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