Automated Nightly Commit - Sat Dec 13 00:00:55 EST 2014

compare/qsub_draw_series.sh

@@ -8,6 +8,7 @@ NCL_DIR=`pwd`
 mkdir log
 
 for PREFIX in olr mjo er kelvin low high
+#for PREFIX in td
 do
   for SERIES_TYPE in pct total anom
   do

compare/qsub_write_series.sh

@@ -7,7 +7,7 @@ echo $JOB_NAME `date`
 NCL_DIR=`pwd`
 mkdir log
 
-for VAR_NAME in high
+for VAR_NAME in td
 do
   for SENSOR in avhrr hirs
   do

compare/sub.draw_series.ncl

@@ -104,7 +104,7 @@ begin ; draw_series
     end do
 
     yData = 100. * yData / totalData
-    yData@units = "% of total variance"
+    yData@units = "% of high variance"
   end if
 
 

compare/sub.draw_wk99_trend.ncl

@@ -61,6 +61,15 @@ begin ; draw_wk99_trend
   copy_VarMeta( earlyPower, diffPower )
   printMinMax( diffPower, True )
 
+  ; smooth the difference
+  ; note that we ignore freq(0) because it's the powerful time-mean
+  do i = 1, 5
+    diffPower(freq|1:,k|:) = wgt_runave( diffPower(freq|1:,k|:), \\
+                                          (/ 0.25, 0.5, 0.25 /), 1 )
+    diffPower(k|:,freq|1:) = wgt_runave( diffPower(k|:,freq|1:), \\
+                                          (/ 0.25, 0.5, 0.25 /), 1 )
+  end do
+
   retVal = kf_draw( io_wks, diffPower(freq|:,k|:), "None", drawSymCurves, \
       drawAntiCurves, i_res )
 

compare/time_series.ncl

@@ -29,6 +29,7 @@ begin ; main
   if(.not.isvar("seriesType") )then
     seriesType = "anom"
   end if
+; sensor = (/ "hirs", "avhrr", "claus" /)
   sensor = (/ "hirs", "avhrr" /)
 
   plotType = "png"

compare/wk99_trend.ncl

@@ -24,7 +24,7 @@ begin ; main
 
   ; These are some parameters that could be useful to have up top
   if(.not.isvar( "sensor" ) ) then
-    sensor = "hirs"
+    sensor = "avhrr"
   end if
 
   plotType = "png"
@@ -49,7 +49,7 @@ begin ; main
   ; Customize base plot
   res = True
   res@cnFillPalette = "CBR_coldhot"
-  res@cnLevels = ispan( -5, 5, 1 ) * 0.1
+  res@cnLevels = ispan( -5, 5, 1 ) * 0.04
   res@cnLinesOn = False
   res@cnFillMode = "CellFill"