anchor.diff:  
   DOUBLE PRECISION, INTENT(IN) :: t
   
  -getcpc = 252.54 + 0.11474*t
  +getcpc = 252.54d0 + 0.11474d0*t
   
   END FUNCTION getcpc
   
  
positive.code: 
  DOUBLE PRECISION, INTENT(IN) :: t
  
  getkc = 7.51 + 2.09E-2*t - 1.45E-5*t**2 + 7.67E-9*t**3
  
  END FUNCTION getkc
  
  
negative.code: 
  DOUBLE PRECISION, INTENT(IN) :: t
  
  getcpf = 162.3d0 + 0.3038d0*t - 2.391d-4*t**2 + 6.404d-8*t**3
  
  END FUNCTION getcpf
  
  