anchor.diff:                   * ( pre(k)/PRE00 )**KAPPA
             tem(k) = max( 200.D0, tem(k) )
   
  -          if( abs(pre_save-pre(k)) <= eps ) exit
  +          if( abs(pre_save-pre(k)) <= eps_hs ) exit
          enddo
   
          if ( itr > itrmax ) then
  
positive.code: 
      real(8) :: dT, f, df
      real(8) :: lat, lon
  
      integer :: n, k, l, itr
      !---------------------------------------------------------------------------
  
negative.code:                     * ( pre(k)/PRE00 )**KAPPA
               tem(k) = max( 200.D0, tem(k) )
  
               if( abs(pre_save-pre(k)) <= eps_hs ) exit
  
            enddo
  
  