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

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

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

Revision 1.4.6.1 - (show annotations) (download)
Mon Apr 8 20:10:39 2002 UTC (22 years, 1 month ago) by heimbach
Branch: release1
CVS Tags: release1_p13_pre, release1_p13, release1_p8, release1_p9, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, release1_chkpt44d_post, release1_p12, release1_p10, release1_p11, release1_p16, release1_p17, release1_p14, release1_p15, release1_p12_pre
Branch point for: release1_50yr
Changes since 1.4: +16 -3 lines 
Changes encapsulated by checkpoint43a-release1mods and chkpt44d_post
on the main trunk.
These are:

 o added missing EXCLUDE_MONITOR flags
 o changed "e" to "_d" in gmredi_slope_limit, gmredi_slope_psi
   (incompatible typ in MIN/MAX expressions caused problems
   on IBM SP3)
 o in genmake added variable MAKEDEPEND
   plus resetting for case SunOS
 o added timer_stats.c routine for IBM SP3
 o removed variables in dynamics
 o real fresh water flux implemented with non-linear free-surface.
 o few fix (mask in shap_s2, EmPmR in external_field_load,
   USE_NATURAL_BCS in solve_for_P);
 o add arguments myIter & myTime to S/R obcs_calc & solve_for_P
 o merge of relevant stuff from the ecco-branch:
   - genmake: removed $S64 overwrite for case SunOS
   - pkg/exf: update and corrections for field swapping and obcs
   - pkg/ecco: parameter lists for the_model_main, the_main_loop
               harmonized between ECCO and MITgcm
   - pkg/autodiff: added flow directives for obcs, mdsio_gl_slice
                   updated checkpointing_lev... lists for obcs
   - model/src: minor changes in forward_step, plot_field
                added directive for divided adjoint in the_main_loop
   - pkg/mdsio: added mdsio_gl_slice
 o check parameters & config (chkpt44a_pre,post)
 o OBC and NonLin_FrSurf.
 o fix bug in mom_vi_del2uv
 o select when filters are applied ; add options to zonal_filter (data.zonfilt)
 o gmredi: fix Pb in the adiabatic form ; add options (.e.g. Bolus advection)
 o update AIM experiments (NCEP input files)
 o improve and extend diagnostics (Monitor, TimeAve with NonLin-FrSurf)
 o added some stuff for AD

These were merged with
cvs co -r release1 -P MITgcm
cd MITgcm
cvs update -kk
cvs update -j checkpoint43a-release1mods -j chkpt44d_post -d -P -kk
1 C $Header$
2 C $Name$
3
4 #include "OBCS_OPTIONS.h"
5
6 SUBROUTINE OBCS_CALC( bi, bj, futureTime, futureIter,
7 & uVel, vVel, wVel, theta, salt,
8 & myThid )
9 C /==========================================================\
10 C | SUBROUTINE OBCS_CALC |
11 C | o Calculate future boundary data at open boundaries |
12 C | at time = futureTime |
13 C |==========================================================|
14 C | |
15 C \==========================================================/
16 IMPLICIT NONE
17
18 C === Global variables ===
19 #include "SIZE.h"
20 #include "EEPARAMS.h"
21 #include "PARAMS.h"
22 #include "OBCS.h"
23
24 C == Routine arguments ==
25 INTEGER bi, bj
26 INTEGER futureIter
27 _RL futureTime
28 _RL uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
29 _RL vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
30 _RL wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
31 _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
32 _RL salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
33 INTEGER myThid
34
35 #ifdef ALLOW_OBCS
36
37 C == Local variables ==
38 INTEGER I, J , K, I_obc, J_obc
39
40 C Eastern OB
41 IF (useOrlanskiEast) THEN
42 CALL ORLANSKI_EAST(
43 & bi, bj, futureTime,
44 & uVel, vVel, wVel, theta, salt,
45 & myThid )
46 ELSE
47 DO K=1,Nr
48 DO J=1-Oly,sNy+Oly
49 I_obc=OB_Ie(J,bi,bj)
50 IF (I_obc.ne.0) THEN
51 OBEu(J,K,bi,bj)=0.
52 OBEv(J,K,bi,bj)=0.
53 OBEt(J,K,bi,bj)=tRef(K)
54 OBEs(J,K,bi,bj)=sRef(K)
55 #ifdef ALLOW_NONHYDROSTATIC
56 OBEw(J,K,bi,bj)=0.
57 #endif
58 #ifdef NONLIN_FRSURF
59 OBEeta(J,bi,bj)=0.
60 #endif
61 ENDIF
62 ENDDO
63 ENDDO
64 ENDIF
65
66 C Western OB
67 IF (useOrlanskiWest) THEN
68 CALL ORLANSKI_WEST(
69 & bi, bj, futureTime,
70 & uVel, vVel, wVel, theta, salt,
71 & myThid )
72 ELSE
73 DO K=1,Nr
74 DO J=1-Oly,sNy+Oly
75 I_obc=OB_Iw(J,bi,bj)
76 IF (I_obc.ne.0) THEN
77 OBWu(J,K,bi,bj)=0.
78 OBWv(J,K,bi,bj)=0.
79 OBWt(J,K,bi,bj)=tRef(K)
80 OBWs(J,K,bi,bj)=sRef(K)
81 #ifdef ALLOW_NONHYDROSTATIC
82 OBWw(J,K,bi,bj)=0.
83 #endif
84 #ifdef NONLIN_FRSURF
85 OBWeta(J,bi,bj)=0.
86 #endif
87 ENDIF
88 ENDDO
89 ENDDO
90 ENDIF
91
92 C Northern OB
93 IF (useOrlanskiNorth) THEN
94 CALL ORLANSKI_NORTH(
95 & bi, bj, futureTime,
96 & uVel, vVel, wVel, theta, salt,
97 & myThid )
98 ELSE
99 DO K=1,Nr
100 DO I=1-Olx,sNx+Olx
101 J_obc=OB_Jn(I,bi,bj)
102 IF (J_obc.ne.0) THEN
103 OBNv(I,K,bi,bj)=0.
104 OBNu(I,K,bi,bj)=0.
105 OBNt(I,K,bi,bj)=tRef(K)
106 OBNs(I,K,bi,bj)=sRef(K)
107 #ifdef ALLOW_NONHYDROSTATIC
108 OBNw(I,K,bi,bj)=0.
109 #endif
110 #ifdef NONLIN_FRSURF
111 OBNeta(I,bi,bj)=0.
112 #endif
113 ENDIF
114 ENDDO
115 ENDDO
116 ENDIF
117
118 C Southern OB
119 IF (useOrlanskiSouth) THEN
120 CALL ORLANSKI_SOUTH(
121 & bi, bj, futureTime,
122 & uVel, vVel, wVel, theta, salt,
123 & myThid )
124 ELSE
125 DO K=1,Nr
126 DO I=1-Olx,sNx+Olx
127 J_obc=OB_Js(I,bi,bj)
128 IF (J_obc.ne.0) THEN
129 OBSu(I,K,bi,bj)=0.
130 OBSv(I,K,bi,bj)=0.
131 OBSt(I,K,bi,bj)=tRef(K)
132 OBSs(I,K,bi,bj)=sRef(K)
133 #ifdef ALLOW_NONHYDROSTATIC
134 OBSw(I,K,bi,bj)=0.
135 #endif
136 #ifdef NONLIN_FRSURF
137 OBSeta(I,bi,bj)=0.
138 #endif
139 ENDIF
140 ENDDO
141 ENDDO
142 ENDIF
143
144 #endif /* ALLOW_OBCS */
145 RETURN
146 END
  ViewVC Help
Powered by ViewVC 1.1.22