Fixed a bug where the polytropic index wasn't being updated during the switch to protostellar stage 4.
Added more logic to EvolveProtostellar.F90 to contrain the initialization of the polytropic index to be between 1.5 and 3.0