/[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.6 by adcroft, Tue Jul 6 20:23:32 2004 UTC revision 1.7 by adcroft, Thu Jul 8 17:03:29 2004 UTC
# Line 129  C                saved data...) Line 129  C                saved data...)
129                 ENDIF                 ENDIF
130  C              update OBC to next timestep  C              update OBC to next timestep
131                 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)+
132       &          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))*
133       &          (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))+
134       &          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)))
135  C              vVel  C              vVel
# Line 155  C                saved data...) Line 155  C                saved data...)
155                 ENDIF                 ENDIF
156  C              update OBC to next timestep  C              update OBC to next timestep
157                 OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+                 OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+
158       &           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))*
159       &           (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))+
160       &           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)))
161  C              Temperature  C              Temperature
# Line 181  C                saved data...) Line 181  C                saved data...)
181                 ENDIF                 ENDIF
182  C              update OBC to next timestep  C              update OBC to next timestep
183                 OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+                 OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+
184       &          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))*
185       &          (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))+
186       &          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)))
187  C              Salinity  C              Salinity
# Line 207  C                saved data...) Line 207  C                saved data...)
207                 ENDIF                 ENDIF
208  C              update OBC to next timestep  C              update OBC to next timestep
209                 OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+                 OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+
210       &           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))*
211       &           (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))+
212       &           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)))
213  C              wVel  C              wVel
# Line 234  C                saved data...) Line 234  C                saved data...)
234                 ENDIF                 ENDIF
235  C                 update OBC to next timestep  C                 update OBC to next timestep
236                    OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+                    OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+
237       &             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))*
238       &            (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))+
239       &            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)))
240  #endif  #endif
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22