<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
Hi Jason,<br>
<br>
the y-axis ranges should be set for all plots. In the script below
(I threw out everything that is not needed for readability) trYMinF,
trYMaxF for the left and right y-axis is set for all plots. It
should go like this.<br>
<br>
-Karin<br>
<br>
<font face="monospace" color="blue">begin<br>
<br>
f0 = "County001-SO2-0012-2013-18full-line.txt"<br>
f1 = "County013-SO2-0006-2013-18full-line.txt" <br>
f2 = "County027-SO2-0007-2013-18full-line.txt"<br>
f3 = "County029-SO2-0002-2013-18full-line.txt"<br>
f4 = "County031-SO2-0003-2013-18full-line.txt"<br>
f5 = "County043-SO2-0005-2013-18full-line.txt"<br>
f6 = "County055-SO2-1007-2013-18full-line.txt"<br>
f7 = "County067-SO2-1015-2013-18full-line.txt"<br>
f8 = "County081-SO2-0124-2013-18full-line.txt"<br>
f9 = "County101-SO2-0003-2013-18full-line.txt"<br>
f10 = "County103-SO2-0009-2013-18full-line.txt"<br>
<br>
data0 = asciiread(f0,(/2191,3/),"double")<br>
data1 = asciiread(f1,(/2191,3/),"double")<br>
data2 = asciiread(f2,(/2191,3/),"double")<br>
data3 = asciiread(f3,(/2191,3/),"double")<br>
data4 = asciiread(f4,(/2191,3/),"double")<br>
data5 = asciiread(f5,(/2191,3/),"double")<br>
data6 = asciiread(f6,(/2191,3/),"double")<br>
data7 = asciiread(f7,(/2191,3/),"double")<br>
data8 = asciiread(f8,(/2191,3/),"double")<br>
data9 = asciiread(f9,(/2191,3/),"double")<br>
data10 = asciiread(f10,(/2191,3/),"double")<br>
<br>
xlabelc = (/"2013","2014","2015","2016","2017","2018"/)<br>
<br>
<br>
; DEC-GEOS dates in<br>
x0 = data0(:,0)<br>
x1 = data1(:,0)<br>
x2 = data2(:,0)<br>
x3 = data3(:,0)<br>
x4 = data4(:,0)<br>
x5 = data5(:,0)<br>
x6 = data6(:,0)<br>
x7 = data7(:,0)<br>
x8 = data8(:,0)<br>
x9 = data9(:,0)<br>
x10 = data10(:,0)<br>
<br>
; DEC data in<br>
y0 = data0(:,1)<br>
y0@_FillValue=integertoshort(-999)<br>
y0=where(y0.lt.-5.0,y0@FillValue,y0) <br>
y1 = data1(:,1)<br>
y1@_FillValue=integertoshort(-999)<br>
y1=where(y1.lt.-5.0,y1@FillValue,y1)<br>
y2 = data2(:,1)<br>
y2@_FillValue=integertoshort(-999)<br>
y2=where(y2.lt.-5.0,y2@FillValue,y2)<br>
y3 = data3(:,1)<br>
y3@_FillValue=integertoshort(-999)<br>
y3=where(y3.lt.-5.0,y3@FillValue,y3)<br>
y4 = data4(:,1)<br>
y4@_FillValue=integertoshort(-999)<br>
y4=where(y4.lt.-5.0,y4@FillValue,y4)<br>
y5 = data5(:,1)<br>
y5@_FillValue=integertoshort(-999)<br>
y5=where(y5.lt.-5.0,y5@FillValue,y5)<br>
y6 = data6(:,1)<br>
y6@_FillValue=integertoshort(-999)<br>
y6=where(y6.lt.-5.0,y6@FillValue,y6)<br>
y7 = data7(:,1)<br>
y7@_FillValue=integertoshort(-999)<br>
y7=where(y7.lt.-5.0,y7@FillValue,y7)<br>
y8 = data8(:,1)<br>
y8@_FillValue=integertoshort(-999)<br>
y8=where(y8.lt.-5.0,y8@FillValue,y8)<br>
y9 = data9(:,1)<br>
y9@_FillValue=integertoshort(-999)<br>
y9=where(y9.lt.-5.0,y9@FillValue,y9)<br>
y10 = data10(:,1)<br>
y10@_FillValue=integertoshort(-999)<br>
y10=where(y10.lt.-5.0,y10@FillValue,y10)<br>
<br>
; GEOS-Chem data in <br>
<br>
GCD0 = data0(:,2)<br>
GCD0@_FillValue=integertoshort(-999)<br>
GCD0=where(GCD0.lt.-5.0,GCD0@FillValue,GCD0)<br>
GCD1 = data1(:,2)<br>
GCD1@_FillValue=integertoshort(-999)<br>
GCD1=where(GCD1.lt.-5.0,GCD1@FillValue,GCD1)<br>
GCD2 = data2(:,2)<br>
GCD2@_FillValue=integertoshort(-999)<br>
GCD2=where(GCD2.lt.-5.0,GCD2@FillValue,GCD2)<br>
GCD3 = data3(:,2)<br>
GCD3@_FillValue=integertoshort(-999)<br>
GCD3=where(GCD3.lt.-5.0,GCD3@FillValue,GCD3)<br>
GCD4 = data4(:,2)<br>
GCD4@_FillValue=integertoshort(-999)<br>
GCD4=where(GCD4.lt.-5.0,GCD4@FillValue,GCD4)<br>
GCD5 = data5(:,2)<br>
GCD5@_FillValue=integertoshort(-999)<br>
GCD5=where(GCD5.lt.-5.0,GCD5@FillValue,GCD5)<br>
GCD6 = data6(:,2)<br>
GCD6@_FillValue=integertoshort(-999)<br>
GCD6=where(GCD6.lt.-5.0,GCD6@FillValue,GCD6)<br>
GCD7 = data7(:,2)<br>
GCD7@_FillValue=integertoshort(-999)<br>
GCD7=where(GCD7.lt.-5.0,GCD7@FillValue,GCD7)<br>
GCD8 = data8(:,2)<br>
GCD8@_FillValue=integertoshort(-999)<br>
GCD8=where(GCD8.lt.-5.0,GCD8@FillValue,GCD8)<br>
GCD9 = data9(:,2)<br>
GCD9@_FillValue=integertoshort(-999)<br>
GCD9=where(GCD9.lt.-5.0,GCD9@FillValue,GCD9)<br>
GCD10 = data10(:,2)<br>
GCD10@_FillValue=integertoshort(-999)<br>
GCD10=where(GCD10.lt.-5.0,GCD10@FillValue,GCD10)<br>
<br>
;************************************************<br>
; create plots <br>
;************************************************<br>
plot = new(11,graphic)<br>
<br>
wks = gsn_open_wks("png","SO2_Multisite_plot.png") ;
ps,pdf,x11,png<br>
<br>
trxmin = 2013 ; min/max values for X axis <br>
trxmax = 2018<br>
<br>
trymin = -2.<br>
trymax = 8.<br>
trymin2 = 0.<br>
trymax2 = 8.<br>
<br>
res = True ; plot mods desired<br>
res@gsnDraw = False<br>
res@gsnFrame = False<br>
<br>
res@vpWidthF = 0.70 ; change aspect ratio of
plot<br>
res@vpHeightF = 0.25<br>
<br>
res@trXMinF = trxmin<br>
res@trXMaxF = trxmax<br>
res@trXReverse = False ; Reverse X axis values<br>
<br>
res@trYMinF = trymin2<br>
res@trYMaxF = trymax2<br>
<br>
res@xyLineThicknessF = 2.0<br>
res@xyLineColor = (/"Red"/)<br>
<br>
res1 = res<br>
res2 = res<br>
res3 = res<br>
res4 = res<br>
res5 = res<br>
res6 = res<br>
res7 = res<br>
res8 = res<br>
res9 = res<br>
res10 = res<br>
res11 = res<br>
res12 = res<br>
<br>
res@xyLineColor = (/"Black"/)<br>
res@trYMinF = trymin<br>
res@trYMaxF = trymax<br>
<br>
res12@xyLineColor = (/"Black"/)<br>
<br>
;**************************************************************************<br>
<br>
res@tiYAxisString = "SO2 (ppb)"<br>
res@tiYAxisFontColor = res@xyLineColor<br>
res@tiYAxisFontHeightF = 0.05<br>
<br>
res@tmXBOn = False<br>
<br>
plot(0) = gsn_csm_xy2 (wks,x0,y0,GCD0,res,res1) ; create plot<br>
plot(1) = gsn_csm_xy2 (wks,x1,y1,GCD1,res,res3) ; create plot<br>
plot(2) = gsn_csm_xy2 (wks,x2,y2,GCD2,res,res4) ; create plot<br>
plot(3) = gsn_csm_xy2 (wks,x3,y3,GCD3,res,res12) ; create plot<br>
plot(4) = gsn_csm_xy2 (wks,x4,y4,GCD4,res,res5) ; create plot<br>
plot(5) = gsn_csm_xy2 (wks,x5,y5,GCD5,res,res6) ; create plot<br>
plot(6) = gsn_csm_xy2 (wks,x6,y6,GCD6,res,res7) ; create plot<br>
plot(7) = gsn_csm_xy2 (wks,x7,y7,GCD7,res,res8) ; create plot<br>
plot(8) = gsn_csm_xy2 (wks,x8,y8,GCD8,res,res9) ; create plot<br>
plot(9) = gsn_csm_xy2 (wks,x9,y9,GCD9,res,res10) ; create plot<br>
plot(10) = gsn_csm_xy2 (wks,x10,y10,GCD10,res,res11) ; create
plot<br>
<br>
<br>
gsn_panel(wks,plot,(/6,2/),False)<br>
<br>
end</font><br>
<br>
<br>
<br>
<div class="moz-cite-prefix">Am 08.12.21 um 01:48 schrieb Herb,
Jason via ncl-talk:<br>
</div>
<blockquote type="cite"
cite="mid:BL3PR04MB807527B1E0161A92714AB401CE6F9@BL3PR04MB8075.namprd04.prod.outlook.com">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style type="text/css" style="display:none;">P {margin-top:0;margin-bottom:0;}</style>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
Hello All,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
I am attempting to do a series of plots together; Black lines
are surface observations and red lines are model simulation. The
issue I am having is that each plot separate box plot should
have the same scale as the upper right most plot. However, NCL
is putting varying scale ranges as can be seen in the attached
plot. Is there a simpler way for this to be done automatically
by NCL than doing the ranges by hand? As I have other plots will
be of similar style. I have attached the script as well.<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
Thanks,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif;
font-size: 12pt; color: rgb(0, 0, 0);">
Jason<br>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
ncl-talk mailing list
<a class="moz-txt-link-abbreviated" href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>
List instructions, subscriber options, unsubscribe:
<a class="moz-txt-link-freetext" href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Dipl. Geophys. Karin Meier-Fleischer
Visualization Group - NCL, CDO, Python
Application Support
Deutsches Klimarechenzentrum GmbH (DKRZ)
Bundesstrasse 45a - D-20146 Hamburg - Germany
Phone: +49 (0)40 460094 126
Fax: +49 (0)40 460094 270
E-Mail: <a class="moz-txt-link-abbreviated" href="mailto:meier-fleischer@dkrz.de">meier-fleischer@dkrz.de</a>
URL: <a class="moz-txt-link-abbreviated" href="http://www.dkrz.de">www.dkrz.de</a>
Geschäftsführer: Prof. Dr. Thomas Ludwig
Sitz der Gesellschaft: Hamburg
Amtsgericht Hamburg HRB 39784</pre>
</body>
</html>