<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class=""><br class=""></div>Hi, looks to me like you need to add the line<div class=""><br class=""></div><div class="">draw(plots)</div><div class=""><br class=""></div><div class="">before </div><div class=""><br class=""></div><div class="">end</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div class="">
<div class="">-- <br class="">Alessandra Giannini<br class="">IRI for Climate and Society - The Earth Institute at Columbia University<br class="">P.O. Box 1000, Palisades NY 10964-8000<br class="">phone/fax: +1 845 680-4473/4864 - email: <a href="mailto:alesall@iri.columbia.edu" class="">alesall@iri.columbia.edu</a></div>
</div>
<div><br class=""><blockquote type="cite" class=""><div class="">On Feb 1, 2018, at 8:43 PM, dale zuri via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class="">Hi Ncl user,<br class=""></div>The program is not producing any error. But the figure is not plotted. <br class=""></div><div class="">I would appreciate any help and suggestion.</div><div class="">Thanks</div><div class="">AkilaS<br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div>begin<br class="">;***********************************************************<br class="">Model="t2m";;;;;;;change this line to include may forecast/march forecast<br class="">Model1="sst";;;;;;;change this line to include may forecast/march forecast<br class="">Model2="q2m";;;;;;;change this line to include may forecast/march forecast<br class="">Model3="t2mL";;;;;;;change this line to include may forecast/march forecast<br class="">;Box="23" ;;;;Both Temp/Prec<br class="">Box="12" ;;;;Both Temp/Prec<br class="">yrstrt="1980"<br class="">yrend="2017"<br class="">unit="Temperature(~F34~0~F~ C)"<br class="">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class="">f = addfile("<a href="http://era1-1980.nc/" class="">era1-1980.nc</a>","r")<br class="">u = f->sst ; units m/s<br class="">printVarSummary(u)<br class="">diri = "/import/c1/NCARCCM/asampath/ERA-Interim/ERA-Alaska-Fire-Summer"<br class="">fils = systemfunc ("csh -c 'cd " + diri + " ; ls era1*'")<br class="">print(fils) ; relative path to files<br class="">f1 = addfiles(fils,"r")<br class="">ListSetType (f1, "cat") ; concatenate (=default)`<br class="">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<br class="">t2air = short2flt(f1[:]->t2m(:,{55:45},{200:210})) ;air temp Ocean ; read T from all files<br class="">t2lan = short2flt(f1[:]->t2m(:,{70:60},{200:210})) ;air temp land ; read T from all files<br class="">tsst = short2flt(f1[:]->sst(:,{55:45},{200:210})) ;sea box ; read T from all files<br class="">DEW2= (dew-273.15);;;;dew point temperature<br class="">T2Air= (t2air-273.15);;;;2mtemperature Ocean<br class="">TSST= (tsst-273.15);;;;SeaSurface temperature<br class="">T2LAN= (t2lan-273.15);;;;air temp land<br class="">AvDEW2 = dim_avg_n_Wrap(DEW2, (/1,2/)) ; ==> zAvg(nlat,nlon)<br class="">AvT2Air = dim_avg_n_Wrap(T2Air, (/1,2/)) ; ==> zAvg(nlat,nlon)<br class=""><br class="">AvTSST = dim_avg_n_Wrap(TSST, (/1,2/)) ; ==> zAvg(nlat,nlon)<br class="">AvT2LAN = dim_avg_n_Wrap(T2LAN, (/1,2/)) ; ==> zAvg(nlat,nlon)<br class="">print(AvT2Air)<br class="">print(AvT2LAN)<br class=""> rc = regline_stats(AvT2Air,AvTSST) ; linear regression coef<br class=""> print(rc)<br class="">;;;;;;;;;;;Relative humidity calculation;;;;surface pressure in mb;;;specific humidity in kg/kg;;dew point temp in degree C;;;;;;;;;<br class="">e =6.112*(exp((17.625*AvDEW2)/(243.04+AvDEW2))) ;;;;;vapor pressure in mb<br class="">q = ((0.622 * e)/(1013.25 - (0.378 * e)))*1000; <br class="">print(q)<br class=""> rc1 = regline_stats(q,AvTSST) ; qOCean<br class="">print(rc1)<br class=""> rc2 = regline_stats(AvT2LAN,AvT2Air) ; land-Air<br class="">;************************************<br class="">qx = ispan (1,228,1)<br class="">;printVarSummary(qx)<br class="">;print(dimsizes(liness1))<br class="">;print(dimsizes(lines1))<br class="">;print(dimsizes(t2))<br class="">;;;;;;;Array adjustment;;;;;;;;;;;;;;;;<br class="">;gg=dimsizes(qx)-dimsizes(liness1)<br class="">;xgg = new(gg,float)<br class="">;print(xgg)<br class="">;lin1= array_append_record (xgg, liness1, 0)<br class="">;print(lin1)<br class="">;t resources;;;;;;;;;;;;;;;;;;O;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br class="">asciiwrite ("temporal_Clim_"+Model+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt", sprintf("%9.3f", AvT2Air))<br class="">asciiwrite ("temporal_Clim_"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt", sprintf("%9.3f", AvTSST))<br class="">asciiwrite ("temporal_Clim_"+Model2+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt", sprintf("%9.3f", q))<br class="">asciiwrite ("temporal_Clim_"+Model3+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+".txt", sprintf("%9.3f", AvT2LAN))<br class="">;wks = gsn_open_wks ("pdf", ""+Model+"-"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")<br class="">;wks = gsn_open_wks ("pdf", ""+Model2+"-"+Model1+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")<br class="">wks = gsn_open_wks ("x11", ""+Model3+"-"+Model+"_PSA_"+yrstrt+"-"+yrend+"_"+Box+"_timeseries_")<br class="">plots = new(1, graphic)<br class=""><br class="">res = True<br class="">res@gsnDraw = False<br class="">res@gsnFrame = False<br class="">res@vpHeightF= 0.25 ; change aspect ratio of plot<br class="">res@vpWidthF = 0.4<br class="">res@xyLineThicknessF = 4.<br class="">res@xyDashPattern = 0<br class="">res@xyMarker = 16<br class=""><br class="">;res@tmYLMinorOn = False ; Turn off left minor tickmarks<br class="">res@tmYLOn = True ; Turn off left tickmarks<br class="">res@tmYLMinorOn = False ; Turn off left minor tickmarks<br class="">res@tmYROn = False<br class="">res@tmXTOn = False<br class=""> ; now change the size of the tickmark labels<br class="">res@tmXBLabelFont = 22 ; resize tick labels<br class="">res@tmXBLabelFontHeightF = 0.009 ; resize tick labels<br class="">res@tmYLLabelFont = 22 ; resize tick labels<br class="">res@tmYLLabelFontHeightF = 0.013<br class="">res@tiMainFont= 22<br class="">res@tiMainFontHeightF = 0.013<br class="">res@tiXAxisFontHeightF = 0.013<br class="">res@tiYAxisFontHeightF = 0.013<br class="">res@tiXAxisFont = 22<br class="">res@tiYAxisFont = 22<br class="">res@tiXAxisString = "Time [years]"<br class="">;res@tiYAxisString = "mm/d"<br class="">res@tiYAxisString ="" <br class="">;*****************************<br class="">; panel plot only resources<br class="">resP = True ; modify the panel plot<br class="">resP@gsnMaximize = True ; large format<br class="">resP@gsnStringFont = "helvetica-bold"<br class="">;resP@gsnPanelFigureStrings = (/"a)","b)","c)","d)","e)","f)","g)","h)"/)<br class="">;resP@gsnPanelFigureStringsFontHeightF = 0.015<br class="">;resP@gsnPanelFigureStringsFont = "helvetica-bold"<br class=""><br class=""> txres = True<br class=""> txres@txFontHeightF = 0.02<br class=""> txres@txJusti = "CenterCenter"<br class=""> txres@txFontThicknessF = 2.0 ; default=1.00<br class="">; txres@txFontHeightF = 0.025 ; default=0.05<br class=""><br class="">;res@xyLineColors = (/ "Red", "Green","Red","Green"/)<br class="">;res@xyMarkLineMode = "MarkLines" ; Markers *and* lines<br class="">;res@xyMarkers = (/4,4/) ; 3 different markers<br class="">;res@xyMarkerColor = (/"Red","blue"/)<br class="">res@xyLineColors = (/"Red","blue"/)<br class="">;;;;;;;;;;;;;;;;;;Legend Label;;;;;;;;;;;;;;;;;;<br class="">;res@xyExplicitLegendLabels = (/"T2m(Air)","SST"/)<br class="">res@xyExplicitLegendLabels = (/"q2m(g/kg)","SST"/)<br class="">;res@xyExplicitLegendLabels = (/"T2m(Air)","T2m(Land)"/)<br class="">res@pmLegendSide = "Top"<br class=""> res@pmLegendDisplayMode = "Always" ; Display a legend.<br class=""> res@pmLegendWidthF = 0.1 ; Make it smaller<br class=""> res@pmLegendHeightF = 0.05 ; Make it smaller<br class=""> res@pmLegendOrthogonalPosF = -0.4 ; Move into to plot<br class=""> res@pmLegendParallelPosF = 0.5 ; Move to right<br class=""> res@lgPerimOn = False ; No legend perimeter.<br class=""> res@lgLabelFont = 22 <br class=""> res@lgLabelFontHeightF = 0.013<br class="">;;;;;;;;;;;;Tick mark line label;;;;;;;;;;;;;;;<br class="">;res@tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc)+") -"+yrstrt+"-"+yrend+""<br class="">res@tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc1)+") -"+yrstrt+"-"+yrend+""<br class="">;res@tiMainString = ""+Box+"-Regcoeff("+sprintf("%1.0g",rc2)+") -"+yrstrt+"-"+yrend+""<br class="">res@tmXBMode = "Explicit"<br class="">;res@trYMinF = -(max(zAvg))<br class="">;res@trYMaxF = (max(zAvg1))<br class="">res@trXMinF = 0<br class="">res@trXMaxF =228 <br class="">res@tmXBValues = ispan(0,227,6) <br class="">res@tmXBLabels = (/"1980","1981","1982","1983","1984","1985","1986", \<br class=""> "1987","1988","1989","1990","1991","1992","1993","1994","1995", \<br class=""> "1996","1997","1998","1999","2000","2001","2002","2003","2004", \<br class=""> "2005","2006","2007","2008","2009","2010","2011","2012","2013", \<br class=""> "2014","2015","2016","2017"/)<br class="">res@tmXBLabelAngleF = 80 <br class="">plots(0) = gsn_csm_xy(wks,qx,(/AvT2Air,AvTSST/),res) ; create plot<br class="">;plots(0) = gsn_csm_xy(wks,qx,(/q,AvTSST/),res) ; create plot<br class="">;plots(0) = gsn_csm_xy(wks,qx,(/AvT2LAN,AvT2Air/),res) ; create plot<br class="">end<br class=""><br class=""></div>
_______________________________________________<br class="">ncl-talk mailing list<br class=""><a href="mailto:ncl-talk@ucar.edu" class="">ncl-talk@ucar.edu</a><br class="">List instructions, subscriber options, unsubscribe:<br class="">http://mailman.ucar.edu/mailman/listinfo/ncl-talk<br class=""></div></blockquote></div><br class=""></div></body></html>