31 parameter(mu = 2.0e-7)
49 parms(1) = shellradius(1)
51 parms(3) = shellmass(1)
54 parms(6) = (mu*temp(1))*log(parms(1)/parms(2))
61 parms(13) = temp(3) + parms(1)
62 parms(14) = parms(5)*log(parms(13)/parms(1))
71 shellcurr(max_shells+1)=0.0
subroutine nshellparm(parms, nparms)