void printSpecies(Model_t *m, FILE *f)
{
int i, j;
Species_t *s;
Compartment_t *c;
fprintf(f, "\n");
fprintf(f, "# Initial Conditions for Species and Compartments:\n");
for ( i=0; i<Model_getNumCompartments(m); i++ ) {
if ( i== 0 ) fprintf(f, "# Compartments:\n");
c = Model_getCompartment(m,i);
if(Compartment_isSetId(c))
fprintf(f, "%s ", Compartment_getId(c));
if(Compartment_isSetName(c))
fprintf(f, "(%s) ", Compartment_getName(c));
if ( Compartment_isSetVolume(c) )
fprintf(f, "= %g; ", Compartment_getSize(c));
fprintf(f, "%s", Compartment_getConstant(c) ? "" : "variable; ");
if(Compartment_isSetOutside(c))
fprintf(f, "outside %s; ", Compartment_getOutside(c));
/* fprintf(f, "\n"); */
fprintf(f, "dimensions %d; ", Compartment_getSpatialDimensions(c));
if(Compartment_isSetUnits(c))
fprintf(f, "[%s]; ", Compartment_getUnits(c));
fprintf(f, "\n");
fprintf(f, "# Species concentrations in `compartment' %s\n",
Compartment_getId(c));
for(j=0;j<Model_getNumSpecies(m);j++){
s = Model_getSpecies(m,j);
if(strcmp(Species_getCompartment(s), Compartment_getId(c))==0){
fprintf(f, "%s ", Species_getId(s));
if(Species_isSetName(s))
fprintf(f, "(%s) ", Species_getName(s));
if ( Species_isSetInitialAmount(s) )
fprintf(f, "= %g/%g; ",
Species_getInitialAmount(s),
Compartment_getSize(c));
else if ( Species_isSetInitialConcentration(s) )
fprintf(f, "= %g; ", Species_getInitialConcentration(s));
else
fprintf(f, "# no initial value;");
fprintf(f, "%s", Species_getBoundaryCondition(s) ? "boundary;" : "");
fprintf(f, "%s", Species_getConstant(s) ? "constant;" : "");
if(Species_isSetCharge(s))
fprintf(f, "charge = %d; ", Species_getCharge(s));
fprintf(f, "\n");
}
}
fprintf(f, "\n");
}
}
END_TEST
START_TEST (test_Species_setCharge1)
{
int i = Species_setCharge(C, 2);
fail_unless( i == LIBSBML_OPERATION_SUCCESS );
fail_unless( Species_isSetCharge(C) );
fail_unless( Species_getCharge(C) == 2);
i = Species_unsetCharge(C);
fail_unless( i == LIBSBML_OPERATION_SUCCESS );
fail_unless( !Species_isSetCharge(C) );
}
/*
setValues: the user can enter a species name and
change its initial condition (amount or concentration)
*/
static void setValues(Model_t *m) {
char *species;
char *newIA;
char *newIC;
Species_t *s;
printf("Please enter the id of the species to change: ");
species = get_line(stdin);
species = util_trim(species);
if ( (s = Model_getSpeciesById(m,species) ) ) {
printf("\n");
printf("Id: %s\n", Species_getId(s));
if ( Species_isSetName(s) ) {
printf("Name: %s\n", Species_getName(s));
}
if ( Species_isSetInitialAmount(s) ) {
printf("Initial Amount: %g", Species_getInitialAmount(s));
}
else if (Species_isSetInitialConcentration(s) ) {
printf("Initial Concentration: %g", Species_getInitialConcentration(s));
}
if ( Species_getHasOnlySubstanceUnits(s) ) {
if ( Species_isSetSubstanceUnits(s) ) {
printf("%s ", Species_getSubstanceUnits(s));
}
} else {
if ( Species_isSetSubstanceUnits(s) ) {
printf("%s ", Species_getSubstanceUnits(s));
}
if ( Species_isSetSpatialSizeUnits(s) ) {
printf("%s%s", "/", Species_getSpatialSizeUnits(s));
}
}
if ( Species_getHasOnlySubstanceUnits(s) ) {
printf(" (has only substance units)");
}
printf("\n");
if ( Species_isSetCharge(s) ) {
printf("Charge: %-10d", Species_getCharge(s));
}
printf("\n");
printf("%s ", Species_getBoundaryCondition(s) ?
"Species is a Boundary\n" : "\n");
printf("%s ", Species_getConstant(s) ?
"Species is set constant" : "\n");
printf("\n");
if ( Species_isSetInitialAmount(s) ) {
printf("Please enter new initial Amount: ");
newIA = get_line(stdin);
newIA = util_trim(newIA);
Species_setInitialAmount(s, (float) atof(newIA));
}
else if ( Species_isSetInitialConcentration(s) ) {
printf("Please enter new initial Concentration: ");
newIC = get_line(stdin);
newIC = util_trim(newIC);
Species_setInitialConcentration(s, (float) atof(newIC));
}
}
else {
printf("%s not found.\n", species);
}
}
