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

Diff of /MITgcm/pkg/obcs/orlanski_south.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_US(I,K,bi,bj)       &                f2*CVEL_US(I,K,bi,bj)
113  C              update OBC to next timestep  C              update OBC to next timestep
114                 OBSu(I,K,bi,bj)=uVel(I,J_obc,K,bi,bj)+                 OBSu(I,K,bi,bj)=uVel(I,J_obc,K,bi,bj)+
115       &          CVEL_US(I,K,bi,bj)*(deltaT/dyU(I,J_obc+1,bi,bj))*       &          CVEL_US(I,K,bi,bj)*deltaT*recip_dyU(I,J_obc+1,bi,bj)*
116       &          (ab1*(uVel(I,J_obc+1,K,bi,bj)-uVel(I,J_obc,K,bi,bj)) +       &          (ab1*(uVel(I,J_obc+1,K,bi,bj)-uVel(I,J_obc,K,bi,bj)) +
117       &         ab2*(US_STORE_1(I,K,bi,bj)-US_STORE_4(I,K,bi,bj)))       &         ab2*(US_STORE_1(I,K,bi,bj)-US_STORE_4(I,K,bi,bj)))
118  C              vVel (to be applied at J_obc+1)  C              vVel (to be applied at J_obc+1)
# Line 132  C              vVel (to be applied at J_ Line 132  C              vVel (to be applied at J_
132       &                f2*CVEL_VS(I,K,bi,bj)       &                f2*CVEL_VS(I,K,bi,bj)
133  C              update OBC to next timestep  C              update OBC to next timestep
134                 OBSv(I,K,bi,bj)=vVel(I,J_obc+1,K,bi,bj)+                 OBSv(I,K,bi,bj)=vVel(I,J_obc+1,K,bi,bj)+
135       &          CVEL_VS(I,K,bi,bj)*(deltaT/dyF(I,J_obc+1,bi,bj))*       &          CVEL_VS(I,K,bi,bj)*deltaT*recip_dyF(I,J_obc+1,bi,bj)*
136       &          (ab1*(vVel(I,J_obc+2,K,bi,bj)-vVel(I,J_obc+1,K,bi,bj))+       &          (ab1*(vVel(I,J_obc+2,K,bi,bj)-vVel(I,J_obc+1,K,bi,bj))+
137       &          ab2*(VS_STORE_1(I,K,bi,bj)-VS_STORE_4(I,K,bi,bj)))       &          ab2*(VS_STORE_1(I,K,bi,bj)-VS_STORE_4(I,K,bi,bj)))
138  C              Temperature  C              Temperature
# Line 152  C              Temperature Line 152  C              Temperature
152       &                f2*CVEL_TS(I,K,bi,bj)       &                f2*CVEL_TS(I,K,bi,bj)
153  C              update OBC to next timestep  C              update OBC to next timestep
154                 OBSt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)+                 OBSt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)+
155       &          CVEL_TS(I,K,bi,bj)*(deltaT/dyC(I,J_obc+1,bi,bj))*       &          CVEL_TS(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc+1,bi,bj)*
156       &          (ab1*(theta(I,J_obc+1,K,bi,bj)-theta(I,J_obc,K,bi,bj))+       &          (ab1*(theta(I,J_obc+1,K,bi,bj)-theta(I,J_obc,K,bi,bj))+
157       &          ab2*(TS_STORE_1(I,K,bi,bj)-TS_STORE_4(I,K,bi,bj)))       &          ab2*(TS_STORE_1(I,K,bi,bj)-TS_STORE_4(I,K,bi,bj)))
158  C              Salinity  C              Salinity
# Line 172  C              Salinity Line 172  C              Salinity
172       &                f2*CVEL_SS(I,K,bi,bj)       &                f2*CVEL_SS(I,K,bi,bj)
173  C              update OBC to next timestep  C              update OBC to next timestep
174                 OBSs(I,K,bi,bj)=salt(I,J_obc,K,bi,bj)+                 OBSs(I,K,bi,bj)=salt(I,J_obc,K,bi,bj)+
175       &          CVEL_SS(I,K,bi,bj)*(deltaT/dyC(I,J_obc+1,bi,bj))*       &          CVEL_SS(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc+1,bi,bj)*
176       &          (ab1*(salt(I,J_obc+1,K,bi,bj)-salt(I,J_obc,K,bi,bj)) +       &          (ab1*(salt(I,J_obc+1,K,bi,bj)-salt(I,J_obc,K,bi,bj)) +
177       &         ab2*(SS_STORE_1(I,K,bi,bj)-SS_STORE_4(I,K,bi,bj)))       &         ab2*(SS_STORE_1(I,K,bi,bj)-SS_STORE_4(I,K,bi,bj)))
178  C              wVel  C              wVel
# Line 193  C              wVel Line 193  C              wVel
193       &                   + f2*CVEL_WS(I,K,bi,bj)       &                   + f2*CVEL_WS(I,K,bi,bj)
194  C                 update OBC to next timestep  C                 update OBC to next timestep
195                    OBSw(I,K,bi,bj)=wVel(I,J_obc,K,bi,bj)+                    OBSw(I,K,bi,bj)=wVel(I,J_obc,K,bi,bj)+
196       &             CVEL_WS(I,K,bi,bj)*(deltaT/dyC(I,J_obc+1,bi,bj))*       &            CVEL_WS(I,K,bi,bj)*deltaT*recip_dyC(I,J_obc+1,bi,bj)*
197       &            (ab1*(wVel(I,J_obc+1,K,bi,bj)-wVel(I,J_obc,K,bi,bj))+       &            (ab1*(wVel(I,J_obc+1,K,bi,bj)-wVel(I,J_obc,K,bi,bj))+
198       &            ab2*(WS_STORE_1(I,K,bi,bj)-WS_STORE_4(I,K,bi,bj)))       &            ab2*(WS_STORE_1(I,K,bi,bj)-WS_STORE_4(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