A problem of riverplume test case

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
xiaom6515
Posts: 3
Joined: Fri Mar 10, 2017 11:22 pm
Location: SYSU

A problem of riverplume test case

#1 Unread post by xiaom6515 »

Dear all,
I want to add the sediment in riverplume.So I modified the test case which is called riverplume1.
At first I modified the ana_psource.h file just like this

Code: Select all

# ifdef SEDIMENT
 
         DO k=1,N(ng)
           DO is=1,Nsrc(ng)
            DO ised=1,NST
             SOURCES(ng)%Tsrc(is,k,idsed(ised))=5.0_r8
	    END DO
           END DO
	  END DO
# endif
But when I run this case,I found there exist a error that the highest sediment concentration is

not 5,even, is less than half of 5.
After I try to modify some of the parameters,I found that if I modified the salinity, the highest
sediment concentration would change too.
I modified the ana_psource.h file just like this(SOURCES(ng)%Tsrc

(is,k,isalt)=0.0_r8→SOURCES(ng)%Tsrc(is,k,isalt)=28.0_r8):

Code: Select all

!  Set-up tracer (tracer units) point Sources/Sinks.
!
# if defined RIVERPLUME1
        IF (DOMAIN(ng)%NorthEast_Test(tile)) THEN
          DO k=1,N(ng)
            DO is=1,Nsrc(ng)
              SOURCES(ng)%Tsrc(is,k,itemp)=T0(ng)    ! T0 is in ocean.in
              SOURCES(ng)%Tsrc(is,k,isalt)=28.0_r8
            END DO
          END DO
        END IF
The highest sediment concentration would be 4,but it was not 5.
But I don't known what the problem is,could anyone tell me how could this be?
Thanks in advance.
Attachments
result2.png
result1.png
ana_psource.h
(16.68 KiB) Downloaded 310 times

Post Reply