model/src/calc_surf_dr.F

C $Header: /u/gcmpack/MITgcm/model/src/calc_surf_dr.F,v 1.6 2002/02/10 20:04:11 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: CALC_SURF_DR
C     !INTERFACE:
      SUBROUTINE CALC_SURF_DR(etaFld,
     I                        myTime, myIter, myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE CALC_SURF_DR                                   
C     | o Calculate the new surface level thickness according to  
C     |   the surface r-position  (Non-Linear Free-Surf)          
C     | o take decision if grid box becomes too thin or too thick 
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SURFACE.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myTime :: Current time in simulation
C     myIter :: Current iteration number in simulation
C     myThid :: Thread number for this instance of the routine.
C     etaFld :: current eta field used to update the hFactor
      _RL myTime
      INTEGER myIter
      INTEGER myThid
      _RL etaFld(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)

#ifdef NONLIN_FRSURF

C     !LOCAL VARIABLES:
C     Local variables in common block
C     Rmin_surf :: minimum r_value of the free surface position 
C                  that satisfy  the hFacInf criteria
      COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
      _RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
C     Local variables
C     i,j,k,bi,bj  :: loop counter
C     rSurftmp     :: free surface r-position that is used to compute hFac_surf
C     adjust_nb_pt :: Nb of grid points where rSurf is adjusted (hFactInf)
C     adjust_volum :: adjustment effect on the volume (domain size)
C     numbWrite    :: count the Number of warning written on STD-ERR file
C     numbWrMax    ::  maximum  Number of warning written on STD-ERR file
      INTEGER i,j,bi,bj
      INTEGER ks, numbWrite, numbWrMax
      _RL hFacInfMOM, Rmin_tmp, hFactmp, adjust_nb_pt, adjust_volum
      _RL rSurftmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RS hhm, hhp
      CHARACTER*(MAX_LEN_MBUF) suff
CEOP
      DATA numbWrite / 0 /
      numbWrMax = Nx*Ny

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      IF (myIter.EQ.-1) THEN

       hFacInfMOM = hFacInf

       DO bj=myByLo(myThid), myByHi(myThid)
        DO bi=myBxLo(myThid), myBxHi(myThid)

C-- Initialise arrays :
         DO j=1-Oly,sNy+Oly
          DO i=1-Olx,sNx+Olx 
           hFac_surfC(i,j,bi,bj) = 0.
           hFac_surfW(i,j,bi,bj) = 0.
           hFac_surfS(i,j,bi,bj) = 0.
           PmEpR(i,j,bi,bj) = 0.
           Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj)
          ENDDO
         ENDDO

C-- Compute the mimimum value of r_surf (used for computing hFac_surfC)
         DO j=1,sNy
          DO i=1,sNx
           ks = ksurfC(i,j,bi,bj)
           IF (ks.LE.Nr) THEN
             Rmin_tmp = rF(ks+1)
             IF ( ks.EQ.ksurfW(i,j,bi,bj))
     &          Rmin_tmp = MAX(Rmin_tmp, R_low(i-1,j,bi,bj))
             IF ( ks.EQ.ksurfW(i+1,j,bi,bj)) 
     &          Rmin_tmp = MAX(Rmin_tmp, R_low(i+1,j,bi,bj))
             IF ( ks.EQ.ksurfS(i,j,bi,bj))
     &          Rmin_tmp = MAX(Rmin_tmp, R_low(i,j-1,bi,bj))
             IF ( ks.EQ.ksurfS(i,j+1,bi,bj))
     &          Rmin_tmp = MAX(Rmin_tmp, R_low(i,j+1,bi,bj))

             Rmin_surf(i,j,bi,bj) =
     &        MAX( MAX(rF(ks+1),R_low(i,j,bi,bj)) + hFacInf*drF(ks),  
     &                                Rmin_tmp + hFacInfMOM*drF(ks)
     &           )
           ENDIF
          ENDDO
         ENDDO

C-    end bi,bj loop.
        ENDDO
       ENDDO

       _EXCH_XY_R8( Rmin_surf, myThid ) 

C-    end of initialization block
      ENDIF

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      adjust_nb_pt = 0.
      adjust_volum = 0.

      DO bj=myByLo(myThid), myByHi(myThid)
       DO bi=myBxLo(myThid), myBxHi(myThid)

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C-- Compute the new fractional thickness of surface level (ksurfC):

        DO j=0,sNy+1
         DO i=0,sNx+1
          rSurftmp(i,j) = Ro_surf(i,j,bi,bj)+etaFld(i,j,bi,bj)
          ks = ksurfC(i,j,bi,bj)
          IF (ks.LE.Nr) THEN
           IF (rSurftmp(i,j) .LT. Rmin_surf(i,j,bi,bj)) THEN
C-- Needs to do something :
            IF (numbWrite.LE.numbWrMax) THEN
             numbWrite = numbWrite + 1 
             hFactmp = ( rSurftmp(i,j)-MAX(rF(ks+1),R_low(i,j,bi,bj))
     &                 )*recip_drF(ks)
             IF (hFactmp.LT.hFacInf) THEN
              WRITE(errorMessageUnit,'(2A,6I4,I10)') 
     &         'WARNING: hFacC < hFacInf at:',
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
             ELSE
              WRITE(errorMessageUnit,'(2A,6I4,I10)') 
     &         'WARNING: hFac < hFacInf near:',
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
             ENDIF
              WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)')
     &         'hFac_n-1,hFac_n,eta =',
     &          hfacC(i,j,ks,bi,bj), hFactmp, etaFld(i,j,bi,bj)
            ENDIF
C-- Decide to STOP :
c             WRITE(errorMessageUnit,'(A)') 
c    &        'STOP in CALC_SURF_DR : too SMALL hFac !'
c             STOP 'ABNORMAL END: S/R CALC_SURF_DR' 
C----------

C-- Continue with Rmin_surf:
             IF ( i.GE.1.AND.i.LE.sNx .AND.
     &            j.GE.1.AND.j.LE.sNy ) THEN
               adjust_nb_pt = adjust_nb_pt + 1.
               adjust_volum = adjust_volum 
     &          + rA(i,j,bi,bj)*(Rmin_surf(i,j,bi,bj)-rSurftmp(i,j))
             ENDIF
             rSurftmp(i,j) = Rmin_surf(i,j,bi,bj)
C----------
           ENDIF

C-- Set hFac_surfC :
           hFac_surfC(i,j,bi,bj) = 
     &         ( rSurftmp(i,j) - MAX(rF(ks+1), R_low(i,j,bi,bj))
     &         )*recip_drF(ks)*maskC(i,j,ks,bi,bj)

           IF (hFac_surfC(i,j,bi,bj).GT.hFacSup) THEN
C-- Usefull warning when hFac becomes very large:
              WRITE(errorMessageUnit,'(2A,6I4,I10)') 
     &         'WARNING: hFacC > hFacSup at:',
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
              WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)') 
     &         'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj),
     &          hFac_surfC(i,j,bi,bj), etaFld(i,j,bi,bj)
C-- Decide to STOP :
c             WRITE(errorMessageUnit,'(A)') 
c    &         'STOP in CALC_SURF_DR : too LARGE hFac !'
c             STOP 'ABNORMAL END: S/R CALC_SURF_DR' 
C----------
           ENDIF
          ENDIF

         ENDDO
        ENDDO

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C-- Compute fractional thickness of surface level, for U & V point:

        DO j=1,sNy
         DO i=1,sNx+1
          ks = ksurfW(i,j,bi,bj)
          IF (ks.LE.Nr) THEN
            hhm = rF(ks)
            IF(ks.EQ.ksurfC(i-1,j,bi,bj)) hhm = rSurftmp(i-1,j)
            hhp = rF(ks)
            IF(ks.EQ.ksurfC(i,j,bi,bj))   hhp = rSurftmp(i,j)  
            hFac_surfW(i,j,bi,bj) = 
     &         ( MIN(hhm,hhp) 
     &          - MAX(rF(ks+1),R_low(i-1,j,bi,bj),R_low(i,j,bi,bj)) 
     &         )*recip_drF(ks)*maskW(i,j,ks,bi,bj)
          ENDIF
         ENDDO
        ENDDO

        DO j=1,sNy+1
         DO i=1,sNx
          ks = ksurfS(i,j,bi,bj)
          IF (ks.LE.Nr) THEN
            hhm = rF(ks)
            IF(ks.EQ.ksurfC(i,j-1,bi,bj)) hhm = rSurftmp(i,j-1)
            hhp = rF(ks)
            IF(ks.EQ.ksurfC(i,j,bi,bj))   hhp = rSurftmp(i,j)
            hFac_surfS(i,j,bi,bj) = 
     &         ( MIN(hhm,hhp) 
     &          - MAX(rF(ks+1),R_low(i,j-1,bi,bj),R_low(i,j,bi,bj)) 
     &         )*recip_drF(ks)*maskS(i,j,ks,bi,bj)
          ENDIF
         ENDDO
        ENDDO

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

C-    end bi,bj loop.
       ENDDO
      ENDDO

C-- Global diagnostic :
      _GLOBAL_SUM_R8( adjust_nb_pt , myThid ) 
      _GLOBAL_SUM_R8( adjust_volum , myThid ) 
      IF (adjust_nb_pt .GE.1.) THEN
        _BEGIN_MASTER( myThid ) 
        WRITE(standardMessageUnit,'(2(A,I10),1PE16.8)') 
     &    ' SURF_ADJUSTMENT: Iter=', myIter,
     &    ' Nb_pts,Vol=', nint(adjust_nb_pt), adjust_volum
        _END_MASTER( myThid )
      ENDIF

      _EXCH_XY_R4(hFac_surfC, myThid ) 
      CALL EXCH_UV_XY_RS(hFac_surfW,hFac_surfS,.FALSE.,myThid)

      IF (useRealFreshWaterFlux .AND. myTime.EQ.startTime)
     & _EXCH_XY_R4( PmEpR, myThid )

C-----
C Note: testing ksurfW,S is equivalent to a full height mask 
C   ==> no need for applying the mask here.
C and with "partial thin wall" ==> mask could be applied in S/R UPDATE_SURF_DR
C-----

      WRITE(suff,'(I10.10)') myIter
c     CALL WRITE_FLD_XY_RS('hFac_surfC.',suff,hFac_surfC,
c    &                     myIter,myThid)                    

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
#endif /* NONLIN_FRSURF */

      RETURN
      END
1	jmc	1.7	C $Header: /u/gcmpack/MITgcm/model/src/calc_surf_dr.F,v 1.6 2002/02/10 20:04:11 jmc Exp $
2	jmc	1.1	C $Name: $
3
4			#include "CPP_OPTIONS.h"
5
6	cnh	1.3	CBOP
7			C !ROUTINE: CALC_SURF_DR
8			C !INTERFACE:
9	jmc	1.1	SUBROUTINE CALC_SURF_DR(etaFld,
10			I myTime, myIter, myThid )
11	cnh	1.3	C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| SUBROUTINE CALC_SURF_DR
14			C \| o Calculate the new surface level thickness according to
15			C \| the surface r-position (Non-Linear Free-Surf)
16			C \| o take decision if grid box becomes too thin or too thick
17			C ==========================================================
18			C \ev
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			#include "SURFACE.h"
28
29	cnh	1.3	C !INPUT/OUTPUT PARAMETERS:
30	jmc	1.1	C == Routine arguments ==
31	cnh	1.3	C myTime :: Current time in simulation
32			C myIter :: Current iteration number in simulation
33			C myThid :: Thread number for this instance of the routine.
34			C etaFld :: current eta field used to update the hFactor
35	jmc	1.1	_RL myTime
36			INTEGER myIter
37			INTEGER myThid
38			_RL etaFld(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
39
40			#ifdef NONLIN_FRSURF
41
42	cnh	1.3	C !LOCAL VARIABLES:
43	jmc	1.2	C Local variables in common block
44	cnh	1.3	C Rmin_surf :: minimum r_value of the free surface position
45	jmc	1.2	C that satisfy the hFacInf criteria
46			COMMON /LOCAL_CALC_SURF_DR/ Rmin_surf
47			_RL Rmin_surf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48	jmc	1.1	C Local variables
49	cnh	1.3	C i,j,k,bi,bj :: loop counter
50			C rSurftmp :: free surface r-position that is used to compute hFac_surf
51			C adjust_nb_pt :: Nb of grid points where rSurf is adjusted (hFactInf)
52			C adjust_volum :: adjustment effect on the volume (domain size)
53	jmc	1.7	C numbWrite :: count the Number of warning written on STD-ERR file
54			C numbWrMax :: maximum Number of warning written on STD-ERR file
55	jmc	1.6	INTEGER i,j,bi,bj
56	jmc	1.7	INTEGER ks, numbWrite, numbWrMax
57	jmc	1.2	_RL hFacInfMOM, Rmin_tmp, hFactmp, adjust_nb_pt, adjust_volum
58			_RL rSurftmp(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
59	jmc	1.1	_RS hhm, hhp
60			CHARACTER*(MAX_LEN_MBUF) suff
61	cnh	1.3	CEOP
62	jmc	1.7	DATA numbWrite / 0 /
63			numbWrMax = Nx*Ny
64	jmc	1.1
65			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
66
67	jmc	1.4	IF (myIter.EQ.-1) THEN
68	jmc	1.2
69			hFacInfMOM = hFacInf
70
71	jmc	1.1	DO bj=myByLo(myThid), myByHi(myThid)
72			DO bi=myBxLo(myThid), myBxHi(myThid)
73	jmc	1.2
74			C-- Initialise arrays :
75	jmc	1.1	DO j=1-Oly,sNy+Oly
76			DO i=1-Olx,sNx+Olx
77			hFac_surfC(i,j,bi,bj) = 0.
78			hFac_surfW(i,j,bi,bj) = 0.
79			hFac_surfS(i,j,bi,bj) = 0.
80	jmc	1.6	PmEpR(i,j,bi,bj) = 0.
81	jmc	1.2	Rmin_surf(i,j,bi,bj) = Ro_surf(i,j,bi,bj)
82	jmc	1.1	ENDDO
83			ENDDO
84	jmc	1.2
85			C-- Compute the mimimum value of r_surf (used for computing hFac_surfC)
86			DO j=1,sNy
87			DO i=1,sNx
88			ks = ksurfC(i,j,bi,bj)
89			IF (ks.LE.Nr) THEN
90			Rmin_tmp = rF(ks+1)
91			IF ( ks.EQ.ksurfW(i,j,bi,bj))
92			& Rmin_tmp = MAX(Rmin_tmp, R_low(i-1,j,bi,bj))
93			IF ( ks.EQ.ksurfW(i+1,j,bi,bj))
94			& Rmin_tmp = MAX(Rmin_tmp, R_low(i+1,j,bi,bj))
95			IF ( ks.EQ.ksurfS(i,j,bi,bj))
96			& Rmin_tmp = MAX(Rmin_tmp, R_low(i,j-1,bi,bj))
97			IF ( ks.EQ.ksurfS(i,j+1,bi,bj))
98			& Rmin_tmp = MAX(Rmin_tmp, R_low(i,j+1,bi,bj))
99
100			Rmin_surf(i,j,bi,bj) =
101			& MAX( MAX(rF(ks+1),R_low(i,j,bi,bj)) + hFacInf*drF(ks),
102			& Rmin_tmp + hFacInfMOM*drF(ks)
103			& )
104			ENDIF
105			ENDDO
106			ENDDO
107
108			C- end bi,bj loop.
109	jmc	1.1	ENDDO
110			ENDDO
111	jmc	1.2
112			_EXCH_XY_R8( Rmin_surf, myThid )
113
114			C- end of initialization block
115	jmc	1.1	ENDIF
116
117			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
118
119	jmc	1.2	adjust_nb_pt = 0.
120			adjust_volum = 0.
121
122	jmc	1.1	DO bj=myByLo(myThid), myByHi(myThid)
123			DO bi=myBxLo(myThid), myBxHi(myThid)
124	jmc	1.2
125			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
126			C-- Compute the new fractional thickness of surface level (ksurfC):
127
128			DO j=0,sNy+1
129			DO i=0,sNx+1
130			rSurftmp(i,j) = Ro_surf(i,j,bi,bj)+etaFld(i,j,bi,bj)
131	jmc	1.1	ks = ksurfC(i,j,bi,bj)
132			IF (ks.LE.Nr) THEN
133	jmc	1.2	IF (rSurftmp(i,j) .LT. Rmin_surf(i,j,bi,bj)) THEN
134			C-- Needs to do something :
135	jmc	1.7	IF (numbWrite.LE.numbWrMax) THEN
136			numbWrite = numbWrite + 1
137	jmc	1.2	hFactmp = ( rSurftmp(i,j)-MAX(rF(ks+1),R_low(i,j,bi,bj))
138			& )*recip_drF(ks)
139	jmc	1.1	IF (hFactmp.LT.hFacInf) THEN
140	jmc	1.7	WRITE(errorMessageUnit,'(2A,6I4,I10)')
141			& 'WARNING: hFacC < hFacInf at:',
142	jmc	1.1	& ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
143	jmc	1.2	ELSE
144	jmc	1.7	WRITE(errorMessageUnit,'(2A,6I4,I10)')
145			& 'WARNING: hFac < hFacInf near:',
146	jmc	1.2	& ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
147	jmc	1.1	ENDIF
148	jmc	1.7	WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)')
149			& 'hFac_n-1,hFac_n,eta =',
150	jmc	1.2	& hfacC(i,j,ks,bi,bj), hFactmp, etaFld(i,j,bi,bj)
151	jmc	1.7	ENDIF
152	jmc	1.2	C-- Decide to STOP :
153	jmc	1.7	c WRITE(errorMessageUnit,'(A)')
154			c & 'STOP in CALC_SURF_DR : too SMALL hFac !'
155	jmc	1.2	c STOP 'ABNORMAL END: S/R CALC_SURF_DR'
156			C----------
157
158			C-- Continue with Rmin_surf:
159			IF ( i.GE.1.AND.i.LE.sNx .AND.
160			& j.GE.1.AND.j.LE.sNy ) THEN
161			adjust_nb_pt = adjust_nb_pt + 1.
162			adjust_volum = adjust_volum
163			& + rA(i,j,bi,bj)*(Rmin_surf(i,j,bi,bj)-rSurftmp(i,j))
164			ENDIF
165			rSurftmp(i,j) = Rmin_surf(i,j,bi,bj)
166			C----------
167			ENDIF
168
169			C-- Set hFac_surfC :
170			hFac_surfC(i,j,bi,bj) =
171			& ( rSurftmp(i,j) - MAX(rF(ks+1), R_low(i,j,bi,bj))
172			& )recip_drF(ks)maskC(i,j,ks,bi,bj)
173
174			IF (hFac_surfC(i,j,bi,bj).GT.hFacSup) THEN
175			C-- Usefull warning when hFac becomes very large:
176	jmc	1.7	WRITE(errorMessageUnit,'(2A,6I4,I10)')
177			& 'WARNING: hFacC > hFacSup at:',
178	jmc	1.1	& ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
179	jmc	1.7	WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)')
180			& 'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj),
181			& hFac_surfC(i,j,bi,bj), etaFld(i,j,bi,bj)
182	jmc	1.2	C-- Decide to STOP :
183	jmc	1.7	c WRITE(errorMessageUnit,'(A)')
184			c & 'STOP in CALC_SURF_DR : too LARGE hFac !'
185	jmc	1.2	c STOP 'ABNORMAL END: S/R CALC_SURF_DR'
186	jmc	1.1	C----------
187			ENDIF
188			ENDIF
189	jmc	1.2
190	jmc	1.1	ENDDO
191			ENDDO
192
193			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
194			C-- Compute fractional thickness of surface level, for U & V point:
195
196			DO j=1,sNy
197			DO i=1,sNx+1
198			ks = ksurfW(i,j,bi,bj)
199			IF (ks.LE.Nr) THEN
200	jmc	1.2	hhm = rF(ks)
201			IF(ks.EQ.ksurfC(i-1,j,bi,bj)) hhm = rSurftmp(i-1,j)
202			hhp = rF(ks)
203			IF(ks.EQ.ksurfC(i,j,bi,bj)) hhp = rSurftmp(i,j)
204			hFac_surfW(i,j,bi,bj) =
205			& ( MIN(hhm,hhp)
206			& - MAX(rF(ks+1),R_low(i-1,j,bi,bj),R_low(i,j,bi,bj))
207			& )recip_drF(ks)maskW(i,j,ks,bi,bj)
208	jmc	1.1	ENDIF
209			ENDDO
210			ENDDO
211	jmc	1.2
212	jmc	1.1	DO j=1,sNy+1
213			DO i=1,sNx
214			ks = ksurfS(i,j,bi,bj)
215			IF (ks.LE.Nr) THEN
216	jmc	1.2	hhm = rF(ks)
217			IF(ks.EQ.ksurfC(i,j-1,bi,bj)) hhm = rSurftmp(i,j-1)
218			hhp = rF(ks)
219			IF(ks.EQ.ksurfC(i,j,bi,bj)) hhp = rSurftmp(i,j)
220			hFac_surfS(i,j,bi,bj) =
221			& ( MIN(hhm,hhp)
222			& - MAX(rF(ks+1),R_low(i,j-1,bi,bj),R_low(i,j,bi,bj))
223			& )recip_drF(ks)maskS(i,j,ks,bi,bj)
224	jmc	1.1	ENDIF
225			ENDDO
226			ENDDO
227	jmc	1.2
228			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
229
230			C- end bi,bj loop.
231	jmc	1.1	ENDDO
232			ENDDO
233
234	jmc	1.2	C-- Global diagnostic :
235			_GLOBAL_SUM_R8( adjust_nb_pt , myThid )
236			_GLOBAL_SUM_R8( adjust_volum , myThid )
237			IF (adjust_nb_pt .GE.1.) THEN
238			_BEGIN_MASTER( myThid )
239	jmc	1.7	WRITE(standardMessageUnit,'(2(A,I10),1PE16.8)')
240			& ' SURF_ADJUSTMENT: Iter=', myIter,
241			& ' Nb_pts,Vol=', nint(adjust_nb_pt), adjust_volum
242			_END_MASTER( myThid )
243	jmc	1.2	ENDIF
244
245			_EXCH_XY_R4(hFac_surfC, myThid )
246	jmc	1.1	CALL EXCH_UV_XY_RS(hFac_surfW,hFac_surfS,.FALSE.,myThid)
247	jmc	1.5
248			IF (useRealFreshWaterFlux .AND. myTime.EQ.startTime)
249			& _EXCH_XY_R4( PmEpR, myThid )
250	jmc	1.1
251	jmc	1.2	C-----
252			C Note: testing ksurfW,S is equivalent to a full height mask
253			C ==> no need for applying the mask here.
254			C and with "partial thin wall" ==> mask could be applied in S/R UPDATE_SURF_DR
255			C-----
256	jmc	1.1
257			WRITE(suff,'(I10.10)') myIter
258			c CALL WRITE_FLD_XY_RS('hFac_surfC.',suff,hFac_surfC,
259			c & myIter,myThid)
260
261			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
262			#endif /* NONLIN_FRSURF */
263
264			RETURN
265			END