Automated Nightly Commit - Sat Dec 20 00:03:03 EST 2014

compare/map.ncl

@@ -18,7 +18,7 @@ load "$CJS_NCL_LIB/lib.cjs_graphics.ncl"
 load "$CJS_NCL_LIB/lib.array.ncl"
 load "$CJS_NCL_LIB/lib.percentiles.ncl"
 load "$CJS_NCL_LIB/print_clock.ncl"
-load "/home/carl/projects/olr_cdr/sub.calc_meanvar.ncl"
+load "/home/carl/projects/olr_cdr/compare/sub.calc_meanvar.ncl"
 
 
 begin ; main
@@ -29,6 +29,7 @@ begin ; main
   plotDpi  = 200
 
   totalNotAnom = False
+  showDiff = False
   if( .not.isvar("varName") ) then
     varName = "olr"
   end if
@@ -41,6 +42,8 @@ begin ; main
   end if
   minLat = -30
   maxLat =  30
+  minLon = 0
+  maxLon = 180
 
   timeUnits = "days since 1800-01-01 00:00:00"
   plotTime = cd_inv_calendar( 2011, 11, 27, 00, 0, 0, timeUnits, 0 )
@@ -78,6 +81,7 @@ begin ; main
   ; Customize base plot
   res = True
   res@lbOrientation = "Horizontal"
+  res@lbLabelBarOn = showDiff
   res@cnLevelSelectionMode = "ManualLevels"
   res@gsnMajorLonSpacing = 45
   res@gsnMinorLonSpacing = 15
@@ -85,6 +89,8 @@ begin ; main
   res@gsnRightString = cd_string( plotTime, "%d %c %Y" )
   res@mpMinLatF = minLat
   res@mpMaxLatF = maxLat
+  res@mpMinLonF = minLon
+  res@mpMaxLonF = maxLon
 
   totalRes = res
   if( totalNotAnom ) then
@@ -113,7 +119,12 @@ begin ; main
   panRes = True
   panRes@gsnPanelRowSpec = False
   panRes@gsnPanelCenter = True
+  panRes@gsnPanelBottom = 0.05
   panRes@gsnPanelYWhiteSpacePercent = 5
+  panRes@gsnPanelLabelBar   = .not.showDiff
+  panRes@lbTitleString      = hirsData@units
+  panRes@lbTitlePosition    = "Right"
+  panRes@lbTitleDirection   = "Across"
 
   print_clock( "Drawing the plot" )
 
@@ -166,7 +177,11 @@ begin ; main
     end do
   end do
 
-  gsn_panel( wks, plots, (/ nPlots, 1 /), panRes )
+  if( showDiff ) then
+    gsn_panel( wks, plots, (/ nPlots, 1 /), panRes )
+  else
+    gsn_panel( wks, plots(:nPlots-2), (/ nPlots-1, 1 /), panRes )
+  end if
 
   print_clock( "Convert the image, if necessary" )
   delete(wks)