<div dir="ltr">Hi all,<div><br></div><div>I am a new grads user and I have never worked with binary files before. I have a large binary file I created in grads (code shown below). I want to split up the binary file by latitude and longitude degree into smaller binary files. While this is possible in GrADS, it is too slow for my needs, so I am trying to do it in Fortran 90. </div><div><br></div><div>However, I am having problems reading the GrADS binary file in Fortran. It gives me the missing value for each of my variables at every time step and every location.</div><div><br></div><div>The binary file contains NLDAS 2 forcing data, with 464 x steps, 224 y steps, and 13688 daily time steps for 4 variables (only 1 Z level). It is a large file at 22GB, and I can successfully display the variables in GrADS from this binary file. This binary file was created from data in Grib files.</div><div><br></div><div>I have tried many various ways of reading in the binary data, using x and y loops, or time loops, or looping through the 4 variables, but nothing has been successful. </div><div><br></div><div>I am working on a Mac (El Capitan).</div><div><br></div><div>Here is the GrADS script I used to create the binary file:</div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="color:rgb(206,121,36)">function</span><span style="color:rgb(0,0,0)"> </span><span style="">main</span><span style="color:rgb(0,0,0)">(</span><span style="">args</span><span style="color:rgb(0,0,0)">) </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">counthr</span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(206,121,36)">subwrd</span><span style="color:rgb(0,0,0)">(</span><span style="">args</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(195,55,32)">1</span><span style="color:rgb(0,0,0)">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'reinit'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'open /Volumes/WebData/NLDAS2/ProcessingFiles/24hours.ctl'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'set undef 9.999e20'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'set gxout fwrite'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'set fwrite -ap /Volumes/WebData/NLDAS2/DailyNLDAS/daily.bin'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span style=""></span><br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'set t 1'</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(52,187,199)">t1</span><span style="color:rgb(0,0,0)">=</span><span style="">1+4</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">t1p</span><span style="color:rgb(0,0,0)">=</span><span style="color:rgb(195,55,32)">1+5</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="color:rgb(206,121,36)">while</span><span style="color:rgb(0,0,0)">(</span><span style="">t1</span><span style="color:rgb(0,0,0)"><=</span><span style="">counthr</span><span style="color:rgb(195,55,32)">+5-23</span><span style="color:rgb(0,0,0)">) </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">t2</span><span style="color:rgb(0,0,0)">=</span><span style="">t1</span><span style="color:rgb(195,55,32)">+23</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">t2p</span><span style="color:rgb(0,0,0)">=</span><span style="">t1p</span><span style="color:rgb(195,55,32)">+23</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'d (max(TMP2m,t='</span><span style="color:rgb(52,187,199)">t1</span><span style="">',t='</span><span style="color:rgb(52,187,199)">t2</span><span style="">')-273.15)'</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'d (min(TMP2m,t='</span><span style="color:rgb(52,187,199)">t1</span><span style="">',t='</span><span style="color:rgb(52,187,199)">t2</span><span style="">')-273.15)'</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'d sum(APCPsfc,t='</span><span style="color:rgb(52,187,199)">t1p</span><span style="">',t='</span><span style="color:rgb(52,187,199)">t2p</span><span style="">')'</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'d (ave(DSWRFsfc,t='</span><span style="color:rgb(52,187,199)">t1</span><span style="">',t='</span><span style="color:rgb(52,187,199)">t2</span><span style="">')*24*3600/1000000)'</span><span style="color:rgb(0,0,0)"> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">t1</span><span style="color:rgb(0,0,0)">=</span><span style="">t1</span><span style="color:rgb(195,55,32)">+24</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,187,199)"><span style="">t1p</span><span style="color:rgb(0,0,0)">=</span><span style="">t1p</span><span style="color:rgb(195,55,32)">+24</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(206,121,36)"><span style="">endwhile</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="">'quit'</span></p></div><div><br></div><div>Here is the Fortran script to read the binary file</div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="color:rgb(213,59,211)">Program</span><span style=""> Readbin</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(52,189,38)"><span style="">implicit</span><span style="color:rgb(0,0,0)"> </span><span style="">none</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span style="color:rgb(52,189,38)">integer</span> irec,i,j,t,days<br><span style=""></span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="color:rgb(52,189,38)">real</span> tmin(<span style="color:rgb(195,55,32)">464</span>,<span style="color:rgb(195,55,32)">224</span>),tmax(<span style="color:rgb(195,55,32)">464</span>,<span style="color:rgb(195,55,32)">224</span>),prcp(<span style="color:rgb(195,55,32)">464</span>,<span style="color:rgb(195,55,32)">224</span>),srad(<span style="color:rgb(195,55,32)">464</span>,<span style="color:rgb(195,55,32)">224</span>)<br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(83,48,225)"><span style="color:rgb(206,121,36)">OPEN</span><span style="color:rgb(0,0,0)">(</span><span style="color:rgb(195,55,32)">17</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">file=</span><span style="color:rgb(195,55,32)">'/</span><span style="color:rgb(195,55,32)">Volumes/WebData/NLDAS2/DailyNLDAS</span><span style="color:rgb(195,55,32)">/daily.bin'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">&</span><br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(206,121,36)">&STATUS=</span><span style="">'Old'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">FORM=</span><span style="">'UNFORMATTED'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">ACCESS=</span><span style="">'DIRECT'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">RECL=</span><span style="">464</span><span style="color:rgb(206,121,36)">*</span><span style="">224</span><span style="color:rgb(206,121,36)">*</span><span style="">4</span><span style="color:rgb(0,0,0)">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;min-height:13px"><span style=""></span><br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(206,121,36)">OPEN</span><span style="color:rgb(0,0,0)">(</span><span style="">18</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">file=</span><span style="">'</span>/Volumes/WebData/NLDAS2/DailyNLDAS/<span style="">testdaily.bin'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">&</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(206,121,36)">&STATUS=</span><span style="">'UNKNOWN'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">FORM=</span><span style="">'UNFORMATTED'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">ACCESS=</span><span style="">'DIRECT'</span><span style="color:rgb(0,0,0)">,</span><span style="color:rgb(206,121,36)">RECL=</span><span style="">464</span><span style="color:rgb(206,121,36)">*</span><span style="">224</span><span style="color:rgb(206,121,36)">*</span><span style="">4</span><span style="color:rgb(0,0,0)">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">days</span><span style="color:rgb(206,121,36)">=1</span><span style="color:rgb(195,55,32)">0</span><span style=""> </span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">irec</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span></p><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="color:rgb(206,121,36)">do </span><span style="">t</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,days</span><br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"> <span style="color:rgb(206,121,36)">read</span>(<span style="color:rgb(195,55,32)">17</span>,<span style="color:rgb(206,121,36)">REC=</span>irec) ((tmax(i,j),i<span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span>,<span style="color:rgb(195,55,32)">8</span>),j<span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span>,<span style="color:rgb(195,55,32)">8</span>)<br></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">read</span><span style="">(</span><span style="color:rgb(195,55,32)">17</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">1</span><span style="">) ((tmin(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">read</span><span style="">(</span><span style="color:rgb(195,55,32)">17</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">2</span><span style="">) ((prcp(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">read</span><span style="">(</span><span style="color:rgb(195,55,32)">17</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">3</span><span style="">) ((srad(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">write</span><span style="">(</span><span style="color:rgb(195,55,32)">18</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec) ((tmin(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">write</span><span style="">(</span><span style="color:rgb(195,55,32)">18</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">1</span><span style="">) ((tmin(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">write</span><span style="">(</span><span style="color:rgb(195,55,32)">18</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">2</span><span style="">) ((prcp(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> </span><span style="color:rgb(206,121,36)">write</span><span style="">(</span><span style="color:rgb(195,55,32)">18</span><span style="">,</span><span style="color:rgb(206,121,36)">REC=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">3</span><span style="">) ((srad(i,j),i</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">),j</span><span style="color:rgb(206,121,36)">=</span><span style="color:rgb(195,55,32)">1</span><span style="">,</span><span style="color:rgb(195,55,32)">8</span><span style="">)</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style=""> irec</span><span style="color:rgb(206,121,36)">=</span><span style="">irec</span><span style="color:rgb(206,121,36)">+</span><span style="color:rgb(195,55,32)">4</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(206,121,36)"><span style="">enddo</span></p></div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(213,59,211)"><span style="">end program</span></p></div><div><br></div><div>The corresponding control file is:</div><div><p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">set ^testdaily</span><span style="color:rgb(206,121,36)">.</span><span style="">bin</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(0,0,0)">undef </span><span style="">9.999E+20</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">title Daily NLDAS2 Data </span><span style="color:rgb(206,121,36)">for</span><span style=""> 00Z through 23Z</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(195,55,32)"><span style="color:rgb(0,0,0)">xdef </span><span style="">8</span><span style="color:rgb(0,0,0)"> linear </span><span style="">-124.9375</span><span style="color:rgb(0,0,0)"> </span><span style="">0.125</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">ydef </span><span style="color:rgb(195,55,32)">8</span><span style=""> linear </span><span style="color:rgb(195,55,32)">25.0625</span><span style=""> </span><span style="color:rgb(195,55,32)">0.125</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">tdef </span><span style="color:rgb(195,55,32)">10</span><span style=""> linear 00Z02Jan1979 1dy</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">zdef </span><span style="color:rgb(195,55,32)">1</span><span style=""> linear </span><span style="color:rgb(195,55,32)">1</span><span style=""> </span><span style="color:rgb(195,55,32)">1</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">vars </span><span style="color:rgb(195,55,32)">4</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">TMAX </span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(195,55,32)">12</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">105</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">2</span><span style=""> </span><span style="color:rgb(206,121,36)">**</span><span style=""> maximum temperature at 2m above surface [C]</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">TMIN </span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(195,55,32)">11</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">105</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">2</span><span style=""> </span><span style="color:rgb(206,121,36)">**</span><span style=""> minimum temperature at 2m above surface [C]</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">PRCP </span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(195,55,32)">61</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">1</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(206,121,36)">**</span><span style=""> total precipitation backward accumulated [kg</span><span style="color:rgb(206,121,36)">/</span><span style="">m^</span><span style="color:rgb(195,55,32)">2</span><span style=""> </span><span style="color:rgb(206,121,36)">or</span><span style=""> mm]</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">SRAD </span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(195,55,32)">204</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">1</span><span style="color:rgb(206,121,36)">,</span><span style="color:rgb(195,55,32)">0</span><span style=""> </span><span style="color:rgb(206,121,36)">**</span><span style=""> total surface downward shortwave radiation flux [mJ</span><span style="color:rgb(206,121,36)">/</span><span style="">m^</span><span style="color:rgb(195,55,32)">2</span><span style="">]</span></p>
<p style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span style="">ENDVARS</span></p></div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Lydia</div><div><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr" style="font-size:12.8000001907349px"><font color="#999999">Lydia D. Rill</font></div><div style="font-size:12.8000001907349px"><font color="#999999">M.S., Michigan State University 2016 </font></div><div dir="ltr" style="font-size:12.8000001907349px"><font color="#999999">B.S., Valparaiso University 2014</font><br style="color:rgb(136,136,136)"><span title="Call with Google Voice" style="color:rgb(136,136,136)"><font color="#999999"><a value="+12244065130" style="color:rgb(17,85,204)">(224) 406-5130</a></font></span><br style="color:rgb(136,136,136)"><font color="#000000"><a href="mailto:Lydia.D.Rill@gmail.com" style="color:rgb(17,85,204)" target="_blank">Lydia.D.Rill@gmail.com</a></font></div></div></div>
</div></div>