anchor.diff:  
   DOUBLE PRECISION :: cvelo, Nu, Re
   
  -cvelo = cflow / (farea * xden * 1000.)        ! Calculate flow velocity (m/s)
  +cvelo = cflow / (farea * xden * 1000.d0)        ! Calculate flow velocity (m/s)
   Re = cvelo * dh / (kv * 1.d-6)                 ! Calculate Reynolds Number
  -Nu = 0.023*(Pr**0.4)*(Re**0.8)                ! Calculate Nusselt Number
  +Nu = 0.023d0*(Pr**0.4d0)*(Re**0.8d0)                ! Calculate Nusselt Number
   geths = (tc / dh) * Nu                        ! Calculate heat transfer coefficient
   
   
  
positive.code: 
  DOUBLE PRECISION, INTENT(IN) :: t
  
  getcpc = 252.54 + 0.11474*t
  
  END FUNCTION getcpc
  
  
negative.code: 
  DOUBLE PRECISION, INTENT(IN) :: t
  
  getkc = 7.51d0 + 2.09d-2*t - 1.45d-5*t**2 + 7.67d-9*t**3
  
  END FUNCTION getkc
  
  