*HEADING
 test marlow model with s4r, incompressible
*preprint,model=yes
*NODE,NSET=ALL
1,
2,1.
3,1.,1.,
4,0.,1.,
*NSET,NSET=FACE1
1,2,3,4
*NSET,NSET=FACE3
1,2
*NSET,NSET=FACE4
3
*NSET,NSET=FACE5
3,4
*NSET,NSET=FACE6
4,1
*NSET, NSET=NDISP
2
*ELEMENT,TYPE=S4R,ELSET=ONE
1,1,2,3,4
*SHELL SECTION,ELSET=ONE,MATERIAL=HYPER
1.0
*MATERIAL,NAME=HYPER
*HYPERELASTIC, MARLOW
*UNIAXIAL TEST DATA,smooth
  0.0, 0.0
  1.5506E5,  0.1338
  2.4367E5,  0.2675
  3.1013E5,  0.3567
  4.2089E5,  0.6242
  5.3165E5,  0.8917
  5.9810E5,  1.1592
  6.8671E5,  1.4268
  8.8608E5,  2.0510
 10.6329E5,  2.5860
 12.4051E5,  3.0318
 16.1709E5,  3.7898
 19.9367E5,  4.3694
 23.4810E5,  4.8153
 27.4684E5,  5.1720
 31.0127E5,  5.4395
 34.5570E5,  5.7070
 38.3228E5,  5.9299
 42.0886E5,  6.0637
 45.6329E5,  6.1975
 49.3987E5,  6.3312
 53.1646E5,  6.4650
 56.9304E5,  6.5541
 64.2405E5,  6.6433
*EQUATION
 2,
 FACE4,1,1,NDISP,1,-1
*BOUNDARY
ALL,4,6
1,1,3
2,2,3
3,2,3
4,1
4,3
*STEP,NLGEOM
*STATIC
 0.01, 1.0, , 0.05
*CLOAD
 NDISP,1,65.0E5
**BOUNDARY,OP=NEW
**FACE4,1,1,5.
*EL PRINT,F=10
S, 
E, 
*NODE PRINT,F=10
U,
RF
*OUTPUT,FIELD,FREQ=10
*ELEMENT OUTPUT
S,
E
*OUTPUT, HISTORY
*NODE OUTPUT,NSET=NDISP
CF1,U1
*END STEP