anchor.diff:                sinal=-1.0D0
              endif
    
  -          !Determinando os valores para o node i 
  +           !Determinando os valores para o node i 
              xx = mesh%v(i)%p(1)
              yy = mesh%v(i)%p(2)
              zz = mesh%v(i)%p(3)
  
positive.code: !          phi%name=trim(adjustl(trim(transpname))//"_phi_t_"//trim(adjustl(trim(atime))))    
         end do
      end do
  !    print*, erro_Linf!, 'erro_Linf'
  !    print*, dsqrt (temp/contador)! , 'erro_L2'
  !    print*, contador 
  
negative.code:            yy = mesh%v(w)%p(2)
             zz = mesh%v(w)%p(3)
             call aplyr(xx,yy,zz,cx,sx,cy,sy,xp,yp,zp) 
  
             !cosseno=(xp/dsqrt(xp**2+yp**2))
             !seno=(yp/dsqrt(xp**2+yp**2))
             !derivada_i = 2.0D0*node(i)%coef(4)*(cosseno**2) + & 
             !2.0D0*node(i)%coef(5)*cosseno*seno + 2.0D0*node(i)%coef(6)*(seno**2)
   
             phi_i = node(i)%coef(1)
  
             ! Compute the second derivative in normal direction
             dir_i = node(i)%G(1)%lvn(n,1:3)
             dir_i = proj_vec_sphere(dir_i, mesh%v(i)%p)
             call aplyr(dir_i(1),dir_i(2),dir_i(3),cx,sx,cy,sy,rot_dir_i(1),rot_dir_i(2),rot_dir_i(3))
             rot_dir_i = rot_dir_i/norm(rot_dir_i)
             derivada_i = 2.d0*node(i)%coef(4)*rot_dir_i(1)*rot_dir_i(1) + &
                          2.d0*node(i)%coef(5)*rot_dir_i(1)*rot_dir_i(2) + &
                          2.d0*node(i)%coef(6)*rot_dir_i(2)*rot_dir_i(2)
  
             !Determinando os valores dos vizinhos do node i 
             xx = mesh%v(w)%p(1)
             yy = mesh%v(w)%p(2)
             zz = mesh%v(w)%p(3)
             call constr(xx,yy,zz,cx,sx,cy,sy)
  
             xx = mesh%v(i)%p(1)
             yy = mesh%v(i)%p(2)
             zz = mesh%v(i)%p(3)
             call aplyr(xx,yy,zz,cx,sx,cy,sy,xp,yp,zp) 
  
             !cosseno=(xp/dsqrt(xp**2+yp**2))
             !seno=(yp/dsqrt(xp**2+yp**2))
             !derivada_j = 2.0D0*node(w)%coef(4)*(cosseno**2) + &
             !2.0D0*node(w)%coef(5)*cosseno*seno + 2.0D0*node(w)%coef(6)*(seno**2)
  
             phi_j = node(w)%coef(1)
  
             ! Compute second derivate in the normal direction to the edge
             dir_j = -node(i)%G(1)%lvn(n,1:3)
             dir_j = proj_vec_sphere(dir_j, mesh%v(w)%p)
             call aplyr(dir_j(1),dir_j(2),dir_j(3),cx,sx,cy,sy,rot_dir_j(1),rot_dir_j(2),rot_dir_j(3))
             rot_dir_j = rot_dir_j/norm(rot_dir_j)
             derivada_j = 2.d0*node(w)%coef(4)*rot_dir_j(1)*rot_dir_j(1) + &
                          2.d0*node(w)%coef(5)*rot_dir_j(1)*rot_dir_j(2) + &
                          2.d0*node(w)%coef(6)*rot_dir_j(2)*rot_dir_j(2)
  
             !Distancia entre o node i e seu respectivo vizinho   
             dist=arclen(mesh%v(i)%p,mesh%v(w)%p)           
  
  !           p1 = mesh%v(i)%p
  !           call cart2sph (p1(1), p1(2), p1(3),lon,lat)
  