open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,
I run coupling T2d-Sisyphe and use subroutine "INIT_SEDIMENT" to define grain size vary in space.
My code:
DO J=1,NPOIN
IF (Y(I)<1152228.D0) THEN
FDM(1)=0.0015D0
FDM(2)=0.002D0
ELSE
FDM(1)=0.0003D0
FDM(2)=0.00005D0
ENDIF
ENDDO
But, error show: "Function y at(1)has no Implicit type"
Help me.
Regards.
Ngoc Anh

Hello,

Y is not in the arguments of the subroutine. It will work if you add at the beginning of the subroutine :

USE DECLARATIONS_SISYPHE, ONLY : Y

note that in this case Y is itself a pointer declared in point_sisyphe.f as :

Y => MESH%Y%R

Regards,

JMH

Thanks you for your responds quickly!
I'll try. Thanks for your helps.

Regards;
NgocAnh.

Hi
I don't think that what you are trying to do will work since the FDM array is independent of space and time - only the mean diameter FDM and the percent of each class AVAIL are allowed to vary spatially -
Instead of changing the FDM,I would recommend to change the avail (percent of each class) array instead and the number of classes.
Hope this helps, Cath

PS for example, for fine sediments on one side and coarse material on the other, put two classes fdm(1) = 0.001 fdm(2)= 0.002

!
DO I=1,NPOIN
DO J= 1,NOMBLAY
IF (Y(J)<1152228.D0) THEN
AVAIL(I,J,1)= 0.
AVAIL(I,J,2)= 1.0D0
ELSE
AVAIL(I,J,1)= 1.0D0
AVAIL(I,J,2)= 0.0
ENDIF
ENDDO
ENDDO

Hi
I don't think that what you are trying to do will work since the FDM array is independent of space and time - only the mean diameter FDM and the percent of each class AVAIL are allowed to vary spatially -
Instead of changing the FDM,I would recommend to change the avail (percent of each class) array instead and the number of classes.
Hope this helps, Cath

PS for example, for fine sediments on one side and coarse material on the other, put two classes fdm(1) = 0.001 fdm(2)= 0.002

!
DO I=1,NPOIN
DO J= 1,NOMBLAY
IF (Y(J)<1152228.D0) THEN
AVAIL(I,J,1)= 0.
AVAIL(I,J,2)= 1.0D0
ELSE
AVAIL(I,J,1)= 1.0D0
AVAIL(I,J,2)= 0.0
ENDIF
ENDDO
ENDDO

Hi, cath
Thanhs for your advices. It is very helpful to me.
Regards,
Ngoc Anh.