anchor.diff:      !Converte para coordenadas esféricas
       call cart2sph (p(1), p(2), p(3), lon, lat)
   
  -    !TESTE 1 
  +    !TESTE 1
       !lap_ext= (((3*dsin(lon)+2*dcos(lon))*(1-dcos(lat)**2+dsin(lat)**2))/6._r8*dcos(lat)) &
       !- 4._r8*dsin(lat)/3._r8
   
  
positive.code:   !==========================================================================
    ! Subroutines for relaxation, interpolation and transfer residue
    !==========================================================================
  
    subroutine lap_cell(mesh,lapc,fap)
      !------------------------------------------------------------
  
negative.code:   end function lap_ext
  
  
  
  
    subroutine relaxac(mesh)
  