Automated Nightly Commit - Sat Jul 9 00:00:20 EDT 2016

compare/wk99_diff.ncl

@@ -30,61 +30,44 @@ begin ; main
   print_clock( "Reading hirs" )
   inPath = "/home/carl/data/olr/compare/hirs/olr.wk99.nc"
   inFile = addfile( inPath, "r" )
-  hirsSymLog = inFile->symLogPower({-20:20},:)
-  hirsAntiLog = inFile->antiLogPower({-20:20},:)
-  hirsSymNorm = inFile->symNormPower({-20:20},:)
-  hirsAntiNorm = inFile->antiNormPower({-20:20},:)
-
-  hirsSymAll = inFile->symAllPower(:,:,{-20:20},:)
-  hirsAntiAll = inFile->antiAllPower(:,:,{-20:20},:)
-  hirsBack = conform( hirsSymAll, inFile->background({-20:20},:), (/ 2, 3 /) )
-  hirsSymAllLog = log10(hirsSymAll)
-  hirsAntiAllLog = log10(hirsAntiAll)
-  hirsSymAllNorm = hirsSymAll / hirsBack
-  hirsAntiAllNorm = hirsAntiAll / hirsBack
-  copy_VarMeta( hirsSymAll, hirsSymAllLog )
-  copy_VarMeta( hirsAntiAll, hirsAntiAllLog )
-  copy_VarMeta( hirsSymAll, hirsSymAllNorm )
-  copy_VarMeta( hirsAntiAll, hirsAntiAllNorm )
+  hirsSymLog = inFile->symLogPower(freq|:,{k|-20:20})
+  hirsAntiLog = inFile->antiLogPower(freq|:,{k|-20:20})
+  hirsSymNorm = inFile->symNormPower(freq|:,{k|-20:20})
+  hirsAntiNorm = inFile->antiNormPower(freq|:,{k|-20:20})
+
+  hirsSymAll = inFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
+  hirsAntiAll = inFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
+  hirsBack = conform( hirsSymAll, inFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
 
 
   print_clock( "Reading avhrr" )
   inPath = "/home/carl/data/olr/compare/avhrr/olr.wk99.nc"
   inFile = addfile( inPath, "r" )
-  avhrrSymLog = inFile->symLogPower({-20:20},:)
-  avhrrAntiLog = inFile->antiLogPower({-20:20},:)
-  avhrrSymNorm = inFile->symNormPower({-20:20},:)
-  avhrrAntiNorm = inFile->antiNormPower({-20:20},:)
-
-  avhrrSymAll = inFile->symAllPower(:,:,{-20:20},:)
-  avhrrAntiAll = inFile->antiAllPower(:,:,{-20:20},:)
-  avhrrBack = conform( avhrrSymAll, inFile->background({-20:20},:), (/ 2, 3 /) )
-  avhrrSymAllLog = log10(avhrrSymAll)
-  avhrrAntiAllLog = log10(avhrrAntiAll)
-  avhrrSymAllNorm = avhrrSymAll / avhrrBack
-  avhrrAntiAllNorm = avhrrAntiAll / avhrrBack
-  copy_VarMeta( avhrrSymAll, avhrrSymAllLog )
-  copy_VarMeta( avhrrAntiAll, avhrrAntiAllLog )
-  copy_VarMeta( avhrrSymAll, avhrrSymAllNorm )
-  copy_VarMeta( avhrrAntiAll, avhrrAntiAllNorm )
-
-  diffSymLog = bootstrap_diff( hirsSymAllLog, avhrrSymAllLog,  \
-      nTests, pThresh, calcDims )
-  printMinMax( diffSymLog, True )
+  avhrrSymLog = inFile->symLogPower(freq|:,{k|-20:20})
+  avhrrAntiLog = inFile->antiLogPower(freq|:,{k|-20:20})
+  avhrrSymNorm = inFile->symNormPower(freq|:,{k|-20:20})
+  avhrrAntiNorm = inFile->antiNormPower(freq|:,{k|-20:20})
 
-  diffAntiLog = bootstrap_diff( hirsAntiAllLog, avhrrAntiAllLog,  \
-      nTests, pThresh, calcDims )
-  printMinMax( diffAntiLog, True )
+  avhrrSymAll = inFile->symAllPower(lat|:,window|:,freq|:,{k|-20:20})
+  avhrrAntiAll = inFile->antiAllPower(lat|:,window|:,freq|:,{k|-20:20})
+  avhrrBack = conform( avhrrSymAll, inFile->background(freq|:,{k|-20:20}), (/ 2, 3 /) )
 
-
-  diffSymNorm = bootstrap_diff( hirsSymAllNorm, avhrrSymAllNorm,  \
+  diffSym = bootstrap_diff( hirsSymAll, avhrrSymAll,  \
       nTests, pThresh, calcDims )
-  printMinMax( diffSymNorm, True )
+  printMinMax( diffSym, True )
 
-  diffAntiNorm = bootstrap_diff( hirsAntiAllNorm, avhrrAntiAllNorm,  \
+  diffAnti = bootstrap_diff( hirsAntiAll, avhrrAntiAll,  \
       nTests, pThresh, calcDims )
-  printMinMax( diffAntiNorm, True )
+  printMinMax( diffAnti, True )
 
+  diffSymLog = hirsSymLog - avhrrSymLog
+  copy_VarMeta( diffSym, diffSymLog )
+  diffAntiLog = hirsAntiLog - avhrrAntiLog
+  copy_VarMeta( diffAnti, diffAntiLog )
+  diffSymNorm = hirsSymNorm - avhrrSymNorm
+  copy_VarMeta( diffSym, diffSymNorm )
+  diffAntiNorm = hirsAntiNorm - avhrrAntiNorm
+  copy_VarMeta( diffAnti, diffAntiNorm )
 
   print_clock( "Opening the workstation" )
   ; ...allows png or gif to work
@@ -117,6 +100,8 @@ begin ; main
   cntRes@cnMaxLevelCount = 8
   cntRes@cnLineLabelFontHeightF   = 0.01
 
+  pvalRes = True
+  
   res@cnLevels = ispan( -5, 5, 1 ) * 0.02
   res@gsnLeftString = "Log Anti-symmetric"
   print_clock( res@gsnLeftString )
@@ -124,12 +109,14 @@ begin ; main
   plotAntiLog = kf_draw( wks, diffAntiLog(freq|:,k|:),  \
       "None", False, True, res )
   cjs_add_contours( wks, plotAntiLog, hirsAntiLog(freq|:,k|:), cntRes )
+  cjs_add_signif( wks, plotAntiLog, diffAntiLog(freq|:,k|:), pvalRes )
 
   res@gsnLeftString = "Log Symmetric"
   print_clock( res@gsnLeftString )
   plotSymLog = kf_draw( wks, diffSymLog(freq|:,k|:),  \
       "None", True, False, res )
   cjs_add_contours( wks, plotSymLog, hirsSymLog(freq|:,k|:), cntRes )
+  cjs_add_signif( wks, plotSymLog, diffSymLog(freq|:,k|:), pvalRes )
 
   res@cnLevels = ispan( -5, 5, 1 ) * 0.05
   cntRes@cnLevels = cspan( 1.0, 2.0, 0.2 )
@@ -138,12 +125,14 @@ begin ; main
   plotAntiNorm = kf_draw( wks, diffAntiNorm(freq|:,k|:),  \
       "None", False, True, res )
   cjs_add_contours( wks, plotAntiNorm, hirsAntiNorm(freq|:,k|:), cntRes )
+  cjs_add_signif( wks, plotAntiNorm, diffAntiNorm(freq|:,k|:), pvalRes )
 
   res@gsnLeftString = "Normalized Symmetric"
   print_clock( res@gsnLeftString )
   plotSymNorm = kf_draw( wks, diffSymNorm(freq|:,k|:),  \
       "None", True, False, res )
   cjs_add_contours( wks, plotSymNorm, hirsSymNorm(freq|:,k|:), cntRes )
+  cjs_add_signif( wks, plotSymNorm, diffSymNorm(freq|:,k|:), pvalRes )