112 |
& f2*CVEL_UW(J,K,bi,bj) |
& f2*CVEL_UW(J,K,bi,bj) |
113 |
C update OBC to next timestep |
C update OBC to next timestep |
114 |
OBWu(J,K,bi,bj)=uVel(I_obc+1,J,K,bi,bj)+ |
OBWu(J,K,bi,bj)=uVel(I_obc+1,J,K,bi,bj)+ |
115 |
& CVEL_UW(J,K,bi,bj)*(deltaT/dxF(I_obc+1,J,bi,bj))* |
& CVEL_UW(J,K,bi,bj)*deltaT*recip_dxF(I_obc+1,J,bi,bj)* |
116 |
& (ab1*(uVel(I_obc+2,J,K,bi,bj)-uVel(I_obc+1,J,K,bi,bj))+ |
& (ab1*(uVel(I_obc+2,J,K,bi,bj)-uVel(I_obc+1,J,K,bi,bj))+ |
117 |
& ab2*(UW_STORE_1(J,K,bi,bj)-UW_STORE_4(J,K,bi,bj))) |
& ab2*(UW_STORE_1(J,K,bi,bj)-UW_STORE_4(J,K,bi,bj))) |
118 |
C vVel |
C vVel |
132 |
& f2*CVEL_VW(J,K,bi,bj) |
& f2*CVEL_VW(J,K,bi,bj) |
133 |
C update OBC to next timestep |
C update OBC to next timestep |
134 |
OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+ |
OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+ |
135 |
& CVEL_VW(J,K,bi,bj)*(deltaT/dxV(I_obc+1,J,bi,bj))* |
& CVEL_VW(J,K,bi,bj)*deltaT*recip_dxV(I_obc+1,J,bi,bj)* |
136 |
& (ab1*(vVel(I_obc+1,J,K,bi,bj)-vVel(I_obc,J,K,bi,bj))+ |
& (ab1*(vVel(I_obc+1,J,K,bi,bj)-vVel(I_obc,J,K,bi,bj))+ |
137 |
& ab2*(VW_STORE_1(J,K,bi,bj)-VW_STORE_4(J,K,bi,bj))) |
& ab2*(VW_STORE_1(J,K,bi,bj)-VW_STORE_4(J,K,bi,bj))) |
138 |
C Temperature |
C Temperature |
152 |
& f2*CVEL_TW(J,K,bi,bj) |
& f2*CVEL_TW(J,K,bi,bj) |
153 |
C update OBC to next timestep |
C update OBC to next timestep |
154 |
OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+ |
OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+ |
155 |
& CVEL_TW(J,K,bi,bj)*(deltaT/dxC(I_obc+1,J,bi,bj))* |
& CVEL_TW(J,K,bi,bj)*deltaT*recip_dxC(I_obc+1,J,bi,bj)* |
156 |
& (ab1*(theta(I_obc+1,J,K,bi,bj)-theta(I_obc,J,K,bi,bj))+ |
& (ab1*(theta(I_obc+1,J,K,bi,bj)-theta(I_obc,J,K,bi,bj))+ |
157 |
& ab2*(TW_STORE_1(J,K,bi,bj)-TW_STORE_4(J,K,bi,bj))) |
& ab2*(TW_STORE_1(J,K,bi,bj)-TW_STORE_4(J,K,bi,bj))) |
158 |
C Salinity |
C Salinity |
172 |
& f2*CVEL_SW(J,K,bi,bj) |
& f2*CVEL_SW(J,K,bi,bj) |
173 |
C update OBC to next timestep |
C update OBC to next timestep |
174 |
OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+ |
OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+ |
175 |
& CVEL_SW(J,K,bi,bj)*(deltaT/dxC(I_obc+1,J,bi,bj))* |
& CVEL_SW(J,K,bi,bj)*deltaT*recip_dxC(I_obc+1,J,bi,bj)* |
176 |
& (ab1*(salt(I_obc+1,J,K,bi,bj)-salt(I_obc,J,K,bi,bj))+ |
& (ab1*(salt(I_obc+1,J,K,bi,bj)-salt(I_obc,J,K,bi,bj))+ |
177 |
& ab2*(SW_STORE_1(J,K,bi,bj)-SW_STORE_4(J,K,bi,bj))) |
& ab2*(SW_STORE_1(J,K,bi,bj)-SW_STORE_4(J,K,bi,bj))) |
178 |
C wVel |
C wVel |
193 |
& + f2*CVEL_WW(J,K,bi,bj) |
& + f2*CVEL_WW(J,K,bi,bj) |
194 |
C update OBC to next timestep |
C update OBC to next timestep |
195 |
OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+ |
OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+ |
196 |
& CVEL_WW(J,K,bi,bj)*(deltaT/dxC(I_obc+1,J,bi,bj))* |
& CVEL_WW(J,K,bi,bj)*deltaT*recip_dxC(I_obc+1,J,bi,bj)* |
197 |
& (ab1*(wVel(I_obc+1,J,K,bi,bj)-wVel(I_obc,J,K,bi,bj))+ |
& (ab1*(wVel(I_obc+1,J,K,bi,bj)-wVel(I_obc,J,K,bi,bj))+ |
198 |
& ab2*(WW_STORE_1(J,K,bi,bj)-WW_STORE_4(J,K,bi,bj))) |
& ab2*(WW_STORE_1(J,K,bi,bj)-WW_STORE_4(J,K,bi,bj))) |
199 |
#endif |
#endif |