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

Automated Nightly Commit - Tue Dec 9 00:02:32 EST 2014

parent a8e7f087
......@@ -24,7 +24,7 @@ function draw_series( \
io_wks [1] : graphic, \
i_prefix [1] : string, \
i_sensor [*] : string, \
i_varName [1] : string, \
i_seriesType [1] : string, \
i_minTime [1] : numeric, \
i_maxTime [1] : numeric, \
i_res [1] : logical \
......@@ -35,7 +35,6 @@ begin ; draw_series
buffer = 0
else
buffer = 365
buffer = 1095
end if
shortSmooth = 365
longSmooth = shortSmooth*5
......@@ -43,6 +42,7 @@ begin ; draw_series
+ i_sensor + ".20S-20N.nc"
res = i_res
res@gsnLeftString = i_prefix
res@trXMinF = i_minTime
res@trXMaxF = i_maxTime
res@xyExplicitLegendLabels = i_sensor
......@@ -78,7 +78,7 @@ begin ; draw_series
inTime := floor(inFile->time)
startInd = buffer
endInd = dimsizes(inTime)-1-buffer
inData := inFile->$i_varName$(startInd:endInd)
inData := inFile->$i_seriesType$(startInd:endInd)
yData(sensorCounter,{inTime(startInd):inTime(endInd)}) = inData
end do
......
......@@ -24,25 +24,25 @@ begin ; main
; These are some parameters that could be useful to have up top
if(.not.isvar("prefix") )then
prefix = "kelvin"
prefix = "olr"
end if
if(.not.isvar("varName") )then
varName = "anom"
if(.not.isvar("seriesType") )then
seriesType = "anom"
end if
sensor = (/ "hirs", "avhrr", "20th", "claus" /)
sensor = (/ "hirs", "avhrr", "claus" /)
plotType = "png"
plotName = "figures/" + prefix + "." + varName
plotName = "figures/" + prefix + "." + seriesType
plotDpi = 200
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( 2014, 01, 01, 00, 0, 0, timeUnits, 0 )
minTime = cd_inv_calendar( 1981, 01, 01, 00, 0, 0, timeUnits, 0 )
maxTime = cd_inv_calendar( 2011, 01, 01, 00, 0, 0, timeUnits, 0 )
; Customize base plot
res = True
res@trYMinF = -1
res@trYMaxF = 1
; res@trYMinF = -1
; res@trYMaxF = 1
print_clock( "Drawing the plot" )
......@@ -56,7 +56,7 @@ begin ; main
; ...open the workstation
wks = gsn_open_wks( plotTypeLocal, plotName )
plot = draw_series( wks, prefix, sensor, varName, minTime, maxTime, res )
plot = draw_series( wks, prefix, sensor, seriesType, minTime, maxTime, res )
draw( plot )
frame( wks )
......
......@@ -25,7 +25,7 @@ begin ; main
; These are some parameters that could be useful to have up top
if(.not.isvar( "sensor" ) ) then
sensor = "avhrr"
sensor = "claus"
end if
if(.not.isvar("varName") ) then
varName = "olr"
......@@ -62,7 +62,11 @@ begin ; main
minLon, maxLon )
inFile = addfile( inPath, "r" )
inData = lon_subset( inFile->$varName$(:,{minLat:maxLat},:), minLon, maxLon )
if( varName.eq."olr" ).and.( sensor.eq."claus" ) then
inData = lon_subset( inFile->anom(:,{minLat:maxLat},:), minLon, maxLon )
else
inData = lon_subset( inFile->$varName$(:,{minLat:maxLat},:), minLon, maxLon )
end if
time = inFile->time
print_clock( "Calculating series" )
......
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