finish
/clear,nostart
!
resume,new-tf_pf_mcwf\tf_pf_wg_mc,db	! metric system 
!
/prep7
mp,rsvx,7,0.1673e-7 
!
! enter currents for coils
!
itf=  4424	! tf current 2T HB at t=0.197 sec
ipf1= -7522	! pf1 current 2T HB at t=0.197 sec
ipf2= -7522	! pf2 current 2T HB at t=0.197 sec
ipf3= -4305	! pf3 current 2T HB at t=0.197 sec
ipf4= -7589	! pf4 current 2T HB at t=0.197 sec
ipf5= 2349	! pf5 current 2T HB at t=0.197 sec
ipf6= 7300	! pf6 current 2T HB at t=0.197 sec
!
imda=35075    	! 2T high bata at t=0.197
imdb=34852   
imdc=31783
!
ipla=-209732		! plasma current
!
/input,new-tf_pf_mcwf\itg_model-1r,txt
!
! ******************* form model ********************
!
! resume,new-tf_pf_mcwf\itg_file_new_itg-1,db !(from itg_model-1r.txt)
! resume,new-tf_pf_mcwf\itg_file_new_itg-1HB197,db  !(name changed)
!
! revise base supports, G10 alpha and modify TF bracket
!
/input,new-tf_pf_mcwf\itg-rbe3-support,txt
mp,ALPX,114,1.50e-05
mp,ALPX,109,1.50e-05
/input,new-tf_pf_mcwf\itg-modify-1,txt
!
!  enter trim load cases from iota 0.19s
!
/input,new-tf_pf_mcwf\itg_tca_iota019s,txt 	! Fx, Fy and Fz
! /input,new-tf_pf_mcwf\itg_tcam_iota019s,txt 	! reverse loads
!
!
/sol		 For DL, EM and cool-down
!
! fdele,all,all
!
wf=1.2427
acel,,,9.806*wf
!
csw=-3865*4.448/3/12		! centerstack/3/12
cmsel,s,csn
f,all,fz,csw
!
cmsel,s,vvn
vvwt=6800+12000			! vacuum vessel wt.
f,all,fz,-vvwt*4.448/32
!
asel,s,loc,z,1.435+z1,1.436+z1	! cryostat
asel,u,area,,16	
! pcro=-3500*4.448/0.7044
pcro=3500*4.448/0.6161
sfa,all,,pres,pcro
!
! /input,new-tf_pf_mcwf\itg_tc_iota019s,txt  ! Fz only
!
allsel
eqslv,sparse
! eqslv,pcg,1.0e-6	! 1.0e-4
nsubst,2,3,1		! nsbstp, nsbmx, nsbmn
solcon,on
nlgeom,off
!
tref,0
esel,s,mat,,102			! Modular coils
nsle
bf,all,temp,-23.256
!
esel,s,mat,,109		! PF shims(G10/G11)
nsle
bf,all,temp,-27.7
!
esel,s,mat,,114		! TF shims(G10/G11)
nsle
bf,all,temp,-27.7
!
esel,s,mat,,103,104		! TF and PF coils
nsle
bf,all,temp,-42.719
!
! esel,s,mat,,105		! TF wedges
! nsle
! bf,all,temp,0
!
esel,s,mat,,11,12		! wing bags
esel,a,mat,,15,17,2
esel,a,mat,,18,21,3
esel,a,mat,,69,70
esel,u,type,,78,79
nsle
bf,all,temp,-117.87
!
allsel
solve
save,new-tf_pf_mcwf\h9-emdlco-2T-HB197s-tca,db
!
finish
!
!
/post1
!
rsys,0
prrsol
plnsol,s,eqv,0