In my opinion, after the codes, as follow, of calculating the Potential Tides,

Code: Select all

```
DO j=JstrR,JendR
DO i=IstrR,IendR
# ifdef TIDES_ASTRO
! Convert Vu_sat from cycles to radians
Ptide(i,j)=Ptide(i,j)+ &
& ramp*POT_Tamp(i,j,itide)* &
& COS(-POT_Tphase(i,j,itide)+ &
& 2.0_r8*pi*Vu_sat(i,j,itide))
# else
Ptide(i,j)=Ptide(i,j)+ &
& ramp*POT_Tamp(i,j,itide)* &
& COS(omega-POT_Tphase(i,j,itide))
# endif
# ifdef MASKING
Ptide(i,j)=Ptide(i,j)*rmask(i,j)
# endif
END DO
END DO
```

Code: Select all

`OCEAN(ng)%zeta = OCEAN(ng)%zeta + Ptide`

But for the index i of

**OCEAN(ng)%zeta(:,:,i)**, it has three options,

**krhs**,

**kstp**,

**knew**. Which one should I use?

Thanks.