open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi,

Couldn't add the other files in the above post.
Please use this link to access the files:

drive.google.com/drive/folders/1utKzmP3f...7UG0E4vk?usp=sharing

Hello Nitesh,

I adjusted this python script based on yours.
I checked the new output selafin on bluekenue and it looks ok

Hope this helps

Regards,

Taoan

Thank you so much Taoan!
This is exactly what I needed!

hi taoan
I tried your script only for the reading part. I changed the path to the pputils-master directory but I get an error on importing the module. i tried using:
set PYTHONPATH=%PYTHONPATH%;C:\opentelemac-mascaret\v8p0r0\scripts\python27\pputils-master
but its not working either
any idea what should I do?

Hello,

import errors are related with path location for the module
is it your first experience with python?

My suggestion is to change for something like this line:

sys.path.append(path.join( path.dirname(sys.argv[0]),'C:\opentelemac-mascaret\\v8p0r0\\scripts\\python27\\pputils-master'))# put your pputils adress here

On windows double \\ is necessary.

Can you post your script file for checking?

With best regards,

Taoan

Hello Taoan,

Thanks for sharing the python script, based on which I created the binary file for the time and space varying variable ZSD (Secchi depth). While checking the output Selafin file with BlueKenue I found the variable name remained unchanged, it seems line 50:didn't work in my case. Could you please have a look at it?

Best regards,
Yunhao

Hello again,

I slightly modified your python script by placing 'slf2.writeHeader()' after line 50 'vnames.append' and problem solved!

Cheers,
Yunhao

Hi all,

I thought I'd contribute a similar version of the read-write script above. My goal was to merge values from two different files, namely copying BOTTOM from a "bathymetry" file to a "previous computation" file (from GAIA) so that sediment information would be preserved, but a new bottom would be applied. I also take the opportunity to reset CUMUL BED EVOL to 0 over the whole domain. If necessary, it would be possible to also use data from the first time step of the "previous computation" file and apply it to the last time step, but I didn't need to do that in this case.

The output file contains the last time step of the "previous computation" file, with the modified variables. This, of course, assumes that the two files have the same number of points. I also modified the script to be slightly more in keeping with Python best practices; as well, I use it from the command line with the following syntax:As in the other script, you will have to modify the path to your pputils install. There is probably a way to add pputils to $PYTHONPATH and have it work automatically, but I couldn't be bothered to figure that out

I hope this can be useful to someone else, even if only as a starting point to your own scripting. Maybe someday, it will be easy to read/write the selafin format without external libraries