hdiff output

r31575/Cv.BS.f90 2016-12-03 16:30:09.943329643 +0000 r31574/Cv.BS.f90 2016-12-03 16:30:10.195332897 +0000
 53:       Z2=Z2+DUMMY*(ENERGY(J2)-ENERGY(1))**2 53:       Z2=Z2+DUMMY*(ENERGY(J2)-ENERGY(1))**2
 54:    ENDDO 54:    ENDDO
 55:    MEANE(J1)=Z1/Z0 55:    MEANE(J1)=Z1/Z0
 56:    ZSAVE(J1)=Z0 56:    ZSAVE(J1)=Z0
 57:    IF (DELTA/TEMPERATURE.LT.1.0D-7) THEN 57:    IF (DELTA/TEMPERATURE.LT.1.0D-7) THEN
 58:       ONEMEXP=-DELTA/TEMPERATURE 58:       ONEMEXP=-DELTA/TEMPERATURE
 59:    ELSE 59:    ELSE
 60:       ONEMEXP= 1.0D0-EXP(DELTA/TEMPERATURE) 60:       ONEMEXP= 1.0D0-EXP(DELTA/TEMPERATURE)
 61:    ENDIF 61:    ENDIF
 62:    WRITE(1,'(7G20.10)') TEMPERATURE, Z0, Z1, Z2, & 62:    WRITE(1,'(7G20.10)') TEMPERATURE, Z0, Z1, Z2, &
 63:   &                     KAPPA/2.0D0 + & 63:   &                     KAPPA - (Z1/(Z0*TEMPERATURE))**2 + Z2/(Z0*TEMPERATURE**2), MEANE(J1), ZSAVE(J1)
 64:   &                      1.0D0*(1.0D0 - DELTA**2*EXP(DELTA/TEMPERATURE)/(ONEMEXP**2*TEMPERATURE**2)) & 
 65:   &                     - (Z1/(Z0*TEMPERATURE))**2 + Z2/(Z0*TEMPERATURE**2), MEANE(J1), ZSAVE(J1) 
 66: ENDDO 64: ENDDO
 67: CLOSE(1) 65: CLOSE(1)
 68:  66: 
 69: OPEN(UNIT=1,FILE='prob.out.BS',STATUS='UNKNOWN') 67: OPEN(UNIT=1,FILE='prob.out.BS',STATUS='UNKNOWN')
 70: DO J1=1,NTEMP 68: DO J1=1,NTEMP
 71:    TEMPERATURE=TMIN+(J1-1)*TINT 69:    TEMPERATURE=TMIN+(J1-1)*TINT
 72:    PSUM=0.0D0 70:    PSUM=0.0D0
 73:    PSUM2=0.0D0 71:    PSUM2=0.0D0
 74:    DO J2=1,NBINS 72:    DO J2=1,NBINS
 75:       DUMMY=EXP(-(ENERGY(J2)-ENERGY(1))/TEMPERATURE+LNWEIGHT(J2)) 73:       DUMMY=EXP(-(ENERGY(J2)-ENERGY(1))/TEMPERATURE+LNWEIGHT(J2))
 76:       DP=DUMMY*(ENERGY(J2)-ENERGY(1)-MEANE(J1))/(ZSAVE(J1)*TEMPERATURE**2) 74:       DP=DUMMY*(ENERGY(J2)-ENERGY(1)-MEANE(J1))/(ZSAVE(J1)*TEMPERATURE**2)
 77:       PSUM=PSUM+DP*(ENERGY(J2)-ENERGY(1)-MEANE(J1)) 75:       PSUM=PSUM+DP*(ENERGY(J2)-ENERGY(1)-MEANE(J1))
 78: !     PRINT '(A,2I6,4G20.10)','J1,J2,T,ZSAVE,DUMMY,DP=',J1,J2,TEMPERATURE,ZSAVE(J1),DUMMY,DP 76: !     PRINT '(A,2I6,4G20.10)','J1,J2,T,ZSAVE,DUMMY,DP=',J1,J2,TEMPERATURE,ZSAVE(J1),DUMMY,DP
 79:       PSUM2=PSUM2+(ZSAVE(J1)/DUMMY)*(TEMPERATURE*DP)**2 77:       PSUM2=PSUM2+(ZSAVE(J1)/DUMMY)*(TEMPERATURE*DP)**2
 80:    ENDDO 78:    ENDDO
 81:    IF (DELTA/TEMPERATURE.LT.1.0D-7) THEN 79:    WRITE(1,'(5G20.10)') TEMPERATURE, KAPPA+1.0D0 - DELTA**2*EXP(DELTA/TEMPERATURE)/(ONEMEXP**2*TEMPERATURE**2)+PSUM, &
 82:       ONEMEXP=-DELTA/TEMPERATURE 80:   &                                  KAPPA+1.0D0 - DELTA**2*EXP(DELTA/TEMPERATURE)/(ONEMEXP**2*TEMPERATURE**2)+PSUM2
 83:    ELSE 
 84:       ONEMEXP= 1.0D0-EXP(DELTA/TEMPERATURE) 
 85:    ENDIF 
 86:    WRITE(1,'(5G20.10)') TEMPERATURE, KAPPA/2.0D0+1.0D0 - DELTA**2*EXP(DELTA/TEMPERATURE)/(ONEMEXP**2*TEMPERATURE**2)+PSUM, & 
 87:   &                                  KAPPA/2.0D0+1.0D0 - DELTA**2*EXP(DELTA/TEMPERATURE)/(ONEMEXP**2*TEMPERATURE**2)+PSUM2 
 88: ENDDO 81: ENDDO
 89: CLOSE(1) 82: CLOSE(1)
 90:  83: 
 91: INQUIRE(FILE='Tanal.BS',EXIST=YESNO) 84: INQUIRE(FILE='Tanal.BS',EXIST=YESNO)
 92:  85: 
 93: IF (YESNO) THEN 86: IF (YESNO) THEN
 94:    NPLUS=0.0D0 87:    NPLUS=0.0D0
 95:    NMINUS=0.0D0 88:    NMINUS=0.0D0
 96:    DPLUS=0.0D0 89:    DPLUS=0.0D0
 97:    DMINUS=0.0D0 90:    DMINUS=0.0D0


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0