Automated Nightly Commit - Thu Mar 29 00:00:15 EDT 2018

compare/bsub_ncl.sh

@@ -6,7 +6,7 @@ else
 fi
 echo $QUEUE
 
-NCL_SCRIPT=wk99_diff
+NCL_SCRIPT=histo2d
 NCL_DIR=`pwd`
 LOG_DIR=$NCL_DIR/log
 mkdir -p $LOG_DIR
@@ -23,7 +23,7 @@ bsub \
   -o $LOG_FILE \
   -e $ERR_FILE \
   -q $QUEUE \
-  -n 1 -W 3:00 \
+  -n 1 -W 12:00 \
   -sp 75 \
   ~/template/run_ncl.sh $NCL_DIR $NCL_SCRIPT "$NCL_OPTION"
 

compare/histo2d.ncl

@@ -24,10 +24,12 @@ begin
   fontHeightF = 0.02
 
 ; shdLevels = ispan(-7,7,1)
-  shdLevels = -999
+; shdLevels = -999
+; shdLevels = 10^ispan(1,5,1)
+  shdLevels = 500 * 2^ispan(1,8,1)
 
-  minLon = 000
-  maxLon = 360
+  minLon = -180
+  maxLon = 180
   minLat = -30
   maxLat =  30
   centerLon = 180
@@ -36,20 +38,22 @@ begin
   minTime = cd_inv_calendar( 1979, 01, 01, 00, 0, 0, timeUnits, 0 )
   maxTime = cd_inv_calendar( 2013, 12, 31, 18, 0, 0, timeUnits, 0 )
 
-  hirsBin  = ispan( 210, 310, 2 )
-  avhrrBin  = ispan( 210, 310, 2 )
+  hirsBin  = ispan( 100, 350, 2 )
+  avhrrBin  = ispan( 100, 350, 2 )
 
   print_clock( "Reading the data" )
 
   print_clock( "Reading AVHRR" )
   inPath = "/home/carl/data/olr/compare/avhrr/olr.total.nc"
   inFile = addfile( inPath, "r" )
-  avhrrData = inFile->olr({minTime:maxTime},{minLat:maxLat},:)
+  avhrrData = lon_subset( inFile->olr({minTime:maxTime},{minLat:maxLat},:),  \
+      minLon, maxLon )
 
   print_clock( "Reading HIRS" )
   inPath = "/home/carl/data/olr/compare/hirs/olr.total.nc"
   inFile = addfile( inPath, "r" )
-  hirsData = inFile->olr({minTime:maxTime},{minLat:maxLat},:)
+  hirsData = lon_subset( inFile->olr({minTime:maxTime},{minLat:maxLat},:),  \
+      minLon, maxLon )
 
   print_clock( "Regridding" )
   hirsRegrid = copy_Gridding( avhrrData, hirsData )
@@ -137,6 +141,9 @@ begin
 
   plot = cjs_draw_hist2d( wks, hirs1d(goodInd), avhrr1d(goodInd),  \
       hirsBin, avhrrBin, res )
+
+  regRes = True
+  regRes@maxNumData = 1000
   cjs_add_regline( wks, plot, hirs1d(goodInd), avhrr1d(goodInd), True )
 
   draw( plot )