mccg.cpp

pic 模拟程序！面向对象

                                     discardDumpFileNGDDataFlag.getValue() ) );

	else if (gas.getValue() == (ostring)"Xe") 
          pMCC->addPackage( new XeMCC(pressure.getValue(), 
		temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
		*SRG->advisorp->get_speciesList(), SRG->get_dtBase(), ionzFraction_i.getValue(),
		ecxFactor.getValue(), icxFactor.getValue(), region, x1MinMKS.getValue(),
                x1MaxMKS.getValue(), x2MinMKS.getValue(), x2MaxMKS.getValue(),
                delayTime.getValue(), stopTime.getValue(), analyticF.getValue(),
                                      discardDumpFileNGDDataFlag.getValue() ) );

	else if (gas.getValue() == (ostring)"H") 
          pMCC->addPackage(new HMCC(pressure.getValue(), 
		temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
		*SRG->advisorp->get_speciesList(), SRG->get_dtBase(), ionzFraction_i.getValue(),
		ecxFactor.getValue(), icxFactor.getValue(), region, x1MinMKS.getValue(),
                x1MaxMKS.getValue(), x2MinMKS.getValue(), x2MaxMKS.getValue(),
                delayTime.getValue(), stopTime.getValue(), analyticF.getValue(),
                                    discardDumpFileNGDDataFlag.getValue() ) );

	else if (gas.getValue() == (ostring)"Li") 
          pMCC->addPackage(new LiMCC(pressure.getValue(), 
		temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
		*SRG->advisorp->get_speciesList(), SRG->get_dtBase(), ionzFraction_i.getValue(),
		ecxFactor.getValue(), icxFactor.getValue(), relativisticMCC.getValue(), region,
                x1MinMKS.getValue(), x1MaxMKS.getValue(), x2MinMKS.getValue(), x2MaxMKS.getValue(),
                delayTime.getValue(), stopTime.getValue(), analyticF.getValue(),
                                     discardDumpFileNGDDataFlag.getValue() ) );

	else if (gas.getValue() == (ostring)"N") 
          pMCC->addPackage(new NMCC(pressure.getValue(), 
		temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
		*SRG->advisorp->get_speciesList(), SRG->get_dtBase(), ionzFraction_i.getValue(),
		ecxFactor.getValue(), icxFactor.getValue(), relativisticMCC.getValue(), region,
                x1MinMKS.getValue(), x1MaxMKS.getValue(), x2MinMKS.getValue(), x2MaxMKS.getValue(),
                delayTime.getValue(), stopTime.getValue(), analyticF.getValue(),
                                    discardDumpFileNGDDataFlag.getValue() ) );
        else if (gas.getValue() == (ostring)"He") {
	  pMCC->addPackage(new HeMCC(pressure.getValue(), 
		temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
		*SRG->advisorp->get_speciesList(), SRG->get_dtBase(), ionzFraction_i.getValue(),
		ecxFactor.getValue(), icxFactor.getValue(), region, x1MinMKS.getValue(),
                x1MaxMKS.getValue(), x2MinMKS.getValue(), x2MaxMKS.getValue(),
                delayTime.getValue(), stopTime.getValue(), analyticF.getValue(),
                                    discardDumpFileNGDDataFlag.getValue() ) );
        }
	else cout << "Unrecognized MCC package requested: " << gas.getValue();
  }

  if (tunnelingIonizationFlag.getValue()) {
        if (gas.getValue() ==(ostring) "H") { 
          try{
            pMCTI->addPackage(new HMCTI(pressure.getValue(), 
                                      temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
                                      *SRG->advisorp->get_speciesList(), SRG->get_dtBase(), 
                                      ionzFraction_i.getValue(), ecxFactor.getValue(), 
                                      icxFactor.getValue(), region, x1MinMKS.getValue(),
                                      x1MaxMKS.getValue(), x2MinMKS.getValue(), 
                                      x2MaxMKS.getValue(), delayTime.getValue(), 
                                      stopTime.getValue(), analyticF.getValue(),
                                      TI_numMacroParticlesPerCell.getValue(), 
                                      ETIPolarizationFlag.getValue(),
                                      EfieldFrequency.getValue(),
                                      discardDumpFileNGDDataFlag.getValue()));
          }
          catch(Oops& oops){
            oops.prepend("MCCParams::CreateCounterPart: Error: \n");//SpatialRegionGroup::CreateCounterPart
            throw oops;
          }

        }
        else if (gas.getValue() == (ostring)"He") {
          try {
            pMCTI->addPackage(new HeMCTI(pressure.getValue(), 
                                       temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
                                       iSpeciesPlusPlusPtr, 
                                       *SRG->advisorp->get_speciesList(), SRG->get_dtBase(), 
                                       ionzFraction_i.getValue(), ecxFactor.getValue(), 
                                       icxFactor.getValue(), region, x1MinMKS.getValue(),
                                       x1MaxMKS.getValue(), x2MinMKS.getValue(), 
                                       x2MaxMKS.getValue(), delayTime.getValue(), 
                                       stopTime.getValue(), analyticF.getValue(),
                                       TI_numMacroParticlesPerCell.getValue(), 
                                       ETIPolarizationFlag.getValue(),
                                       EfieldFrequency.getValue(),
                                       discardDumpFileNGDDataFlag.getValue()));
          }
          catch(Oops& oops){
            oops.prepend("MCCParams::CreateCounterPart: Error: \n"); //OK
            throw oops;
          }
        }
        else if (gas.getValue() == (ostring)"Li") {
          try {
          cout << "Creating a LiMCTI instance" << endl; 
          pMCTI->addPackage(new LiMCTI(pressure.getValue(), 
                                       temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
                                       *SRG->advisorp->get_speciesList(), SRG->get_dtBase(), 
                                       ionzFraction_i.getValue(), ecxFactor.getValue(), 
                                       icxFactor.getValue(), region, x1MinMKS.getValue(),
                                       x1MaxMKS.getValue(), x2MinMKS.getValue(), 
                                       x2MaxMKS.getValue(), delayTime.getValue(), 
                                       stopTime.getValue(), analyticF.getValue(),
                                       TI_numMacroParticlesPerCell.getValue(), 
                                       ETIPolarizationFlag.getValue(),
                                       EfieldFrequency.getValue(),
                                       discardDumpFileNGDDataFlag.getValue()));
          }
          catch(Oops& oops){
            oops.prepend("MCCParams::CreateCounterPart: Error: \n"); //OK
            throw oops;
          }
        }
        else if (gas.getValue() == (ostring)"Cs") {
          cout << "Creating a CsMCTI instance" << endl;  
          try{
            pMCTI->addPackage(new CsMCTI(pressure.getValue(), 
                                       temperature.getValue(), eSpeciesPtr, iSpeciesPtr, 
                                       iSpeciesPlusPlusPtr, 
                                       *SRG->advisorp->get_speciesList(), SRG->get_dtBase(), 
                                       ionzFraction_i.getValue(), ecxFactor.getValue(), 
                                       icxFactor.getValue(), region, x1MinMKS.getValue(),
                                       x1MaxMKS.getValue(), x2MinMKS.getValue(), 
                                       x2MaxMKS.getValue(), delayTime.getValue(), 
                                       stopTime.getValue(), analyticF.getValue(),
                                       TI_numMacroParticlesPerCell.getValue(), 
                                       ETIPolarizationFlag.getValue(),
                                       EfieldFrequency.getValue(),
                                       discardDumpFileNGDDataFlag.getValue()));
          }
          catch(Oops& oops){
            oops.prepend("MCCParams::CreateCounterPart: Error:\n"); //OK
            throw oops; 
          }
        }
	else if (gas.getValue() ==(ostring) "Ar") 
          // pMCTI->addPackage( new ArMCTI(...) );
          ;
	else if (gas.getValue() ==(ostring) "Ne") 
          // pMCTI->addPackage( new NeMCTI(...) );
          ;
	else if (gas.getValue() == (ostring)"Xe") 
          // pMCTI->addPackage( new XeMCTI(...) );
          ;
	else if (gas.getValue() == (ostring)"N") 
          // pMCTI->addPackage( new NMCTI(...) );
          ;
	else cout << "Unrecognized MCTI package requested: " << gas.getValue();
  }
  return;
}

void MCCParams::checkRules()
{
	char buffer[512];
   BOOL gasTest = FALSE; 
   if (gas.getValue() ==(ostring) "Ar") gasTest = TRUE;
	else if (gas.getValue() ==(ostring) "Ne") gasTest = TRUE;
	else if (gas.getValue() == (ostring)"Xe") gasTest = TRUE;
	else if (gas.getValue() == (ostring)"H") gasTest = TRUE;
	else if (gas.getValue() == (ostring)"He") gasTest = TRUE;
	else if (gas.getValue() == (ostring)"Li") gasTest = TRUE;
	else if (gas.getValue() == (ostring)"N") gasTest = TRUE;
        else if (gas.getValue() == (ostring)"Cs") gasTest = TRUE;
	if (!gasTest) {
	  sprintf(buffer, "unrecognized gas package: %s", gas.getValue().c_str());
	  ruleMessages.add(new ostring(buffer)); 
	}
	ostring name = eSpecies.getValue();
	if (name == (ostring)"NULL") ruleMessages.add(new ostring("eSpecies is required!"));
	else if (!SRG->KnownSpecies(name))
	{
		sprintf(buffer, "Unrecognized speciesName: %s", name());
		ruleMessages.add(new ostring(buffer));
	}
	name = iSpecies.getValue();
	if (name == (ostring)"NULL") ruleMessages.add(new ostring("iSpecies is required!"));
	else if (!SRG->KnownSpecies(name))
	{
		sprintf(buffer, "Unrecognized speciesName: %s", name());
		ruleMessages.add(new ostring(buffer));
	}
}