Fortran - Apertura de *.NAV - La rutina no corre

 
Vista:

Apertura de *.NAV - La rutina no corre

Publicado por Roiter Egoavil (1 intervención) el 17/11/2006 02:41:09
Hola amigos!, me gustaria que ayudaran con esta pequeña rutina que sirver para poder leer un archivo con extensión NAV que se baja del Goes (Imagenes satelitales), lo que sucede es que al momento de compilar me sale errores en las lineas;

character file_name*80

type *,' width =',xsize, ' height =',ysize

el archivo que debe de ser leido se llama Goes.nav
la rutina es la siguiente:

subroutine getnav(lat,lon)
!c************************************************************************
!c* This routine reads goes.nav files with a maximum size of *
!c* 2000 pixels wide and 2000 lines/records and stores *
!c* latitude/longitude data into lat lon *
!c* the format of the goes.nav data is as follows *
!c* *
!c* i*2 i*2 *
!c* | | *
!c* V V *
!c* record 1 x size ysize *
!c* record 2 lat *100 lon*100 *
!c* record 3 lat *100 lon*100 *
!c* record 4 etc etc *
!c* *
!c************************************************************************ *
character file_name*80
real lat(2000,2000),lon(2000,2000)
!c
!c
integer*2 tempstor(2)
integer*4 buffout,record_count,xsize,ysize,yyyy,xxxx
equivalence (buffout,tempstor(1))
!C
!C
!C PLEASE NOTE: recl=1 in the open statment below means
!C that the record length is 32 bits ( 4 bytes)
!C not 8 bits (1 byte)
open(23,file='goes.nav',access='direct',recl=1,status='old')
READ(23,REC=1)buffout
!C
!C
xsize=tempstor(1)
!C
!C
ysize=tempstor(2)
type *,' width =',xsize, ' height =',ysize
!C
!C

record_count=1
do yyyy=1,ysize
do xxxx=1,xsize
record_count = record_count + 1
read(23,rec=record_count)buffout
lat(xxxx,yyyy)=tempstor(1) * .01
lon(xxxx,yyyy)=tempstor(2) * .01
!c ...............................................
!c NOTE: LATITUDE VALUES OFF THE EARTH ARE 99.00
!C AND LONGITUDE VALUES ARE -256.35999 AT THIS
!C POINT IN THE SUBROUTINE
!C ...............................................
end do
end do
RETURN
END
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder