anchor.diff:          IF(ISYS.EQ.0) NSYS=1
           IF(ISYS.EQ.1) NSYS=2
           IF(ISYS.EQ.2) NSYS=4
  -On branch master
  -Your branch is up-to-date with 'origin/master'.
  -Changes to be committed:
  -  (use "git reset HEAD <file>..." to unstage)
  -
  -	new file:   .mesh.f90.swp
  -	modified:   mesh.f90
   
           READ(3,*)   NELEMF
       
  -        !    mvarÖÐNELEM: number of total elements
  +        ! =========================================
   
           NELEM=NELEMB+NELEMF
   
  -        !WRITE(11,*) ' NELEM=',NELEM,'  IOPL=',IPOL
           ALLOCATE (NCONB(NELEMB,8),NCONDB(NELEMB,8))!body node/normal list
           ALLOCATE (XYZB(3,NNB),DXYZB(3,NNBD))    !body node/normal
           
  
positive.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)
  !    print *, M, IPOLAR(I), XOFSET(I),YOFSET(I),ZOFSET(I)
  5       CONTINUE
  !
          DO 10 INODE=1, NNB
  
negative.code:         allocate(NNORMN(8*NELEM) )
  
          ALLOCATE( XYZE(3,8,NELEM),DXYZE(3,8,NELEM))
          !allocate(DAMPE(8,NELEM))
          ALLOCATE( TXYZE(3,8,NELEM))
          ALLOCATE( XYZTP(3,8*NELEM),DXYZTP(3,8*NELEM))!,DAMPTP(8*NELEM))
  
          !ALLOCATE(SAMB(NELEM,16,0:8),SAMBXY(NELEM,16,3),
          !&         DSAMB(NELEM,16,6))
  
  
          call MESHFS4()! Read in data on free surface mesh
          call MESHBD(IPOL) ! Read in data on body mesh
  
  
  
  
  
      end subroutine
  
  !       MESHFS4 + MESHBD 
  