17 fprintf(
ioQQQ,
" PROBLEM DISASTER More than one density command was entered.\n" );
25 fprintf(
ioQQQ,
" DISASTER The density MUST be entered with this command. STOP\n" );
37 " DISASTER - the log of the entered hydrogen density is %.3f - too extreme for this processor.\n",
41 " DISASTER - the log of the largest hydrogen density this processor can do is %.3f.\n",
45 " DISASTER - the log of the smallest hydrogen density this processor can do is %.3f.\n",
47 fprintf(
ioQQQ,
" Sorry.\n" );
57 fprintf(
ioQQQ,
"This density is too high. This version of Cloudy does not permit densities greater than %e cm-3.\n",
MAX_DENSITY );
63 fprintf(
ioQQQ,
" PROBLEM DISASTER Hydrogen density must be > 0.\n" );
85 fprintf(
ioQQQ,
" The column density MUST be set if the col den option is to be used.\n" );
88 strcpy(
dense.chDenseLaw,
"POWC" );
90 else if( p.
nMatch(
"DEPT") )
96 fprintf(
ioQQQ,
" The scale depth MUST be set if the depth option is to be used.\n" );
99 strcpy(
dense.chDenseLaw,
"POWD" );
104 strcpy(
dense.chDenseLaw,
"POWR" );
118 if( strcmp(
dense.chDenseLaw ,
"CDEN") == 0 ||
119 strcmp(
dense.chDenseLaw ,
"CPRE") == 0 ||
120 strcmp(
dense.chDenseLaw ,
"WIND") == 0 ||
121 strcmp(
dense.chDenseLaw ,
"DYNA") == 0
129 else if( strcmp(
dense.chDenseLaw,
"POWR") == 0 )
137 else if( strcmp(
dense.chDenseLaw,
"POWC") == 0 )
139 strcpy(
optimize.chVarFmt[
optimize.nparm],
"HDEN=%f LOG, power=%f, column=%f" );
146 else if( strcmp(
dense.chDenseLaw,
"POWD") == 0 )
148 strcpy(
optimize.chVarFmt[
optimize.nparm],
"HDEN=%f LOG, power=%f, depth=%f" );
157 fprintf(
ioQQQ,
" Internal error in HDEN\n" );
#define DEBUG_ENTRY(funcname)
bool nMatch(const char *chKey) const
void ParseHDEN(Parser &p)