diff --git a/src/Fvcom.fs b/src/Fvcom.fs index 060bc4d..d756d86 100644 --- a/src/Fvcom.fs +++ b/src/Fvcom.fs @@ -165,6 +165,17 @@ let readOmega (ds: DataSet) t l = Log.Error $"{e}" Array.empty +let readOmegaBlock (ds: DataSet) t = + try + let n = ds.Dimensions["node"].Length + let l = ds.Dimensions["siglev"].Length + let u = ds[ "omega" ].GetData([| t; 0; 0 |], [| 1; l; n |]) :?> single [,,] + u[0, *, *] + with + | e -> + Log.Error $"{e}" + Array2D.zeroCreate 0 0 + let readOmegas (ds: DataSet) t l es = try let n = ds.Dimensions["node"].Length