Commit 3a19b1c9 authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Thu Oct 9 00:01:29 EDT 2014

parent 508c4cd9
......@@ -59,8 +59,9 @@ begin ; main
res = True
res@cnFillPalette = "colorbrewer_wetdry"
shdPath = "/home/carl/data/olr/cdr/" + sensor + ".anom.nc"
cntPath = "/home/carl/data/olr/cdr/" + sensor + ".anom.waves.nc"
res@gsnLeftString = str_upper(sensor)
shdPath = "/home/carl/data/olr/compare/" + sensor + "/olr.anom.nc"
cntPath = "/home/carl/data/olr/compare/" + sensor + "/olr.waves.nc"
shdLevels = ispan(-5,5,1) * 10
cntLevels = ispan( -10,10,1) * 10
plotA = draw_hov( wks, shdPath, cntPath, shdLevels, cntLevels, \
......
......@@ -27,37 +27,28 @@ begin ; main
window = 96 * 2
overlap = 64 * 2
sensor = "hirs prelim"
sensor = "avhrr"
plotType = "png"
plotName = "figures/prelim_" + sensor
plotName = "figures/spectra_" + sensor
plotDpi = 200
fontHeightF = 0.02
print( "Reading the data..." )
; fileName = "/home/carl/data/olr/cdr/" + sensor + ".anom.nc"
fileName = "/home/carl/data/olr/cdr/daily/v01r02_test/olr.anom.nc"
fileName = "/home/carl/data/olr/compare/" + sensor + "/olr.anom.nc"
fin = addfile( fileName, "r" )
timeUnits = fin->time@units
print( (/ timeUnits /) )
minTime = cd_inv_calendar( 2013, 01, 01, 00, 00, 00, timeUnits, 0 )
data = fin->olr({tofloat(minTime):},{-15:15},:)
data = fin->olr(:,{-15:15},:)
obsPerDay = 1
; data@_FillValue = 0
printVarSummary( data )
printMinMax( data, True )
print( num( ismissing( data ) ) )
if( any( ismissing(data) ) ) then
print_clock( "WARNING: Setting missing data to zero" )
data = where( ismissing(data), 0, data )
end if
maxMiss = 0
do t = 0, dimsizes(data&time)-1
miss = num( ismissing( data(t,:,:) ) )
if( miss.ne.0 ) then
print( ut_string( data&time(t), "" ) + ": " + miss )
end if
end do
; Find north and south
north = data(time|:,{lat|0:},lon|:)
south = data(time|:,{lat|:0:-1},lon|:)
......
......@@ -158,7 +158,7 @@ begin
res@tmYRLabels = (/ 48, 20, 10, 6.67, 5, 4, 3.33, 2.86, 2.5, 2.22, 2, 1.82, 1.67, 1.54, 1.43, 1.33, 1.25 /)
res@tmYRValues = 1 / res@tmYRLabels
res@gsnLeftString = inttochar(97+i_plotNumber) + ") Anti-symmetric " + str_upper(i_sensor)
res@gsnLeftString = inttochar(97+i_plotNumber) + ") " + str_upper(i_sensor) + ": Anti-symmetric"
; res@tiMainString = "Anti-symmetric"
res@trYMinF = 0.0
......@@ -179,7 +179,7 @@ begin
res@tiYAxisSide = "Left"
plotA = gsn_csm_contour( io_wks, antiPower(freq|:,k|:), res )
res@gsnLeftString = inttochar(98+i_plotNumber) + ") Symmetric " + str_upper(i_sensor)
res@gsnLeftString = inttochar(98+i_plotNumber) + ") " + str_upper(i_sensor) + ": Symmetric"
; res@tiMainString = "Symmetric"
print( min( symPower({freq|:.51},{k|-20:20}) ) + " " + max( symPower({freq|:.51},{k|-20:20}) ) )
......@@ -331,10 +331,10 @@ begin
plotType = "png"
plotName = "spectra"
plotName = "figures/spectra"
plotRes = 300
fontHeightF = 0.02
fontHeightF = 0.03
if( ( plotType.eq."png" ).or.( plotType.eq."gif" ) ) then
plotTypeLocal = "eps"
......@@ -353,28 +353,25 @@ begin
panRes@lbTitleDirection = "Across"
panRes@lbOrientation = "Horizontal"
panRes@lbLabelAutoStride = True
panRes@lbTitleFontHeightF = fontHeightF * 0.8
panRes@lbLabelFontHeightF = fontHeightF * 0.8
panRes@lbTitleFontHeightF = fontHeightF * 0.5
panRes@lbLabelFontHeightF = fontHeightF * 0.5
panRes@pmLabelBarOrthogonalPosF = fontHeightF * 0
panRes@pmLabelBarHeightF = fontHeightF * 2
wks = gsn_open_wks( plotTypeLocal, plotName )
gsn_merge_colormaps( wks, "posneg_carl", "default" )
sensor = "avhrr"
avhrrPlots = plotOne( wks, sensor, 0, fontHeightF )
sensor = "hirs"
hirsPlots = plotOne( wks, sensor, 2, fontHeightF )
hirsPlots = plotOne( wks, sensor, 0, fontHeightF )
sensor = "20th"
centPlots = plotOne( wks, sensor, 2, fontHeightF )
sensor = "avhrr"
avhrrPlots = plotOne( wks, sensor, 2, fontHeightF )
gsn_panel( wks, ndtooned( (/ hirsPlots, avhrrPlots, centPlots /) ), (/ 3, 2 /), panRes )
gsn_panel( wks, ndtooned( (/ hirsPlots, avhrrPlots /) ), (/ 3, 2 /), panRes )
delete(wks)
if( ( plotType.eq."png" ).or.( plotType.eq."gif" ) ) then
system( "convert -trim +repage -density " + plotRes + " " \\
+ plotName + ".eps " + plotName + "." + plotType )
; system( "rm -f " + plotName + ".eps" )
system( "rm -f " + plotName + ".eps" )
end if
end
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