CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C C THOMAS J. 03/2008 C C Subroutine to check if all information from Sander C C are correctly put in Divcon C C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC subroutine checking(file1) implicit double precision (a-h,o-z) #include "divcon.dim" #include "divcon.h" integer ierror, i, j character*80 file1 call opnfil(1,ierror) C open(unit=115,file=file1,action='write',iostat=ierror) if (ierror .ne. 0) return i = 1 j = 1 do while (keywrd(i+70:j+139) .ne. '') if(keywrd(j+69:j+69) .eq. ' ') then write(inpt,'(A,A)') keywrd(i:j+69), ' &' i = j + 70 j = i else j = j - 1 endif enddo write(inpt,'(A,A)') keywrd(i:i+69), ' &' if (pbc .or. resclust) then if (pif .or. mais) then write(inpt,'(A)') 'CARTESIAN DIRECT GRADIENT CENTRAL & RESIDUE PBC PM3' else write(inpt,'(A)') 'CARTESIAN DIRECT GRADIENT CENTRAL & RESIDUE PBC' endif else if (pif .or. mais) then write(inpt,'(A)') 'CARTESIAN DIRECT GRADIENT CENTRAL & PM3 RESIDUE' else write(inpt,'(A)') 'CARTESIAN DIRECT GRADIENT CENTRAL' endif endif write(inpt,*) j = 1 do i=1,natoms if (irpnt(j) .eq. i) then write(inpt,'(I7,A3,F14.7,F14.7,F14.7,A5)') i, & symbol(iatnum(i)), xyz(1,i), xyz(2,i), xyz(3,i), ' RES' j = j + 1 else write(inpt,'(I7,A3,F14.7,F14.7,F14.7)') i, & symbol(iatnum(i)), xyz(1,i), xyz(2,i), xyz(3,i) endif enddo write(inpt,'(A)') 'END_COORD' if (pbc) then write(inpt,'(A)') 'BOX' write(inpt,'(A,E14.8,A,E14.8,A,E14.8)') ' XBOX= ',dbox(1), & ' YBOX= ', dbox(2), ' ZBOX= ', dbox(3) write(inpt,'(A)') 'END_BOX' endif i = 1 do while (options2(i) .ne. ' ') write(inpt,'(A80)') options2(i) i = i + 1 enddo close(inpt) end