Changeset 14129 for NEMO/branches/2020/r14118_ticket2377/src/OCE/stpctl.F90
- Timestamp:
- 2020-12-08T11:30:38+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r14118_ticket2377/src/OCE/stpctl.F90
r14053 r14129 26 26 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 27 27 USE lib_mpp ! distributed memory computing 28 USE eosbn2 29 ! 28 30 USE netcdf ! NetCDF library 29 30 31 IMPLICIT NONE 31 32 PRIVATE … … 65 66 INTEGER , DIMENSION(9) :: iareasum, iareamin, iareamax 66 67 INTEGER , DIMENSION(3,4) :: iloc ! min/max loc indices 67 REAL(wp) :: zzz 68 REAL(wp) :: zzz, zminsal, zmaxsal ! local real 68 69 REAL(wp), DIMENSION(9) :: zmax, zmaxlocal 69 70 LOGICAL :: ll_wrtstp, ll_colruns, ll_wrtruns, ll_0oce … … 71 72 CHARACTER(len=20) :: clname 72 73 !!---------------------------------------------------------------------- 73 !74 74 IF( nstop > 0 .AND. ngrdstop > -1 ) RETURN ! stpctl was already called by a child grid 75 75 ! … … 180 180 END DO 181 181 IF( kt == nitend ) istatus = NF90_CLOSE(nrunid) 182 END IF182 END IF 183 183 ! !== error handling ==! 184 184 ! !== done by all processes at every time step ==! 185 185 ! 186 IF ( ln_seos.AND.(rn_lambda2==0._wp) ) THEN ! Discard checks on salinity 187 zmaxsal = +1.e38 ! if not used in eos 188 zminsal = -1.e38 189 ELSE 190 zmaxsal = 100._wp 191 zminsal = 0._wp 192 ENDIF 193 ! 186 194 IF( zmax(1) > 20._wp .OR. & ! too large sea surface height ( > 20 m ) 187 195 & zmax(2) > 10._wp .OR. & ! too large velocity ( > 10 m/s) 188 & zmax(3) <= 0._wp.OR. & ! negative or zero sea surface salinity189 & zmax(4) >= 100._wp.OR. & ! too large sea surface salinity ( > 100 )190 & zmax(4) < 0._wp.OR. & ! too large sea surface salinity (keep this line for sea-ice)196 & zmax(3) <= zminsal .OR. & ! negative or zero sea surface salinity 197 & zmax(4) >= zmaxsal .OR. & ! too large sea surface salinity ( > 100 ) 198 & zmax(4) < zminsal .OR. & ! too large sea surface salinity (keep this line for sea-ice) 191 199 & ISNAN( zmax(1) + zmax(2) + zmax(3) ) .OR. & ! NaN encounter in the tests 192 200 & ABS( zmax(1) + zmax(2) + zmax(3) ) > HUGE(1._wp) ) THEN ! Infinity encounter in the tests
Note: See TracChangeset
for help on using the changeset viewer.