Commit b4054626 authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Sat Sep 27 00:01:23 EDT 2014

parent a7f57e13
......@@ -37,28 +37,73 @@ begin ; main
fileName = "olr.anom.nc"
end if
sensor = (/ "hirs", "avhrr", "20th", "20th" /)
minYear = (/ 1981, 1981, 1981, 1905 /)
maxYear = (/ 2010, 2010, 2010, 1934 /)
dataList = NewList( "fifo" )
sensor = "hirs"
minYear = 1981
maxYear = 2010
inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
long_name = sensor + " (" + minYear + " - " + maxYear + ")"
print_clock( "Calculating " + long_name )
hirsData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
hirsData@long_name = long_name
ListPush( dataList, hirsData )
sensor = "avhrr"
minYear = 1981
maxYear = 2010
inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
long_name = sensor + " (" + minYear + " - " + maxYear + ")"
print_clock( "Calculating " + long_name )
avhrrData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
avhrrData@long_name = long_name
ListPush( dataList, avhrrData )
sensor = "20th"
minYear = 1981
maxYear = 2010
inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
nData = dimsizes(sensor)
long_name = sensor + " (" + minYear + " - " + maxYear + ")"
print_clock( "Calculating " + long_name )
modernData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
modernData@long_name = long_name
ListPush( dataList, modernData )
sensor = "20th"
minYear = 1911
maxYear = 1940
inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
long_name = sensor + " (" + minYear + " - " + maxYear + ")"
print_clock( "Calculating " + long_name )
ancientData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
ancientData@long_name = long_name
ListPush( dataList, ancientData )
dataList = NewList( "fifo" )
do plotCounter = 0, nData-1
long_name = sensor(plotCounter) + " (" \
+ minYear(plotCounter) + " - " + maxYear(plotCounter) + ")"
print_clock( "Calculating " + long_name )
data := calc_meanvar( inPath(plotCounter), varName, \
minYear(plotCounter), maxYear(plotCounter), meanNotVar )
data@long_name = long_name
ListPush( dataList, data )
end do
nData = ListCount( dataList )
; Customize base plot
res = True
totalRes = True
totalRes@cnFillPalette = "CBR_wet"
totalRes@cnLevelSelectionMode = "ExplicitLevels"
if( meanNotVar ) then
totalRes@cnLevels = nice_mnmxintvl( 200, 300, 10, False )
else
totalRes@cnLevels = nice_mnmxintvl( 100, 1800, 10, False )
end if
diffRes = True
diffRes@cnFillPalette = "CBR_drywet"
diffRes@cnLevelSelectionMode = "ExplicitLevels"
if( meanNotVar ) then
diffRes@cnLevels = nice_mnmxintvl( -200, 200, 10, False )
else
diffRes@cnLevels = nice_mnmxintvl( -1200, 1200, 10, False )
end if
; Customize panel
panRes = True
panRes@gsnPanelRowSpec = True
panRes@gsnPanelCenter = False
print_clock( "Drawing the plot" )
......@@ -75,15 +120,17 @@ begin ; main
nPlots = nData * ( nData + 1 ) / 2
plots = new( nPlots, graphic )
meanNotDiff = new( nPlots, logical )
plotCounter = 0
do xCounter = 0, nData-1
xData := dataList[xCounter]
res@gsnLeftString = inttochar(97+plotCounter) + ") " \
totalRes@gsnLeftString = inttochar(97+plotCounter) + ") " \
+ xData@long_name
print_clock( res@gsnLeftString )
print_clock( totalRes@gsnLeftString )
printMinMax( xData, True )
plots(plotCounter) = cjs_draw_shaded_map( wks, xData, \
fontHeightF, res )
fontHeightF, totalRes )
plotCounter = plotCounter + 1
do yCounter = 0, xCounter-1
......@@ -91,14 +138,15 @@ begin ; main
yRegrid := copy_Gridding( xData, yData )
xRegrid := copy_Gridding( yData, xData )
diffData = yData - xRegrid
diffData := yData - xRegrid
copy_VarMeta( yData, diffData )
res@gsnLeftString = inttochar(97+plotCounter) + ") " \
diffRes@gsnLeftString = inttochar(97+plotCounter) + ") " \
+ yData@long_name + " minus " + xData@long_name
print_clock( res@gsnLeftString )
print_clock( diffRes@gsnLeftString )
printMinMax( diffData, True )
plots(plotCounter) = cjs_draw_shaded_map( wks, diffData, \
fontHeightF, res )
fontHeightF, diffRes )
plotCounter = plotCounter + 1
end do
end do
......
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