Wave direction using nesting approach

[Quan T. Nguyen]

Dear FUNWAVE Community,

I have a question related to using the nesting tool. I would like to simulate the wave propagation for a large domain and apply the single-grid nesting tool for the smaller domain, but the wave direction in the child grid did not match the parent grid. 

Detailed information on my work is as:

• The size of the parent grid is 18km x 13km with 5m-resolution and that of the child grid is 2km x 2km with 1m-resolution.
• I've used the make-station-single-block Matlab file in the "nesting_tools" folder to create the stations for nesting and the coupling-file-single-block file to generate the coupling.txt file. In which, boundaries of nested domain are defined in the order E, W, S, then N, with E and W starting at the south ends and ending at the north and N and S starting at the west sides and ending at the east. 
• For the bathymetry data of the nested domain, I have transposed the dataset to match the defined boundaries above. But I'm not sure whether we need to transpose the bathymetry of the nested domain or not.
• About the wavemaker, I have applied the new irregular wave with 0% coherence (please check the attachment for my input files). The direction of wave propagation is from the left to the right of the study domain.

The result of the nested domain showed that the waves propagated in both left-to-right and south-to-north directions. So, my question is how to avoid this problem? 

Thank you in advance.

Regards,

Quan.

[Quan T. Nguyen]

Dear Prof. Fengyan Shi,

Thank you for the suggestion of using the automatic nesting program FUNWAVE-MGN.

I have tried the model and got some problems as listed below:

• Compiler error: when I compiled the model, there was an error of mismatch in argument "mask" for the subroutine "CONSTRUCTION_HO_Minmod" in fluxes.F file. Therefore, I have checked the file and modified the function CONSTRUCTION_HO to CONSTRUCTION_HO_minmod for consistency (I have attached the fixed fluxes.F file in this email for your revision). Now the model could be compiled normally. But the calculation is still not correct as explained below.
• Incorrect simulation results: When I used the compiled model above for the benchmark case of Osu_conical_island, besides the initial solitary wave as defined by the wavemaker, the model also generated minor waves surrounding the conical island (red circle in attached images), which therefore interacted with the solitary waves and impacted the simulation results.

My question is how to avoid this problem and is the modification of fluxes.F file correct? 

I have also attached the Makefile and Input file for your reference. Thank you in advance and I am looking forward to your response.

Quan.

On Sat, May 14, 2022 at 5:25 AM Fengyan Shi <fy...@udel.edu> wrote:

Check the coupling file. I don’t see wave signors coming from the southern and northern boundaries in the small domain. Or you can try the automatic nesting program in 

https://github.com/fengyanshi/FUNWAVE-MGN

-Fengyan

--
You received this message because you are subscribed to the Google Groups "funwave-tvd" group.
To unsubscribe from this group and stop receiving emails from it, send an email to funwave-tvd...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/funwave-tvd/CACfbzFbeh4yCGuga%2Bxbun3x%2BXsKuP6ZBDrq3NwkbvPW48f73Qw%40mail.gmail.com.
<input_coupling_child_grid.txt><input_parent_grid.txt><2D_surf_G0_R5m.png><2D_eta_G3_R1m.png>