*heading
TEST EXTRAPOLATION=PARABOLIC. ROTATE SHELL ON EDGE BY 45 DEG. 
*restart,write
*node
1,0,0,0
2,1,0,0
3,1,1,0
4,0,1,0
*element,type=s4r,elset=eall
1,1,2,3,4
*shell section,elset=eall,mat=elas
0.1, 
*material,name=elas
*elastic
1.e9, 
*boundary
1,1,3
1,5,6
2,1,3
*step,nlgeom,extrap=parabolic
**step,nlgeom
*static
0.1,1.0
*boundary
1,4,4,0.78539812
*node file
u,rf
*controls,parameter=field
** adjust the intial time average flux
 ,,0.1
*end step