void buildsystemlist() {
  if (Online) {NSYSTEMS = 5;} else {NSYSTEMS = 7;}
  systems = new dynamicalSystem[NSYSTEMS];
  initialSystem = 0;

  int i=0;
  systems[i] = new NPZsystem();
  systems[i].showVarnames = true;
  systems[i].name = "NPZ (nutrients + phytoplankton + zooplankton)";
  systems[i].shortname = "simple";

  i++;
  systems[i] = new allo_upwelling_system1D(6,6,1,20);
  systems[i].name = "allometric, six classes of P and Z";
  systems[i].shortname = "size spectrum";
  systems[i].showFluxes = true;
  systems[i].params[systems[i].find("xpmax")].initial = 35;
  systems[i].showHistory = true;
  systems[i].showVarnames = true;
  
  i++;
  systems[i] = new allo_upwelling_system1D(20,20,1,10);
  systems[i].name = "allometric, 20 classes of P and Z";
  systems[i].shortname = "smooth size spectrum";
  systems[i].showFluxes = true;
  systems[i].params[systems[i].find("xpmax")].initial = 35;
  systems[i].showHistory = false;
  systems[i].showVarnames = true;

  i++;
  systems[i] = new allo_upwelling_system1D(6,6,4,30);
  systems[i].name = "six classes of P and Z; cross-shelf upwelling circulation";
  systems[i].shortname = "cross-shelf conveyor";
  systems[i].showFluxes = true;
  systems[i].params[systems[i].find("b")].initial = 40;
  systems[i].params[systems[i].find("supp_intensity")].initial = 3;
  systems[i].showVarnames = true;
  
  i++;
  systems[i] = new myNPZD();
  systems[i].showVarnames = true;
  systems[i].showHistory = false;
  systems[i].name = "NPZD from Banas et al. 2009 (JGR, submitted)";
  systems[i].shortname = "NPZD, Banas et al 09";
  statevar s = systems[i].vars[systems[i].find("1")];
  s.display = true;
  s.xc = -0.8;
  s.yc = -0.8;
  s.shortname = "1 µM nitrogen";
  systems[i].showVarDiagnostics = true;
  systems[i].params[systems[i].find("attSW")].frozen=true;
  systems[i].params[systems[i].find("attChl")].frozen=true;

  if (Online) return;
  
  i++;
  systems[i] = new romsBiophysicalSlice(new myNPZD(), "analytical");
  systems[i].name = "NPZD (dinoflagellate) + idealized upwelling";
  systems[i].shortname = "dino, 2-D";
  systems[i].dot_stdsize *= 5;
  
  i++;
  dynamicalSystem DS = new myNPZD();
  DS.params[DS.find("Imax")].initial = 0.6;
  DS.params[DS.find("mz")].initial = 0.1;
  DS.shortname += " slow";
  systems[i] = new romsBiophysicalSlice(DS, "analytical");
  systems[i].name = "NPZD (copepod) + idealized upwelling";
  systems[i].shortname = "copepod, 2-D";
  systems[i].dot_stdsize *= 5; 
}




/* system list for Coastal Stories lab spr 2006 (v1.2)
void buildsystemlist() {
  NSYSTEMS = 3;
  systems = new dynamicalSystem[NSYSTEMS];
  initialSystem = 0;
  
  systems[0] = new NPZsystem();
  systems[0].name = "NPZ (nutrients + phytoplankton + zooplankton)";
  systems[0].shortname = "simple";
  
  systems[1] = new allo_upwelling_system1D(6,6,1,20);
  systems[1].name = "six classes of P, Z (# = size in microns)";
  systems[1].shortname = "less simple";
  systems[1].showFluxes = true;
  systems[1].params[systems[1].find("xpmax")].initial = 35;
  systems[1].showHistory = true;
    
  systems[2] = new allo_upwelling_system1D(6,6,4,30);
  systems[2].name = "six classes of P, Z; upwelling circulation";
  systems[2].shortname = "complicated";
  systems[2].showFluxes = true;
  systems[2].params[systems[2].find("b")].initial = 40;
  systems[2].params[systems[2].find("supp_intensity")].initial = 3;
}
*/




/* system list for Ocean Sciences 2006
void buildsystemlist() {
  NSYSTEMS = 6;
  systems = new dynamicalSystem[NSYSTEMS];
  
  systems[0] = new NNPZDsystem();
  
  systems[1] = new allo_upwelling_system(2,2,5,30);
  systems[1].name = "two size classes / upwelling";
  systems[1].shortname = "upwelling, 2 P,Z";
  systems[1].params[systems[1].find("xpmin")].initial = 1;
  systems[1].params[systems[1].find("xpmax")].initial = 10;
  systems[1].params[systems[1].find("xzmin")].initial = 3;
  systems[1].params[systems[1].find("xzmax")].initial = 150;
  
  systems[2] = new allo_upwelling_system(4,4,5,30);
  systems[2].name = "four size classes / upwelling";
  systems[2].shortname = "upwelling, 4 P,Z";
  
  systems[3] = new allo_upwelling_system(3,3,1,30);
  systems[3].shortname = "one box, 3 P,Z";

  systems[4] = new allo_upwelling_system(10,10,1,30);
  systems[4].shortname = "one box, 10 P,Z";

  systems[5] = new allo_upwelling_system(20,20,1,30);
  systems[5].shortname = "one box, 20 P,Z";

    
  initialSystem = 2;
} */