Welcome, Guest
Username: Password: Remember me
Forgot your password? Forgot your username? Create an account
Forum
System Modules
GAIA / SISYPHE
SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE
Flat Threaded

TOPIC: SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 9 months ago #39691

  • nshettigar
  • nshettigar's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 16
  • Thank you received: 1
I am trying to understand the correct use of the keyword in GAIA, SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE. When used with telemac keyword, GLOBAL NUMBERS OF SOURCE NODES and WATER DISCHARGE OF SOURCES, does not give any sediment flux at the given node number. I want to just impose a time-independent discharge and sediment flux at 4 node numbers. Does anyone here have some experience to share on this? (I am using version v8p1r2)
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 9 months ago #39692

  • pavans
  • pavans's Avatar
  • OFFLINE
  • Moderator
  • Posts: 46
  • Thank you received: 17
Hello,

I'm sorry but sediments at sources are available only in 3D and not still in 2D. A warning message will be add soon to warn users.

Best regards,

Sara
The administrator has disabled public write access.
The following user(s) said Thank You: nshettigar

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 9 months ago #39693

  • nshettigar
  • nshettigar's Avatar
  • OFFLINE
  • Fresh Boarder
  • Posts: 16
  • Thank you received: 1
Thank you for informing. However is there any way it could be still achieved without relying on SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE? Or is it not possible altogether to impose a sediment concentration at the source point?
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 9 months ago #39695

  • pavans
  • pavans's Avatar
  • OFFLINE
  • Moderator
  • Posts: 46
  • Thank you received: 17
No, up to now it is not possible at all since a modification to the call of the subroutine cvdftr, in telemac2d.F is mandatory to have sediments at sources.
This feature is in the “to do list” but still not realized.

Hope this answer to your question.

Regards

Sara
The administrator has disabled public write access.
The following user(s) said Thank You: nshettigar

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41232

  • toby.jhnsn
  • toby.jhnsn's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 161
  • Thank you received: 8
Hi,

Has this feature been implemented yet?

In the middle of a project deliverable here and just came accross this as I was experiencing the same issue. If not can this be accomplished with SISYPHE?

Regards,
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41237

  • toby.jhnsn
  • toby.jhnsn's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 161
  • Thank you received: 8
So after more testing - it seems that this feature has not yet been implemented which is a shame - would be a very useful feature.

I need to release suspended sediments as point sources in my T2D Model. Can anybody advise which part of the GAIA source I could directly modify to change the depth-averaged suspended sediment concentrations at each timestep to get around this.

Thanks,
Toby
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41239

  • pavans
  • pavans's Avatar
  • OFFLINE
  • Moderator
  • Posts: 46
  • Thank you received: 17
Hello Toby,

I'm sorry but this feature has not been implemented as you say.

You can try the following modification:
go to the subroutine telemac2d.F and replace the FIRST call to cvdftr (the one for the case of coupling with GAIA) with these new lines:
          IF(DEBUG.GT.0) WRITE(LU,*) 'CALLING CVDFTR SOLSYS=',SOLSYS_GAI
          CALL CVDFTR(T%ADR(ITRAC)%P,TTILD%ADR(ITRAC)%P,
       &              TN%ADR(ITRAC)%P,TSCEXP%ADR(ITRAC)%P,
       &              DIFT,ICONVFT(ITRAC),.TRUE.,HN_GAI,HOLD,HPROP,
       &              UCONV_GAI,VCONV_GAI,DM1,ZCONV,SOLSYS_GAI,
       &              VISCT%ADR(ITRAC)%P,VISC_S,FLUER_ADV,SMH,YASMH,
       &              FLUDPT_ADV,.TRUE.,AM1,AM2,ZF,
       &              TBOR%ADR(ITRAC)%P,ATBOR%ADR(ITRAC)%P,
       &              BTBOR%ADR(ITRAC)%P,IT1,MASKTR,MESH,TB,
       &              T1,T2,T4,T10,TE1,TE2,TE3,
       &              KDIR,KDDL,
  !                                              BILAN
       &              DT,ENTET,TETAT,1.D0,INFOGT,.TRUE.,                                                 
       &              OPTADV_TR(ITRAC),1,OPDTRA,OPTBAN,                                                  
       &              MSK,MASKEL,MASKPT,MBOR,S,MASSOU(ITRAC),
  !                   OPTSOU=1 FOR SEDIMENTS
       &              1,SLVTRA(ITRAC),FLBOR,VOLU2D,V2DPAR,UNSV2D,
       &              2,FLBORTRA,FLULIM_GAI,YAFLULIM_GAI,                                                
       &              FLULIM_GAI,.FALSE.,                                                                
       &              DIRFLU,RAIN,PLUIE,TRAIN(ITRAC),
       &              FLULIM_GAI,.FALSE.,MAXADV,TB2,NCO_DIST,NSP_DIST)
          IF(DEBUG.GT.0) WRITE(LU,*) 'BACK FROM CVDFTR'


What is different: TSCEXP%ADR(ITRAC)%P,SMH and YASMH (variables necessary for sources) are now included in the call.
This should be enough.

Hope it helps,

Sara
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41241

  • toby.jhnsn
  • toby.jhnsn's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 161
  • Thank you received: 8
Hi Sara,

Thanks very much - you're a lifesaver. Have implemented the code and successfully compiled. Just running some tests now to confirm everything is working as expected.

Can you please also confirm for me the units of the 'SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCES'. Unable to find any mention of this term in the user manual and the units are not mentioned in the dictionary file either.

Is this the standard g/L, kg/m3? Only reason I'm asking is that the suspended sediment concentration output in the GAIA results is in the units of KG/L.

Many thanks,
Toby
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41252

  • toby.jhnsn
  • toby.jhnsn's Avatar
  • OFFLINE
  • Expert Boarder
  • Posts: 161
  • Thank you received: 8
So the model is working perfect - sediments are entering the domain via the sources now - on first look everything looks correct. I just have one question regarding the GAIA printout 'QUANTITY CREATED BY SOURCE TERM'.

I assume this should be the total quantity of sediment entering via the sources per computational timestep. Strangely I am seeing fluctuating values that are sometimes positive and sometimes negative. Given that my sources discharges are generally steady flows with constant sediment concentations at the sources I am not sure why the printout would be producing such results. Any ideas?

Regards,
The administrator has disabled public write access.

SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE 2 years 1 month ago #41256

  • pavans
  • pavans's Avatar
  • OFFLINE
  • Moderator
  • Posts: 46
  • Thank you received: 17
Hi Toby,

some remarks and hints:
- printout QUANTITY CREATED BY SOURCE TERM for sediments implies : quantity created by sediment release + deposition&erosion terms (which are also considered as source terms).
- if your test is not complex you could compare this value with the value of CUMULATED EROSION and CUMULATED DEPOSITION of the Gaia listing for the bed. Maybe you can try just a few iterations.
- the proposed solution in my previous post was rapidly checked in a simple test case with LIPS advection scheme (SCHEME FOR ADVECTION OF SUSPENDED SEDIMENTS = 4 & SCHEME OPTION FOR ADVECTION OF SUSPENDED SEDIMENTS = 4) and it works but it is not working for NERDS scheme (SCHEME FOR ADVECTION OF SUSPENDED SEDIMENTS = 13). Sorry for that, I think that several changes are necessary for this. But if you do not have tidal flats, it also works for N and PSI schemes (SCHEME FOR ADVECTION OF SUSPENDED SEDIMENTS = 4 or 5). So, check also your advection scheme.
- units for suspended concentration: g/L, kg/m3. I think you were looking at the output for layer concentration and not the one for suspended concentration.

Hope this helps,

Sara
The administrator has disabled public write access.
Flat Threaded
Moderators: Pablo, pavans
Forum
System Modules
GAIA / SISYPHE
SUSPENDED SEDIMENTS CONCENTRATION VALUES AT THE SOURCE

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