anchor.diff:    integer itr
     
   
  -  tx0 = 200 ! T ! get initial guess
  +  tx0 = T ! get initial guess
    if(tx0>330.D0) then
       tx0=300.D0  
    endif
  
positive.code:   logical succ
    
  
    tx0 = 200 ! T ! get initial guess
   if(tx0 > 330.D0) THEN 
      tx0=300.D0  
   endif
  
  write(*,*) tx0
    
    if (tx0 <1D0) then ! obtain initial guess
      tx0 = 100D0 ! or use constant cp 
    endif 
  
  
negative.code:   tx=tx0
    dt=1D-7
    
  
  write(*,*) tx0
  
    do 
    succ=supst(tx, p, dg, hg)
     hx = hg * energyscale
  