[A fortran program: convert the LP ODF data to ASCII records]

This fortran code will convert the LP ODF data on PDS node to ASCII records so that it would be much easier to convert it to another data format.

 

This program is designed to run on a PC with window OS and Compaq Visual Fortran (6.6 or heigher virsion (a lower vision donnot have data type integer*8), you don't have to change it too much to transfer it to another platform or compile it with other compliers.

 

If you have any questions, please donnot hesitate to contact me at sunyuwoody@126.com

 

---------------------- the following is the original code -----------------------------

 

 

C Program ReadODF

Program ReadODF

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

C C

C This program will transform the LP ODF format (descriped by the interface C

C file trk-2-18) to ASCII records C

C C

C Note* This Program need to be run in a 6.6 or higher virsion.Since there C

C is no data type integer*8 in Compaq Visual Fortran C

C C

C SUN YU.( 2010 ) Shandong University of Science and Technology, Qingdao 266510, Shanddong Province, China C

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

 

 

C ------------ Begin of Data Type Definition------------------

Implicit None

Integer*4 i,R,Iramp

Integer*4 Idat_four(10),IRampG_four(10),IRampGH_four(10)

Integer*4 IOrbitG_four(10),IOrbitGH_four(10)

Real*8 RRange_M,RM_M,Rmout(19)

Integer*8 IOrbitG_eight(2)

Character*50 PrefixInputfile,Inputfile,TrimInputfile

Character*8 Char_eight(2)

Character*20 Char_twenty(2)

Logical Flag,Alive,exist

 

C -------------End of Data Type Definition--------------------

 

 

C ---------------Format-----------------

105 Format(F15.3,F15.10,F20.10,3I6,I4,2I5,6I3,F17.3,I5,2I9)

c108 Format(I15,f6.3,f22.10,2I4,I20,f20.5,F22.10,I4)

c109 Format(I10,I10,3I5,2F25.10)

c110 Format(2I10,2I5,2D25.14)

C ---------------------------------------

C Begin of executable codes

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

C Open binary files C

C LP ODF file has eight independent data files, all of them are C

C in the form of binary,the format description is included in C

C the interface file "trk-2-18" which is also provided together C

C with the data. C

C C

C * Note C

C To read this file, you need to notice what kind of work station C

C or computer you are using in oder to figure out how the binary C

C files are handled on your machine, it is either Big_endian(MSB) or C

C Little_endian(LSB). For example, if you are using a PC, it would be C

C Little_endian, but if it is a SUN work station, there is a good C

C chance that it is Big_endian. C

C C

C In this case, we are running the program on a PC and we need to be C

C aware that it is Little_endian. However, The LP ODF files cannot C

C be recognized as you can see in the interface file that the data C

C type is MSB_INTEGER. If you are using a Big_endian machine, there C

C should be no problem. C

C C

C * How to solve this problem C

C This problem can be easily solved by adding the OPTION C

C "conver='big_endian'" in the OPEN STATMENT C

C C

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

 

C Inputfile='./LPODF/lpodfcase1'

Write(6,*) ' ***********************'

Write(6,*)

Write(6,*) ' PROGRAM READ LPODF '

Write(6,*)

Write(6,*) ' ***********************'

Write(6,*)

Write(6,*)

Write(6,*) ' -------------- INPUTFILE-------------'

Write(6,*) ' '

Write(6,*) ' Plese Select Inputfiles '

Write(6,*) ' From lpodfcase1 to 8 '

Write(6,*) ' '

Write(6,*) ' -------------------------------------'

Write(6,*)

Write(6,*)

Write(6,*) '-----------------------'

Write(6,*) ' The Input file is:'

C Read (*,*) Inputfile

Inputfile='lpodfcase1'

PrefixInputfile='./LPODF/'//Inputfile

 

TrimInputfile=PrefixInputfile(1:Len(Trim(PrefixInputfile)))

Write(6,*) '------------------------'

Write(6,*) 'The ODF file selected is'

write(*,*) TrimInputfile

Inquire(file=TrimInputfile,exist=Alive)

If (alive) Then

Write(6,*) ' ----------------READIND--------------'

Write(6,*)

Write(6,*) ' Start to read file ','Inputfile',

1 '…………'

Write(6,*)

Write(6,*) ' -------------------------------------'

Flag=Alive

Else

Write(6,*) ' -----------------ERROR---------------'

Write(6,*)

Write(6,*) ' The Inputfile dose not exist '

Write(6,*) ' Please check the Inputfile '

Write(6,*)

Write(6,*) ' -------------------------------------'

Stop

Endif

 

If (Flag) Then

Open(10,file=TrimInputfile,convert='big_endian',form='binary',

1 access='direct',recl=36)

 

C Open the OUTPUT files

Open(20,file='LPODF_ALLData_ASCII')

C First, the first 4 records we be read and write out

C !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

C Record 1:ODF lable group header

Read(10,rec=1) (Idat_four(i),i=1,9)

 

Write(20,*) 'LPODF ASCII OUTPUT:'

Write(20,*) '=====================================

1=============='

Write(20,*) '----------Label Group: Header Record-

1--------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (Idat_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' Label Group: Header Record'

Endif

 

C Record 2: ODF lable group

Read(10,rec=2) (Char_eight(i),i=1,2),(Idat_four(i),i=1,5)

Write(20,*) 'SystemID ProgramID SCID CrDate CrTime

1 FRefDate FRefTime'

Write(20,'(A10,A9,I6,I8,X,I8,3X,I8,3X,I6.6)')

1 (Char_eight(i),i=1,2),(Idat_four(i),i=1,5)

Write(20,*) '======================================

1=============='

Write(6,*)

Write(6,*) ' Label Group Data'

 

C Record 3:ODF identifier goup header

Read(10,rec=3) (Idat_four(i),i=1,9)

Write(20,*) '--------Identifier Group: Hearder Record-

1--------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (Idat_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' Identifier Group:Header Record'

 

C Record 4: ODF identifier goup data

Read(10,rec=4) (Char_eight(i),i=1,2),Char_twenty(1)

Write(20,*) ' Item_1 Item_2 Item_3'

Write(20,'(4X,A8,3X,A8,5X,A20)')

1 (Char_eight(i),i=1,2),Char_twenty(1)

Write(20,*) '======================================

1=============='

Write(6,*)

Write(6,*) ' Identifier Group Data'

 

C !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

R=5

Do While (.True.) !This program will never stop

Read(10,rec=R) Idat_four(1) !unless you tell them to

If (Idat_four(1).EQ.-1) Then !-1 means the EOF

Read(10,rec=R) (Idat_four(i),i=1,9) !we always have to figure whether

!it's the end of the data file

Write(20,*)'--------End of File Group: Hearder Record

1------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (Idat_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' End of the ODF file'

Goto 9999 !if it is the end of the data file

!then read the End of File Group

!and write it done, after that,exit

 

C!!!!!!!!!!!!!!!!!!!!!!! Reading the Orbit Data Group !!!!!!!!!!!!!!!!!!!!!!

 

Elseif (Idat_four(1).EQ.109) Then !find out if it is the Orbit Data

!Group.

 

3333 Read(10,rec=R) (IOrbitGH_four(i),i=1,9)

Write(20,*) '----------Orbit Data Group: Hearder Record

1------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (IOrbitGH_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' Orbit Data Group: Hearder'

Write(6,*)

Write(6,*) ' Orbit Data Group'

Write(20,'(A12,A15,A20,A10,2A6,A4,A6,7A3,A17,3A8)')

1'Time Tag', 'DL_delay','Observable','Fmt','DSSr','DSSt','Net',

2'D_Typ','DL','UL','Ex','V','2S','SC','17',

3'Reference Freq','Item_20','Item_21','Item_22'

R=R+1 !if it is the Orbit Data Group,

!write out the Header first. after

!the Header is the data group,

!read these data and write out

!till the EOF

Do While (.True.)

Read(10,rec=R) Idat_four(1)

If (Idat_four(1).EQ.-1) Then

Read(10,rec=R) (Idat_four(i),i=1,9)

Write(20,*)'--------End of file Group: Hearder Record

1------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (Idat_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' End of the ODF file'

Goto 9999

Elseif (Idat_four(1).EQ.109) Then !this is in case there is more

!than one Orbit Data Group,

!though, there shoulb be only one

Goto 3333

Elseif (Idat_four(1).EQ.2030) Then

Goto 4444

Else

 

Read(10,rec=R) (IOrbitG_four(i),i=1,5),(IOrbitG_eight(i),i=1,2)

 

Write(20,105)

1IOrbitG_four(1)+ibits(IOrbitG_four(2),22,10)*1.D-3, !Record Time Tag

3ibits(IOrbitG_four(2),0,22)*1.D-9, !Downlink_delay

4IOrbitG_four(3)+IOrbitG_four(4)*1.D-9, !Observable

5ibits(IOrbitG_four(5),29,3), !Format ID=2

5ibits(IOrbitG_four(5),22,7), !Receiving Station ID NO.--DSSr

5ibits(IOrbitG_four(5),15,7), !Transmittiong Station ID NO.--DSSt

5ibits(IOrbitG_four(5),13,2), !Network ID for Transmitting Station

5ibits(IOrbitG_four(5),7,6), !Data Type ID.--D_Typ

5ibits(IOrbitG_four(5),5,2), !Downlink band ID.--DL

5ibits(IOrbitG_four(5),3,2), !Uplink band ID.--UL

5ibits(IOrbitG_four(5),1,2), !Exciter band ID.--Ex

5ibits(IOrbitG_four(5),0,1), !Data Validity Indicator (0 for good)

6ibits(IOrbitG_eight(1),57,7), !2nd Receiving station ID

6ibits(IOrbitG_eight(1),47,10), !Spacecraft ID

6ibits(IOrbitG_eight(1),46,1), !Receiver/Exciter Independet flag(1:yes)

6ibits(IOrbitG_eight(1),24,22)*1.D4+

6ibits(IOrbitG_eight(1),0,24), !

7ibits(IOrbitG_eight(2),44,20), !

7ibits(IOrbitG_eight(2),22,22), !

7ibits(IOrbitG_eight(2),0,22) !

Endif

 

 

R=R+1

Enddo

!the most important data needed in

!the ODF file is the Orbit Data Groups

!and there should be only one of them

!in each ODF file.

4444 Continue

 

 

C!!!!!!!!!!!!!!!!!!! Reading the Ramp Group !!!!!!!!!!!!!!!!!!!!!!!!!!

Elseif (Idat_four(1).EQ.2030) Then !if it is not the EOF, then find

!out if it is the Ramp Group.

Iramp=1

 

2222 Read(10,rec=R) (IRampGH_four(i),i=1,9)

Write(20,*) '-----------Ramp Group: Hearder Record-

1--------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (IRampGH_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' Ramp Group',Iramp

Iramp=Iramp+1

Write(20,'(A5,A23,A18,2A21)')

1'DSS', 'Ramp Start Time','Ramp Rate',

2'Start Frequency','Ramp End Time'

R=R+1 !if it is the Ramp Group, then

!write out the Header first. after

!the Header is the data group,

!read these ramp data and write out

!till the EOF

Do While (.True.)

Read(10,rec=R) Idat_four(1)

If (Idat_four(1).EQ.-1) Then

Read(10,rec=R) (Idat_four(i),i=1,9)

Write(20,*)'--------End of file Group: Hearder Record

1------------'

Write(20,*) 'Prime_Key Second_Key Log_Rec_Len

1 Gp_St_Pkt_No'

Write(20,'(I8,2X,I8,6X,I8,8X,I8)') (Idat_four(i),i=1,4)

Write(20,*) '----------------------Data-----------

1--------------'

Write(6,*)

Write(6,*) ' End of the ODF file'

Goto 9999

Elseif (Idat_four(1).EQ.2030) Then !there is more than one ramp group

!in an ODF file.

Goto 2222

Else

Read(10,rec=R) (IRampG_four(i),i=1,9)

 

Write(20,'(I4,F25.10,F18.10,2F23.10)')

1ibits(IrampG_four(5),0,10),

2(IRampG_four(1)+IRampG_four(2)*1.D-9),

3IRampG_four(3)+IRampG_four(4)*1.D-9,

4ibits(IRampG_four(5),10,22)*1.D9+IRampG_four(6)+

4IRampG_four(7)*1.D-9,

5(IRampG_four(8)+IRampG_four(9)*1.D-9)

Endif

R=R+1

Enddo

 

Endif

R=R+1

Enddo

 

Close(10)

Close(20)

 

9999 Continue

End Program