/[MITgcm]/MITgcm/model/src/update_etah.F
ViewVC logotype

Diff of /MITgcm/model/src/update_etah.F

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

revision 1.11 by jmc, Wed Apr 27 22:20:22 2011 UTC revision 1.12 by jmc, Fri May 20 16:29:49 2011 UTC
# Line 11  C     !INTERFACE: Line 11  C     !INTERFACE:
11  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
12  C     *==========================================================*  C     *==========================================================*
13  C     | SUBROUTINE UPDATE_ETAH  C     | SUBROUTINE UPDATE_ETAH
14  C     | o Update etaH at the begining of the time step.  C     | o Update etaH after mom-correction-step/integr_continuity
15  C     |  (required with NLFS to derive surface layer thickness)  C     |  (required with NLFS to derive surface layer thickness)
16    C     | o Also derive SSH tendency at grid-cell Western and
17    C     |   Southern edges (for hybrid sigma-coordinate)
18  C     *==========================================================*  C     *==========================================================*
19  C     \ev  C     \ev
20    
# Line 79  C     Incorporate the Explicit part of - Line 81  C     Incorporate the Explicit part of -
81          ENDIF          ENDIF
82    
83  #ifdef ALLOW_OBCS  #ifdef ALLOW_OBCS
84  C- note (with Non-Lin Free-Surface):  C--    Apply OBC to etaH (NonLin-FreeSurf): needed since viscous terms
85  C        1) needs to apply OBC to etaH since viscous terms depend on hFacZ.  C       depend on hFacZ which is not only function of boundary hFac values.
86  C           that is not only function of boundaries hFac values.         IF ( useOBCS.AND.nonlinFreeSurf.GT.0 )
87  C        2) has to be done before calc_surf_dr; but since obcs_calc is       &    CALL OBCS_APPLY_ETA( bi, bj, etaH, myThid )
 C           called later, hFacZ will lag 1 time step behind OBC update.  
 C        3) avoid also unrealistic value of etaH in OB regions that  
 C           might produce many "WARNING" message from calc_surf_dr.  
 C-------  
 C--    Apply OBC to etaH if NonLin-FreeSurf, reset to zero otherwise:  
        IF ( useOBCS ) CALL OBCS_APPLY_ETA( bi, bj, etaH, myThid )  
88  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
89    
90  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
# Line 97  C- end bi,bj loop. Line 93  C- end bi,bj loop.
93         ENDDO         ENDDO
94        ENDDO        ENDDO
95    
96        IF (implicDiv2Dflow .NE. 1. _d 0 .OR. useOBCS )        IF ( implicDiv2Dflow .NE. 1. _d 0 .OR.
97         &    ( useOBCS.AND.nonlinFreeSurf.GT.0 ) )
98       &    CALL EXCH_XY_RL( etaH, myThid )       &    CALL EXCH_XY_RL( etaH, myThid )
99    
100  c     IF (useRealFreshWaterFlux .AND. myTime.EQ.startTime)  c     IF (useRealFreshWaterFlux .AND. myTime.EQ.startTime)
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22