FunctionSpace { 

  // Magnetic vector potential a (b = curl a)
  { Name Hcurl_a_3D ; Type Form1 ;
    BasisFunction {// a = a_e * s_e
      { Name se ; NameOfCoef ae ; Function BF_Edge ;
        Support Region[{conducting, source, non_conducting}] ; Entity EdgesOf[ All ] ; }
    }
    Constraint {
      { NameOfCoef ae  ; EntityType EdgesOf ; NameOfConstraint MVP_3D ; }
	  { NameOfCoef ae  ; EntityType EdgesOfTreeIn ; EntitySubType StartingOn ;
	NameOfConstraint Gauge ; }
    }
  }

  // Electric scalar potential (3D)
  { Name Hregion_u_3D ; Type Form1 ;
    BasisFunction {
      { Name sr ; NameOfCoef e ; Function BF_GradNode ;
        Support Region[{conducting, source}] ; Entity NodesOf[ All , Not terms ] ; }
	  { Name sc1; NameOfCoef V1; Function BF_GroupOfEdges;
	Support Region[{conducting, source}]; Entity GroupsOfEdgesOf[{Sur_CURRENT_CH, Sur_VOLTAGE_CH}]; }
    }
	
    GlobalQuantity {
      { Name U ; Type AliasOf        ; NameOfCoef V1 ; }
      { Name I ; Type AssociatedWith ; NameOfCoef V1 ; }
    }
    Constraint {
      { NameOfCoef U ; EntityType GroupsOfEdgesOf ; NameOfConstraint V_3D ; }
      { NameOfCoef I ; EntityType GroupsOfEdgesOf ; NameOfConstraint I_3D ; }
    }
	
  }
  
  { Name Hgrad_T; Type Form0; 
    BasisFunction { 
      { Name sn; NameOfCoef Tn; Function BF_Node; Support Tot_The; 
        Entity NodesOf[All]; } 
    } 
    Constraint { 
      { NameOfCoef Tn; EntityType NodesOf ; NameOfConstraint temperature; } 
    } 
  } 
}