<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.0.8">
</HEAD>
<BODY>
Hi Quin Hao,<BR>
<BR>
Now it works. Many Thanks.<BR>
<BR>
Paolo<BR>
<BR>
On Fri, 2011-04-01 at 14:30, Qian Hao wrote:
<BLOCKQUOTE TYPE=CITE>
<FONT COLOR="#737373"><I>Well, hello Paolo,<BR>
I found that 'set dfile' could not change the dimension configuration of grads.When you open a file for the first time, grads will set up its dimensional coordinates accordingly.When you open a second file, the dimensional coordinates do not change. You may check this with the 'q dims' command.<BR>
And that's the problem. So you have to close the first opened file first if you want to display another variable from a second data file with a different coordinate settings. You my change you .gs file like this:<BR>
<BR>
'open suolo.ctl'<BR>
...<BR>
'close 1'<BR>
'open msg.ctl'<BR>
...<BR>
<BR>
I think this would work.<BR>
2011/4/1 Paolo Boi <</FONT><A HREF="mailto:pboi@arpa.sardegna.it"><U>pboi@arpa.sardegna.it</U></A><FONT COLOR="#737373">>
<BLOCKQUOTE>
Thanks Quian Hao,<BR>
<BR>
The .gs script contains the lines 'd file1' and 'd file2', but it plots only the first data file, no matter it is. I tried to exchange the two files, but it plots only the first one.<BR>
Here is the script .gs:<BR>
<BR>
'reinit'<BR>
<BR>
* Palette colori<BR>
* --------------<BR>
<BR>
'run /opt/grads-1.8sl11/lib/colors/BYR-02'<BR>
<BR>
* Apertura file .ctl<BR>
* ------------------<BR>
<BR>
'open suolo.ctl'<BR>
'open msg.ctl'<BR>
<BR>
<BR>
* Mappa<BR>
* -----<BR>
<BR>
'set mpdraw on'<BR>
'set display color black'<BR>
'set mproj nps'<BR>
'set mpdset hires'<BR>
'set map 13'<BR>
'set lon 0 20'<BR>
'set lat 35 45'<BR>
<BR>
*<BR>
i=1<BR>
*ntime=1<BR>
*while (i<=ntime)<BR>
'clear'<BR>
'set t 'i<BR>
*<BR>
* Non riesce a passare da un grafico all'altro.<BR>
* Grafica il primo file (sia esso mslp o msg) ma il secondo no<BR>
* se i due grib hanno la stessa risoluzione (vedi geop e temp 850)<BR>
* allora li grafica<BR>
* perch??<BR>
*<BR>
* Traccia le isolinee di mslp<BR>
* -------------------<BR>
'set dfile 1'<BR>
<BR>
'set gxout contour'<BR>
'set cint 4 '<BR>
'set ccolor 0'<BR>
'set csmooth on'<BR>
'set clopts 1 1 0.08'<BR>
'set cthick 6'<BR>
'd MSLsfc/100'<BR>
<BR>
<BR>
'set dfile 2'<BR>
<BR>
* Plot shading<BR>
* ------------<BR>
'set gxout shaded'<BR>
*'set ccols 13 5 10 3 2 6 9 8 16 18 20 22 24 26 28 30 40 42 44 46 48 50 52 54 56 58 60 62 64 66 68 70 72 74 76 '<BR>
*'set clevs -80 -70 -65 -60 -58 -56 -54 -52 -50 -48 -46 -44 -42 -40 -38 -36 -34 -32 -30 -28 -26 -24 -22 -20 -15 -10 -5 0 2 4 6 8 10 12 14 16 18 20 22 24 26 30 40 '<BR>
'set clevs 0 3 9 15 21 27 30 36 41 47 50 60 66 70 75 80 86 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 205 210 215 220 225 230 235 240 245 250 255'<BR>
'set ccols 5 13 10 3 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 54 52 50 48 46 44 42 40 38 36 34 32 30 28 26 24 22 20 18 16'<BR>
'd IMGD/4 '<BR>
<BR>
<BR>
* legenda<BR>
<BR>
'run cbarn 1 1'<BR>
<BR>
'q time'<BR>
res = subwrd(result,3)<BR>
data=substr(res,1,15)<BR>
<BR>
'set string 1 tc 8'<BR>
'set strsiz 0.13 0.16'<BR>
'draw title IR T(C) 'data''<BR>
*'draw title IR T(C) 'day' 'month' 'year' 'hour''<BR>
'draw string 8.2 0.5 Dati EUMETSAT Elaborazioni SAR'<BR>
*'draw string 1.2 0.5 Analisi'<BR>
<BR>
<BR>
'enable print </FONT><A HREF="http://msg.gm"><U>msg.gm</U></A><FONT COLOR="#737373">'<BR>
'print'<BR>
'disable print'<BR>
*i=i+1<BR>
*endwhile<BR>
'quit'<BR>
<BR>
<BR>
<BR>
On Fri, 2011-04-01 at 02:28, Qian Hao wrote:
<BLOCKQUOTE TYPE=CITE>
Please also provide your .gs script file. I suggest you 'set dfile' to the corresponding file you are going to dispaly. Hope this may help.<BR>
<BR>
2011/3/30 Paolo Boi <</FONT><A HREF="mailto:pboi@arpa.sardegna.it"><U>pboi@arpa.sardegna.it</U></A><FONT COLOR="#737373">>
<BLOCKQUOTE>
Hi,<BR>
<BR>
I tried to plot two different set of data in GRIB format. Here are the two .ctl files, with different step in lat long:<BR>
<BR>
dset ^12140800.w73.grb<BR>
index ^12140800.w73.grb.idx<BR>
undef 9.999E+20<BR>
title 12140800.w73.grb<BR>
* produced by grib2ctl v0.9.12.5p16<BR>
dtype grib 255<BR>
options yrev<BR>
ydef 1024 linear 28.108000 0.0269494<BR>
xdef 1280 linear -4.311000 0.0269494<BR>
tdef 1 linear 08Z14dec2009 1mo<BR>
zdef 1 linear 1 1<BR>
vars 1<BR>
IMGD 0 127,122,22 ** Image data []<BR>
ENDVARS<BR>
<BR>
dset ^suolo.grb<BR>
index ^suolo.grb.idx<BR>
undef 9.999E+20<BR>
title suolo.grb<BR>
* produced by grib2ctl v0.9.12.5p16<BR>
dtype grib 255<BR>
options yrev<BR>
ydef 25 linear 20.000000 2.5<BR>
xdef 65 linear -90.000000 2.5<BR>
tdef 14 linear 12Z21mar2011 12hr<BR>
zdef 1 linear 1 1<BR>
vars 2<BR>
MSLsfc 0 151,1,0 ** Mean Sea Level (MSL) Pressure [Pa] Pa<BR>
Tprs 0 130,100,850 ** Temperature [K]<BR>
ENDVARS<BR>
<BR>
The two data set have a completely different step in lat long. <BR>
Here is the output error. I' using grads-1.8sl11.<BR>
<BR>
Data Request Error: Invalid grid coordinates<BR>
World coordinates convert to non-integer grid coordinates<BR>
Variable = imgd Dimension = 3<BR>
Error ocurred at column 1<BR>
DISPLAY error: Invalid expression<BR>
Expression = imgd/4<BR>
Cannot plot color bar: No shading information<BR>
Hardcopy output file is closed<BR>
No hardcopy metafile open<BR>
GX package terminated<BR>
<BR>
<BR>
<BR>
Many thanks<BR>
Paolo<BR>
<BR>
</BLOCKQUOTE>
</I></FONT><BR>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
</TR>
<TR>
</TR>
<TR>
</TR>
<TR>
</TR>
<TR>
<TD>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
</TR>
<TR>
</TR>
<TR>
</TR>
<TR>
</TR>
<TR>
<TD>
---------------------------------------------------<BR>
<B>Paolo Boi</B><BR>
<BR>
ARPAS- Agenzia Regionale per la Protezione dell'Ambiente della Sardegna<BR>
Dipartimento Specialistico Regionale Idrometeoclimatico<BR>
Servizio Idrometeoclimatico<BR>
<BR>
<FONT SIZE="1">Viale Porto Torres 119 - 07100 Sassari<BR>
Tel + 39 079 258600<BR>
Fax + 39 079 262681</FONT><BR>
<BR>
</TD>
</TR>
</TABLE>
</TD>
</TR>
</TABLE>
<BR>
<FONT COLOR="#737373"><I>_______________________________________________<BR>
gradsusr mailing list</FONT><BR>
<A HREF="mailto:gradsusr@gradsusr.org"><U>gradsusr@gradsusr.org</A><BR>
<A HREF="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</U></A><BR>
<FONT COLOR="#737373"><BR>
</BLOCKQUOTE>
<BR>
<BR>
<BR>
-- <BR>
此致<BR>
敬礼！<BR>
<BR>
<BR>
<HR>
<PRE>_______________________________________________
gradsusr mailing list</FONT>
<A HREF="mailto:gradsusr@gradsusr.org"><U>gradsusr@gradsusr.org</A>
<A HREF="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</U></I></A>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
</TR>
<TR>
</TR>
<TR>
<TD>
<BR>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
</TR>
<TR>
</TR>
<TR>
<TD>
<BR>
---------------------------------------------------<BR>
<BR>
<B>Paolo Boi</B><BR>
<BR>
<BR>
<BR>
ARPAS- Agenzia Regionale per la Protezione dell'Ambiente della Sardegna<BR>
<BR>
Dipartimento Specialistico Regionale Idrometeoclimatico<BR>
<BR>
Servizio Idrometeoclimatico<BR>
<BR>
<BR>
<BR>
<FONT SIZE="1">Viale Porto Torres 119 - 07100 Sassari<BR>
<BR>
Tel + 39 079 258600<BR>
<BR>
Fax + 39 079 262681</FONT><BR>
<BR>
<BR>
</TD>
</TR>
</TABLE>
</TD>
</TR>
</TABLE>
</PRE>
<FONT COLOR="#737373"><I>_______________________________________________<BR>
gradsusr mailing list</FONT><BR>
<A HREF="mailto:gradsusr@gradsusr.org"><U>gradsusr@gradsusr.org</A><BR>
<A HREF="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</U></A><BR>
<FONT COLOR="#737373"><BR>
</BLOCKQUOTE>
<BR>
<BR>
<BR>
-- <BR>
此致<BR>
敬礼！<BR>
<BR>
<HR>
<PRE>_______________________________________________
gradsusr mailing list
gradsusr@gradsusr.org</FONT>
<A HREF="http://gradsusr.org/mailman/listinfo/gradsusr"><U>http://gradsusr.org/mailman/listinfo/gradsusr</U></I></A></PRE>
</BLOCKQUOTE>
<PRE><TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
<TD>
<TABLE CELLSPACING="0" CELLPADDING="0" WIDTH="100%">
<TR>
<TD>
---------------------------------------------------<BR>
<B>Paolo Boi</B><BR>
<BR>
ARPAS- Agenzia Regionale per la Protezione dell'Ambiente della Sardegna<BR>
Dipartimento Specialistico Regionale Idrometeoclimatico<BR>
Servizio Idrometeoclimatico<BR>
<BR>
<FONT SIZE="1">Viale Porto Torres 119 - 07100 Sassari<BR>
Tel + 39 079 258600<BR>
Fax + 39 079 262681</FONT><BR>
<BR>
</TD>
</TR>
</TABLE>
</TD>
</TR>
</TABLE>
</PRE>
</BODY>
</HTML>