No subject


Wed Apr 14 09:37:02 EDT 2010


getting the whole globe.


> > One very simple solution is to add a GRIB-2 to GRIB-1 filter along the
> > way, to convert the GRIB-2 to GRIB-1 format. Grib-2 is a good format for
> > data transmission but is far less less optimal for day to day use. For
> > example, the performance of your meteograms scripts will go up by a lot
> > with grib-1. Of course, you have to factor in the cost of conversion.
>
> That sounds rather complex. I did some modifications and I wait for the
> feedback from my provider. If there still is a problem, I consider to
> take the meteograms out.
> I was thinking, if there is another option.
>
> grads can write data to disk, right?
> Would it be possible to change the procedure like that
>
> while timestep <= 61
>  Dowload timestep data
>  start grads
>  render charts
>  write the relevant data for the cities for this timestep to disk
> done
> start another grads script
> read the data for the meteograms and render those
>

If I recall well your meteograms are for a relatively confined region. Use
lats4d to pull the data you want

% lats4d -i filename -o subset -lon $lon1 $lon2 -lat $lat1 $lat2 -format
grads_grib \
              -func 'g2(@)' -v

then run the meteogram script on the subset; ($lon1,$lat1,$lon2,$lat2) is
your bounding box.
(You can also run lats4d from within your chart rendering script; in this
case omit "-i filename")


> What would be an appropiate format to write the data to disk so another
> script can easily read them?
>

If using lats4d I'd use grib-1 (-format grads_grib).

    Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu

------=_Part_22732_8864895.1221702277219
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: 7bit
Content-Disposition: inline

<div dir="ltr"><div class="gmail_quote">On Wed, Sep 17, 2008 at 5:59 PM, Stefan Gofferje <span dir="ltr">&lt;<a href="mailto:gradsusers at gofferje.homelinux.org">gradsusers at gofferje.homelinux.org</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
<br>
Arlindo da Silva schrieb:<br>
<div class="Ih2E3d">&gt; In your case, you could use lats4d with the g2 extension to create a<br>
&gt; scandinaian subset for your meteograms.<br>
<br>
</div>I already do some limiting in the download sript:<br></blockquote><div><br></div><div>From what I can tell you are selecting levels and variables, but still getting the whole globe.&nbsp;</div><div>&nbsp;</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="Ih2E3d">
&gt; One very simple solution is to add a GRIB-2 to GRIB-1 filter along the<br>
&gt; way, to convert the GRIB-2 to GRIB-1 format. Grib-2 is a good format for<br>
&gt; data transmission but is far less less optimal for day to day use. For<br>
&gt; example, the performance of your meteograms scripts will go up by a lot<br>
&gt; with grib-1. Of course, you have to factor in the cost of conversion.<br>
<br>
</div>That sounds rather complex. I did some modifications and I wait for the<br>
feedback from my provider. If there still is a problem, I consider to<br>
take the meteograms out.<br>
I was thinking, if there is another option.<br>
<br>
grads can write data to disk, right?<br>
Would it be possible to change the procedure like that<br>
<br>
while timestep &lt;= 61<br>
 &nbsp;Dowload timestep data<br>
 &nbsp;start grads<br>
 &nbsp;render charts<br>
 &nbsp;write the relevant data for the cities for this timestep to disk<br>
done<br>
start another grads script<br>
read the data for the meteograms and render those<br></blockquote><div><br></div><div>If I recall well your meteograms are for a relatively confined region. Use lats4d to pull the data you want</div><div><br></div><div>% lats4d -i filename -o subset -lon $lon1 $lon2 -lat $lat1 $lat2 -format grads_grib \</div>
<div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;-func &#39;g2(@)&#39; -v</div><div>&nbsp;</div><div>then run the meteogram script on the subset; ($lon1,$lat1,$lon2,$lat2) is your bounding box.</div><div>(You can also run lats4d from within your chart rendering script; in this case omit &quot;-i filename&quot;)</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><br>
What would be an appropiate format to write the data to disk so another<br>
script can easily read them?<br>
<div><div></div><div class="Wj3C7c"></div></div></blockquote><div><br></div><div>If using lats4d I&#39;d use grib-1 (-format grads_grib).</div><div><br></div><div>&nbsp;&nbsp; &nbsp;Arlindo</div><div><br></div></div>-- <br>Arlindo da Silva<br>
<a href="mailto:dasilva at alum.mit.edu">dasilva at alum.mit.edu</a><br>
</div>

------=_Part_22732_8864895.1221702277219--



More information about the gradsusr mailing list