79 |
ENDIF |
ENDIF |
80 |
|
|
81 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
82 |
#ifdef NONLIN_FRSURF |
C- note (with Non-Lin Free-Surface): |
83 |
C- note: 1) needs to apply OBC to etaH since viscous terms depend on hFacZ. |
C 1) needs to apply OBC to etaH since viscous terms depend on hFacZ. |
84 |
C that is not only function of boundaries hFac values. |
C that is not only function of boundaries hFac values. |
85 |
C 2) has to be done before calc_surf_dr; but since obcs_calc is |
C 2) has to be done before calc_surf_dr; but since obcs_calc is |
86 |
C called later, hFacZ will lag 1 time step behind OBC update. |
C called later, hFacZ will lag 1 time step behind OBC update. |
87 |
C 3) avoid also unrealistic value of etaH in OB regions that |
C 3) avoid also unrealistic value of etaH in OB regions that |
88 |
C might produce many "WARNING" message from calc_surf_dr. |
C might produce many "WARNING" message from calc_surf_dr. |
89 |
C------- |
C------- |
90 |
IF ( useOBCS .AND. nonlinFreeSurf.GT.0 ) |
C-- Apply OBC to etaH if NonLin-FreeSurf, reset to zero otherwise: |
91 |
& CALL OBCS_APPLY_ETA( bi, bj, etaH, myThid ) |
IF ( useOBCS ) CALL OBCS_APPLY_ETA( bi, bj, etaH, myThid ) |
|
#endif /* NONLIN_FRSURF */ |
|
92 |
#endif /* ALLOW_OBCS */ |
#endif /* ALLOW_OBCS */ |
93 |
|
|
94 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |