/*
 * test.java
 */

import com.comsol.model.*;
import com.comsol.model.util.*;
import java.io.*;
import java.util.*;


/** Model exported on Aug 14 2014, 23:58 by COMSOL 4.3.2.189. */
public class add_inlets {

  public static void main(String[] args) {
    run();
  }

  public static Model run() {
      Model model;
      try
	  {
	      model=ModelUtil.load("Model","/data/fopiroot/fopiroot_dev/padresponse_COMSOL/test2.mph");
	  }
      catch(IOException e)
	  {
	      model=ModelUtil.create("Model"); 
	      return model;
	  }

      double timesigma=1.24e-9;
      int ntimesteps=10;
      double timestep=timesigma*4/ntimesteps;
      int maxnElec=1000;

      int irel=0;
      String inlet_name;
      double reltime=0;


      int nElec=0;
      for (irel=0;irel<ntimesteps;irel++)
	  {
	      inlet_name="inlet_"+Integer.toString(irel);

	      nElec=(int) ( gauss(reltime,timesigma*4/2,timesigma) * maxnElec );
	      System.out.println("adding an inlet: time="+Double.toString(reltime)+" and nElec="+Integer.toString(nElec));

	      model.physics("cpt").feature().create(inlet_name, "Inlet", 1);
	      model.physics("cpt").feature(inlet_name).selection().set(new int[]{20});
	      model.physics("cpt").feature(inlet_name).set("rt", reltime);
	      model.physics("cpt").feature(inlet_name).set("InitialPosition", "Density");
	      model.physics("cpt").feature(inlet_name).set("N", nElec);
	      model.physics("cpt").feature(inlet_name).set("dpro", "gp1(x)*rn3(x,t)");
	      model.physics("cpt").feature(inlet_name)
		  .set("v0", new String[][]{{"rn_vxy(particleindex,t)"},
					    {"rn_vz(particleindex,t)"}, 
					    {"0"}});
	      reltime+=timestep;
	  }

    return model;
  }
    public static double gauss(double x, double mean, double sig)
    {
	double retval=0;
	double norm=1./(sig*Math.sqrt(2*Math.PI));
	double exponent=-0.5*( Math.pow( (x-mean)/sig , 2 ) );
	return Math.exp(exponent);
    }

}