code: public static String getSKSAT(String[] soilParas) {
          if (soilParas != null && soilParas.length >= 3) {
              if (soilParas.length >= 4) {
                  return divide(calcSatBulk(soilParas[0], soilParas[1], soilParas[2], divide(soilParas[3], "100")), "10", 3);
              } else {
                  return divide(calcSatMatric(soilParas[0], soilParas[1], soilParas[2]), "10", 3);
              }
          } else {
              return null;
          }
      }
docs: For calculating SKSAT
  
  @param soilParas should include 1. Sand weight percentage by layer
  ([0,100]%), 2. Clay weight percentage by layer ([0,100]%), 3. Organic
  matter weight percentage by layer ([0,100]%), (= SLOC * 1.72), 4. Grave
  weight percentage by layer ([0,100]%)
  @return Saturated hydraulic conductivity, cm/h