/[MITgcm]/MITgcm/pkg/obcs/orlanski_north.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/orlanski_north.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by adcroft, Thu Sep 27 18:13:13 2001 UTC revision 1.4 by jmc, Thu Jul 11 16:22:30 2002 UTC
# Line 112  C              uVel Line 112  C              uVel
112       &                f2*CVEL_UN(I,K,bi,bj)       &                f2*CVEL_UN(I,K,bi,bj)
113  C              update OBC to next timestep  C              update OBC to next timestep
114                 OBNu(I,K,bi,bj)=uVel(I,J_obc,K,bi,bj)-                 OBNu(I,K,bi,bj)=uVel(I,J_obc,K,bi,bj)-
115       &           CVEL_UN(I,K,bi,bj)*(deltaT/dyU(I,J_obc,bi,bj))*       &           CVEL_UN(I,K,bi,bj)*deltaT*recip_dyU(I,J_obc,bi,bj)*
116       &           (ab1*(uVel(I,J_obc,K,bi,bj)-uVel(I,J_obc-1,K,bi,bj)) +       &           (ab1*(uVel(I,J_obc,K,bi,bj)-uVel(I,J_obc-1,K,bi,bj)) +
117       &           ab2*(UN_STORE_4(I,K,bi,bj)-UN_STORE_1(I,K,bi,bj)))       &           ab2*(UN_STORE_4(I,K,bi,bj)-UN_STORE_1(I,K,bi,bj)))
118  C              vVel  C              vVel
# Line 132  C              vVel Line 132  C              vVel
132       &                f2*CVEL_VN(I,K,bi,bj)       &                f2*CVEL_VN(I,K,bi,bj)
133  C              update OBC to next timestep  C              update OBC to next timestep
134                 OBNv(I,K,bi,bj)=vVel(I,J_obc,K,bi,bj)-                 OBNv(I,K,bi,bj)=vVel(I,J_obc,K,bi,bj)-
135       &           CVEL_VN(I,K,bi,bj)*(deltaT/dyF(I,J_obc-1,bi,bj))*       &           CVEL_VN(I,K,bi,bj)*deltaT*recip_dyF(I,J_obc-1,bi,bj)*
136       &           (ab1*(vVel(I,J_obc,K,bi,bj)-vVel(I,J_obc-1,K,bi,bj)) +       &           (ab1*(vVel(I,J_obc,K,bi,bj)-vVel(I,J_obc-1,K,bi,bj)) +
137       &           ab2*(VN_STORE_4(I,K,bi,bj)-VN_STORE_1(I,K,bi,bj)))             &           ab2*(VN_STORE_4(I,K,bi,bj)-VN_STORE_1(I,K,bi,bj)))      
138  C              Temperature  C              Temperature
# Line 152  C              Temperature Line 152  C              Temperature
152       &                f2*CVEL_TN(I,K,bi,bj)       &                f2*CVEL_TN(I,K,bi,bj)
153  C              update OBC to next timestep  C              update OBC to next timestep
154                 OBNt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)-                 OBNt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)-
155       &          CVEL_TN(I,K,bi,bj)*(deltaT/dyC(I,J_obc,bi,bj))*       &          CVEL_TN(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc,bi,bj)*
156       &          (ab1*(theta(I,J_obc,K,bi,bj)-theta(I,J_obc-1,K,bi,bj))+       &          (ab1*(theta(I,J_obc,K,bi,bj)-theta(I,J_obc-1,K,bi,bj))+
157       &          ab2*(TN_STORE_4(I,K,bi,bj)-TN_STORE_1(I,K,bi,bj)))       &          ab2*(TN_STORE_4(I,K,bi,bj)-TN_STORE_1(I,K,bi,bj)))
158  C              Salinity  C              Salinity
# Line 172  C              Salinity Line 172  C              Salinity
172       &                f2*CVEL_SN(I,K,bi,bj)       &                f2*CVEL_SN(I,K,bi,bj)
173  C              update OBC to next timestep  C              update OBC to next timestep
174                 OBNs(I,K,bi,bj)=salt(I,J_obc,K,bi,bj)-                 OBNs(I,K,bi,bj)=salt(I,J_obc,K,bi,bj)-
175       &          CVEL_SN(I,K,bi,bj)*(deltaT/dyC(I,J_obc,bi,bj))*       &          CVEL_SN(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc,bi,bj)*
176       &          (ab1*(salt(I,J_obc,K,bi,bj)-salt(I,J_obc-1,K,bi,bj)) +       &          (ab1*(salt(I,J_obc,K,bi,bj)-salt(I,J_obc-1,K,bi,bj)) +
177       &          ab2*(SN_STORE_4(I,K,bi,bj)-SN_STORE_1(I,K,bi,bj)))       &          ab2*(SN_STORE_4(I,K,bi,bj)-SN_STORE_1(I,K,bi,bj)))
178  C              wVel  C              wVel
# Line 193  C              wVel Line 193  C              wVel
193       &                   + f2*CVEL_WN(I,K,bi,bj)       &                   + f2*CVEL_WN(I,K,bi,bj)
194  C                 update OBC to next timestep  C                 update OBC to next timestep
195                    OBNw(I,K,bi,bj)=wVel(I,J_obc,K,bi,bj)-                    OBNw(I,K,bi,bj)=wVel(I,J_obc,K,bi,bj)-
196       &             CVEL_WN(I,K,bi,bj)*(deltaT/dyC(I,J_obc,bi,bj))*       &             CVEL_WN(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc,bi,bj)*
197       &             (ab1*(wVel(I,J_obc,K,bi,bj)-wVel(I,J_obc-1,K,bi,bj))+       &             (ab1*(wVel(I,J_obc,K,bi,bj)-wVel(I,J_obc-1,K,bi,bj))+
198       &             ab2*(WN_STORE_4(I,K,bi,bj)-WN_STORE_1(I,K,bi,bj)))       &             ab2*(WN_STORE_4(I,K,bi,bj)-WN_STORE_1(I,K,bi,bj)))
199  #endif  #endif
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22