oceanbox/Oceanbox.FvcomKit
6
0
Fork 0
Code Issues Pull Requests 2 Actions Packages Projects Releases 80 Wiki Activity

Merge branch 'feat/omega-block' into 'main'

Browse Source 
feat: add read omega block

See merge request oceanbox/Oceanbox.FvcomKit!31
This commit was merged in pull request #32.
This commit is contained in:
Stig Rune Jensen
2025-02-12 15:27:54 +00:00
parent 02de20d50b 9abd9d45bc
commit 86e1307d79
1 changed files with 11 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/Fvcom.fs
View File

@@ -165,6 +165,17 @@ let readOmega (ds: DataSet) t l =
Log.Error $"{e}" Log.Error $"{e}"
Array.empty 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 = let readOmegas (ds: DataSet) t l es =
try try
let n = ds.Dimensions["node"].Length let n = ds.Dimensions["node"].Length