feat: read bed stress data

src/Fvcom.fs

@@ -99,6 +99,16 @@ let readUV (ds: DataSet) t l =
         Log.Error $"{e}"
         Array.empty
 
+let readBedStress (ds: DataSet) t =
+    try
+        let n = ds.Dimensions["nele"].Length
+        let tau = ds[ "tauc" ].GetData([| t; 0 |], [| 1; n |]) :?> single [,]
+        tau[0, *]
+    with
+    | e ->
+        Log.Error $"{e}"
+        Array.empty
+
 let readUV' (ds: DataSet) t l =
     readUV ds t l
     |> Array.collect (fun (x, y) -> [| x; y |])
@@ -448,10 +458,10 @@ let getGrid (ds: DataSet) : FvcomGrid =
             |> Array.map float
         let nv = ds[ "nv" ].GetData() :?> int [,]
 
-        let h = readBathymetry ds
-        let siglay = Siglay.readSiglay ds
-        let siglay_c = Siglay.readSiglayAtCenter ds
-        let siglev = ds[ "siglev" ].GetData() :?> single [,]
+        let h = Array.zeroCreate 1 //readBathymetry ds
+        let siglay = Array2D.zeroCreate 1 1 //Siglay.readSiglay ds
+        let siglay_c = Array2D.zeroCreate 1 1 //Siglay.readSiglayAtCenter ds
+        let siglev = Array2D.zeroCreate 1 1 //ds[ "siglev" ].GetData() :?> single [,]
         let elem =
             Array.zip3 nv[0, *] nv[1, *] nv[2, *]
             |> Array.map (fun (a, b, c) -> a - 1, b - 1, c - 1)