Hi everyone,<div><br></div><div>       Now I want to draw a figure that contains many variables from different files under Grads.  The script is as follows.<br><br>&#39;reinit&#39;<br>&#39;set grid off&#39;<br>&#39;set grads off&#39;<br>
<br>*** for Windows<br>*&#39;sdfopen C:\daily4times\<a href="http://hgt.2012.nc">hgt.2012.nc</a>&#39; <br>*&#39;sdfopen C:\daily4times\<a href="http://hgt.meanrunmean.1981-2010.nc">hgt.meanrunmean.1981-2010.nc</a>&#39;<br>
*&#39;sdfopen C:\daily4times\<a href="http://hgt.meanrunstd.1981-2010.nc">hgt.meanrunstd.1981-2010.nc</a>&#39;<br><br>*** for Linux<br>&#39;sdfopen /host/daily4times/<a href="http://hgt.2012.nc">hgt.2012.nc</a>&#39; <br>&#39;sdfopen /host/daily4times/<a href="http://hgt.meanrunmean.1981-2010.nc">hgt.meanrunmean.1981-2010.nc</a>&#39;<br>
&#39;sdfopen /host/daily4times/<a href="http://hgt.meanrunstd.1981-2010.nc">hgt.meanrunstd.1981-2010.nc</a>&#39;<br><br><br><br>*** for output figure (type: eps)<br>*&#39;enable print D:\Writing\disaster_management\hgtgap.gmf&#39;<br>
*&#39;enable print dif_hgt.eps&#39;<br>&#39;set parea 1 9 1 8&#39;<br>&#39;set lat 10 60&#39;<br>&#39;set xlint 10&#39;<br>&#39;set lon 100 160&#39;<br>&#39;set ylint 10&#39;<br>&#39;set cint 0.5&#39;<br>&#39;set clskip 1&#39;<br>
&#39;set rbcols 14 9 4 11 5 13 0 10 7 12 8 2 6&#39;<br>&#39;set strmden 1&#39;<br>&#39;set clevs -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6&#39;<br>&#39;set xlopts 1 4 0.15&#39;<br>&#39;set ylopts 1 4 0.15&#39;<br>&#39;define stdanomalies=((hgt.1(time=12Z14JUL2012,LEV=500)-hgt.2(time=06Z29DEC2010,LEV=500))/hgt.3(time=06Z29DEC2010,LEV=500))&#39;<br>
&#39;set gxout shaded&#39;<br>&#39;d stdanomalies&#39;<br>&#39;cbarn 1 0.5 9.3 4.5&#39;<br><span style="color:rgb(255,0,0)">&#39;open /host/GFS/20120714/gfs_3_20120714_1200_fff.ctl&#39;<br>&#39;set dfile 2&#39;<br>&#39;set lat 10 60&#39;<br>
&#39;set xlint 10&#39;<br>&#39;set lon 100 160&#39;<br>&#39;set ylint 10&#39;<br>&#39;set gxout contour&#39;<br>&#39;set t 1&#39;<br>&#39;set lev 500&#39;<br>&#39;d hgt.4&#39;<br>&#39;draw title a.GFS 500 hPa hgtprs 12Z14JUL2012&#39;</span><br>
<br>*** for output figure (type: eps)<br>*&#39;print<br>*&#39;disable print&#39;<br>*** for outut figure (type: png)<br>&#39;printim \host\daily4times\hgtstdanomalies_500_20120714.png white&#39;<br></div><div><br></div><div>
    As you see, it contains two parts. One is to display &#39;d stdanomalies&#39;, the other is to <span style="color:rgb(255,0,0)">&#39;d hgt.4&#39;. <span style="color:rgb(0,0,0)">I can display the first one, but cannot display the second one. Because the second variable is come from GFS forecast data(</span></span><span style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)"><span style="color:rgb(255,0,0)">GFS/20120714/gfs_3_20120714_1200_fff.ctl&#39;).<span style="color:rgb(0,0,0)"> But I think Grads can display different variables from different files.<br>
    So, anyone can help me how to work it.  Thank you very much!!!!<br><br>   Best wishes,<br><br>Duan<br><br>  <br></span></span> <img alt="Inline image 1" src="cid:ii_13cd9125854edbb5"></span><br></span><div><br></div>-- <br>
Weili Duan<br>Disaster Prevention Research Institute (DPRI),
Kyoto University, Uji, Kyoto 611-0011, Japan<br>TEL: <a href="tel:%2B81-774-38-4131" value="+81774384131" target="_blank">+81-774-38-4131</a> FAX: <a href="tel:%2B81-774-38-4130" value="+81774384130" target="_blank">+81-774-38-4130</a><br>
E-mail: <a href="mailto:duan.scut.cn@gmail.com" target="_blank">duan.scut.cn@gmail.com</a><br>           <a href="mailto:dyuan.weiri.88r@st.kyoto-u.ac.jp" target="_blank">duanwl1986@flood.dpri.kyoto-u.ac.jp</a><a href="mailto:dyuan.weiri.88r@st.kyoto-u.ac.jp" target="_blank"></a>
</div>