class WaterColumn extends Slice {

  WaterColumn(Ecosystem eco, int K, float Htot) {
    name = "1D water column";
    shortname = "watercol";
    createDomain(1, K, 1, Htot);
    linkToEcosystem(eco);
    
    internalTimestep = 0.001;

    Statevar mass = addVar("mass","mass","",Ncells);
    getVar("mass").setInitial(1);
    mass.advects = true;
    mass.mixes = true; 
    
    allocateStorageForPhysics();

    addParam("surface PAR","PAR0","W/m2", 150, 0, 300);
    addParam("vertical mixing","logkappa","m2/s", -5, -6, -4);
    
    addParam("D sinking rate", "wsinkD", "m/day", 10, 0, 100);
    addSinking("D", getParamVal("wsinkD"), 0);

    startArchive(1);
  }
  

  void updatePrivateParams() {
    float kappa0 = pow(10,getParamVal("logkappa")) * 86400;
    for (int m=0; m<M; m++) for (int k=0; k<K; k++) kappa[k][m] = kappa0;
    sinkingRate[getVarIndex("D")] = getParamVal("wsinkD");
    super.updatePrivateParams();
  }

}