27 if(
thermal.ConstGrainTemp <= 10. )
33 p.
NoNumb(
"grain temperature");
37 else if( p.
nMatch(
"TEMP") )
40 thermal.lgTemperatureConstant =
true;
41 thermal.lgTemperatureConstantCommandParsed =
true;
50 else if( p.
nMatch(
" KEV") )
67 fprintf(
ioQQQ,
" PROBLEM temperature entered as a log but is too large "\
68 "for this processor. I am interpreting it as the linear temperature.\n");
74 thermal.ConstTemp *= convert_to_Kelvin;
80 fprintf(
ioQQQ,
" PROBLEM Te too low, reset to %g K.\n",
86 fprintf(
ioQQQ,
" PROBLEM Te too high, reset to %g K.\n",
113 else if( p.
nMatch(
"DENS") )
116 strcpy(
dense.chDenseLaw,
"CDEN" );
118 pressure.lgPres_radiation_ON =
false;
119 pressure.lgPres_magnetic_ON =
false;
123 else if( p.
nMatch(
"PRES") )
126 strcpy(
dense.chDenseLaw,
"CPRE" );
135 dense.lgDenseInitConstant =
false;
141 dense.lgDenseInitConstant =
true;
149 dense.lgDenseInitConstant =
false;
150 dense.lgPressureVaryTime =
true;
157 fprintf(
ioQQQ,
" PROBLEM the constant pressure time command requires"
158 " a positive timescale.\n");
168 fprintf(
ioQQQ,
" PROBLEM the constant pressure time command requires"
169 " two numbers, the timescale for the variation and an index.\n");
178 pressure.lgPres_radiation_ON =
false;
184 pressure.lgPres_magnetic_ON =
false;
195 pressure.lgPres_radiation_ON =
true;
225 pressure.lgPressureInitialSpecified =
true;
229 p.
NoNumb(
"initial pressure" );
232 pressure.PressureInitialSpecified = pow(10.,
pressure.PressureInitialSpecified)*
236 pressure.lgPressureInitialSpecified =
false;
245 strcpy(
optimize.chVarFmt[
optimize.nparm],
"CONSTANT PRESSURE SET %f LOG" );
246 if( !
dense.lgDenseInitConstant )
266 fprintf(
ioQQQ,
" The keyword should be TEMPerature, DENSity, GAS or PRESsure, sorry.\n" );