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

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

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

Revision 1.9 - (hide annotations) (download)
Thu Nov 15 15:55:42 2012 UTC (11 years, 6 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 dimitri 1.9 C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_apply_eta.F,v 1.8 2012/09/18 20:09:17 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "OBCS_OPTIONS.h"
5    
6 jmc 1.6 CBOP
7     C !ROUTINE: OBCS_APPLY_ETA
8     C !INTERFACE:
9 jmc 1.1 SUBROUTINE OBCS_APPLY_ETA( bi, bj,
10 jmc 1.6 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 jmc 1.1 IMPLICIT NONE
22     C == Global variables ==
23     #include "SIZE.h"
24     #include "EEPARAMS.h"
25     #include "PARAMS.h"
26     #include "GRID.h"
27 jmc 1.4 #include "SURFACE.h"
28 dimitri 1.9 #include "OBCS_PARAMS.h"
29 jmc 1.7 #include "OBCS_GRID.h"
30     #include "OBCS_FIELDS.h"
31 jmc 1.1
32 jmc 1.6 C !INPUT/OUTPUT PARAMETERS:
33 jmc 1.1 C == Routine Arguments ==
34 jmc 1.6 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 jmc 1.1 _RL etaFld(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39     INTEGER myThid
40 jmc 1.6 CEOP
41 jmc 1.1
42 jmc 1.6 #ifdef NONLIN_FRSURF
43     C !LOCAL VARIABLES:
44 jmc 1.3 C == Local variables ==
45 jmc 1.6 INTEGER i, j
46     INTEGER Iobc, Jobc
47 jmc 1.1
48 jmc 1.6 c IF ( nonlinFreeSurf.GT.0 ) THEN
49 jmc 1.1
50 jmc 1.3 C Set model variable to OB values on North/South Boundaries
51 jmc 1.6 IF ( tileHasOBN(bi,bj) ) THEN
52 jmc 1.1 C Northern boundary
53 jmc 1.8 DO i=1-OLx,sNx+OLx
54 jmc 1.6 Jobc = OB_Jn(i,bi,bj)
55 jmc 1.8 IF ( Jobc.NE.OB_indexNone ) THEN
56 jmc 1.6 IF ( kSurfS(i,Jobc,bi,bj).LE.Nr )
57     & etaFld(i,Jobc,bi,bj) = OBNeta(i,bi,bj)
58 jmc 1.3 ENDIF
59 jmc 1.6 ENDDO
60     ENDIF
61     IF ( tileHasOBS(bi,bj) ) THEN
62 jmc 1.1 C Southern boundary
63 jmc 1.8 DO i=1-OLx,sNx+OLx
64 jmc 1.6 Jobc = OB_Js(i,bi,bj)
65 jmc 1.8 IF ( Jobc.NE.OB_indexNone ) THEN
66 jmc 1.6 IF ( kSurfS(i,Jobc+1,bi,bj).LE.Nr )
67     & etaFld(i,Jobc,bi,bj) = OBSeta(i,bi,bj)
68 jmc 1.3 ENDIF
69 jmc 1.6 ENDDO
70     ENDIF
71 jmc 1.1
72     C Set model variables to OB values on East/West Boundaries
73 jmc 1.6 IF ( tileHasOBE(bi,bj) ) THEN
74 jmc 1.1 C Eastern boundary
75 jmc 1.8 DO j=1-OLy,sNy+OLy
76 jmc 1.6 Iobc = OB_Ie(j,bi,bj)
77 jmc 1.8 IF ( Iobc.NE.OB_indexNone ) THEN
78 jmc 1.6 IF ( kSurfW(Iobc,j,bi,bj).LE.Nr )
79     & etaFld(Iobc,j,bi,bj) = OBEeta(j,bi,bj)
80 jmc 1.3 ENDIF
81 jmc 1.6 ENDDO
82     ENDIF
83     IF ( tileHasOBW(bi,bj) ) THEN
84 jmc 1.1 C Western boundary
85 jmc 1.8 DO j=1-OLy,sNy+OLy
86 jmc 1.6 Iobc = OB_Iw(j,bi,bj)
87 jmc 1.8 IF ( Iobc.NE.OB_indexNone ) THEN
88 jmc 1.6 IF ( kSurfW(Iobc+1,j,bi,bj).LE.Nr )
89     & etaFld(Iobc,j,bi,bj) = OBWeta(j,bi,bj)
90 jmc 1.3 ENDIF
91 jmc 1.6 ENDDO
92     ENDIF
93 jmc 1.1
94 jmc 1.6 c ENDIF
95 jmc 1.1 #endif /* NONLIN_FRSURF */
96 jmc 1.3
97 jmc 1.1 RETURN
98     END
  ViewVC Help
Powered by ViewVC 1.1.22