Ask Question

Name:
Title:
Your Question:

Answer Question

Name:
Your Answer:
User Submitted Source Code!


Description:
  readhkl
Language: FORTRAN
Code:
           REAL*4 ESD
      OPEN(UNIT=1,FILE='hklIsig.txt',FORM='FORMATTED',STATUS='OLD')
      OPEN(UNIT=2,FILE='Factheo.hkl',FORM='FORMATTED')
c
      XP=0.
      IH = 0
      II=1
      ESD2 = 0.01
   10 READ(1,*,END=20)JH,JK,JL,Y,ESD1,II
      if(Y.NE.0.0) then
      WRITE(2,100) JH,JK,JL,Y/10.,ESD2,II
      endif
c
      GO TO 10
c
   20 WRITE(2,100)IH,IH,IH,XP,XP,II
c
 100  format(3I4,2X,F8.2,2X,F8.2,I4)
      Print*, (' FINI ')
      END     
          
Comments: