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

Annotation of /MITgcm/pkg/obcs/obcs_apply_uv.F

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

Revision 1.3 - (hide annotations) (download)
Mon Sep 20 23:22:57 2004 UTC (19 years, 8 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint62c, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint57, checkpoint56, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62m, checkpoint62l, checkpoint58y_post, checkpoint58t_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint57t_post, checkpoint55c_post, checkpoint57v_post, checkpoint57f_post, checkpoint60, checkpoint61, checkpoint62, checkpoint57a_post, checkpoint57h_pre, checkpoint58w_post, checkpoint57h_post, checkpoint57y_pre, checkpoint55g_post, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint55d_post, checkpoint58e_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint55h_post, checkpoint58n_post, checkpoint57e_post, checkpoint55b_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint56a_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint57a_pre, checkpoint55a_post, checkpoint57o_post, checkpoint61q, checkpoint57k_post, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint55e_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.2: +14 -3 lines 
o merged code to
  * prescribe/read time-dependent open boundaries
    (works in conjunction with exf, cal)
  * sponge layer code for open boundaries
  * each boundary N/S/E/W now has its own CPP option
    (healthy for the adjoint)
1 heimbach 1.3 C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_apply_uv.F,v 1.2.8.1 2002/02/05 20:23:59 heimbach Exp $
2     C $Name: $
3 adcroft 1.2
4     #include "OBCS_OPTIONS.h"
5    
6     SUBROUTINE OBCS_APPLY_UV( bi, bj, K,
7     U uFld, vFld,
8     I myThid )
9     C /==========================================================\
10     C | S/R OBCS_APPLY_UV |
11     C \==========================================================/
12     IMPLICIT NONE
13     C == Global variables ==
14     #include "SIZE.h"
15     #include "EEPARAMS.h"
16     #include "PARAMS.h"
17     #include "GRID.h"
18     #include "OBCS.h"
19    
20     C == Routine Arguments ==
21     INTEGER bi,bj,K
22     _RL uFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
23     _RL vFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24     INTEGER myThid
25    
26     #ifdef ALLOW_OBCS
27    
28     C == Local variables ==
29     INTEGER i,j
30    
31     C Set model variables to OB values on North/South Boundaries
32 heimbach 1.3 #ifdef ALLOW_OBCS_NORTH
33 adcroft 1.2 DO I=1-Olx,sNx+Olx
34     C Northern boundary
35     IF (OB_Jn(I,bi,bj).NE.0) THEN
36     vFld(I,OB_Jn(I,bi,bj),K,bi,bj)=OBNv(I,K,bi,bj)
37     & *_maskS(I,OB_Jn(I,bi,bj),K,bi,bj)
38     uFld(I,OB_Jn(I,bi,bj),K,bi,bj)=OBNu(I,K,bi,bj)
39     & *_maskW(I,OB_Jn(I,bi,bj),K,bi,bj)
40     vFld(I,OB_Jn(I,bi,bj)+1,K,bi,bj)=OBNv(I,K,bi,bj)
41     & *_maskS(I,OB_Jn(I,bi,bj),K,bi,bj)
42     ENDIF
43 heimbach 1.3 ENDDO
44     #endif
45     #ifdef ALLOW_OBCS_SOUTH
46     DO I=1-Olx,sNx+Olx
47 adcroft 1.2 C Southern boundary
48     IF (OB_Js(I,bi,bj).NE.0) THEN
49     vFld(I,OB_Js(I,bi,bj)+1,K,bi,bj)=OBSv(I,K,bi,bj)
50     & *_maskS(I,OB_Js(I,bi,bj)+1,K,bi,bj)
51     uFld(I,OB_Js(I,bi,bj),K,bi,bj)=OBSu(I,K,bi,bj)
52     & *_maskW(I,OB_Js(I,bi,bj),K,bi,bj)
53     vFld(I,OB_Js(I,bi,bj),K,bi,bj)=OBSv(I,K,bi,bj)
54     & *_maskS(I,OB_Js(I,bi,bj)+1,K,bi,bj)
55     ENDIF
56     ENDDO
57 heimbach 1.3 #endif
58 adcroft 1.2
59    
60     C Set model variables to OB values on East/West Boundaries
61 heimbach 1.3 #ifdef ALLOW_OBCS_EAST
62 adcroft 1.2 DO J=1-Oly,sNy+Oly
63     C Eastern boundary
64     IF (OB_Ie(J,bi,bj).NE.0) THEN
65     uFld(OB_Ie(J,bi,bj),J,K,bi,bj)=OBEu(J,K,bi,bj)
66     & *_maskW(OB_Ie(J,bi,bj),J,K,bi,bj)
67     vFld(OB_Ie(J,bi,bj),J,K,bi,bj)=OBEv(J,K,bi,bj)
68     & *_maskS(OB_Ie(J,bi,bj),J,K,bi,bj)
69     uFld(OB_Ie(J,bi,bj)+1,J,K,bi,bj)=OBEu(J,K,bi,bj)
70     & *_maskW(OB_Ie(J,bi,bj),J,K,bi,bj)
71     ENDIF
72 heimbach 1.3 ENDDO
73     #endif
74     #ifdef ALLOW_OBCS_WEST
75     DO J=1-Oly,sNy+Oly
76 adcroft 1.2 C Western boundary
77     IF (OB_Iw(J,bi,bj).NE.0) THEN
78     uFld(OB_Iw(J,bi,bj)+1,J,K,bi,bj)=OBWu(J,K,bi,bj)
79     & *_maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj)
80     vFld(OB_Iw(J,bi,bj),J,K,bi,bj)=OBWv(J,K,bi,bj)
81     & *_maskS(OB_Iw(J,bi,bj),J,K,bi,bj)
82     uFld(OB_Iw(J,bi,bj),J,K,bi,bj)=OBWu(J,K,bi,bj)
83     & *_maskW(OB_Iw(J,bi,bj)+1,J,K,bi,bj)
84     ENDIF
85     ENDDO
86 heimbach 1.3 #endif
87 adcroft 1.2
88     #endif
89     RETURN
90     END
  ViewVC Help
Powered by ViewVC 1.1.22