I am trying to run a simple analytical model (a rectangular-shape domain with uniform depth). I put M2 tidal elevation and current along the eastern boundary and radiation OBC along western boundary so that tides are radiated out. The northern and southern boundaries are closed. The grid size is 100 x 50. Temperature and salinity are constant throughout the domain.

I want to to put passive tracers (dyes) for a square-area on the center of the domain (11x11 grids) and the following are what I did.

1. I defined T_PASSIVE and ANA_PASSIVE in cppdefs file

Code: Select all

```
#defined T_PASSIVE
#defined ANA_PASSIVE
```

Code: Select all

```
NPT=1
Hout(inert) == T
```

Code: Select all

```
DO ip=1,NPT
itrc=inert(ip)
DO k=1,N(ng)
DO j=20,30
DO i=40,50
t(i,j,k,1,itrc)=20.0_r8
! t(i,j,k,1,1)=20.0_r8
! t(i,j,k,2,itrc)=t(i,j,k,1,itrc)
END DO
END DO
END DO
END DO
```

My cppdef file is also attached. I use COAWST, but the cppdef is pretty much the same as that for a recent Rutgers ROMS.

Any comments would be appreciated.

Thanks in advance.

DJ Kobashi

TAMU