[gradsusr] Station Data issues

[Jennifer Adams]

On Sep 29, 2013, at 7:40 PM, Mike Manning wrote:

> For anyone else chasing problems like this it was a problem in Fortran when writing out the data – this has been rectified now.
Very good, glad to hear it is working now. Thanks for closing the book on this question. 
--Jennifer
>  
> Cheers, Mike
>  
> From: gradsusr-bounces at gradsusr.org [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Mike Manning
> Sent: Friday, 27 September 2013 4:31 PM
> To: gradsusr at gradsusr.org
> Subject: [gradsusr] Station Data issues
>  
> 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
>  
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
111 Research Hall, Mail Stop 2B3
George Mason University
4400 University Drive
Fairfax, VA 22030 





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130929/e063fe68/attachment-0001.html