<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>hello,</DIV>
<DIV>it seems you can't type " rc = gsfallow("on") " directly in the command line, but instead you can write it in a script file (.gs). here is an example for using panels() function to draw multi-panel plot, hope it can help you.</DIV>
<DIV> </DIV>
<DIV>* draw multi-panel plot, the example is 4 panels<BR>*Modified By Jian Zeng from beijing meteorological bureau<BR>*ensure panels.gsf exists or download it from <A href="ftp://grads.iges.org/grads/scripts/panels.gsf">ftp://grads.iges.org/grads/scripts/panels.gsf</A><BR>*<BR>'reinit'<BR>rc = gsfallow("on")<BR>nrows=2<BR>ncols=2 <BR>args=nrows%' '%ncols<BR>'sdfopen E:\some\place\somedata.dat'</DIV>
<DIV>'set mproj scaled'<BR>panels(args)<BR>*p is the panel number, p=1 means the first panel<BR>p = 1<BR>*ptot indicate the total numbers of panels (or virtual pages)<BR>ptot = nrows * ncols<BR>*here is the loop, draw panel one by one<BR>while (p <= ptot)<BR>*_vpg.p seting the virtual page, NO.p<BR> _vpg.p<BR> 'set parea 1 9.8 1 7.5'<BR> 'set t 'p<BR> 'set grads off'<BR> 'set gxout shaded'<BR> 'set font 2'<BR> 'set xlopts 1 4 0.2';'set ylopts 1 4 0.2' <BR> 'd skt'<BR> 'set font 0'<BR> 'cbar_interp 1 1 1'<BR> 'zb 30 30 0.1'<BR> 'set font 5';'set strsiz 0.3 0.3'<BR> 'draw string 3.5 8 skt(t='p',panel='p')' <BR> p = p + 1<BR>endwhile<BR> 'set vpage off'<BR>'printim e:\some\place\somepic.png png white x800 y600'<BR>;<BR><BR>Jian Zeng<BR><BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV><PRE><BR>At 2012-09-18 18:41:51,"TYTECA Sophie" <sophie.tyteca@meteo.fr> wrote:
>Hello,
>
>I want to use GrADS script function "panels.gsf" to draw multi-panel plot
>all scripts are here: /opt/GrADS/lib and GASCRP=/opt/GrADS/lib
>I see in panels_demo.gs that I have to use "gsfallow" first, right?
>See below my try :
>
>Grid Analysis and Display System (GrADS) Version 2.0.1
>Copyright (c) 1988-2011 by Brian Doty and the
>Institute for Global Environment and Society (IGES)
>GrADS comes with ABSOLUTELY NO WARRANTY
>See file COPYRIGHT for more information
>
>Config: v2.0.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5
>opendap-grids,stn geotiff shapefile
>Issue 'q config' command for more detailed configuration information
>GX Package Initialization: Size = 8.5 11
>ga-> open PLGEO500.ctl
>Scanning description file: PLGEO500.ctl
>Data file PLGEO500.SG3.1959.dat is open as file 1
>LON set to 0 360
>LAT set to -87.864 87.864
>LEV set to 500 500
>Time values set: 1959:1:1:0 1959:1:1:0
>E set to 1 1
>ga-> rc = gsfallow("on")
>Syntax Error: Invalid Operand
> 'gsfallow' not a variable or function name
> Error ocurred at column 1
>DEFINE error: Invalid expression.
>
>thank you for your help
>_______________________________________________
>gradsusr mailing list
>gradsusr@gradsusr.org
>http://gradsusr.org/mailman/listinfo/gradsusr
</PRE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>