anchor.diff:  END FUNCTION geths
   
   
  -SUBROUTINE th_trans3(xpline, h)
  +SUBROUTINE th_trans4(xpline, h)
   
   !
   ! Purpose:
  
positive.code: 	bc2 = bc
  	ke1 = ke2
  	ke2 = ke
      WRITE(ounit,'(I5, F15.1, F23.5, ES15.5, ES17.5)') n, bc, Ke, ser, fer
  	IF ((rKe == 100000) .AND. (ser < 1.e-2) .AND. (fer < 1.e-2)) EXIT
  	n = n + 1
  	IF (bc > 2999. .AND. bc < 3000.) THEN
  
negative.code: REAL :: t1, rho1, ent1
  REAL :: t2, rho2, ent2
  
  IF ((t < 473.15) .OR. (t > 617.91)) THEN
      WRITE(ounit,*) '  Coolant temp. : ', t
      WRITE(ounit,*) '  ERROR : MODERATOR TEMP. IS OUT OF THE RANGE OF DATA IN THE STEAM TABLE'
      WRITE(ounit,*) '  CHECK INPUT MASS FLOW RATE OR POWER'
  