Commit 5d5d5dea authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Tue Oct 7 00:01:39 EDT 2014

parent 5bc0a9d4
......@@ -34,13 +34,35 @@ begin ; main
if( meanNotVar ) then
fileName = "olr.total.nc"
plotName = "figures/map.mean"
units = "W m-2"
factor = 1
else
if( varName.eq."olr" ) then
fileName = "olr.anom.nc"
plotName = "figures/map.variance"
units = "(W m-2)^2"
factor = 1
else
fileName = "olr.std.waves.nc"
plotName = "figures/map." + varName
units = "%"
factor = 100
end if
end if
if( varName.eq."olr" ) then
minLat = -60
maxLat = 60
pmLabelBarOrthogonalPosF = 0.3
else
if( varName.eq."kelvin" ) then
minLat = -20
maxLat = 20
pmLabelBarOrthogonalPosF = 0.6
else
minLat = -30
maxLat = 30
pmLabelBarOrthogonalPosF = 0.5
end if
end if
......@@ -58,8 +80,9 @@ begin ; main
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},:) ) /) )
ancientData = factor * ancientData
minTotal = min( (/ minTotal, min( ancientData({minLat:maxLat},:) ) /) )
maxTotal = max( (/ maxTotal, max( ancientData({minLat:maxLat},:) ) /) )
ListPush( dataList, ancientData )
sensor = "20th"
......@@ -70,8 +93,9 @@ begin ; main
print_clock( "Calculating " + long_name )
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},:) ) /) )
modernData = factor * modernData
minTotal = min( (/ minTotal, min( modernData({minLat:maxLat},:) ) /) )
maxTotal = max( (/ maxTotal, max( modernData({minLat:maxLat},:) ) /) )
ListPush( dataList, modernData )
end if
......@@ -83,8 +107,9 @@ 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},:) ) /) )
avhrrData = factor * avhrrData
minTotal = min( (/ minTotal, min( avhrrData({minLat:maxLat},:) ) /) )
maxTotal = max( (/ maxTotal, max( avhrrData({minLat:maxLat},:) ) /) )
ListPush( dataList, avhrrData )
sensor = "hirs"
......@@ -95,18 +120,27 @@ begin ; main
print_clock( "Calculating " + long_name )
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},:) ) /) )
hirsData = factor * hirsData
minTotal = min( (/ minTotal, min( hirsData({minLat:maxLat},:) ) /) )
maxTotal = max( (/ maxTotal, max( hirsData({minLat:maxLat},:) ) /) )
ListPush( dataList, hirsData )
nData = ListCount( dataList )
print( dataList )
; Customize base plot
totalRes = True
totalRes@lbOrientation = "Vertical"
res = True
res@lbOrientation = "Horizontal"
res@cnLevelSelectionMode = "ManualLevels"
res@lbTitleString = units
res@mpMinLatF = minLat
res@mpMaxLatF = maxLat
res@gsnMajorLonSpacing = 45
res@gsnMinorLonSpacing = 15
res@pmLabelBarOrthogonalPosF = pmLabelBarOrthogonalPosF
totalRes = res
totalRes@cnFillPalette = "CBR_wet"
totalRes@cnLevelSelectionMode = "ManualLevels"
if( meanNotVar ) then
totalRes@gsnSpreadColorStart = 10
totalRes@gsnSpreadColorEnd = 0
......@@ -115,15 +149,9 @@ begin ; main
totalRes@cnMinLevelValF = mnmxintvl(0)
totalRes@cnMaxLevelValF = mnmxintvl(1)
totalRes@cnLevelSpacingF = mnmxintvl(2)
totalRes@pmLabelBarOrthogonalPosF = 0.30
totalRes@tmXBLabelsOn = False
diffRes = True
diffRes@lbOrientation = "Vertical"
diffRes = res
diffRes@cnFillPalette = "CBR_drywet"
diffRes@cnLevelSelectionMode = "ManualLevels"
diffRes@pmLabelBarOrthogonalPosF = 0.30
diffRes@tmXBLabelsOn = True
; Customize panel
panRes = True
......@@ -169,7 +197,7 @@ begin ; main
+ yData@long_name + " minus " + xData@long_name
print_clock( diffRes@gsnLeftString )
currMax = 0.9 * max( abs(diffData({-60:60},:)) )
currMax = 0.9 * max( abs(diffData({minLat:maxLat},:)) )
mnmxintvl := nice_mnmxintvl( -currMax, currMax, 10, False )
diffRes@cnMinLevelValF = mnmxintvl(0)
diffRes@cnMaxLevelValF = mnmxintvl(1)
......
......@@ -9,16 +9,16 @@ mkdir log
for VAR_NAME in olr mjo er kelvin low
do
NCL_OPTION="varName=\"$VAR_NAME\""
JOB_NAME="map.$VAR_NAME"
qsub \
-N $JOB_NAME \
-o $NCL_DIR/log/$JOB_NAME.log -j oe \
-q san \
-l walltime=12:00:00,nodes=1:ppn=1 \
-v NCL_DIR=$NCL_DIR,NCL_SCRIPT=$NCL_SCRIPT,NCL_OPTION="$NCL_OPTION" \
$NCL_DIR/run_ncl.pbs
NCL_OPTION="varName=\"$VAR_NAME\""
JOB_NAME="map.$VAR_NAME"
qsub \
-N $JOB_NAME \
-o $NCL_DIR/log/$JOB_NAME.log -j oe \
-q san \
-l walltime=12:00:00,nodes=1:ppn=1 \
-v NCL_DIR=$NCL_DIR,NCL_SCRIPT=$NCL_SCRIPT,NCL_OPTION="$NCL_OPTION" \
$NCL_DIR/run_ncl.pbs
done
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment