// begin of Neumann selection blocks For 2)Load Surfaces
*Set Var IntervalNum=0
*loop intervals
*Set Var IntervalNum=Operation(IntervalNum+1)
*#
*If(ndime==3)
*Set Cond Neumann_Temperature_for_Surface *elems *CanRepeat 
*Add Cond Neumann_Humidity_for_Surface *elems *CanRepeat 
*EndIf
*#
*If(CondNumEntities(int)>0)
// selection has *CondNumEntities entities
*loop elems *OnlyInCond
*# ElemsNum *ElemsType *ElemsNnode SurfaceNum = *Cond(Entity_number,Int)
*If(ElemsType==4 && ElemsNnode == 4)
SELECTION *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum" LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) ;
*ElseIf(ElemsType==4 && ElemsNnode==10)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum" LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) *GlobalNodes(5) *GlobalNodes(6) ;
*ElseIf(ElemsType==5 && ElemsNnode==8)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) ;
*ElseIf(ElemsType==5 && ElemsNnode==20)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) *GlobalNodes(5) *GlobalNodes(6) *GlobalNodes(7) *GlobalNodes(8) ;
*ElseIf(ElemsType==5 && ElemsNnode==27)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) *GlobalNodes(5) *GlobalNodes(6) *GlobalNodes(7) *GlobalNodes(8) *GlobalNodes(9) ;
*ElseIf(ElemsType==6 && ElemsNnode==6)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum" LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) ; 
*ElseIf(ElemsType==3 && ElemsNnode == 4)
SELECTION *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum" LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) ;
*ElseIf(ElemsType==3 && ElemsNnode==8)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum" LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) *GlobalNodes(5) *GlobalNodes(6) *GlobalNodes(7) *GlobalNodes(8) ;
*ElseIf(ElemsType==2 && ElemsNnode==3)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) ;
*ElseIf(ElemsType==2 && ElemsNnode==6)
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum *ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes(1) *GlobalNodes(2) *GlobalNodes(3) *GlobalNodes(4) *GlobalNodes(5) *GlobalNodes(6) ;
*Else
SELECTION  *\
*Format "%i%i%i%s%i%i"
	   "$N$I*IntervalNum Else*ElemsType *ElemsNnode Neumann*Cond(Boundary)Num*Cond(Entity_number,Int) ElemNum*elemsnum"  LIST *GlobalNodes ; 
*EndIf
*End elems
*EndIf
*# 1=Linear,2=Triangle,3=Quadrilateral,4=Tetrahedra,5=Hexahedra,6=Prism,7=OnlyPoints
*End intervals

