Merge pull request #467 from FESOM/refactoring_bugfix_allowsalinitybl…

…owup synchronize the salinity cutoff and blowup criteria for salinity…
FESOM · Jun 22, 2023 · b802a4f · b802a4f
2 parents 56acfe1 + a239529
commit b802a4f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/write_step_info.F90 b/src/write_step_info.F90
@@ -469,11 +469,11 @@ subroutine check_blowup(istep, ice, dynamics, tracers, partit, mesh)
           !_______________________________________________________________
           ! check salt
           if ( (tracers%data(2)%values(nz, n) /= tracers%data(2)%values(nz, n)) .or.  &
-             tracers%data(2)%values(nz, n) < 0 .or. tracers%data(2)%values(nz, n)>50 ) then
+             tracers%data(2)%values(nz, n) <=3.0_WP .or. tracers%data(2)%values(nz, n)>=45.0_WP ) then
 !$OMP CRITICAL
              found_blowup_loc=1
              write(*,*) '___CHECK FOR BLOW UP___________ --> mstep=',istep
-             write(*,*) ' --STOP--> found salinity becomes NaN or <0, >50'
+             write(*,*) ' --STOP--> found salinity becomes NaN or <=3.0, >=45.0'
              write(*,*) 'mype     = ',mype
              write(*,*) 'mstep    = ',istep
              write(*,*) 'node     = ',n