[gradsusr] Station Data issues

[Mike Manning]

Hi all,

 

I'm writing station data to a binary file using Fortran for an 8 day period.
When I create the binary file and run the stnmap util it creates the map
file fine with all 8 time series:

 

Starting scan of station data binary file.

Binary data file open: ocf.bin

 

Processing time step 1

     Time = 1 has stn count = 762

 

Processing time step 2

     Time = 2 has stn count = 762

 

Processing time step 3

     Time = 3 has stn count = 762

 

Processing time step 4

     Time = 4 has stn count = 762

 

Processing time step 5

     Time = 5 has stn count = 762

 

Processing time step 6

     Time = 6 has stn count = 762

 

Processing time step 7

     Time = 7 has stn count = 762

 

Processing time step 8

     Time = 8 has stn count = 762

  Max reports per time:  762 reports at t = 1

  Max data elements in largest report: 3

 

Version 2 Station map file created: /grads2/ocf/ocf.map

 

 

 

When I run a gradscript to display the data, it works perfect for the first
timestep (t=1) however when I try and set it to t=2, I get the following
error:

 

Low Level I/O Error:  Read error on data file

  Data file name = ocf.bin

  Error reading 28 bytes at location 139536556

Data Request Error:  Variable is 'amax'

  Error ocurred at column 1

DISPLAY error:  Invalid expression

  Expression = amax(t=2)

You must display a variable before using basemap

No hardcopy metafile open

GX package terminated

 

 

For my fortran script to write the binary data, I'm using this write clause:
open(21, file=outputfile, form='unformatted',
status='unknown',access='stream')

 

Anyone got any ideas on why it'd be struggling reading anything more than
the first time series?

 

Cheers, Mike

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130927/cf9cd471/attachment.html