<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 14 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><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;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
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;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="color:#1F497D">Before you “quit” you have to display the variable, like “d t2m” or whatever it’s called in the ERA5 file.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D">Another tip -- instead of using world coordinates for bounds (“set lat” and “set lon”), use “set x” and “set y” and the exact grid dimensions (I think they’re 1440 for nx and 721 for ny).<o:p></o:p></span></p>
<p class="MsoNormal"><a name="_MailEndCompose"><span style="color:#1F497D"><o:p> </o:p></span></a></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> gradsusr [mailto:gradsusr-bounces@gradsusr.org]
<b>On Behalf Of </b>Stephan Steyn<br>
<b>Sent:</b> Tuesday, May 7, 2019 7:58 AM<br>
<b>To:</b> gradsusr@gradsusr.org<br>
<b>Subject:</b> [gradsusr] converting GRIB to .dat for FORTRAN<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="EN-ZA">Dear GrADS users<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">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></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">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></span></p>
<p class="MsoNormal"><span lang="EN-ZA">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></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">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></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">I am pretty much stuck now... Please help!<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">Stephan<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">*** convert.gs ***<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">*---------------------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">function main()<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">* Please edit the following line to match your computer environment:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">mydir='/home/stephan/ERA5'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'open ERA5_sfcvar12_197901.ctl'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'set lon 0.00 359.75'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'set lat -90.00 90.00'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'set time jan1979'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'set gxout fwrite'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'set fwrite 'mydir'/ERA5_sfcvar12_197901.dat'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">'quit'<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">*----------------------------------------------------------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">Below the content of "ERA_sfcvar12_197901.ctl": -----------------------<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">dset ^ERA5_sfcvar12_197901.grib<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">index ^ERA5_sfcvar12_197901.grib.idx<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">undef 9.999E+20<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">title ERA5_sfcvar12_197901.grib<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">* produced by grib2ctl v0.9.15<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">* command line options: ERA5_sfcvar12_197901.grib<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">dtype grib 255<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">options yrev<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">ydef 721 linear -90.000000 0.25<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">xdef 1440 linear 0.000000 0.250000<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">tdef 31 linear 12Z01jan1979 1dy<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">zdef 1 linear 1 1<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">vars 4<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">no10Usfc 0 165,1,0 ** surface 10 metre U wind component [m s**-1]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">no10Vsfc 0 166,1,0 ** surface 10 metre V wind component [m s**-1]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">no2Dsfc 0 168,1,0 ** surface 2 metre dewpoint temperature [K]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">no2Tsfc 0 167,1,0 ** surface 2 metre temperature [K]<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA">ENDVARS<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-ZA"><o:p> </o:p></span></p>
<div class="MsoNormal" align="center" style="text-align:center"><span lang="EN-ZA" style="font-size:12.0pt;font-family:"Times New Roman","serif"">
<hr size="3" width="100%" align="center">
</span></div>
<p class="MsoNormal"><span lang="EN-ZA" style="font-size:12.0pt;font-family:"Times New Roman","serif"">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.<o:p></o:p></span></p>
</div>
</body>
</html>