Automated Nightly Commit - Sat Dec 6 00:02:25 EST 2014

qsub_write_series.sh

@@ -9,7 +9,7 @@ mkdir log
 
 for VAR_NAME in olr mjo er kelvin low
 do
-  for SENSOR in avhrr hirs 20th
+  for SENSOR in claus
   do
 
     NCL_OPTION="varName=\"$VAR_NAME\" sensor=\"$SENSOR\""

sub.draw_series.ncl

@@ -35,12 +35,13 @@ begin ; draw_series
     buffer = 0
   else
     buffer = 365
-    buffer = 365
+    buffer = 1095
   end if
   shortSmooth = 365
   longSmooth = shortSmooth*5
   inPath = "/home/carl/data/olr/compare/time_series/" + i_prefix + "."  \
       + i_sensor + ".20S-20N.nc"
+
   res = i_res
   res@trXMinF = i_minTime
   res@trXMaxF = i_maxTime
@@ -64,6 +65,8 @@ begin ; draw_series
   end do
   xData = remove_duplicates( allTimes )
   xData@units = inTime@units
+  xData!0 = "time"
+  xData&time = xData
 
   print_clock("Reading data")
   yData = new( (/ dimsizes(i_sensor), dimsizes(xData) /), float )
@@ -79,6 +82,14 @@ begin ; draw_series
     yData(sensorCounter,{inTime(startInd):inTime(endInd)}) = inData
   end do
 
+  print_clock("Trending")
+  trendData = new( dimsizes(yData), double )
+  do sensorCounter = 0, dimsizes(i_sensor)-1
+    rc = regline( xData({i_minTime:i_maxTime}),  \
+        yData(sensorCounter,{i_minTime:i_maxTime}) )
+    trendData(sensorCounter,:) = rc * xData + rc@yintercept
+  end do
+
   print_clock("Smoothing")
   yData = runave_n_Wrap( yData, shortSmooth, 0, 1 )
   smoothData = runave_n_Wrap( yData, longSmooth, 0, 1 )
@@ -93,6 +104,11 @@ begin ; draw_series
   res@xyBoldColors = False
   rawPlot = cjs_draw_xy( io_wks, xData, yData, res )
 
+  res@xyLineThicknessF = 1
+  res@xyBoldColors = False
+  trendPlot = cjs_draw_xy( io_wks, xData, trendData, res )
+  overlay( rawPlot, trendPlot )
+
   res@xyLineThicknessF = 6
   res@xyBoldColors = True
   smoothPlot = cjs_draw_xy( io_wks, xData, smoothData, res )

time_series.ncl

@@ -29,7 +29,7 @@ begin ; main
   if(.not.isvar("varName") )then
     varName = "anom"
   end if
-  sensor = (/ "hirs", "avhrr", "20th" /)
+  sensor = (/ "hirs", "avhrr", "20th", "claus" /)
 
   plotType = "png"
   plotName = "figures/" + prefix + "." + varName

write_time_series.ncl

@@ -35,6 +35,11 @@ begin ; main
   else
     fileName = "olr.waves.nc"
   end if
+  if( sensor.eq."claus" ) then
+    rawVar = "brtmp"
+  else
+    rawVar = "olr"
+  end if
 
   minLon =  00
   maxLon = 360
@@ -53,7 +58,7 @@ begin ; main
 
   print_clock( "Reading" )
   rawFile = addfile( rawPath, "r" )
-  rawData = lon_subset( rawFile->olr(:,{minLat:maxLat},:),  \
+  rawData = lon_subset( rawFile->$rawVar$(:,{minLat:maxLat},:),  \
       minLon, maxLon )
 
   inFile = addfile( inPath, "r" )
@@ -83,9 +88,15 @@ begin ; main
   print_clock( "Writing" )
   system("rm " + outPath )
   outFile = addfile( outPath, "c" )
-  outFile->time = time
-  outFile->total = totalSeries
-  outFile->anom = anomSeries
+  if( sensor.eq."claus" ) then
+    outFile->time = time(::4)
+    outFile->total = totalSeries(::4)
+    outFile->anom = anomSeries(::4)
+  else
+    outFile->time = time
+    outFile->total = totalSeries
+    outFile->anom = anomSeries
+  end if
 
 
   print_clock( "Thank you, come again." )