Commit d94705b1 authored by Carl Schreck's avatar Carl Schreck
Browse files

Automated Nightly Commit - Thu Jun 21 00:00:17 EDT 2018

parent 5531996b
......@@ -7,7 +7,10 @@ fi
echo $QUEUE
#NCL_SCRIPT=histo2d
NCL_SCRIPT=wave_meanvar
#NCL_SCRIPT=wave_meanvar
#NCL_SCRIPT=total_meanvar
#NCL_SCRIPT=wk99_diff
NCL_SCRIPT=graphic_abstract
NCL_DIR=`pwd`
LOG_DIR=$NCL_DIR/log
mkdir -p $LOG_DIR
......
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; meanvar.ncl
; Carl Schreck (cjschrec@ncsu.edu)
; September 2014
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Description: Compare maps of variance or mean
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/time_axis_labels.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
load "$CJS_NCL_LIB/lib.cjs_graphics.ncl"
load "$CJS_NCL_LIB/lib.array.ncl"
load "$CJS_NCL_LIB/print_clock.ncl"
load "/home/carl/projects/olr_cdr/compare/sub.calc_meanvar.ncl"
load "/home/carl/projects/olr_cdr/compare/sub.draw_meanvar.ncl"
begin ; main
print_clock( "Here we go!" )
; These are some parameters that could be useful to have up top
plotType = "png"
plotDpi = 200
plotName = "figures/graphic_abstract"
; ...allows png or gif to work
if( ( plotType.eq."png" ).or.( plotType.eq."gif" ) ) then
plotTypeLocal = "eps"
else
plotTypeLocal = plotType
end if
; ...open the workstation
wks = gsn_open_wks( plotTypeLocal, plotName )
plotVar = "low"
res = True
cntRes = True
res@fontHeightF = 0.013
res@cnLevels = ispan(-5,5,1) * 1
res@mpMinLatF = -30
res@mpMaxLatF = 30
res@gsnLeftString = "MJO signal in HIRS (cont) and Difference with AVHRR (shaded)"
res@gsnRightString = ""
plotVar = "mjo"
mjoPlot = draw_meanvar( wks, plotVar, res, cntRes )
draw(mjoPlot)
frame(wks)
print_clock( "Convert the image, if necessary" )
delete(wks)
if( isStrSubset( plotType, "png" ).or.isStrSubset( plotType, "gif" ) ) then
system( "convert -trim -border 5x5 -bordercolor white " \
+ "+repage -density " + plotDpi + " " \\
+ plotName + ".eps " + plotName + "." + plotType )
end if
print_clock( "Thank you, come again." )
end; main
......@@ -44,11 +44,13 @@ begin ; main
res = True
res@fontHeightF = 0.015
res@cnLevels := ispan(-5,5,1) * 3
res@gsnLeftString = "a) Mean"
cntRes = True
meanPlot = draw_meanvar( wks, plotVar, res, cntRes )
plotVar = "variance"
res@cnLevels := ispan(-5,5,1) * 75
res@gsnLeftString = "b) Variance"
cntRes = True
varPlot = draw_meanvar( wks, plotVar, res, cntRes )
......
......@@ -49,12 +49,14 @@ begin ; main
res@mpMaxLatF = 30
res@lbLabelBarOn = False
res@tmXBLabelsOn = False
res@gsnLeftString = "a) Low frequency (> 120 days)"
lowPlot = draw_meanvar( wks, plotVar, res, cntRes )
plotVar = "mjo"
res = True
res@fontHeightF = 0.013
cntRes = True
res@gsnLeftString = "b) MJO"
mjoPlot = draw_meanvar( wks, plotVar, res, cntRes )
plotVar = "kelvin"
......@@ -62,6 +64,7 @@ begin ; main
res@fontHeightF = 0.013
cntRes = True
cntRes@cnLevels = ispan( 6, 20, 2 )
res@gsnLeftString = "c) Kelvin waves"
kelvinPlot = draw_meanvar( wks, plotVar, res, cntRes )
plotVar = "er"
......@@ -69,6 +72,7 @@ begin ; main
res@tmXBLabelsOn = True
res@fontHeightF = 0.013
cntRes = True
res@gsnLeftString = "c) ER waves"
erPlot = draw_meanvar( wks, plotVar, res, cntRes )
......@@ -78,7 +82,7 @@ begin ; main
panRes@gsnPanelBottom = 0.10
panRes@gsnPanelLabelBar = True
panRes@lbTitleString = "%"
panRes@lbTitlePosition = "Right"
panRes@lbTitlePosition = "Bottom"
panRes@lbTitleDirection = "Across"
panRes@lbOrientation = "Horizontal"
panRes@lbTitleFontHeightF = 0.02
......
......@@ -124,7 +124,7 @@ begin ; main
; pvalRes@cnLinesOn = False
res@cnLevels = ispan( -5, 5, 1 ) * 0.02
res@gsnLeftString = "Log Anti-symmetric"
res@gsnLeftString = "a) Log Anti-symmetric"
print_clock( res@gsnLeftString )
printVarSummary(diffAntiLog)
plotAntiLog = kf_draw( wks, diffAntiLog(freq|:,k|:), \
......@@ -132,7 +132,7 @@ begin ; main
cjs_add_contours( wks, plotAntiLog, aAntiLog(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotAntiLog, diffAntiLog(freq|:,k|:), pvalRes )
res@gsnLeftString = "Log Symmetric"
res@gsnLeftString = "b) Log Symmetric"
print_clock( res@gsnLeftString )
plotSymLog = kf_draw( wks, diffSymLog(freq|:,k|:), \
"None", True, False, res )
......@@ -141,14 +141,14 @@ begin ; main
res@cnLevels = ispan( -5, 5, 1 ) * 0.02
cntRes@cnLevels = cspan( 1.0, 2.0, 0.2 )
res@gsnLeftString = "Normalized Anti-symmetric"
res@gsnLeftString = "c) Normalized Anti-symmetric"
print_clock( res@gsnLeftString )
plotAntiNorm = kf_draw( wks, diffAntiNorm(freq|:,k|:), \
"None", False, True, res )
cjs_add_contours( wks, plotAntiNorm, aAntiNorm(freq|:,k|:), cntRes )
cjs_add_signif( wks, plotAntiNorm, diffAntiNorm(freq|:,k|:), pvalRes )
res@gsnLeftString = "Normalized Symmetric"
res@gsnLeftString = "d) Normalized Symmetric"
print_clock( res@gsnLeftString )
plotSymNorm = kf_draw( wks, diffSymNorm(freq|:,k|:), \
"None", True, False, res )
......
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