Automated Nightly Commit - Tue Sep 30 00:01:59 EDT 2014

meanvar.ncl

@@ -25,29 +25,49 @@ begin ; main
 
   ; These are some parameters that could be useful to have up top
   plotType = "png"
-  plotName = "plot"
   plotDpi  = 200
 
   fontHeightF = 0.02
-  meanNotVar = False
+  meanNotVar = True
   varName = "olr"
   if( meanNotVar ) then
     fileName = "olr.total.nc"
+    plotName = "mean"
   else
     fileName = "olr.anom.nc"
+;   fileName = "olr.std.waves.nc"
+    plotName = "variance"
   end if
+; plotName = varName
 
-  dataList = NewList( "fifo" )
+  dataList = NewList( "lifo" )
+  minTotal = -999.
+  minTotal@_FillValue = minTotal
+  maxTotal = minTotal
 
-  sensor = "hirs"
+  sensor = "20th"
+  minYear = 1911
+  maxYear = 1940
+  inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
+  long_name = sensor + " (" + minYear + " - " + maxYear + ")"
+  print_clock( "Calculating " + long_name )
+  ancientData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
+  ancientData@long_name = long_name
+  minTotal = min( (/ minTotal, min( ancientData({-60:60},:) ) /) )
+  maxTotal = max( (/ maxTotal, max( ancientData({-60:60},:) ) /) )
+  ListPush( dataList, ancientData )
+
+  sensor = "20th"
   minYear = 1981
   maxYear = 2010
   inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
   long_name = sensor + " (" + minYear + " - " + maxYear + ")"
   print_clock( "Calculating " + long_name )
-  hirsData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
-  hirsData@long_name = long_name
-  ListPush( dataList, hirsData )
+  modernData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
+  modernData@long_name = long_name
+  minTotal = min( (/ minTotal, min( modernData({-60:60},:) ) /) )
+  maxTotal = max( (/ maxTotal, max( modernData({-60:60},:) ) /) )
+  ListPush( dataList, modernData )
 
   sensor = "avhrr"
   minYear = 1981
@@ -57,48 +77,43 @@ begin ; main
   print_clock( "Calculating " + long_name )
   avhrrData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
   avhrrData@long_name = long_name
+  minTotal = min( (/ minTotal, min( avhrrData({-60:60},:) ) /) )
+  maxTotal = max( (/ maxTotal, max( avhrrData({-60:60},:) ) /) )
   ListPush( dataList, avhrrData )
 
-  sensor = "20th"
+  sensor = "hirs"
   minYear = 1981
   maxYear = 2010
   inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
   long_name = sensor + " (" + minYear + " - " + maxYear + ")"
   print_clock( "Calculating " + long_name )
-  modernData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
-  modernData@long_name = long_name
-  ListPush( dataList, modernData )
-
-  sensor = "20th"
-  minYear = 1911
-  maxYear = 1940
-  inPath = "/home/carl/data/olr/compare/" + sensor + "/" + fileName
-  long_name = sensor + " (" + minYear + " - " + maxYear + ")"
-  print_clock( "Calculating " + long_name )
-  ancientData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
-  ancientData@long_name = long_name
-  ListPush( dataList, ancientData )
+  hirsData := calc_meanvar( inPath, varName, minYear, maxYear, meanNotVar )
+  hirsData@long_name = long_name
+  minTotal = min( (/ minTotal, min( hirsData({-60:60},:) ) /) )
+  maxTotal = max( (/ maxTotal, max( hirsData({-60:60},:) ) /) )
+  ListPush( dataList, hirsData )
 
   nData = ListCount( dataList )
+  print( dataList )
 
   ; Customize base plot
   totalRes = True
   totalRes@cnFillPalette = "CBR_wet"
-  totalRes@cnLevelSelectionMode = "ExplicitLevels"
+  totalRes@cnLevelSelectionMode = "ManualLevels"
   if( meanNotVar ) then
-    totalRes@cnLevels = nice_mnmxintvl( 200, 300, 10, False )
-  else
-    totalRes@cnLevels = nice_mnmxintvl( 100, 1800, 10, False )
+    totalRes@gsnSpreadColorStart = 10
+    totalRes@gsnSpreadColorEnd = 0
   end if
+  mnmxintvl = nice_mnmxintvl( minTotal, maxTotal, 10, False )
+  totalRes@cnMinLevelValF = mnmxintvl(0)
+  totalRes@cnMaxLevelValF = mnmxintvl(1)
+  totalRes@cnLevelSpacingF = mnmxintvl(2)
+  totalRes@pmLabelBarOrthogonalPosF = 0.30
 
   diffRes = True
   diffRes@cnFillPalette = "CBR_drywet"
-  diffRes@cnLevelSelectionMode = "ExplicitLevels"
-  if( meanNotVar ) then
-    diffRes@cnLevels = nice_mnmxintvl( -200, 200, 10, False )
-  else
-    diffRes@cnLevels = nice_mnmxintvl( -1200, 1200, 10, False )
-  end if
+  diffRes@cnLevelSelectionMode = "ManualLevels"
+  diffRes@pmLabelBarOrthogonalPosF = 0.30
 
   ; Customize panel
   panRes = True
@@ -128,7 +143,6 @@ begin ; main
     totalRes@gsnLeftString = inttochar(97+plotCounter) + ") "  \
         + xData@long_name
     print_clock( totalRes@gsnLeftString )
-    printMinMax( xData, True )
     plots(plotCounter) = cjs_draw_shaded_map( wks, xData,  \
         fontHeightF, totalRes )
     plotCounter = plotCounter + 1
@@ -144,7 +158,13 @@ begin ; main
       diffRes@gsnLeftString = inttochar(97+plotCounter) + ") "  \
           + yData@long_name + " minus " + xData@long_name
       print_clock( diffRes@gsnLeftString )
-      printMinMax( diffData, True )
+
+      currMax = 0.9 * max( abs(diffData({-60:60},:)) )
+      mnmxintvl := nice_mnmxintvl( -currMax, currMax, 10, False )
+      diffRes@cnMinLevelValF = mnmxintvl(0)
+      diffRes@cnMaxLevelValF = mnmxintvl(1)
+      diffRes@cnLevelSpacingF = mnmxintvl(2)
+
       plots(plotCounter) = cjs_draw_shaded_map( wks, diffData,  \
           fontHeightF, diffRes )
       plotCounter = plotCounter + 1