New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10415 for NEMO/trunk/src/OCE – NEMO

Changeset 10415 for NEMO/trunk/src/OCE


Ignore:
Timestamp:
2018-12-19T12:28:25+01:00 (5 years ago)
Author:
clem
Message:

correct the advection for ice ponds, make a couple of cosmetic changes and fix the conservation issue for sea ice advection (just a bugged control print)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/stpctl.F90

    r10364 r10415  
    5050      !!              - Print it each 50 time steps 
    5151      !!              - Stop the run IF problem encountered by setting indic=-3 
    52       !!                Problems checked: |ssh| maximum larger than 10 m 
     52      !!                Problems checked: |ssh| maximum larger than 20 m 
    5353      !!                                  |U|   maximum larger than 10 m/s  
    5454      !!                                  negative sea surface salinity 
     
    136136      ENDIF 
    137137      ! 
    138       IF (  zmax(1) >   15._wp .OR.   &                    ! too large sea surface height ( > 15 m ) 
     138      IF (  zmax(1) >   20._wp .OR.   &                    ! too large sea surface height ( > 20 m ) 
    139139         &  zmax(2) >   10._wp .OR.   &                    ! too large velocity ( > 10 m/s) 
    140140         &  zmax(3) >=   0._wp .OR.   &                    ! negative or zero sea surface salinity 
     
    159159         IF(lwp) THEN 
    160160            WRITE(numout,cform_err) 
    161             WRITE(numout,*) ' stp_ctl: |ssh| > 10 m  or  |U| > 10 m/s  or  S <= 0  or  S >= 100  or  NaN encounter in the tests' 
     161            WRITE(numout,*) ' stp_ctl: |ssh| > 20 m  or  |U| > 10 m/s  or  S <= 0  or  S >= 100  or  NaN encounter in the tests' 
    162162            WRITE(numout,*) ' ======= ' 
    163163            WRITE(numout,9100) kt,   zmax(1), iih , ijh 
Note: See TracChangeset for help on using the changeset viewer.