mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-20 02:52:00 -05:00
26 lines
633 B
Fortran
26 lines
633 B
Fortran
subroutine get_crc24(mc,ncrc)
|
|
!
|
|
! 1. To calculate 24-bit CRC, mc(1:50) is the message and mc(51:74) are zero.
|
|
! 2. To check a received CRC, mc(1:74) is the received message plus CRC.
|
|
! ncrc will be zero if the received message/CRC are consistent.
|
|
!
|
|
character c24*24
|
|
|
|
integer*1 mc(74),r(25),p(25)
|
|
integer ncrc
|
|
! polynomial for 24-bit CRC 0x100065b
|
|
data p/1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,1,0,1,1,0,1,1/
|
|
|
|
! divide by polynomial
|
|
r=mc(1:25)
|
|
do i=0,49
|
|
r(25)=mc(i+25)
|
|
r=mod(r+r(1)*p,2)
|
|
r=cshift(r,1)
|
|
enddo
|
|
|
|
write(c24,'(24b1)') r(1:24)
|
|
read(c24,'(b24.24)') ncrc
|
|
|
|
end subroutine get_crc24
|