anchor.diff:  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'
  +    WRITE(ounit,*) '  CHECK INPUT COOLANT MASS FLOW RATE OR CORE POWER'
       WRITE(*,*) '  Coolant temp. : ', t
       WRITE(*,*) '  ERROR : MODERATOR TEMP. IS OUT OF THE RANGE OF DATA IN THE STEAM TABLE'
  -    WRITE(*,*) '  CHECK INPUT MASS FLOW RATE OR POWER'
  +    WRITE(*,*) '  CHECK INPUT COOLANT MASS FLOW RATE OR CORE POWER'
       STOP
   END IF
   
  
positive.code: IMPLICIT NONE
  
  SAVE
  
  ! REAL :: kv   ! Water kinematic viscosity
  ! REAL :: Pr  !Prandtl Number
  ! REAL :: tcon  ! Thermal conductivity (W/mK)
  
  CONTAINS
  
  
negative.code: IF ((ent < 858341.5) .OR. (ent > 1624307.1)) THEN
      WRITE(ounit,*) '  Enthalpy. : ', ent
      WRITE(ounit,*) '  ERROR : ENTHALPY IS OUT OF THE RANGE OF DATA IN THE STEAM TABLE'
      WRITE(ounit,*) '  CHECK INPUT COOLANT MASS FLOW RATE OR CORE POWER'
      WRITE(*,*) '  Enthalpy. : ', ent
      WRITE(*,*) '  ERROR : ENTHALPY IS OUT OF THE RANGE OF DATA IN THE STEAM TABLE'
      WRITE(*,*) '  CHECK INPUT COOLANT MASS FLOW RATE OR CORE POWER'
      STOP
  END IF
  
  