Welcome, Guest
Username: Password: Remember me

TOPIC: two Initial conditions

two Initial conditions 13 years 3 weeks ago #2945

  • m.mahgoob
  • m.mahgoob's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 106
can one define two initial conditions using Telemac2d, for example constant water level and constant velocity.

Or telemac takes only one of them, actually I tried it and every time I define initial velocity however telemac starts with zero velocity and consider only the initial elevation. I did it as follows:

INITIAL CONDITIONS: CONSTANT ELEVATION
INITIAL ELEVATION : 0.37
INITIAL CONDITIONS: 'PARTICULAR'
FORTRAN FILE: condin.f
The administrator has disabled public write access.

Re: two Initial conditions 13 years 3 weeks ago #2948

  • jmhervouet
  • jmhervouet's Avatar
Hello,

It should be OK. INITIAL CONDITIONS : 'PARTICULAR' will treat the depth. After in condin.f you will find that the velocity is set to 0.D0 by lines:

CALL OS('X=C ',X=U,C=0.D0)
CALL OS('X=C ',X=V,C=0.D0)

that come after the section on depth.

If you have initialised your velocity in the 'PARTICULAR' section you have to remove these lines. You can also keep them and change the constants C=0.D0 into your required values.

I understand that the keyword INITIAL CONDITIONS is misleading, it should be INITIAL CONDITIONS FOR DEPTH.

That should be it, otherwise send your file condin.f,

With best regards,

Jean-Michel Hervouet
The administrator has disabled public write access.

Re: two Initial conditions 13 years 3 weeks ago #2950

  • m.mahgoob
  • m.mahgoob's Avatar
  • OFFLINE
  • Senior Boarder
  • Posts: 106
I have already changed the lines you mentioned, however the telemac starting usually from zero velocity.

anyway, attached the condin.f and the steering file, I want to have an initial condition of a velocity U= 0.05, V = 0.02, and initial elevation = 0.37 m

thanks for your help
Attachments:
The administrator has disabled public write access.

Re: two Initial conditions 13 years 3 weeks ago #2955

  • jmhervouet
  • jmhervouet's Avatar
Hello,

Try the attached condin.f. In what you have written, you do not give the depth after giving the velocities, so I added it, plus a message indicating that we do go in the 'PARTICULAR' section. This gives, if for some reason the file is not attached:

C ZONE A MODIFIER
WRITE(LU,*) 'CONSTANT VELOCITIES GIVEN'
CALL OS('X=C ',X=U,C=0.05D0)
CALL OS('X=C ',X=V,C=0.02D0)
! copied from 'CONSTANT ELEVATION' section (I do not use COTINI but
! you can if you want)
CALL OS('X=C ',X=H,C=0.37D0)
CALL OS('X=X-Y ',X=H,Y=ZF)
C FIN DE LA ZONE A MODIFIER

That should do it,

With best regards,

Jean-Michel Hervouet
Attachments:
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.