From c00b2fe4280338f5aceb95f1a3958d8a38018af7 Mon Sep 17 00:00:00 2001 From: Pavel Demin Date: Fri, 20 Apr 2018 21:54:27 +0200 Subject: [PATCH] remove unpackbits --- packjt.f90 | 27 --------------------------- 1 file changed, 27 deletions(-) diff --git a/packjt.f90 b/packjt.f90 index c23e685..e33b955 100644 --- a/packjt.f90 +++ b/packjt.f90 @@ -1,34 +1,7 @@ module packjt -! These variables are accessible from outside via "use packjt": - integer jt_itype,jt_nc1,jt_nc2,jt_ng,jt_k1,jt_k2 - character*6 jt_c1,jt_c2,jt_c3 - contains - subroutine unpackbits(sym,nsymd,m0,dbits) - - ! Unpack bits from sym() into dbits(), one bit per byte. - ! NB: nsymd is the number of input words, and m0 their length. - ! there will be m0*nsymd output bytes, each 0 or 1. - - integer sym(:) - integer*1 dbits(:) - - k=0 - do i=1,nsymd - mask=ishft(1,m0-1) - do j=1,m0 - k=k+1 - dbits(k)=0 - if(iand(mask,sym(i)).ne.0) dbits(k)=1 - mask=ishft(mask,-1) - enddo - enddo - - return - end subroutine unpackbits - subroutine unpackcall(ncall,word,iv2,psfx) parameter (NBASE=37*36*10*27*27*27)