*Heading
 This input deck is run to make sure that *temperature
 and *pressure sress (and similarly *massflowrate, *field)
 are properly processed in pre/impl/ntfile.f when applied to
 different number of nodes (V58_2722).
*Preprint, echo=NO, model=YES, history=YES, contact=NO
**
** units kN and KPa and m
*parameter
 ultI   = 0.3200e+3
 ultII  = 1.0e+10
** ultII  = 0.3200e+3
 Emod   = 8.5E7
** GIc    = 0.33e2
** GIIc   = 0.80e2
 GIc    = 28.
 GIIc   = 28.
** GIIc   = 0.00001
 eta    = 2.284
 thick  = 0.01325e-1
 width  = 1.0
**
** First run coupled_cod.inp, then run massdiff_mod.inp
**
*Node, Nset=Nall
 1, 0,0
 2, 1,0
 3, 1,1
 4, 0,1
 5, 2,0
 6, 2,1
 7, 1.1,0
 8, 1.1,1
*Element, Type=Cpe4t, Elset=Eall
 1, 1,2,3,4
 2, 7,5,6,8
*Element, Type=COH2D4T, Elset=Coh
 3, 3,2,7,8
*Nset,Nset=Nsome
 1,2,5
*Nset, Nset=Fix
 1,4
*Nset, Nset=Pull
 5,6
*Solid Section, Elset=Eall, Material=Steel
*COHESIVE SECTION,ELSET=Coh,MATERIAL=MAT1,RESPONSE=TRACTION SEPARATION, THICKNESS=SPECIFIED
0.1,1.0
*Material, Name=Steel
*Elastic
 2.11e+05,
*Conductivity
 4.349e-02,
**
*MATERIAL,NAME=MAT1
*ELASTIC,TYPE=TRACTION
<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY
<GIc>
*GAP CONDUCTANCE
1.0,0.0
1.1,1.0
*GAP RADIATION
1.0,1.0
1.,0.
0.,1.0
*PHYSICAL CONSTANTS, ABSOLUTE ZERO=-460., STEFAN=1.E-8   
**PHYSICAL CONSTANTS, ABSOLUTE ZERO=-273.16,STEFANB=5.6697E-8
**
*Initial Conditions, Type=Temperature
 Nall, 100.0
**Step, Inc=200, unsymm=yes, nlgeom=yes
*Step, Inc=200, nlgeom=yes
*Coupled Temperature-Displacement,Steady State
 .5,1.
*Boundary
 Fix, 1,2
 Fix, 11,11,200.0
 Pull, 1,1, 0.5
 Pull,11,11,150.0
*Output, Field, time interval=0.25
*Element Output
 S,
*Output, Field, time interval=0.25
*Node Output
 Nt,u
*Element Output, Position=Nodes
 Sinv,sdeg
*Node Print, Freq=0
*El Print, Freq=0
*Node File, Freq=1, Nset=Nsome
 Nt,
*El File, Freq=1, Position=Averaged At Nodes
 Sinv,
*End Step