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

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

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

Revision 1.1 - (show annotations) (download)
Thu Jan 30 00:05:30 2003 UTC (21 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint57m_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint57g_pre, checkpoint50c_post, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint52d_pre, checkpoint57g_post, checkpoint48e_post, checkpoint56b_post, checkpoint50c_pre, checkpoint57y_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint54d_post, checkpoint54e_post, checkpoint51l_post, checkpoint48i_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint50d_pre, checkpoint52k_post, checkpoint58, checkpoint55, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint58f_post, checkpoint52f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint50b_pre, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint51d_post, checkpoint51t_post, checkpoint58t_post, checkpoint51n_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint57t_post, checkpoint55c_post, checkpoint48d_pre, checkpoint51j_post, checkpoint52e_pre, checkpoint57v_post, checkpoint57f_post, checkpoint52e_post, checkpoint51n_pre, checkpoint53d_post, checkpoint57a_post, checkpoint48d_post, checkpoint57h_pre, checkpoint48f_post, checkpoint52b_pre, checkpoint54b_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint57y_pre, checkpoint55g_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint58o_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint58p_post, checkpoint58q_post, checkpoint52f_pre, checkpoint55d_post, checkpoint58e_post, checkpoint54a_pre, checkpoint53c_post, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, branchpoint-genmake2, checkpoint54a_post, checkpoint55h_post, checkpoint58n_post, checkpoint51r_post, checkpoint51i_post, checkpoint57e_post, checkpoint55b_post, checkpoint51b_post, checkpoint51c_post, checkpoint53a_post, checkpoint55f_post, checkpoint52d_post, checkpoint53g_post, checkpoint57p_post, checkpint57u_post, checkpoint50g_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint56a_post, checkpoint58l_post, checkpoint53f_post, checkpoint57h_done, checkpoint52j_post, checkpoint50e_post, checkpoint57j_post, checkpoint57f_pre, checkpoint58g_post, branch-netcdf, checkpoint52l_post, checkpoint52n_post, checkpoint53b_pre, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint51e_post, checkpoint57a_pre, checkpoint55a_post, checkpoint49, checkpoint57o_post, checkpoint51o_post, checkpoint57k_post, checkpoint51f_pre, checkpoint53b_post, checkpoint52a_post, checkpoint57w_post, checkpoint58i_post, checkpoint51g_post, ecco_c52_e35, checkpoint57x_post, checkpoint50b_post, checkpoint58c_post, checkpoint58u_post, checkpoint51m_post, checkpoint53d_pre, checkpoint58s_post, checkpoint55e_post, checkpoint54c_post, checkpoint51a_post, checkpoint51p_post, checkpoint48g_post, checkpoint51u_post
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch 
r* Coord.: needs to apply OBC to rStar_Factor_W,S before updating hFacW,S
1 C $Header: $
2 C $Name: $
3
4 #include "OBCS_OPTIONS.h"
5
6 SUBROUTINE OBCS_APPLY_R_STAR(
7 I bi, bj, ksurfC, ksurfW, ksurfS,
8 U rStarFldC, rStarFldW, rStarFldS,
9 I myTime, myIter, myThid )
10 C *==========================================================*
11 C | S/R OBCS_APPLY_R_STAR
12 C *==========================================================*
13 IMPLICIT NONE
14 C == Global variables ==
15 #include "SIZE.h"
16 #include "EEPARAMS.h"
17 #include "PARAMS.h"
18 #include "GRID.h"
19 #include "OBCS.h"
20
21 C == Routine Arguments ==
22 INTEGER bi,bj
23 INTEGER ksurfC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
24 INTEGER ksurfW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
25 INTEGER ksurfS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
26 _RL rStarFldC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
27 _RL rStarFldW(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
28 _RL rStarFldS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
29 _RL myTime
30 INTEGER myIter, myThid
31
32 #ifdef ALLOW_OBCS
33 #ifdef NONLIN_FRSURF
34
35 C == Local variables ==
36 INTEGER i,j
37
38 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
39
40 C- Set model rStar_Factor to OB values on North/South Boundaries
41 DO i=1-Olx,sNx+Olx
42 C Northern boundary
43 j = OB_Jn(i,bi,bj)
44 IF (j.NE.0) THEN
45 rStarFldS(i,j,bi,bj) = 1. _d 0
46 IF (ksurfS(i,j,bi,bj).LE.Nr) THEN
47 rStarFldS(i,j,bi,bj) = 1. _d 0
48 & + OBNeta(i,bi,bj)
49 & *MAX(recip_Rcol(i,j-1,bi,bj),recip_Rcol(i,j,bi,bj))
50 ENDIF
51 ENDIF
52 C Southern boundary
53 j = OB_Js(i,bi,bj)
54 IF (j.NE.0) THEN
55 rStarFldS(i,j+1,bi,bj) = 1. _d 0
56 IF (ksurfS(i,j+1,bi,bj).LE.Nr) THEN
57 rStarFldS(i,j+1,bi,bj) = 1. _d 0
58 & + OBSeta(i,bi,bj)
59 & *MAX(recip_Rcol(i,j+1,bi,bj),recip_Rcol(i,j,bi,bj))
60 ENDIF
61 ENDIF
62 ENDDO
63
64 C- Set model rStar_Factor to OB values on East/West Boundaries
65 DO j=1-Oly,sNy+Oly
66 C Eastern boundary
67 i = OB_Ie(J,bi,bj)
68 IF (i.NE.0) THEN
69 rStarFldW(i,j,bi,bj) = 1. _d 0
70 IF (ksurfW(i,j,bi,bj).LE.Nr) THEN
71 rStarFldW(i,j,bi,bj) = 1. _d 0
72 & + OBEeta(i,bi,bj)
73 & *MAX(recip_Rcol(i-1,j,bi,bj),recip_Rcol(i,j,bi,bj))
74 ENDIF
75 ENDIF
76 C Western boundary
77 i = OB_Iw(j,bi,bj)
78 IF (i.NE.0) THEN
79 rStarFldW(i+1,j,bi,bj) = 1. _d 0
80 IF (ksurfW(i+1,j,bi,bj).LE.Nr) THEN
81 rStarFldW(i+1,j,bi,bj) = 1. _d 0
82 & + OBWeta(i,bi,bj)
83 & *MAX(recip_Rcol(i+1,j,bi,bj),recip_Rcol(i,j,bi,bj))
84 ENDIF
85 ENDIF
86 ENDDO
87
88 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
89
90 #endif /* NONLIN_FRSURF */
91 #endif /* ALLOW_OBCS */
92 RETURN
93 END
  ViewVC Help
Powered by ViewVC 1.1.22