possible bug in t3dbc_im.F, u3dbc_im.F

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
User avatar
jivica
Posts: 169
Joined: Mon May 05, 2003 2:41 pm
Location: The University of Western Australia, Perth, Australia
Contact:

possible bug in t3dbc_im.F, u3dbc_im.F

#1 Unread post by jivica »

Think that there is a possible bug :?: in routines for applying boundary conditions for tracers and velocity when used NUDGING option. Issue is regarding time step in the last term t(Istr-1,j,k,nstp,itrc) should be t(Istr-1,j,k,nout,itrc) to be consistent.

Code: Select all

#   ifdef WEST_TNUDGING
            t(Istr-1,j,k,nout,itrc)=t(Istr-1,j,k,nout,itrc)+            &
     &                              tau*(BOUNDARY(ng)%t_west(j,k,itrc)- &
     &                                   t(Istr-1,j,k,nstp,itrc))
#   endif
We are taking part t=(1-tau)*t + tau*BRY which is fine but time
records are not synced correctly.

Think it should read:

Code: Select all

#   ifdef WEST_TNUDGING
            t(Istr-1,j,k,nout,itrc)=t(Istr-1,j,k,nout,itrc)+            &
     &                              tau*(BOUNDARY(ng)%t_west(j,k,itrc)- &
     &                                   t(Istr-1,j,k,nout,itrc))
#   endif

This is true for all T_ and U_ and all sides..

Or possibly I am wrong and not getting something correctly.
Cheers,
Ivica

Cheers,
Ivica

Post Reply