<html><head></head><body>using define keyword or not does not affect the the result<div>the problem is that you cannot use the define statement (define x=... or x=... are the same statement because in the latter 'define' is implicit) starting from a station file as default file, you should have a gridded file.</div><div>The solution should be easy: set dfile gridded_file (1, 2, ...)</div><div>eg.: gridded file unit 2, stnfile unit 1:</div><div>set dfile 2</div><div>define x=oacres(t,t.1)</div><div><br></div><div>hoping it helps</div><div>Davide<br><br><div style="margin: 0px; padding: 0px; box-sizing: content-box; color: rgb(34, 34, 34);"><i style="box-sizing: content-box;"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;">Dott. Davide Sacchetti</span></font></i></div><div style="margin: 0px; padding: 0px; box-sizing: content-box; color: rgb(34, 34, 34);"><i style="box-sizing: content-box;"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;"></span></font></i></div><i style="box-sizing: content-box; color: rgb(34, 34, 34);"><font size="2" style="box-sizing: content-box;"><span id="mcntmcntext-gen4583" style="box-sizing: content-box; font-family: Calibri, sans-serif;">Dipartimento Stato dell’Ambiente e Tutela dai Rischi Naturali<br style="box-sizing: content-box;">U.O. Clima Meteo Idro<br style="box-sizing: content-box;">Viale Brigate Partigiane 2, 16129 Genova (I)<br style="box-sizing: content-box;">tel: +39 0106437535<br style="box-sizing: content-box;">mail: <a href="mailto:mauro.damonte@arpal.liguria.it" class="mailto" title="Send email to mauro.damonte@arpal.liguria.it" target="_blank" style="box-sizing: content-box;">davide.sacchetti@arpal.liguria.it</a> web: <a href="http://www.meteoliguria.it/" title="http://www.meteoliguria.it" target="_blank" style="box-sizing: content-box;">https://www.arpal.liguria.it/homepage/meteo.html</a></span></font></i><br><br><br><div><strong>
Da:
</strong>
Isimar de Azevedo Santos <isimar.uenf@gmail.com>
<br>
<strong>
A:
</strong>
<pavel@cyclone.tau.ac.il>, GrADS Users Forum <gradsusr@gradsusr.org>
<br>
<strong>
Inviato:
</strong>
13/09/2021 22:31
<br>
<strong>
Oggetto:
</strong>
Re: [gradsusr] How to store a map, based on station data, in GRADS, in a NetCDF-file?
<br><br><blockquote class="mori" style="margin:0 0 0 .8ex;border-left:1px solid #CCC;padding-left:1ex;"><div>Hi, Pavel<div>Why not use simply: "x= ... ", without "define"</div><div>Isimar from Brazil.</div></div><br><div class="mcntgmail_quote"><div class="mcntgmail_attr">Em seg., 13 de set. de 2021 às 12:43, Pavel Kishcha <<a href="mailto:pavel@cyclone.tau.ac.il" title="Invia e-mail a pavel@cyclone.tau.ac.il" class="mailto" target="_blank">pavel@cyclone.tau.ac.il</a>> escreveu:<br></div><blockquote class="mcntgmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="mcntgmail_quote"><div><p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">Hello to all users of grads,</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:12pt"><br></span></p><p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:12pt">A Cressman objective analysis </span><span style="font-size:12pt">was
conducted on the station data to obtain a map </span>representing<span style="font-size:12pt"> the station data.</span><br></p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:12pt">'set gxout shaded'</span><br></p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">'d oacres (gexpr, sexpr)'</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:12pt">As a result, the map was successfully created.</span></p><p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"><span style="font-size:12pt">To store the obtained map in a NetCDF-file, we tried to use
the following GRADS commands:</span><br></p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">'set sdfwrite <a title="map.nc" href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fmap.nc&e=7312d17f&h=82f565dc&f=y&p=n" target="_blank">map.nc</a>'</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">‘define x = oacres (gexpr, sexpr)'</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">'sdfwrite x'</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"> </p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">Unfortunately, we got the following error info:</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">DEFINE error: Define not yet valid for station data</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"> Default
file is a station data file</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"> </p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">How to solve this problem and to store the obtained map (based
on station data) in GRADS, in a NetCDF-file? Any ideas?</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black"> </p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">Thank you,</p>
<p class="mcntMsoNormal" style="margin:0cm 0cm 10pt;line-height:115%;font-size:12pt;font-family:"Times New Roman",serif;color:black">Pavel Kishcha</p></div><div> </div><blockquote class="mcntgmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="text-align:left;margin-left:40px">_______________________________________________<br></div></div>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank" title="Invia e-mail a gradsusr@gradsusr.org" class="mailto">gradsusr@gradsusr.org</a><br>
<a title="gradsusr.org" href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=7312d17f&h=5ee7e1e3&f=y&p=n" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div></div>
_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank" title="Invia e-mail a gradsusr@gradsusr.org" class="mailto">gradsusr@gradsusr.org</a><br>
<a title="gradsusr.org" href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=7312d17f&h=5ee7e1e3&f=y&p=n" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div>
<br><br>_______________________________________________
<br>gradsusr mailing list
<br>gradsusr@gradsusr.org
<br> https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgradsusr.org%2Fmailman%2Flistinfo%2Fgradsusr&e=7312d17f&h=5ee7e1e3&f=y&p=n
<br></blockquote></div></div></body></html>