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

Automated Nightly Commit - Sat Oct 18 00:01:07 EDT 2014

parent 3e09a29e
......@@ -28,7 +28,7 @@ begin ; main
plotType = "png"
plotDpi = 200
totalNotAnom = True
totalNotAnom = False
if( .not.isvar("varName") ) then
varName = "olr"
end if
......@@ -39,12 +39,12 @@ begin ; main
fileName = "olr.anom.nc"
plotName = "figures/map.anom"
end if
minLat = -60
maxLat = 60
minLat = -30
maxLat = 30
timeUnits = "days since 1800-01-01 00:00:00"
plotTime = cd_inv_calendar( 2011, 11, 27, 00, 0, 0, timeUnits, 0 )
plotTime = cd_inv_calendar( 1987, 05, 15, 00, 0, 0, timeUnits, 0 )
; plotTime = cd_inv_calendar( 1987, 05, 15, 00, 0, 0, timeUnits, 0 )
plotName = plotName + cd_string( plotTime, ".%Y-%N-%D" )
dataList = NewList( "lifo" )
......@@ -81,8 +81,10 @@ begin ; main
res@cnLevelSelectionMode = "ManualLevels"
res@gsnMajorLonSpacing = 45
res@gsnMinorLonSpacing = 15
res@pmLabelBarOrthogonalPosF = 0.3
res@pmLabelBarOrthogonalPosF = 0.4
res@gsnRightString = cd_string( plotTime, "%d %c %Y" )
res@mpMinLatF = minLat
res@mpMaxLatF = maxLat
totalRes = res
if( totalNotAnom ) then
......
......@@ -9,7 +9,7 @@ mkdir log
for PREFIX in olr mjo er kelvin low
do
for VAR_NAME in total anom
for VAR_NAME in anom
do
NCL_OPTION="varName=\"$VAR_NAME\" prefix=\"$PREFIX\""
......
......@@ -31,8 +31,13 @@ function draw_series( \
)
begin ; draw_series
; These are some parameters that could be useful to have up top
if( i_prefix.eq."olr" ) then
buffer = 0
else
buffer = 365
end if
shortSmooth = 365
longSmooth = 365*5
longSmooth = shortSmooth*5
inPath = "/home/carl/data/olr/compare/time_series/" + i_prefix + "." \
+ i_sensor + ".20S-20N.nc"
res = i_res
......@@ -42,6 +47,7 @@ begin ; draw_series
res@tiXAxisString = ""
res@xyMarkLineMode = "Lines"
res@xyMonoLineThickness = True
res@vpHeightF = 0.4
print_clock( "Reading times" )
do sensorCounter = 0, dimsizes(i_sensor)-1
......@@ -66,8 +72,10 @@ begin ; draw_series
do sensorCounter = 0, dimsizes(i_sensor)-1
inFile = addfile( inPath(sensorCounter), "r" )
inTime := floor(inFile->time)
inData := inFile->$i_varName$
yData(sensorCounter,{min(inTime):max(inTime)}) = inData
startInd = buffer
endInd = dimsizes(inTime)-1-buffer
inData := inFile->$i_varName$(startInd:endInd)
yData(sensorCounter,{inTime(startInd):inTime(endInd)}) = inData
end do
print_clock("Smoothing")
......@@ -80,10 +88,12 @@ begin ; draw_series
resTick@ttmAxis = "XB"
month_axis_labels( i_minTime, 60, 12, 1, res, resTick )
res@xyLineThicknessF = 1
res@xyLineThicknessF = 2
res@xyBoldColors = False
rawPlot = cjs_draw_xy( io_wks, xData, yData, res )
res@xyLineThicknessF = 3
res@xyLineThicknessF = 6
res@xyBoldColors = True
smoothPlot = cjs_draw_xy( io_wks, xData, smoothData, res )
overlay( rawPlot, smoothPlot )
......
......@@ -37,7 +37,7 @@ begin ; main
timeUnits = "days since 1800-01-01 00:00:00"
minTime = cd_inv_calendar( 1974, 01, 01, 00, 0, 0, timeUnits, 0 )
maxTime = cd_inv_calendar( 2013, 12, 31, 18, 0, 0, timeUnits, 0 )
maxTime = cd_inv_calendar( 2014, 01, 01, 00, 0, 0, timeUnits, 0 )
; Customize base plot
res = True
......
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