<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-ZA" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Dear GrADS users<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">As described on https://ds.data.jma.go.jp/tcc/tcc/library/MRCS_SV12/sstdata/readme_grib_e.htm I am attempting to use a GrADS script to convert a GRIB file (in the form of ERA5 gridded surface data) to a .dat file that I can read with FORTRAN
(4-byte real floating point format). I already have fortran programs that can ingest this format.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I managed to use "grib2ctl.pl" to create the .ctl file (copied below) and "gribmap" to create a .grib.idx file and successfully display the data in GrADS.<o:p></o:p></p>
<p class="MsoNormal">Upon executing "run convert.gs" (copied below) within GrADS I don't receive any error messages, only to find that the .dat file was never created. I also tried to make subsets of the data with the same non-result.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am using GrADS Version 2.0.2 run on 64-bit Fedora in a Oracle Virtualbox on Lenovo Thinkpad (Core i7).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I am pretty much stuck now... Please help!<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Stephan<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">*** convert.gs ***<o:p></o:p></p>
<p class="MsoNormal">*---------------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">function main()<o:p></o:p></p>
<p class="MsoNormal">* Please edit the following line to match your computer environment:<o:p></o:p></p>
<p class="MsoNormal">mydir='/home/stephan/ERA5'<o:p></o:p></p>
<p class="MsoNormal">'open ERA5_sfcvar12_197901.ctl'<o:p></o:p></p>
<p class="MsoNormal">'set lon 0.00 359.75'<o:p></o:p></p>
<p class="MsoNormal">'set lat -90.00 90.00'<o:p></o:p></p>
<p class="MsoNormal">'set time jan1979'<o:p></o:p></p>
<p class="MsoNormal">'set gxout fwrite'<o:p></o:p></p>
<p class="MsoNormal">'set fwrite 'mydir'/ERA5_sfcvar12_197901.dat'<o:p></o:p></p>
<p class="MsoNormal">'quit'<o:p></o:p></p>
<p class="MsoNormal">*----------------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Below the content of "ERA_sfcvar12_197901.ctl": -----------------------<o:p></o:p></p>
<p class="MsoNormal">dset ^ERA5_sfcvar12_197901.grib<o:p></o:p></p>
<p class="MsoNormal">index ^ERA5_sfcvar12_197901.grib.idx<o:p></o:p></p>
<p class="MsoNormal">undef 9.999E+20<o:p></o:p></p>
<p class="MsoNormal">title ERA5_sfcvar12_197901.grib<o:p></o:p></p>
<p class="MsoNormal">* produced by grib2ctl v0.9.15<o:p></o:p></p>
<p class="MsoNormal">* command line options: ERA5_sfcvar12_197901.grib<o:p></o:p></p>
<p class="MsoNormal">dtype grib 255<o:p></o:p></p>
<p class="MsoNormal">options yrev<o:p></o:p></p>
<p class="MsoNormal">ydef 721 linear -90.000000 0.25<o:p></o:p></p>
<p class="MsoNormal">xdef 1440 linear 0.000000 0.250000<o:p></o:p></p>
<p class="MsoNormal">tdef 31 linear 12Z01jan1979 1dy<o:p></o:p></p>
<p class="MsoNormal">zdef 1 linear 1 1<o:p></o:p></p>
<p class="MsoNormal">vars 4<o:p></o:p></p>
<p class="MsoNormal">no10Usfc 0 165,1,0 ** surface 10 metre U wind component [m s**-1]<o:p></o:p></p>
<p class="MsoNormal">no10Vsfc 0 166,1,0 ** surface 10 metre V wind component [m s**-1]<o:p></o:p></p>
<p class="MsoNormal">no2Dsfc 0 168,1,0 ** surface 2 metre dewpoint temperature [K]<o:p></o:p></p>
<p class="MsoNormal">no2Tsfc 0 167,1,0 ** surface 2 metre temperature [K]<o:p></o:p></p>
<p class="MsoNormal">ENDVARS<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<HR>University of the Free State: This message and its contents are subject to a disclaimer. <BR>
Please refer to http://www.ufs.ac.za/disclaimer for full details.<BR>
</body>
</html>