28.3.3 symmetry

This implements a symmetric boundary condition across the equator. Specifically, the second last row of velocity points must be defined on the equator ( $\phi^U_{jmt-1} = 0.0$). Note that the equator is at the northern end of the domain. The condition applies when j in the memory window corresponds to jrow=jmt-1 and is given by

 

t_i,k,j+1,n = t_i,k,j,n

u_i,k,j+1,1 = u_i,k,j,1

u_i,k,j+1,2 = -u_i,k,j,2

psi_i,jrow+1 = -psi_i,jrow

kmt_i,jrow+1 = kmt_i,jrow (28.30)