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

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

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

Revision 1.9 - (show annotations) (download)
Thu Nov 15 15:55:42 2012 UTC (11 years, 5 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65o, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, HEAD
Changes since 1.8: +2 -1 lines 
adding tidal velocity forcing capability to obcs
 Modified Files:
  model/src/dynamics.F forward_step.F the_main_loop.F
  pkg/obcs/OBCS_FIELDS.h OBCS_OPTIONS.h OBCS_PARAMS.h
  OBCS_SEAICE.h obcs_apply_eta.F obcs_apply_r_star.F
  obcs_apply_surf_dr.F obcs_apply_ts.F obcs_apply_w.F
  obcs_calc.F obcs_check.F obcs_init_variables.F obcs_readparms.F
  verification/seaice_obcs/code/OBCS_OPTIONS.h
 Added Files:
  pkg/obcs/obcs_add_tides.F
  verification/seaice_obcs/input.tides/*
  verification/seaice_obcs/results/output.tides.txt
1 C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_apply_eta.F,v 1.8 2012/09/18 20:09:17 jmc Exp $
2 C $Name: $
3
4 #include "OBCS_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: OBCS_APPLY_ETA
8 C !INTERFACE:
9 SUBROUTINE OBCS_APPLY_ETA( bi, bj,
10 U etaFld,
11 I myThid )
12
13 C !DESCRIPTION:
14 C *==========================================================*
15 C | S/R OBCS_APPLY_ETA
16 C | Apply surface position anomaly (Eta) OB values
17 C | to corresponding field array
18 C *==========================================================*
19
20 C !USES:
21 IMPLICIT NONE
22 C == Global variables ==
23 #include "SIZE.h"
24 #include "EEPARAMS.h"
25 #include "PARAMS.h"
26 #include "GRID.h"
27 #include "SURFACE.h"
28 #include "OBCS_PARAMS.h"
29 #include "OBCS_GRID.h"
30 #include "OBCS_FIELDS.h"
31
32 C !INPUT/OUTPUT PARAMETERS:
33 C == Routine Arguments ==
34 C bi, bj :: indices of current tile
35 C etaFld :: surface r-position anomaly field
36 C myThid :: my Thread Id number
37 INTEGER bi, bj
38 _RL etaFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39 INTEGER myThid
40 CEOP
41
42 #ifdef NONLIN_FRSURF
43 C !LOCAL VARIABLES:
44 C == Local variables ==
45 INTEGER i, j
46 INTEGER Iobc, Jobc
47
48 c IF ( nonlinFreeSurf.GT.0 ) THEN
49
50 C Set model variable to OB values on North/South Boundaries
51 IF ( tileHasOBN(bi,bj) ) THEN
52 C Northern boundary
53 DO i=1-OLx,sNx+OLx
54 Jobc = OB_Jn(i,bi,bj)
55 IF ( Jobc.NE.OB_indexNone ) THEN
56 IF ( kSurfS(i,Jobc,bi,bj).LE.Nr )
57 & etaFld(i,Jobc,bi,bj) = OBNeta(i,bi,bj)
58 ENDIF
59 ENDDO
60 ENDIF
61 IF ( tileHasOBS(bi,bj) ) THEN
62 C Southern boundary
63 DO i=1-OLx,sNx+OLx
64 Jobc = OB_Js(i,bi,bj)
65 IF ( Jobc.NE.OB_indexNone ) THEN
66 IF ( kSurfS(i,Jobc+1,bi,bj).LE.Nr )
67 & etaFld(i,Jobc,bi,bj) = OBSeta(i,bi,bj)
68 ENDIF
69 ENDDO
70 ENDIF
71
72 C Set model variables to OB values on East/West Boundaries
73 IF ( tileHasOBE(bi,bj) ) THEN
74 C Eastern boundary
75 DO j=1-OLy,sNy+OLy
76 Iobc = OB_Ie(j,bi,bj)
77 IF ( Iobc.NE.OB_indexNone ) THEN
78 IF ( kSurfW(Iobc,j,bi,bj).LE.Nr )
79 & etaFld(Iobc,j,bi,bj) = OBEeta(j,bi,bj)
80 ENDIF
81 ENDDO
82 ENDIF
83 IF ( tileHasOBW(bi,bj) ) THEN
84 C Western boundary
85 DO j=1-OLy,sNy+OLy
86 Iobc = OB_Iw(j,bi,bj)
87 IF ( Iobc.NE.OB_indexNone ) THEN
88 IF ( kSurfW(Iobc+1,j,bi,bj).LE.Nr )
89 & etaFld(Iobc,j,bi,bj) = OBWeta(j,bi,bj)
90 ENDIF
91 ENDDO
92 ENDIF
93
94 c ENDIF
95 #endif /* NONLIN_FRSURF */
96
97 RETURN
98 END
  ViewVC Help
Powered by ViewVC 1.1.22