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

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

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

Revision 1.3 - (show annotations) (download)
Mon Sep 20 23:22:57 2004 UTC (19 years, 9 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 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
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 #ifdef ALLOW_OBCS_NORTH
33 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 ENDDO
44 #endif
45 #ifdef ALLOW_OBCS_SOUTH
46 DO I=1-Olx,sNx+Olx
47 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 #endif
58
59
60 C Set model variables to OB values on East/West Boundaries
61 #ifdef ALLOW_OBCS_EAST
62 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 ENDDO
73 #endif
74 #ifdef ALLOW_OBCS_WEST
75 DO J=1-Oly,sNy+Oly
76 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 #endif
87
88 #endif
89 RETURN
90 END
  ViewVC Help
Powered by ViewVC 1.1.22