[ncl-talk] Sub: Figure is disappearing
dale zuri
dalezuri at gmail.com
Thu Feb 1 18:43:58 MST 2018
Hi Ncl user,
The program is not producing any error. But the figure is not plotted.
I would appreciate any help and suggestion.
Thanks
AkilaS
begin
;***********************************************************
Model="t2m";;;;;;;change this line to include may forecast/march forecast
Model1="sst";;;;;;;change this line to include may forecast/march forecast
Model2="q2m";;;;;;;change this line to include may forecast/march forecast
Model3="t2mL";;;;;;;change this line to include may forecast/march forecast
;Box="23" ;;;;Both Temp/Prec
Box="12" ;;;;Both Temp/Prec
yrstrt="1980"
yrend="2017"
unit="Temperature(~F34~0~F~ C)"
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
f = addfile("era1-1980.nc","r")
u = f->sst ; units m/s
printVarSummary(u)
diri = "/import/c1/NCARCCM/asampath/ERA-Interim/ERA-Alaska-Fire-Summer"
fils = systemfunc ("csh -c 'cd " + diri + " ; ls era1*'")
print(fils) ; relative path to files
f1 = addfiles(fils,"r")
ListSetType (f1, "cat") ; concatenate (=default)`
dew = short2flt(f1[:]->d2m(:,{55:45},{200:210})) ; read T
from all filesu = short2flt(f1[:]->t2m(:,{55:45},{200:210})) ;air temp
Ocean ; read T from all files
t2air = short2flt(f1[:]->t2m(:,{55:45},{200:210})) ;air temp
Ocean ; read T from all files
t2lan = short2flt(f1[:]->t2m(:,{70:60},{200:210})) ;air temp
land ; read T from all files
tsst = short2flt(f1[:]->sst(:,{55:45},{200:210})) ;sea box
; read T from all files
DEW2= (dew-273.15);;;;dew point temperature
T2Air= (t2air-273.15);;;;2mtemperature Ocean
TSST= (tsst-273.15);;;;SeaSurface temperature
T2LAN= (t2lan-273.15);;;;air temp land
AvDEW2 = dim_avg_n_Wrap(DEW2, (/1,2/)) ; ==> zAvg(nlat,nlon)
AvT2Air = dim_avg_n_Wrap(T2Air, (/1,2/)) ; ==> zAvg(nlat,nlon)
AvTSST = dim_avg_n_Wrap(TSST, (/1,2/)) ; ==> zAvg(nlat,nlon)
AvT2LAN = dim_avg_n_Wrap(T2LAN, (/1,2/)) ; ==> zAvg(nlat,nlon)
print(AvT2Air)
print(AvT2LAN)
rc = regline_stats(AvT2Air,AvTSST) ; linear regression coef
print(rc)
;;;;;;;;;;;Relative humidity calculation;;;;surface pressure in
mb;;;specific humidity in kg/kg;;dew point temp in degree C;;;;;;;;;
e =6.112*(exp((17.625*AvDEW2)/(243.04+AvDEW2))) ;;;;;vapor pressure in mb
q = ((0.622 * e)/(1013.25 - (0.378 * e)))*1000;
print(q)
rc1 = regline_stats(q,AvTSST) ; qOCean
print(rc1)
rc2 = regline_stats(AvT2LAN,AvT2Air) ; land-Air
;************************************
qx = ispan (1,228,1)
;printVarSummary(qx)
;print(dimsizes(liness1))
;print(dimsizes(lines1))
;print(dimsizes(t2))
;;;;;;;Array adjustment;;;;;;;;;;;;;;;;
;gg=dimsizes(qx)-dimsizes(liness1)
;xgg = new(gg,float)
;print(xgg)
;lin1= array_append_record (xgg, liness1, 0)
;print(lin1)
;t resources;;;;;;;;;;;;;;;;;;O;;;;;;;;;;;;;;;;;;;;;;;;;;;;
asciiwrite ("temporal_Clim_"+Model+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt",
sprintf("%9.3f", AvT2Air))
asciiwrite
("temporal_Clim_"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt",
sprintf("%9.3f", AvTSST))
asciiwrite
("temporal_Clim_"+Model2+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt",
sprintf("%9.3f", q))
asciiwrite
("temporal_Clim_"+Model3+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt",
sprintf("%9.3f", AvT2LAN))
;wks = gsn_open_wks ("pdf",
""+Model+"-"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")
;wks = gsn_open_wks ("pdf",
""+Model2+"-"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")
wks = gsn_open_wks ("x11",
""+Model3+"-"+Model+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")
plots = new(1, graphic)
res = True
res at gsnDraw = False
res at gsnFrame = False
res at vpHeightF= 0.25 ; change aspect ratio of plot
res at vpWidthF = 0.4
res at xyLineThicknessF = 4.
res at xyDashPattern = 0
res at xyMarker = 16
;res at tmYLMinorOn = False ; Turn off left minor tickmarks
res at tmYLOn = True ; Turn off left tickmarks
res at tmYLMinorOn = False ; Turn off left minor tickmarks
res at tmYROn = False
res at tmXTOn = False
; now change the size of the tickmark labels
res at tmXBLabelFont = 22 ; resize tick labels
res at tmXBLabelFontHeightF = 0.009 ; resize tick labels
res at tmYLLabelFont = 22 ; resize tick labels
res at tmYLLabelFontHeightF = 0.013
res at tiMainFont= 22
res at tiMainFontHeightF = 0.013
res at tiXAxisFontHeightF = 0.013
res at tiYAxisFontHeightF = 0.013
res at tiXAxisFont = 22
res at tiYAxisFont = 22
res at tiXAxisString = "Time [years]"
;res at tiYAxisString = "mm/d"
res at tiYAxisString =""
;*****************************
; panel plot only resources
resP = True ; modify the panel plot
resP at gsnMaximize = True ; large format
resP at gsnStringFont = "helvetica-bold"
;resP at gsnPanelFigureStrings = (/"a)","b)","c)","d)","e)","f)","g)","h)"/)
;resP at gsnPanelFigureStringsFontHeightF = 0.015
;resP at gsnPanelFigureStringsFont = "helvetica-bold"
txres = True
txres at txFontHeightF = 0.02
txres at txJusti = "CenterCenter"
txres at txFontThicknessF = 2.0 ; default=1.00
; txres at txFontHeightF = 0.025 ; default=0.05
;res at xyLineColors = (/ "Red", "Green","Red","Green"/)
;res at xyMarkLineMode = "MarkLines" ; Markers *and* lines
;res at xyMarkers = (/4,4/) ; 3 different markers
;res at xyMarkerColor = (/"Red","blue"/)
res at xyLineColors = (/"Red","blue"/)
;;;;;;;;;;;;;;;;;;Legend Label;;;;;;;;;;;;;;;;;;
;res at xyExplicitLegendLabels = (/"T2m(Air)","SST"/)
res at xyExplicitLegendLabels = (/"q2m(g/kg)","SST"/)
;res at xyExplicitLegendLabels = (/"T2m(Air)","T2m(Land)"/)
res at pmLegendSide = "Top"
res at pmLegendDisplayMode = "Always" ; Display a legend.
res at pmLegendWidthF = 0.1 ; Make it smaller
res at pmLegendHeightF = 0.05 ; Make it smaller
res at pmLegendOrthogonalPosF = -0.4 ; Move into to plot
res at pmLegendParallelPosF = 0.5 ; Move to right
res at lgPerimOn = False ; No legend perimeter.
res at lgLabelFont = 22
res at lgLabelFontHeightF = 0.013
;;;;;;;;;;;;Tick mark line label;;;;;;;;;;;;;;;
;res at tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc)+")
-"+yrstrt+"-"+yrend+""
res at tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc1)+")
-"+yrstrt+"-"+yrend+""
;res at tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc2)+")
-"+yrstrt+"-"+yrend+""
res at tmXBMode = "Explicit"
;res at trYMinF = -(max(zAvg))
;res at trYMaxF = (max(zAvg1))
res at trXMinF = 0
res at trXMaxF =228
res at tmXBValues = ispan(0,227,6)
res at tmXBLabels = (/"1980","1981","1982","1983","1984","1985","1986", \
"1987","1988","1989","1990","1991","1992","1993","1994","1995", \
"1996","1997","1998","1999","2000","2001","2002","2003","2004", \
"2005","2006","2007","2008","2009","2010","2011","2012","2013", \
"2014","2015","2016","2017"/)
res at tmXBLabelAngleF = 80
plots(0) = gsn_csm_xy(wks,qx,(/AvT2Air,AvTSST/),res) ; create plot
;plots(0) = gsn_csm_xy(wks,qx,(/q,AvTSST/),res) ; create plot
;plots(0) = gsn_csm_xy(wks,qx,(/AvT2LAN,AvT2Air/),res) ; create plot
end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180201/59832451/attachment.html>
More information about the ncl-talk
mailing list