using System; using System.Runtime.InteropServices; using NUnit.Framework; using General.tests; using System.IO; using System.Linq; using System.Threading; //To compile and run this test you add NUnit to your visual studio solution and // enable the build in configuration manager // todo: enable pre-build events to generate the version number namespace gridgeom.Tests { public class ODugrid { //Constructor loads the library static ODugrid() { string filename = TestHelper.GetLibraryPath(GridGeomLibWrapper.LibDetails.LIB_NAME); _gridgeom_libptr = TestHelper.LoadLibrary(filename); //we should chek the pointer is not null Assert.That(_gridgeom_libptr, Is.Not.Null); //load netcdf for reading in meshgeom TestHelper.SetSharedPath(GridGeomLibWrapper.LibDetails.NETCDF_DEP); filename = TestHelper.GetLibraryPath(GridGeomLibWrapper.LibDetails.NETCDF_LIB_NAME); _netcdf_libptr = TestHelper.LoadLibrary(filename); Assert.That(_netcdf_libptr, Is.Not.Null); } //pointer to the loaded dll public static IntPtr _gridgeom_libptr; //pointer to the loaded dll public static IntPtr _netcdf_libptr; private void getMeshid(int ioncid, ref int meshid, int meshType, ref IoNetcdfLibWrapper wrapper) { // get the number meshes int numMeshes = -1; int ierr = wrapper.ionc_get_number_of_meshes(ref ioncid, ref meshType, ref numMeshes); Assert.That(ierr, Is.EqualTo(0)); Assert.That(numMeshes, Is.EqualTo(1)); // get the mesh id IntPtr c_meshids = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * numMeshes); ierr = wrapper.ionc_ug_get_mesh_ids(ref ioncid, ref meshType, ref c_meshids, ref numMeshes); Assert.That(ierr, Is.EqualTo(0)); int[] rc_meshids = new int[numMeshes]; Marshal.Copy(c_meshids, rc_meshids, 0, numMeshes); meshid = rc_meshids[0]; Marshal.FreeCoTaskMem(c_meshids); } [Test] [NUnit.Framework.Category("ougridTests")] public void convert1DUgridToXY() { //dimension info int nmeshpoints = 10; int nbranches = 3; int ngeopoints = 9; int jsferic = 0; //branches info double[] branchlengths = { 4.0, 3.0, 3.0 }; int[] nbranchgeometrynodes = { 3, 3, 3 }; //geometry info double[] geopointsX = { 1.0, 3.0, 5.0, 5.0, 5.0, 5.0, 5.0, 7.0, 8.0 }; double[] geopointsY = { 4.0, 4.0, 4.0, 1.0, 2.0, 4.0, 4.0, 4.0, 4.0 }; //mesh geometry int[] branchids = { 1, 1, 1, 1, 2, 2, 2, 3, 3, 3 }; double[] branchoffsets = { 0.0, 2.0, 3.0, 4.0, 0.0, 1.5, 3.0, 0.0, 1.5, 3.0 }; // Create the netcdf files double[] meshXCoords = { 1.0, 3.0, 4.0, 5.0, 5.0, 5.0, 5.0, 5.0, 6.5, 8.0 }; double[] meshYCoords = { 4.0, 4.0, 4.0, 4.0, 1.0, 2.5, 4.0, 4.0, 4.0, 4.0 }; IntPtr c_branchids = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nmeshpoints); IntPtr c_branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_geopointsX = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * ngeopoints); IntPtr c_geopointsY = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * ngeopoints); IntPtr c_nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_branchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nbranches); IntPtr c_meshXCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_meshYCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); try { var wrapper = new GridGeomLibWrapper(); Marshal.Copy(branchids, 0, c_branchids, nmeshpoints); Marshal.Copy(branchoffsets, 0, c_branchoffsets, nmeshpoints); Marshal.Copy(geopointsX, 0, c_geopointsX, ngeopoints); Marshal.Copy(geopointsY, 0, c_geopointsY, ngeopoints); Marshal.Copy(nbranchgeometrynodes, 0, c_nbranchgeometrynodes, nbranches); Marshal.Copy(branchlengths, 0, c_branchlengths, nbranches); //call the function and assert for validity int ierr = wrapper.ggeo_get_xy_coordinates(ref c_branchids, ref c_branchoffsets, ref c_geopointsX, ref c_geopointsY, ref c_nbranchgeometrynodes, ref c_branchlengths, ref c_meshXCoords, ref c_meshYCoords, ref jsferic, ref nbranches, ref ngeopoints, ref nmeshpoints); Assert.That(ierr, Is.EqualTo(0)); double[] rc_meshXCoords = new double[nmeshpoints]; double[] rc_meshYCoords = new double[nmeshpoints]; Marshal.Copy(c_meshXCoords, rc_meshXCoords, 0, nmeshpoints); Marshal.Copy(c_meshYCoords, rc_meshYCoords, 0, nmeshpoints); // check conversion is correct for (int i = 0; i < nmeshpoints; i++) { Assert.That(rc_meshXCoords[i], Is.EqualTo(meshXCoords[i])); Assert.That(rc_meshYCoords[i], Is.EqualTo(meshYCoords[i])); } } finally { Marshal.FreeCoTaskMem(c_branchids); Marshal.FreeCoTaskMem(c_branchoffsets); Marshal.FreeCoTaskMem(c_geopointsX); Marshal.FreeCoTaskMem(c_geopointsY); Marshal.FreeCoTaskMem(c_nbranchgeometrynodes); Marshal.FreeCoTaskMem(c_branchlengths); Marshal.FreeCoTaskMem(c_meshXCoords); Marshal.FreeCoTaskMem(c_meshYCoords); } } [Test] [NUnit.Framework.Category("ougridTests")] public void LargeNumberOfPoints() { //dimension info int s_nmeshpoints = 6; int s_nbranches = 2; int s_ngeopoints = 4; int repetition = 100000; int jsferic = 0; //branches info double[] s_branchlengths = { 2.0, 2.0 }; int[] s_nbranchgeometrynodes = { 2, 2 }; //geometry info double[] s_geopointsX = { 0.0, 2.0, 2.0, 2.0 }; double[] s_geopointsY = { 0.0, 0.0, 0.0, -2.0 }; //mesh geometry int[] s_branchids = { 1, 1, 1, 2, 2, 2 }; double[] s_branchoffsets = { 0.0, 1.0, 2.0, 0.0, 1.0, 2.0 }; //mesh coordinates double[] s_meshXCoords = { 0.0, 1.0, 2.0, 2.0, 2.0, 2.0 }; double[] s_meshYCoords = { 0.0, 0.0, 0.0, 0.0, -1.0, -2.0 }; //repeat small structure int nmeshpoints = s_nmeshpoints * repetition; int nbranches = s_nbranches * repetition; int ngeopoints = s_ngeopoints * repetition; double[] branchlengths = new double[nbranches]; int[] nbranchgeometrynodes = new int[nbranches]; double[] geopointsX = new double[ngeopoints]; double[] geopointsY = new double[ngeopoints]; int[] branchids = new int[nmeshpoints]; double[] branchoffsets = new double[nmeshpoints]; double[] meshXCoords = new double[nmeshpoints]; double[] meshYCoords = new double[nmeshpoints]; //generate a large number of point, by repeating the structure above int brid = 0; int geoid = 0; int mid = 0; for (int i = 0; i < repetition; i++) { for (int j = 0; j < s_nbranches; j++) { branchlengths[brid] = s_branchlengths[j]; nbranchgeometrynodes[brid] = s_nbranchgeometrynodes[j]; brid = brid + 1; } for (int j = 0; j < s_ngeopoints; j++) { geopointsX[geoid] = s_geopointsX[j] + i; //add positive offset geopointsY[geoid] = s_geopointsY[j] - i; //add negative offset geoid = geoid + 1; } for (int j = 0; j < s_nmeshpoints; j++) { meshXCoords[mid] = s_meshXCoords[j] + i; //add positive offset meshYCoords[mid] = s_meshYCoords[j] - i; //add negative offset branchids[mid] = s_branchids[j] + i * 2; //add branch ids branchoffsets[mid] = s_branchoffsets[j]; mid = mid + 1; } } IntPtr c_branchids = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nmeshpoints); IntPtr c_branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_geopointsX = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * ngeopoints); IntPtr c_geopointsY = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * ngeopoints); IntPtr c_nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_branchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nbranches); IntPtr c_meshXCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_meshYCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); try { var wrapper = new GridGeomLibWrapper(); Marshal.Copy(branchids, 0, c_branchids, nmeshpoints); Marshal.Copy(branchoffsets, 0, c_branchoffsets, nmeshpoints); Marshal.Copy(geopointsX, 0, c_geopointsX, ngeopoints); Marshal.Copy(geopointsY, 0, c_geopointsY, ngeopoints); Marshal.Copy(nbranchgeometrynodes, 0, c_nbranchgeometrynodes, nbranches); Marshal.Copy(branchlengths, 0, c_branchlengths, nbranches); //call the function int ierr = wrapper.ggeo_get_xy_coordinates(ref c_branchids, ref c_branchoffsets, ref c_geopointsX, ref c_geopointsY, ref c_nbranchgeometrynodes, ref c_branchlengths, ref c_meshXCoords, ref c_meshYCoords, ref jsferic, ref nbranches, ref ngeopoints, ref nmeshpoints); Assert.That(ierr, Is.EqualTo(0)); double[] rc_meshXCoords = new double[nmeshpoints]; double[] rc_meshYCoords = new double[nmeshpoints]; Marshal.Copy(c_meshXCoords, rc_meshXCoords, 0, nmeshpoints); Marshal.Copy(c_meshYCoords, rc_meshYCoords, 0, nmeshpoints); // check conversion is correct for (int i = 0; i < nmeshpoints; i++) { Assert.That(rc_meshXCoords[i], Is.EqualTo(meshXCoords[i])); Assert.That(rc_meshYCoords[i], Is.EqualTo(meshYCoords[i])); } } finally { Marshal.FreeCoTaskMem(c_branchids); Marshal.FreeCoTaskMem(c_branchoffsets); Marshal.FreeCoTaskMem(c_geopointsX); Marshal.FreeCoTaskMem(c_geopointsY); Marshal.FreeCoTaskMem(c_nbranchgeometrynodes); Marshal.FreeCoTaskMem(c_branchlengths); Marshal.FreeCoTaskMem(c_meshXCoords); Marshal.FreeCoTaskMem(c_meshYCoords); } } /// /// Read a netcdf file and populate meshgeom datastructure /// [Test] [NUnit.Framework.Category("ougridTests")] public void createLinksFrom1d2dFile() { //mesh2d int twoddim = 2; int twodnumnode = 452; int twodnumedge = 825; int twodnumface = 375; int twodmaxnumfacenodes = 4; int twodnumlayer = -1; int twodlayertype = -1; int jsferic = 0; //mesh1d int oneddim = 1; int onednumnode = 25; int onednumedge = 24; int onednumface = 0; int onedmaxnumfacenodes = 0; int onednumlayer = -1; int onedlayertype = -1; int nt_nbranches = 1; int nt_ngeometry = 25; int numnodes = 10; int numedge = 1; string c_path = TestHelper.TestFilesDirectoryPath() + @"\river1_full_net.nc"; Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); int meshid = -1; int meshType = 2; int l_networkid = 0; // invalid network, not really intrested in getting the networ var meshtwoddim = new meshgeomdim(); getMeshid(ioncid, ref meshid, meshType, ref wrapperNetcdf); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref l_networkid, ref meshtwoddim); Assert.That(ierr, Is.EqualTo(0)); Assert.That(meshtwoddim.dim, Is.EqualTo(twoddim)); Assert.That(meshtwoddim.numnode, Is.EqualTo(twodnumnode)); Assert.That(meshtwoddim.numedge, Is.EqualTo(twodnumedge)); Assert.That(meshtwoddim.numface, Is.EqualTo(twodnumface)); Assert.That(meshtwoddim.maxnumfacenodes, Is.EqualTo(twodmaxnumfacenodes)); Assert.That(meshtwoddim.numlayer, Is.EqualTo(twodnumlayer)); Assert.That(meshtwoddim.layertype, Is.EqualTo(twodlayertype)); //You need to know in advance the number of mesh points var meshtwod = new meshgeom(); meshtwod.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshtwoddim.numedge * 2); //var gridGeomWrapper = new GridGeomLibWrapper(); bool includeArrays = true; int start_index = 1; //the base index of the arrays ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref l_networkid, ref meshtwod, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); double[] rc_twodnodex = new double[twodnumnode]; double[] rc_twodnodey = new double[twodnumnode]; double[] rc_twodnodez = new double[twodnumnode]; Marshal.Copy(meshtwod.nodex, rc_twodnodex, 0, twodnumnode); Marshal.Copy(meshtwod.nodey, rc_twodnodey, 0, twodnumnode); Marshal.Copy(meshtwod.nodez, rc_twodnodez, 0, twodnumnode); // mesh1d var meshoneddim = new meshgeomdim(); meshType = 1; int meshonedid = -1; getMeshid(ioncid, ref meshonedid, meshType, ref wrapperNetcdf); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshonedid, ref l_networkid, ref meshoneddim); Assert.That(ierr, Is.EqualTo(0)); Assert.That(meshoneddim.dim, Is.EqualTo(oneddim)); Assert.That(meshoneddim.numnode, Is.EqualTo(onednumnode)); Assert.That(meshoneddim.numedge, Is.EqualTo(onednumedge)); Assert.That(meshoneddim.numlayer, Is.EqualTo(onednumlayer)); Assert.That(meshoneddim.ngeometry, Is.EqualTo(nt_ngeometry)); Assert.That(meshoneddim.layertype, Is.EqualTo(onedlayertype)); var meshoned = new meshgeom(); //mesh variables meshoned.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numedge * 2); //network variables meshoned.nnodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.nnodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.branchidx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numnode); meshoned.nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); meshoned.branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.ngeopointx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.ngeopointy = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.nbranchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); meshoned.nedge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); meshoned.nbranchorder = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); //need to produce a file with coordinate_space = string variable ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshonedid, ref meshonedid, ref meshoned, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); double[] rc_onednodex = new double[onednumnode]; double[] rc_onednodey = new double[onednumnode]; double[] rc_onednodez = new double[onednumnode]; Marshal.Copy(meshoned.nodex, rc_onednodex, 0, onednumnode); Marshal.Copy(meshoned.nodey, rc_onednodey, 0, onednumnode); Marshal.Copy(meshoned.nodez, rc_onednodez, 0, onednumnode); //gridwrapper var wrapperGridgeom = new GridGeomLibWrapper(); ierr = wrapperGridgeom.ggeo_get_xy_coordinates(ref meshoned.branchidx, ref meshoned.branchoffsets, ref meshoned.ngeopointx, ref meshoned.ngeopointy, ref meshoned.nbranchgeometrynodes, ref meshoned.nbranchlengths, ref meshoned.nodex, ref meshoned.nodey, ref jsferic, ref meshoneddim.nbranches, ref meshoneddim.ngeometry, ref meshoneddim.numnode); Assert.That(ierr, Is.EqualTo(0)); //ggeo_convert to fill in kn matrix, so we can call make1D2Dinternalnetlinks_dll int startIndex = 0; ierr = wrapperGridgeom.ggeo_convert(ref meshoned, ref meshoneddim, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); //convert ugrid to xy coordinates Assert.That(ierr, Is.EqualTo(0)); //call make1d2dlinks, no argument needed (all in memory) // empty polygon int c_npl = 0; IntPtr c_xpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_ypl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_zpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); int c_nOneDMask = 0; IntPtr c_oneDmask = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * 0); ierr = wrapperGridgeom.ggeo_make1D2Dinternalnetlinks(ref c_npl, ref c_xpl, ref c_ypl, ref c_zpl, ref c_nOneDMask, ref c_oneDmask); Assert.That(ierr, Is.EqualTo(0)); //deallocate memory of gridgeom ierr = wrapperGridgeom.ggeo_deallocate(); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperNetcdf.ionc_close(ref ioncid); Assert.That(ierr, Is.EqualTo(0)); //free arrays //2d Marshal.FreeCoTaskMem(meshtwod.nodex); Marshal.FreeCoTaskMem(meshtwod.nodey); Marshal.FreeCoTaskMem(meshtwod.nodez); Marshal.FreeCoTaskMem(meshtwod.edge_nodes); //1d //mesh variables Marshal.FreeCoTaskMem(meshoned.nodex); Marshal.FreeCoTaskMem(meshoned.nodey); Marshal.FreeCoTaskMem(meshoned.nodez); //Marshal.FreeCoTaskMem(meshoned.edge_nodes); //network variables //Marshal.FreeCoTaskMem(meshoned.nnodex); Marshal.FreeCoTaskMem(meshoned.nnodey); Marshal.FreeCoTaskMem(meshoned.branchidx); Marshal.FreeCoTaskMem(meshoned.nbranchgeometrynodes); Marshal.FreeCoTaskMem(meshoned.branchoffsets); Marshal.FreeCoTaskMem(meshoned.ngeopointx); Marshal.FreeCoTaskMem(meshoned.ngeopointy); Marshal.FreeCoTaskMem(meshoned.nbranchlengths); Marshal.FreeCoTaskMem(meshoned.nedge_nodes); Marshal.FreeCoTaskMem(meshoned.nbranchorder); } /// /// In this test we read a 2d grid from file and we provide the 1d discretization points. /// We emulate the case where Delta Shell has a 2d file opened and the user creates a 1d mesh and wants to generate the links. /// 1D discretization and links are saved in memory and written afterwards. /// [Test] [Category("ougridTests")] public void CreateLinksFrom2DFile() { //mesh2d int twoddim = 2; int twodnumnode = 16; int twodnumedge = 24; int twodnumface = 9; int twodmaxnumfacenodes = 4; int twodnumlayer = -1; int twodlayertype = -1; int startIndex = 1; // the indexes in the array are zero based //mesh1d //discretization points information int nmeshpoints = 4; int nbranches = 1; int[] branchids = { 1, 1, 1, 1 }; double[] meshXCoords = { -6, 5, 23, 34 }; double[] meshYCoords = { 22, 16, 16, 7 }; double[] branchoffset = { 0, 10, 20, 100 }; /// important are the first and last offset double[] branchlength = { 100 }; int[] sourcenodeid = { 1 }; int[] targetnodeid = { 2 }; //links int[] arrayfrom = { 2, 8, 7 }; int[] arrayto = { 2, 3, 4 }; //1. open the file with the 2d mesh string c_path = TestHelper.TestFilesDirectoryPath() + @"\2d_ugrid_net.nc"; Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); //2. get the 2d mesh id int meshid = 1; ierr = wrapperNetcdf.ionc_get_2d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); //3. get the dimensions of the 2d mesh var meshtwoddim = new meshgeomdim(); int l_networkid = 0; ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref l_networkid, ref meshtwoddim); Assert.That(ierr, Is.EqualTo(0)); Assert.That(meshtwoddim.dim, Is.EqualTo(twoddim)); Assert.That(meshtwoddim.numnode, Is.EqualTo(twodnumnode)); Assert.That(meshtwoddim.numedge, Is.EqualTo(twodnumedge)); Assert.That(meshtwoddim.numface, Is.EqualTo(twodnumface)); Assert.That(meshtwoddim.maxnumfacenodes, Is.EqualTo(twodmaxnumfacenodes)); Assert.That(meshtwoddim.numlayer, Is.EqualTo(twodnumlayer)); Assert.That(meshtwoddim.layertype, Is.EqualTo(twodlayertype)); //4. allocate the arrays in meshgeom for storing the 2d mesh coordinates, edge_nodes var meshtwod = new meshgeom(); meshtwod.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshtwoddim.numedge * 2); //5. get the meshgeom arrays bool includeArrays = true; int start_index = 1; int l_network1d = 0; ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref l_network1d, ref meshtwod, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); double[] rc_twodnodex = new double[twodnumnode]; double[] rc_twodnodey = new double[twodnumnode]; double[] rc_twodnodez = new double[twodnumnode]; Marshal.Copy(meshtwod.nodex, rc_twodnodex, 0, twodnumnode); Marshal.Copy(meshtwod.nodey, rc_twodnodey, 0, twodnumnode); Marshal.Copy(meshtwod.nodez, rc_twodnodez, 0, twodnumnode); //6. allocate the 1d arrays for storing the 1d coordinates and edge_nodes IntPtr c_meshXCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_meshYCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_branchids = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nmeshpoints); IntPtr c_sourcenodeid = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_targetnodeid = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_branchoffset = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_branchlength = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nbranches); Marshal.Copy(branchids, 0, c_branchids, nmeshpoints); Marshal.Copy(meshXCoords, 0, c_meshXCoords, nmeshpoints); Marshal.Copy(meshYCoords, 0, c_meshYCoords, nmeshpoints); Marshal.Copy(sourcenodeid, 0, c_sourcenodeid, nbranches); Marshal.Copy(targetnodeid, 0, c_targetnodeid, nbranches); Marshal.Copy(branchoffset, 0, c_branchoffset, nmeshpoints); Marshal.Copy(branchlength, 0, c_branchlength, nbranches); //7. fill kn (Herman datastructure) for creating the links var wrapperGridgeom = new GridGeomLibWrapper(); ierr = wrapperGridgeom.ggeo_convert_1d_arrays(ref c_meshXCoords, ref c_meshYCoords, ref c_branchoffset, ref c_branchlength, ref c_branchids, ref c_sourcenodeid, ref c_targetnodeid, ref nbranches, ref nmeshpoints, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); //9. make the links int c_npl = 0; IntPtr c_xpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_ypl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_zpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); int c_nOneDMask = 0; IntPtr c_oneDmask = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * 0); ierr = wrapperGridgeom.ggeo_make1D2Dinternalnetlinks(ref c_npl, ref c_xpl, ref c_ypl, ref c_zpl, ref c_nOneDMask, ref c_oneDmask); Assert.That(ierr, Is.EqualTo(0)); //10. check if we can call ggeo_get_links_count two times: we need to (1) deallocate the memory stored in network_ggeo_data and (2) reload the 1d and 2d arrays in memory //ierr = wrapperGridgeom.ggeo_deallocate(); //Assert.That(ierr, Is.EqualTo(0)); //ierr = wrapperGridgeom.ggeo_convert_1d_arrays(ref c_meshXCoords, ref c_meshYCoords, ref c_branchoffset, ref c_branchlength, ref c_branchids, ref c_sourcenodeid, ref c_targetnodeid, ref nbranches, ref nmeshpoints, ref startIndex); //Assert.That(ierr, Is.EqualTo(0)); //ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim); //Assert.That(ierr, Is.EqualTo(0)); //ierr = wrapperGridgeom.ggeo_make1D2Dinternalnetlinks(); //Assert.That(ierr, Is.EqualTo(0)); //10. get the number of links int n1d2dlinks = 0; int linkType = 3; ierr = wrapperGridgeom.ggeo_get_links_count(ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); //11. get the links: arrayfrom = 2d cell index, arrayto = 1d node index IntPtr c_arrayfrom = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * n1d2dlinks); //2d cell number IntPtr c_arrayto = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * n1d2dlinks); //1d node ierr = wrapperGridgeom.ggeo_get_links(ref c_arrayfrom, ref c_arrayto, ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); int[] rc_arrayfrom = new int[n1d2dlinks]; int[] rc_arrayto = new int[n1d2dlinks]; Marshal.Copy(c_arrayfrom, rc_arrayfrom, 0, n1d2dlinks); Marshal.Copy(c_arrayto, rc_arrayto, 0, n1d2dlinks); for (int i = 0; i < n1d2dlinks; i++) { Assert.That(rc_arrayfrom[i], Is.EqualTo(arrayfrom[i])); Assert.That(rc_arrayto[i], Is.EqualTo(arrayto[i])); } //for writing the links look io_netcdf ionc_def_mesh_contact, ionc_put_mesh_contact ierr = wrapperNetcdf.ionc_close(ref ioncid); Assert.That(ierr, Is.EqualTo(0)); //Free 2d arrays Marshal.FreeCoTaskMem(meshtwod.nodex); Marshal.FreeCoTaskMem(meshtwod.nodey); Marshal.FreeCoTaskMem(meshtwod.nodez); Marshal.FreeCoTaskMem(meshtwod.edge_nodes); //Free 1d arrays Marshal.FreeCoTaskMem(c_meshXCoords); Marshal.FreeCoTaskMem(c_meshYCoords); Marshal.FreeCoTaskMem(c_branchids); //Free from and to arrays describing the links Marshal.FreeCoTaskMem(c_arrayfrom); Marshal.FreeCoTaskMem(c_arrayto); } /// /// In this test we read a 2d grid from file and we provide the 1d discretization points. /// We emulate the case where Delta Shell has a 2d file opened and the user creates a 1d mesh and wants to generate the links. /// 1D discretization and links are saved in memory and written afterwards. /// [Test] [NUnit.Framework.Category("ougridTests")] public void createLinksFrom2dFileThreeBranches() { //mesh2d int twoddim = 2; int twodnumnode = 16; int twodnumedge = 24; int twodnumface = 9; int twodmaxnumfacenodes = 4; int twodnumlayer = -1; int twodlayertype = -1; int startIndex = 1; //mesh1d //discretization points information int nmeshpoints = 9; int nbranches = 3; int[] branchids = { 1, 1, 1, 1, 2, 2, 2, 3, 3 }; double[] meshXCoords = { 7.5, 12.5, 17.5, 22.5, 22.5, 22.5, 22.5, 17.5, 12.5 }; double[] meshYCoords = { 22.5, 22.5, 22.5, 22.5, 17.5, 12.5, 7.5, 12.5, 17.5 }; double[] branchoffset = { 0, 1, 2, 10, 1, 2, 10, 1, 2 }; /// the actual values of the offset are not important double[] branchlength = { 10, 10, 10 }; int[] sourcenodeid = { 1, 2, 3 }; int[] targetnodeid = { 2, 3, 1 }; //1. open the file with the 2d mesh string c_path = TestHelper.TestFilesDirectoryPath() + @"\2d_ugrid_net.nc"; Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); //2. get the 2d mesh id int meshid = 1; ierr = wrapperNetcdf.ionc_get_2d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); //3. get the dimensions of the 2d mesh var meshtwoddim = new meshgeomdim(); int l_networkid = 0; //sets an invalid network ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref l_networkid, ref meshtwoddim); Assert.That(ierr, Is.EqualTo(0)); Assert.That(meshtwoddim.dim, Is.EqualTo(twoddim)); Assert.That(meshtwoddim.numnode, Is.EqualTo(twodnumnode)); Assert.That(meshtwoddim.numedge, Is.EqualTo(twodnumedge)); Assert.That(meshtwoddim.numface, Is.EqualTo(twodnumface)); Assert.That(meshtwoddim.maxnumfacenodes, Is.EqualTo(twodmaxnumfacenodes)); Assert.That(meshtwoddim.numlayer, Is.EqualTo(twodnumlayer)); Assert.That(meshtwoddim.layertype, Is.EqualTo(twodlayertype)); //4. allocate the arrays in meshgeom for storing the 2d mesh coordinates, edge_nodes var meshtwod = new meshgeom(); meshtwod.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * twodnumnode); meshtwod.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshtwoddim.numedge * 2); //5. get the meshgeom arrays bool includeArrays = true; int start_index = 1; int l_network1d = 0; ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref l_network1d, ref meshtwod, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); double[] rc_twodnodex = new double[twodnumnode]; double[] rc_twodnodey = new double[twodnumnode]; double[] rc_twodnodez = new double[twodnumnode]; Marshal.Copy(meshtwod.nodex, rc_twodnodex, 0, twodnumnode); Marshal.Copy(meshtwod.nodey, rc_twodnodey, 0, twodnumnode); Marshal.Copy(meshtwod.nodez, rc_twodnodez, 0, twodnumnode); //6. allocate the 1d arrays for storing the 1d coordinates and edge_nodes IntPtr c_meshXCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_meshYCoords = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_branchids = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nmeshpoints); IntPtr c_sourcenodeid = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_targetnodeid = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * nbranches); IntPtr c_branchoffset = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nmeshpoints); IntPtr c_branchlength = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * nbranches); Marshal.Copy(branchids, 0, c_branchids, nmeshpoints); Marshal.Copy(meshXCoords, 0, c_meshXCoords, nmeshpoints); Marshal.Copy(meshYCoords, 0, c_meshYCoords, nmeshpoints); Marshal.Copy(sourcenodeid, 0, c_sourcenodeid, nbranches); Marshal.Copy(targetnodeid, 0, c_targetnodeid, nbranches); Marshal.Copy(branchoffset, 0, c_branchoffset, nmeshpoints); Marshal.Copy(branchlength, 0, c_branchlength, nbranches); //7. fill kn (Herman datastructure) for creating the links var wrapperGridgeom = new GridGeomLibWrapper(); ierr = wrapperGridgeom.ggeo_convert_1d_arrays(ref c_meshXCoords, ref c_meshYCoords, ref c_branchoffset, ref c_branchlength, ref c_branchids, ref c_sourcenodeid, ref c_targetnodeid, ref nbranches, ref nmeshpoints, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); //9. make the links int c_npl = 0; IntPtr c_xpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_ypl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); IntPtr c_zpl = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * 0); int c_nOneDMask = 0; IntPtr c_oneDmask = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * 0); ierr = wrapperGridgeom.ggeo_make1D2Dinternalnetlinks(ref c_npl, ref c_xpl, ref c_ypl, ref c_zpl, ref c_nOneDMask, ref c_oneDmask); Assert.That(ierr, Is.EqualTo(0)); //10. get the number of links int n1d2dlinks = 0; int linkType = 3; ierr = wrapperGridgeom.ggeo_get_links_count(ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); //11. get the links: arrayfrom = 2d cell index, arrayto = 1d node index IntPtr c_arrayfrom = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * n1d2dlinks); //2d cell number IntPtr c_arrayto = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * n1d2dlinks); //1d node ierr = wrapperGridgeom.ggeo_get_links(ref c_arrayfrom, ref c_arrayto, ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); int[] rc_arrayfrom = new int[n1d2dlinks]; int[] rc_arrayto = new int[n1d2dlinks]; Marshal.Copy(c_arrayfrom, rc_arrayfrom, 0, n1d2dlinks); Marshal.Copy(c_arrayto, rc_arrayto, 0, n1d2dlinks); //for writing the links look io_netcdf ionc_def_mesh_contact, ionc_put_mesh_contact ierr = wrapperNetcdf.ionc_close(ref ioncid); Assert.That(ierr, Is.EqualTo(0)); //Free 2d arrays Marshal.FreeCoTaskMem(meshtwod.nodex); Marshal.FreeCoTaskMem(meshtwod.nodey); Marshal.FreeCoTaskMem(meshtwod.nodez); Marshal.FreeCoTaskMem(meshtwod.edge_nodes); //Free 1d arrays Marshal.FreeCoTaskMem(c_meshXCoords); Marshal.FreeCoTaskMem(c_meshYCoords); Marshal.FreeCoTaskMem(c_branchids); //Free from and to arrays describing the links Marshal.FreeCoTaskMem(c_arrayfrom); Marshal.FreeCoTaskMem(c_arrayto); } [Test] [TestCase("default")] //[TestCase(@"pathToNetFile")] [NUnit.Framework.Category("findcells")] public void testFindCells(string path) { //int stackSize = 1024 * 1024 * 32; //LC: in C# the default stack size is 1 MB, increase it to something larger for this test! //Thread th = new Thread(() => //{ //1. open the file with the 2d mesh string c_path; if (path == "default") { c_path = TestHelper.TestFilesDirectoryPath() + @"\2d_ugrid_net.nc"; } else { c_path = path; } Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); //2. get the 2d mesh id int meshid = 1; ierr = wrapperNetcdf.ionc_get_2d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); //3. get the dimensions of the 2d mesh var meshDimIn = new meshgeomdim(); int l_networkid = 0; //sets an invalid network ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref l_networkid, ref meshDimIn); Assert.That(ierr, Is.EqualTo(0)); //4. allocate meshgeom to pass to find_cells var meshIn = new meshgeom(); meshIn.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshDimIn.numnode); meshIn.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshDimIn.numnode); meshIn.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshDimIn.numnode); meshIn.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshDimIn.numedge * 2); //5. get the meshgeom arrays from file bool includeArrays = true; int start_index = 1; //retrive as 1 based int l_network1d = -1; //no network present ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref l_network1d, ref meshIn, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); //5. declare but do not allocate meshgeom. it will be allocated by gridgeom (fortran) var meshOut = new meshgeom(); var meshDimOut = new meshgeomdim(); meshOut.face_nodes = IntPtr.Zero; meshOut.facex = IntPtr.Zero; meshOut.facey = IntPtr.Zero; //6. call find cells int startIndex = 1; // provide 1 based (read from netcdf), return 1 based var wrapperGridgeom = new GridGeomLibWrapper(); ierr = wrapperGridgeom.ggeo_find_cells(ref meshDimIn, ref meshIn, ref meshDimOut, ref meshOut, ref startIndex); Assert.That(ierr, Is.EqualTo(0)); //7. copy face_nodes to array int[] face_nodes = new int[meshDimOut.maxnumfacenodes * meshDimOut.numface]; Marshal.Copy(meshOut.face_nodes, face_nodes, 0, meshDimOut.maxnumfacenodes * meshDimOut.numface); //8. deallocate memory allocated by fortran ierr = wrapperGridgeom.ggeo_meshgeom_destructor(ref meshDimOut, ref meshOut); Assert.That(ierr, Is.EqualTo(0)); //9. deallocate memory allocated by c# Marshal.FreeCoTaskMem(meshIn.nodex); Marshal.FreeCoTaskMem(meshIn.nodey); Marshal.FreeCoTaskMem(meshIn.nodez); Marshal.FreeCoTaskMem(meshIn.edge_nodes); //}, stackSize); //th.Start(); //th.Join(); } [Test] [TestCase("default")] //[TestCase(@"pathToNetFile")] [NUnit.Framework.Category("make1D2DEmbeddedLinks")] public void make1D2DEmbeddedLinks(string path) { //1. open the file with the 1d network, 1d mesh and 2d grid. string c_path = TestHelper.TestFilesDirectoryPath() + @"\1d2dForEmbedded.nc"; Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); var wrapperGridgeom = new GridGeomLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); int meshid = -1; int networkid = -1; bool includeArrays = true; int start_index = 1; int jsferic = 0; //--------------------------------------------------------------------------------------// //2a. get network and mesh ids ierr = wrapperNetcdf.ionc_get_1d_network_id(ref ioncid, ref networkid); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperNetcdf.ionc_get_1d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); //2b. get dimensions var meshoneddim = new meshgeomdim(); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref networkid, ref meshoneddim); Assert.That(ierr, Is.EqualTo(0)); //2c. allocate var meshoned = new meshgeom(); //mesh variables meshoned.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); //network variables meshoned.nnodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.nnodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.branchidx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numnode); meshoned.nedge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nnodes * 2); meshoned.nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); meshoned.branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.ngeopointx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.ngeopointy = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.nbranchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); meshoned.nbranchorder = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); //2d. get 1d meshgeom array ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref networkid, ref meshoned, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_count_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoneddim.nbranches, ref meshoneddim.numnode, ref meshoneddim.numedge, ref start_index); Assert.That(ierr, Is.EqualTo(0)); meshoned.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numedge * 2); ierr = wrapperGridgeom.ggeo_create_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoned.edge_nodes, ref meshoneddim.nbranches, ref meshoneddim.numnode, ref meshoneddim.numedge, ref start_index); Assert.That(ierr, Is.EqualTo(0)); //Read the created edge nodes int[] edge_nodes = new int[meshoneddim.numedge * 2]; Marshal.Copy(meshoned.edge_nodes, edge_nodes, 0, meshoneddim.numedge * 2); Assert.That(edge_nodes.Contains(0),Is.EqualTo(false)); //--------------------------------------------------------------------------------------// //3a. get mesh ids ierr = wrapperNetcdf.ionc_get_2d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); networkid = -1; //3b. get dimensions var meshtwoddim = new meshgeomdim(); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref networkid, ref meshtwoddim); Assert.That(ierr, Is.EqualTo(0)); //3c. allocate var meshtwod = new meshgeom(); meshtwod.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshtwoddim.numedge * 2); //3d. get meshgeom ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref networkid, ref meshtwod, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); //8. close file ierr = wrapperNetcdf.ionc_close(ref ioncid); Assert.That(ierr, Is.EqualTo(0)); //--------------------------------------------------------------------------------------// //4. fill kn (Herman datastructure) for creating the links ierr = wrapperGridgeom.ggeo_get_xy_coordinates(ref meshoned.branchidx, ref meshoned.branchoffsets, ref meshoned.ngeopointx, ref meshoned.ngeopointy, ref meshoned.nbranchgeometrynodes, ref meshoned.nbranchlengths, ref meshoned.nodex, ref meshoned.nodey, ref jsferic, ref meshoneddim.nbranches, ref meshoneddim.ngeometry, ref meshoneddim.numnode); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshoned, ref meshoneddim, ref start_index); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim, ref start_index); Assert.That(ierr, Is.EqualTo(0)); //5. make the links int c_npl = 0; int c_nOneDMask = meshoneddim.numnode; IntPtr c_oneDmask = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numnode); int[] oneDmask = new int[meshoneddim.numnode]; for (int i=0;i< meshoneddim.numnode; i++) { oneDmask[i] = 1; } Marshal.Copy(oneDmask, 0, c_oneDmask, meshoneddim.numnode); ierr = wrapperGridgeom.ggeo_make1D2DEmbeddedLinks(ref c_nOneDMask, ref c_oneDmask); Assert.That(ierr, Is.EqualTo(0)); //6. get the number of links int n1d2dlinks = 0; int linkType = 3; ierr = wrapperGridgeom.ggeo_get_links_count(ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); IntPtr c_arrayfrom = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * n1d2dlinks); //2d cell number IntPtr c_arrayto = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * n1d2dlinks); //1d node ierr = wrapperGridgeom.ggeo_get_links(ref c_arrayfrom, ref c_arrayto, ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); int[] rc_arrayfrom = new int[n1d2dlinks]; int[] rc_arrayto = new int[n1d2dlinks]; Marshal.Copy(c_arrayfrom, rc_arrayfrom, 0, n1d2dlinks); Marshal.Copy(c_arrayto, rc_arrayto, 0, n1d2dlinks); var arrayfrom = new int[] { 1, 11, 12, 22, 23, 33, 34, 44, 45, 55, 56, 66, 67, 77, 78, 79, 89 }; var arrayto = new int[] { 1, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5 , 5 }; for (int i = 0; i < n1d2dlinks; i++) { Assert.That(rc_arrayfrom[i], Is.EqualTo(arrayfrom[i])); Assert.That(rc_arrayto[i], Is.EqualTo(arrayto[i])); } // 7.deallocate memory of gridgeom ierr = wrapperGridgeom.ggeo_deallocate(); Assert.That(ierr, Is.EqualTo(0)); } // The network of this test can be seen in packages\test_data\1dnetworkMake1D2DRiverLinksTest.png [Test] [TestCase("default")] [NUnit.Framework.Category("make1D2DRiverLinks")] public void make1D2DRiverLinks(string path) { //1. open the file with the 1d network, 1d mesh and 2d grid. string c_path = TestHelper.TestFilesDirectoryPath() + @"\1d2dForEmbeddedRiverLinks.nc"; Assert.IsTrue(File.Exists(c_path)); int ioncid = 0; //file variable int mode = 0; //create in read mode var wrapperNetcdf = new IoNetcdfLibWrapper(); var wrapperGridgeom = new GridGeomLibWrapper(); int iconvtype = 2; double convversion = 0.0; var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); Assert.That(ierr, Is.EqualTo(0)); int meshid = -1; int networkid = -1; bool includeArrays = true; int start_index = 1; int jsferic = 0; //--------------------------------------------------------------------------------------// //2a. get network and mesh ids ierr = wrapperNetcdf.ionc_get_1d_network_id(ref ioncid, ref networkid); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperNetcdf.ionc_get_1d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); //2b. get dimensions var meshoneddim = new meshgeomdim(); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref networkid, ref meshoneddim); Assert.That(ierr, Is.EqualTo(0)); //2c. allocate var meshoned = new meshgeom(); //mesh variables meshoned.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); //network variables meshoned.nnodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.nnodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); meshoned.branchidx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numnode); meshoned.nedge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nnodes * 2); meshoned.nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); meshoned.branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); meshoned.ngeopointx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.ngeopointy = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); meshoned.nbranchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); meshoned.nbranchorder = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); //2d. get 1d meshgeom array ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref networkid, ref meshoned, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_count_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoneddim.nbranches, ref meshoneddim.nnodes, ref meshoneddim.numedge, ref start_index); Assert.That(ierr, Is.EqualTo(0)); meshoned.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numedge * 2); ierr = wrapperGridgeom.ggeo_create_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoned.edge_nodes, ref meshoneddim.nbranches, ref meshoneddim.nnodes, ref meshoneddim.numedge, ref start_index); Assert.That(ierr, Is.EqualTo(0)); //Read the created edge nodes int[] edge_nodes = new int[meshoneddim.numedge * 2]; Marshal.Copy(meshoned.edge_nodes, edge_nodes, 0, meshoneddim.numedge * 2); Assert.That(edge_nodes.Contains(0), Is.EqualTo(false)); //--------------------------------------------------------------------------------------// //3a. get mesh ids ierr = wrapperNetcdf.ionc_get_2d_mesh_id(ref ioncid, ref meshid); Assert.That(ierr, Is.EqualTo(0)); networkid = -1; //3b. get dimensions var meshtwoddim = new meshgeomdim(); ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref networkid, ref meshtwoddim); Assert.That(ierr, Is.EqualTo(0)); //3c. allocate var meshtwod = new meshgeom(); meshtwod.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshtwoddim.numnode); meshtwod.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshtwoddim.numedge * 2); //3d. get meshgeom ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref networkid, ref meshtwod, ref start_index, ref includeArrays); Assert.That(ierr, Is.EqualTo(0)); //8. close file ierr = wrapperNetcdf.ionc_close(ref ioncid); Assert.That(ierr, Is.EqualTo(0)); //--------------------------------------------------------------------------------------// //4. Herman datastructure for creating the links. ierr = wrapperGridgeom.ggeo_get_xy_coordinates(ref meshoned.branchidx, ref meshoned.branchoffsets, ref meshoned.ngeopointx, ref meshoned.ngeopointy, ref meshoned.nbranchgeometrynodes, ref meshoned.nbranchlengths, ref meshoned.nodex, ref meshoned.nodey, ref jsferic, ref meshoneddim.nbranches, ref meshoneddim.ngeometry, ref meshoneddim.numnode); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshoned, ref meshoneddim, ref start_index); Assert.That(ierr, Is.EqualTo(0)); ierr = wrapperGridgeom.ggeo_convert(ref meshtwod, ref meshtwoddim, ref start_index); Assert.That(ierr, Is.EqualTo(0)); //5. make the links int c_npl = 0; int c_nOneDMask = 0; double c_searchRadius = 100.0; IntPtr c_oneDmask = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * 0); ierr = wrapperGridgeom.ggeo_make1D2DRiverLinks(ref c_nOneDMask, ref c_oneDmask, ref c_searchRadius); Assert.That(ierr, Is.EqualTo(0)); //6. get the number of links int n1d2dlinks = 0; int linkType = 3; ierr = wrapperGridgeom.ggeo_get_links_count(ref n1d2dlinks, ref linkType); Assert.That(ierr, Is.EqualTo(0)); //7. deallocate memory of gridgeom ierr = wrapperGridgeom.ggeo_deallocate(); Assert.That(ierr, Is.EqualTo(0)); } // Test Jan Mooiman ugrid files //[Test] //[TestCase(@"D:\carniato\LUCA\ENGINES\delft3d\ossTest9\src\janm\ugrid_1d_network_mesh_1.nc")] ////[TestCase(@"pathToNetFile")] //[NUnit.Framework.Category("generateMeshEdgeNodes")] //public void generateMeshEdgeNodes(string path) //{ // //1. open file // var c_path = path; // Assert.IsTrue(File.Exists(c_path)); // int ioncid = 0; //file variable // int mode = 0; //create in read mode // var wrapperNetcdf = new IoNetcdfLibWrapper(); // var wrapperGridgeom = new GridGeomLibWrapper(); // int iconvtype = 2; // double convversion = 0.0; // var ierr = wrapperNetcdf.ionc_open(c_path, ref mode, ref ioncid, ref iconvtype, ref convversion); // Assert.That(ierr, Is.EqualTo(0)); // int meshid = -1; // int networkid = -1; // bool includeArrays = true; // int start_index = 1; // int jsferic = 0; // //--------------------------------------------------------------------------------------// // //2a. get network and mesh ids // ierr = wrapperNetcdf.ionc_get_1d_network_id(ref ioncid, ref networkid); // Assert.That(ierr, Is.EqualTo(0)); // ierr = wrapperNetcdf.ionc_get_1d_mesh_id(ref ioncid, ref meshid); // Assert.That(ierr, Is.EqualTo(0)); // //2b. get dimensions // var meshoneddim = new meshgeomdim(); // ierr = wrapperNetcdf.ionc_get_meshgeom_dim(ref ioncid, ref meshid, ref networkid, ref meshoneddim); // Assert.That(ierr, Is.EqualTo(0)); // //2c. allocate // var meshoned = new meshgeom(); // //mesh variables // meshoned.nodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); // meshoned.nodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); // meshoned.nodez = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); // //network variables // meshoned.nnodex = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); // meshoned.nnodey = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nnodes); // meshoned.branchidx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numnode); // meshoned.nedge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nnodes * 2); // meshoned.nbranchgeometrynodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.nbranches); // meshoned.branchoffsets = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.numnode); // meshoned.ngeopointx = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); // meshoned.ngeopointy = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.ngeometry); // meshoned.nbranchlengths = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); // meshoned.nbranchorder = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(double)) * meshoneddim.nbranches); // //2d. get 1d meshgeom array // ierr = wrapperNetcdf.ionc_get_meshgeom(ref ioncid, ref meshid, ref networkid, ref meshoned, ref start_index, // ref includeArrays); // Assert.That(ierr, Is.EqualTo(0)); // ierr = wrapperGridgeom.ggeo_count_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, // ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoneddim.nbranches, ref meshoneddim.nnodes, // ref meshoneddim.numedge, ref start_index); // Assert.That(ierr, Is.EqualTo(0)); // meshoned.edge_nodes = Marshal.AllocCoTaskMem(Marshal.SizeOf(typeof(int)) * meshoneddim.numedge * 2); // ierr = wrapperGridgeom.ggeo_create_edge_nodes(ref meshoned.branchoffsets, ref meshoned.nbranchlengths, // ref meshoned.branchidx, ref meshoned.nedge_nodes, ref meshoned.edge_nodes, ref meshoneddim.nbranches, // ref meshoneddim.nnodes, // ref meshoneddim.numedge, ref start_index); // Assert.That(ierr, Is.EqualTo(0)); // //Read the created edge nodes // int[] edge_nodes = new int[meshoneddim.numedge * 2]; // Marshal.Copy(meshoned.edge_nodes, edge_nodes, 0, meshoneddim.numedge * 2); // Assert.That(edge_nodes.Contains(0), Is.EqualTo(false)); // //8. close file // ierr = wrapperNetcdf.ionc_close(ref ioncid); // Assert.That(ierr, Is.EqualTo(0)); //} } }