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

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

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

Revision 1.3 - (show annotations) (download)
Fri Dec 3 16:47:30 2004 UTC (19 years, 5 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint57m_post, checkpoint62u, checkpoint57g_pre, checkpoint62t, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint57y_post, checkpoint62c, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint57, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint62s, checkpoint58a_post, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint57z_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62x, checkpoint58y_post, checkpoint58t_post, checkpoint57l_post, checkpoint57h_post, checkpoint57t_post, checkpoint57v_post, checkpoint57f_post, checkpoint60, checkpoint61, checkpoint62, checkpoint57a_post, checkpoint57h_pre, checkpoint57x_post, checkpoint58w_post, checkpoint57y_pre, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, checkpoint58m_post, checkpoint57c_pre, checkpoint58r_post, checkpoint58n_post, checkpoint57e_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint62d, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint58j_post, checkpoint57a_pre, checkpoint57o_post, checkpoint61q, checkpoint57k_post, checkpoint61z, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61x, checkpoint61y
Changes since 1.2: +13 -5 lines 
o extend application of tracer open boundary conditions
  to a strip beyond the boundary of width Olx/Oly
1 C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_apply_tloc.F,v 1.2 2004/09/20 23:22:57 heimbach Exp $
2 C $Name: $
3
4 #include "OBCS_OPTIONS.h"
5
6 SUBROUTINE OBCS_APPLY_TLOC( bi, bj, K,
7 U tFld,
8 I myThid )
9 C /==========================================================\
10 C | S/R OBCS_APPLY_TLOC |
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 tFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
23 INTEGER myThid
24
25 #ifdef ALLOW_OBCS
26
27 C == Local variables ==
28 INTEGER I,J,I_obc,J_obc
29 _RL obc_mask
30
31 C Set model variables to OB values on North/South Boundaries
32 DO I=1-Olx,sNx+Olx
33 C Northern boundary
34 #ifdef ALLOW_OBCS_NORTH
35 J_obc = OB_Jn(I,bi,bj)
36 IF (J_obc.NE.0) THEN
37 obc_mask = _maskS(I,J_obc,K,bi,bj)
38 DO J = J_obc, J_obc+Oly
39 tFld(I,J)=OBNt(I,K,bi,bj)*obc_mask
40 ENDDO
41 ENDIF
42 #endif
43 C Southern boundary
44 #ifdef ALLOW_OBCS_SOUTH
45 J_obc = OB_Js(I,bi,bj)
46 IF (J_obc.NE.0) THEN
47 obc_mask = _maskS(I,J_obc+1,K,bi,bj)
48 DO J = J_obc-Oly, J_obc
49 tFld(I,J)=OBSt(I,K,bi,bj)*obc_mask
50 ENDDO
51 ENDIF
52 #endif
53 ENDDO
54
55 C Set model variables to OB values on East/West Boundaries
56 DO J=1-Oly,sNy+Oly
57 C Eastern boundary
58 #ifdef ALLOW_OBCS_EAST
59 I_obc = OB_Ie(J,bi,bj)
60 IF (I_obc.NE.0) THEN
61 obc_mask = _maskW(I_obc,J,K,bi,bj)
62 DO I = I_obc, I_obc+Olx
63 tFld(I,J)=OBEt(J,K,bi,bj)*obc_mask
64 ENDDO
65 ENDIF
66 #endif
67 C Western boundary
68 #ifdef ALLOW_OBCS_WEST
69 I_obc=OB_Iw(J,bi,bj)
70 IF (I_obc.NE.0) THEN
71 obc_mask = _maskW(I_obc+1,J,K,bi,bj)
72 DO I = I_obc-Olx, I_obc
73 tFld(I,J)=OBWt(J,K,bi,bj)*obc_mask
74 ENDDO
75 ENDIF
76 #endif
77 ENDDO
78
79 #endif
80 RETURN
81 END
  ViewVC Help
Powered by ViewVC 1.1.22