38 long int nloop_pres_oscil;
39 double TemperatureInitial;
75 " ConvPresTempEdenIoniz1 entered, will call ConvIoniz to initialize\n");
98 " ConvPresTempEdenIoniz1 ConvIoniz found following converged: Pres;%c, Eden;%c, Temp;%c, Ion:%c Pops:%c\n",
118 "\n ConvPresTempEdenIoniz1 entering main pressure loop.\n");
156 hden_chng_old = hden_chng;
167 enum{DEBUG_LOC=
false};
171 fprintf(
ioQQQ,
"%li\t%.2e\t%.2e\t%.2e\n",
188 if( ( ( hden_chng*hden_chng_old < 0. ) ) &&
196 nloop_pres_oscil = loop;
198 dP_chng_factor =
MAX2(0.1 , dP_chng_factor * 0.75 );
219 " ConvPresTempEdenIoniz1 %.2f l:%3li nH:%.4e ne:%.4e PCurnt:%.4e PCorct:%.4e err:%6.3f%% dP/dn:%.2e Te:%.4e Osc:%c\n",
239 if( loop - nloop_pres_oscil > 4 )
244 if( loop == LoopMax && !lgPresOscil )
246 LoopMax =
MIN2( 100 , LoopMax*2 );
271 # if !defined(NDEBUG)