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

Diff of /MITgcm/pkg/obcs/orlanski_west.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 (to be applied at I_ Line 112  C              uVel (to be applied at I_
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
# Line 132  C              vVel Line 132  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
# Line 152  C              Temperature Line 152  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
# Line 172  C              Salinity Line 172  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
# Line 193  C              wVel Line 193  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
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22