171 |
& ( rSurftmp(i,j) - MAX(rF(ks+1), R_low(i,j,bi,bj)) |
& ( rSurftmp(i,j) - MAX(rF(ks+1), R_low(i,j,bi,bj)) |
172 |
& )*recip_drF(ks)*maskC(i,j,ks,bi,bj) |
& )*recip_drF(ks)*maskC(i,j,ks,bi,bj) |
173 |
|
|
|
IF (hFac_surfC(i,j,bi,bj).GT.hFacSup) THEN |
|
174 |
C-- Usefull warning when hFac becomes very large: |
C-- Usefull warning when hFac becomes very large: |
175 |
|
IF ( numbWrite.LE.numbWrMax .AND. |
176 |
|
& hFac_surfC(i,j,bi,bj).GT.hFacSup ) THEN |
177 |
|
numbWrite = numbWrite + 1 |
178 |
WRITE(errorMessageUnit,'(2A,6I4,I10)') |
WRITE(errorMessageUnit,'(2A,6I4,I10)') |
179 |
& 'WARNING: hFacC > hFacSup at:', |
& 'WARNING: hFacC > hFacSup at:', |
180 |
& ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter |
& ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter |
181 |
WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)') |
WRITE(errorMessageUnit,'(A,2F10.6,1PE14.6)') |
182 |
& 'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj), |
& 'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj), |
183 |
& hFac_surfC(i,j,bi,bj), etaFld(i,j,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---------- |
|
184 |
ENDIF |
ENDIF |
185 |
|
C---------- |
186 |
ENDIF |
ENDIF |
187 |
|
|
188 |
ENDDO |
ENDDO |