anchor.diff:  !C *                                                                 *
   !C *******************************************************************
   !C 
  -       SUBROUTINE MESHFS4()
  -
  -!     USE MVAR_MOD
  -        IMPLICIT   NONE  
  -        
  -
  -        INTEGER IE,J,M
  +    subroutine meshfs4()
  +
  +        implicit none
  +
  +        integer :: IE,J,M
           
           XYZE(3,1:8, 1:NELEMF)=0.0d0
  -!
  +        
           DO 100 IE=1, NELEMF
               IETYPE(IE)=2
               READ(3, *)    M, NCN(IE)
  
positive.code:       !dampe(:,:)=0.0d0
  
  100    CONTINUE
  !
  !
  !  ÊýŸÝÎÄŒþÖÐÎïÃæ·š·œÏòÖžÈëÁ÷ÌåÎªÕý
  !
        DXYZE(1, 1:8, 1:NELEMF)= 0.0d0
        DXYZE(2, 1:8, 1:NELEMF)= 0.0d0
        DXYZE(3, 1:8, 1:NELEMF)= 1.0d0 
  !     DXYZE(3, 1:8, 1:NELEMF)=-1.0d0 
  
  
        RETURN
        END
  
negative.code: !     INTEGER NELEM,NELEMB,NELEMF,NNODE,NNODED,NNB,NNBD,NNF,NNTCH
  
  
          INTEGER NCOR,IPL,IPOLAR(50)
          INTEGER I,IE,M,INODE,NCNN,K,KK
          REAL*8  A1,A2,R1,R2,Z1
          REAL*8  XOFSET(50),YOFSET(50),ZOFSET(50)
          real(8),parameter :: pi = 3.141592653589793
          !
          ! --------------------------------------------
          DO 5 I=1, NCOR
          READ(2,*)  M, IPOLAR(I), XOFSET(I),YOFSET(I),ZOFSET(I)
  5       CONTINUE
  !
          DO 10 INODE=1, NNB
  