/[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.8 by heimbach, Mon Sep 20 23:22:58 2004 UTC
# Line 87  C     == Routine arguments == Line 87  C     == Routine arguments ==
87        INTEGER myThid        INTEGER myThid
88    
89  #ifdef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
90    #ifdef ALLOW_OBCS_WEST
91    
92  C     == Local variables ==  C     == Local variables ==
93        INTEGER J, K, I_obc        INTEGER J, K, I_obc
# Line 129  C                saved data...) Line 130  C                saved data...)
130                 ENDIF                 ENDIF
131  C              update OBC to next timestep  C              update OBC to next timestep
132                 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)+
133       &          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))*
134       &          (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))+
135       &          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)))
136  C              vVel  C              vVel
# Line 155  C                saved data...) Line 156  C                saved data...)
156                 ENDIF                 ENDIF
157  C              update OBC to next timestep  C              update OBC to next timestep
158                 OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+                 OBWv(J,K,bi,bj)=vVel(I_obc,J,K,bi,bj)+
159       &           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))*
160       &           (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))+
161       &           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)))
162  C              Temperature  C              Temperature
# Line 181  C                saved data...) Line 182  C                saved data...)
182                 ENDIF                 ENDIF
183  C              update OBC to next timestep  C              update OBC to next timestep
184                 OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+                 OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+
185       &          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))*
186       &          (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))+
187       &          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)))
188  C              Salinity  C              Salinity
# Line 207  C                saved data...) Line 208  C                saved data...)
208                 ENDIF                 ENDIF
209  C              update OBC to next timestep  C              update OBC to next timestep
210                 OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+                 OBWs(J,K,bi,bj)=salt(I_obc,J,K,bi,bj)+
211       &           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))*
212       &           (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))+
213       &           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)))
214  C              wVel  C              wVel
# Line 234  C                saved data...) Line 235  C                saved data...)
235                 ENDIF                 ENDIF
236  C                 update OBC to next timestep  C                 update OBC to next timestep
237                    OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+                    OBWw(J,K,bi,bj)=wVel(I_obc,J,K,bi,bj)+
238       &             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))*
239       &            (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))+
240       &            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)))
241  #endif  #endif
# Line 285  C              copy (current time) t to Line 286  C              copy (current time) t to
286           ENDDO           ENDDO
287        ENDDO        ENDDO
288    
289    #endif
290  #endif /* ALLOW_ORLANSKI */  #endif /* ALLOW_ORLANSKI */
291        RETURN        RETURN
292        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22