anchor.diff:            real*8, dimension(npro,nflow) :: temp0, temp1 ! local temporary array
             integer :: iel                                                  
   c
  -          factor = one
  +          factor = zero
   c
             u_ref_0(:,1) = 1.0d0 ! hacking, rho_ref, gas phase
             u_ref_0(:,2) = 1.0d0 * 1.0d2 ! hacking, rho_ref*v_ref
  
positive.code: c          
          end subroutine calc_ch
  c
          subroutine e3if_dc_res
  c...To be filled        
          end subroutine e3if_dc_res
  c
          subroutine e3if_dc_egmass
  c...To be filled
          end subroutine e3if_dc_egmass
          
          subroutine e3if_dc
  
negative.code:                                                           ! proj is symetric
            real*8, dimension(npro,nsd,nsd) :: pt_g_p0, pt_g_p1 ! proj_{ik} g^ij proj_{jm}
            integer :: iel
  c... get the tangential projector          
            call calc_projector(proj, nv0) ! get the tangential projector
  c... get the c^h
            call calc_ch(ch0, ch1, f_jump)
  c          
            call e3giju (giju0, dxidx0, npro, nsd, lcsyst) ! get g^{ij}, 0 side
            call e3giju (giju1, dxidx1, npro, nsd, lcsyst) ! get g^{ij}, 0 side
  c.... get the full g^ij matrix
  