<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=utf-8"><meta name=Generator content="Microsoft Word 12 (filtered medium)"><style><!--
/* Font Definitions */
@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";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@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-GB link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Good morning, I am using grads to convert ERA operational analysis data from grads format to grib. The data is separated into monthly folders containing 6 hourly ‘snapshots’ comprising a .ctl file a .idx file and a file with no extension. The data can be opened and visualised using grads so as far as I can tell its integrity is not in question. I have written a script to sequentially process the grads files and convert them to grib (shown below) which produces a grib file output for the first time interval but halts with a segmentation fault (shown below script) when it appears to be dealing with the time range. I am new to programming and grads so please be forgiving.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Many thanks <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>d = 1<o:p></o:p></p><p class=MsoNormal>while (d <= 31)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>a = 0<o:p></o:p></p><p class=MsoNormal>while (a <= 18)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>if (d<10)<o:p></o:p></p><p class=MsoNormal>name = 'lisf96010'd''<o:p></o:p></p><p class=MsoNormal>e = 0d<o:p></o:p></p><p class=MsoNormal>else<o:p></o:p></p><p class=MsoNormal>name = 'lisf9601'd''<o:p></o:p></p><p class=MsoNormal>e = d<o:p></o:p></p><p class=MsoNormal>endif<o:p></o:p></p><p class=MsoNormal>if (a<10)<o:p></o:p></p><p class=MsoNormal>name = ''name'0'a''<o:p></o:p></p><p class=MsoNormal>f = 0a<o:p></o:p></p><p class=MsoNormal>else<o:p></o:p></p><p class=MsoNormal>name = 'name''a'<o:p></o:p></p><p class=MsoNormal>f = a<o:p></o:p></p><p class=MsoNormal>endif<o:p></o:p></p><p class=MsoNormal>'lats4d -v -i 'name'.ctl -o 9601'd''a'_sf -format grib -table .grads.lats.table -time 'f'Z'e'JAN96 'f'Z'e'JAN96'<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>a = a + 06<o:p></o:p></p><p class=MsoNormal>reset<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>endwhile<o:p></o:p></p><p class=MsoNormal>d = d + 01<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>endwhile<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!<o:p></o:p></p><p class=MsoNormal>-bash-3.2$ grads<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1<o:p></o:p></p><p class=MsoNormal>Copyright (c) 1988-2010 by Brian Doty and the<o:p></o:p></p><p class=MsoNormal>Institute for Global Environment and Society (IGES)<o:p></o:p></p><p class=MsoNormal>GrADS comes with ABSOLUTELY NO WARRANTY<o:p></o:p></p><p class=MsoNormal>See file COPYRIGHT for more information<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile<o:p></o:p></p><p class=MsoNormal>Issue 'q config' command for more detailed configuration information<o:p></o:p></p><p class=MsoNormal>Loading User Defined Extensions table </usr/local/WRF-NMM/tar_archives/grads-2.0.a9.oga.1/Classic/bin/gex/udxt> ... ok.<o:p></o:p></p><p class=MsoNormal>Landscape mode? ('n' for portrait):<o:p></o:p></p><p class=MsoNormal>GX Package Initialization: Size = 11 8.5<o:p></o:p></p><p class=MsoNormal>ga-> run do<o:p></o:p></p><p class=MsoNormal>lats4d: Version 2.1.5 of 28 June 2010<o:p></o:p></p><p class=MsoNormal>No hardcopy metafile open<o:p></o:p></p><p class=MsoNormal>All files closed; all defined objects released;<o:p></o:p></p><p class=MsoNormal>All GrADS attributes have been reinitialized<o:p></o:p></p><p class=MsoNormal>lats4d: Data file is<o:p></o:p></p><p class=MsoNormal>File 1 : lisf96010100<o:p></o:p></p><p class=MsoNormal> Descriptor: lisf96010100.ctl<o:p></o:p></p><p class=MsoNormal> Binary: lisf96010100<o:p></o:p></p><p class=MsoNormal> Type = Gridded<o:p></o:p></p><p class=MsoNormal> Xsize = 144 Ysize = 73 Zsize = 1 Tsize = 1 Esize = 1<o:p></o:p></p><p class=MsoNormal> Number of Variables = 21<o:p></o:p></p><p class=MsoNormal> no10u 0 165 ** surface 10 metre u wind component m s**-1<o:p></o:p></p><p class=MsoNormal> no10v 0 166 ** surface 10 metre v wind component m s**-1<o:p></o:p></p><p class=MsoNormal> no2d 0 168 ** surface 2 metre dewpoint temperature K<o:p></o:p></p><p class=MsoNormal> no2t 0 167 ** surface 2 metre temperature K<o:p></o:p></p><p class=MsoNormal> al 0 174 ** surface Albedo -<o:p></o:p></p><p class=MsoNormal> asq 0 233 ** surface Apparent Surface Humidity kg kg**-1<o:p></o:p></p><p class=MsoNormal> cdst 0 183 ** Clim deep soil tmp/soil tmp lev3(930804) K<o:p></o:p></p><p class=MsoNormal> cdsw 0 184 ** Clim deep soil wet/soil wet lev3(930803) m (of water)<o:p></o:p></p><p class=MsoNormal> dst 0 170 ** Deep soil tmp/soil temp lev2(from 930804) K<o:p></o:p></p><p class=MsoNormal> dsw 0 171 ** Deep soil wet/soil wet lev2(from 930803) m (of water)<o:p></o:p></p><p class=MsoNormal> lsrh 0 234 ** surface Log of surface roughness length for heat -<o:p></o:p></p><p class=MsoNormal> msl 0 151 ** surface Mean sea level pressure Pa<o:p></o:p></p><p class=MsoNormal> sd 0 141 ** surface Snow depth m (of water equivalent)<o:p></o:p></p><p class=MsoNormal> skt 0 235 ** surface Skin Temperature K<o:p></o:p></p><p class=MsoNormal> sr 0 173 ** surface Surface roughness m<o:p></o:p></p><p class=MsoNormal> src 0 198 ** surface Skin reservoir content m (of water)<o:p></o:p></p><p class=MsoNormal> ssw 0 140 ** Surf soil wet/soil wet lev1(from 930803) m (of water)<o:p></o:p></p><p class=MsoNormal> st 0 139 ** Surf.temp/soil temp lev 1 (from 930804) K<o:p></o:p></p><p class=MsoNormal> stl4 0 236 ** Soil temperature level 4 K<o:p></o:p></p><p class=MsoNormal> swl4 0 237 ** Soil wetness level 4 m<o:p></o:p></p><p class=MsoNormal> veg 0 199 ** surface Percentage of vegetation %<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>lats4d: Dimension environment file same as data file<o:p></o:p></p><p class=MsoNormal>lats4d: time range: 00Z01JAN1996 00Z01JAN1996 by 1, delta t: 6 hourly<o:p></o:p></p><p class=MsoNormal>lats4d: vertical levels: 1<o:p></o:p></p><p class=MsoNormal>lats4d: using PRESSURE for vertical coordinate<o:p></o:p></p><p class=MsoNormal>lats4d: latitudinal range: -90 90<o:p></o:p></p><p class=MsoNormal>lats4d: longitudinal range: 0 357.5<o:p></o:p></p><p class=MsoNormal>lats4d: surface variables: no10u no10v no2d no2t al asq cdst cdsw dst dsw lsrh msl sd skt sr src ssw st stl4 swl4 veg<o:p></o:p></p><p class=MsoNormal>lats4d: writing to grib_only file 960110_sf.grb on 00Z01JAN1996<o:p></o:p></p><p class=MsoNormal>lats4d: created grib_only file 960110_sf.grb<o:p></o:p></p><p class=MsoNormal>Most GrADS attributes have been reset<o:p></o:p></p><p class=MsoNormal>lats4d: Version 2.1.5 of 28 June 2010<o:p></o:p></p><p class=MsoNormal>No hardcopy metafile open<o:p></o:p></p><p class=MsoNormal>All files closed; all defined objects released;<o:p></o:p></p><p class=MsoNormal>All GrADS attributes have been reinitialized<o:p></o:p></p><p class=MsoNormal>lats4d: Data file is<o:p></o:p></p><p class=MsoNormal>File 1 : lisf96010106<o:p></o:p></p><p class=MsoNormal> Descriptor: lisf96010106.ctl<o:p></o:p></p><p class=MsoNormal> Binary: lisf96010106<o:p></o:p></p><p class=MsoNormal> Type = Gridded<o:p></o:p></p><p class=MsoNormal> Xsize = 144 Ysize = 73 Zsize = 1 Tsize = 1 Esize = 1<o:p></o:p></p><p class=MsoNormal> Number of Variables = 21<o:p></o:p></p><p class=MsoNormal> no10u 0 165 ** surface 10 metre u wind component m s**-1<o:p></o:p></p><p class=MsoNormal> no10v 0 166 ** surface 10 metre v wind component m s**-1<o:p></o:p></p><p class=MsoNormal> no2d 0 168 ** surface 2 metre dewpoint temperature K<o:p></o:p></p><p class=MsoNormal> no2t 0 167 ** surface 2 metre temperature K<o:p></o:p></p><p class=MsoNormal> al 0 174 ** surface Albedo -<o:p></o:p></p><p class=MsoNormal> asq 0 233 ** surface Apparent Surface Humidity kg kg**-1<o:p></o:p></p><p class=MsoNormal> cdst 0 183 ** Clim deep soil tmp/soil tmp lev3(930804) K<o:p></o:p></p><p class=MsoNormal> cdsw 0 184 ** Clim deep soil wet/soil wet lev3(930803) m (of water)<o:p></o:p></p><p class=MsoNormal> dst 0 170 ** Deep soil tmp/soil temp lev2(from 930804) K<o:p></o:p></p><p class=MsoNormal> dsw 0 171 ** Deep soil wet/soil wet lev2(from 930803) m (of water)<o:p></o:p></p><p class=MsoNormal> lsrh 0 234 ** surface Log of surface roughness length for heat -<o:p></o:p></p><p class=MsoNormal> msl 0 151 ** surface Mean sea level pressure Pa<o:p></o:p></p><p class=MsoNormal> sd 0 141 ** surface Snow depth m (of water equivalent)<o:p></o:p></p><p class=MsoNormal> skt 0 235 ** surface Skin Temperature K<o:p></o:p></p><p class=MsoNormal> sr 0 173 ** surface Surface roughness m<o:p></o:p></p><p class=MsoNormal> src 0 198 ** surface Skin reservoir content m (of water)<o:p></o:p></p><p class=MsoNormal> ssw 0 140 ** Surf soil wet/soil wet lev1(from 930803) m (of water)<o:p></o:p></p><p class=MsoNormal> st 0 139 ** Surf.temp/soil temp lev 1 (from 930804) K<o:p></o:p></p><p class=MsoNormal> stl4 0 236 ** Soil temperature level 4 K<o:p></o:p></p><p class=MsoNormal> swl4 0 237 ** Soil wetness level 4 m<o:p></o:p></p><p class=MsoNormal> veg 0 199 ** surface Percentage of vegetation %<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>lats4d: Dimension environment file same as data file<o:p></o:p></p><p class=MsoNormal>Segmentation fault<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='color:black'>James Hughes<o:p></o:p></span></p><p class=MsoNormal><span style='color:black'>Centre for Renewable Energy Systems Technology</span><o:p></o:p></p><p class=MsoNormal><span style='color:black'>Loughborough University</span><o:p></o:p></p><p class=MsoNormal><span style='color:black'>Tel: +44 1509 635312<o:p></o:p></span></p><p class=MsoNormal><span style='color:black'>Email: j.hughes@lboro.ac.uk<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>