6 integer ib, ibk, na, iflgs, ic_type, ic_counter
30 g(n1,ib) = 1.0 / (r1+1.0e-20)
38 g(n2,ib) = 1.0 / (r3+1.0e-20)
48 if (ic_type .eq. vcapacitor1)
then
49 vinitial = value_init(ic_counter)
52 ecapsource = ecapsource + 0.5*
c(n1,ib)*vinitial*vinitial
53 else if (ic_type .eq. cinductor)
then
54 cinitial = value_init(ic_counter)
56 zirn(n1,ib) = cinitial
57 eindsource = eindsource + 0.5*zlr(n1,ib)*cinitial*cinitial
58 else if (ic_type .eq. vcapacitor3)
then
59 vinitial = value_init(ic_counter)
62 ecapsource = ecapsource + 0.5*
c(n2,ib)*vinitial*vinitial
subroutine setup_pisection(ib, ibk, na, iflgs, ic_type, ic_counter)
c *****************************************************************************c Various format statements for read_screamer_data output c To get these into made format to be characters c for each line c