TELEMAC-MASCARET Forum: USER_CONDI3D

Welcome, Guest

TOPIC: USER_CONDI3D_H subroutine

USER_CONDI3D_H subroutine 3 months 4 weeks ago #45236

giuliaperego
OFFLINE
Fresh Boarder
Posts: 22
Thank you received: 1

Hello!
I am trying to do a simulation in a channel.
I would like to start from an initial condition where I have a water depth of 2 meters only in the main channel, and 0 everywhere else in the domain (picture of the domain is attached).
I tried with a fortran file USER_CONDI3D_H.

I added to my .cas file:
FORTRAN FILE : 'USER_CONDI3D_H.f'
INITIAL CONDITIONS = SPECIAL

and in my fortran file I specified the x,y coordinates of the main channel in MODULE CHANNEL_POLYGON, and afterward I set a water depth of 2 m for the points inside those borders.

I tried to run the simulation but a lot of errors appear...
I am a beginner with Fortran programming so my code might be full of errors. I would be very grateful if anybody could help me to correct my mistakes.
I attach fortran file, error screenshot and .cas

Thank you so much
Giulia

Attachments:

The administrator has disabled public write access.

USER_CONDI3D_H subroutine 3 months 4 weeks ago #45237

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	Hi fortran should be in the fixed format, not in the free form! You have to follow the rules by starting at the 7th column ... Hope this helps Otherwise as you only want to fix the initial condition, you probably could use BK to create a selafin file with desired the water depth and use this file as a previous computation file. Simpler and without any request of fortran programming...
	Christophe The administrator has disabled public write access. The following user(s) said Thank You: giuliaperego

USER_CONDI3D_H subroutine 3 months 4 weeks ago #45238

giuliaperego OFFLINE Fresh Boarder Posts: 22 Thank you received: 1	Thank you for your response. How do I set the water table level in Blue Kenue? I didn't know it was possible Giulia
	The administrator has disabled public write access.

USER_CONDI3D_H subroutine 3 months 4 weeks ago #45239

c.coulet OFFLINE Moderator Posts: 3722 Thank you received: 1031	You should create a new selafin file (from the geo file) and add a new variable, then set the values like you want (use mapping of polygon ...) Regards
	Christophe The administrator has disabled public write access. The following user(s) said Thank You: giuliaperego

USER_CONDI3D_H subroutine 3 months 4 weeks ago #45241

giuliaperego OFFLINE Fresh Boarder Posts: 22 Thank you received: 1	Thank you for this precious advice! It worked :)
	The administrator has disabled public write access.