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 7901 for branches/2017/dev_r7832_HPC08_lbclnk_3rd_dim/NEMOGCM/NEMO/OPA_SRC/stpctl.F90 – NEMO

Ignore:
Timestamp:
2017-04-13T05:46:23+02:00 (7 years ago)
Author:
gm
Message:

#1880: (HPC-08) minor corrections

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7832_HPC08_lbclnk_3rd_dim/NEMOGCM/NEMO/OPA_SRC/stpctl.F90

    r7897 r7901  
    1010   !!            2.0  ! 2009-07  (G. Madec)  Add statistic for time-spliting 
    1111   !!            3.7  ! 2016-09  (G. Madec)  Remove solver 
    12    !!            4.0  ! 2017-04  (G. Madec)  regroup  
     12   !!            4.0  ! 2017-04  (G. Madec)  regroup global communications 
    1313   !!---------------------------------------------------------------------- 
    1414 
     
    4444      !!              - Print it each 50 time steps 
    4545      !!              - Stop the run IF problem encountered by setting indic=-3 
    46       !!                Problems checked: U max>10 m/s and SSS min < 0 
     46      !!                Problems checked: |U| and |ssh| maximum larger than 10 m/s  
     47      !!                                  sea surface salinity (SSS) minimum < 0 
    4748      !! 
    4849      !! ** Actions :   'time.step' file containing the last ocean time-step 
     
    5455      INTEGER  ::   ji, jj, jk             ! dummy loop indices 
    5556      INTEGER  ::   ii, ij, ik             ! local integers 
     57      REAL(wp) ::   zzt                    ! local real  
    5658      INTEGER , DIMENSION(3) ::   ilocu    !  
    5759      INTEGER , DIMENSION(2) ::   ilocs    !  
     
    6971      ENDIF 
    7072      ! 
    71       IF(lwp) WRITE ( numstp, '(1x, i8)' )   kt    !==  current time step  ==!   ("time.step" file) 
    72       IF(lwp) REWIND( numstp ) 
     73      IF(lwp) THEN                        !==  current time step  ==!   ("time.step" file) 
     74         WRITE ( numstp, '(1x, i8)' )   kt 
     75         REWIND( numstp ) 
     76      ENDIF 
    7377      ! 
    7478      !                                            !==  test of extrema  ==! 
     
    8084      ! 
    8185      IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 
    82          WRITE(numout,*) ' ==>> time-step= ',kt,' abs(U) max: ', zmax(1), ' SSS min:', - zmax(2) 
     86         WRITE(numout,*) ' ==>> time-step= ',kt,' |U| max: ', zmax(1), ' SSS min:', - zmax(2) 
    8387      ENDIF 
    8488      ! 
     
    106110      IF( -zmax(2) < 0._wp ) THEN                     !* negative salinity 
    107111         IF( lk_mpp ) THEN 
    108             CALL mpp_minloc( tsn(:,:,1,jp_sal),tmask(:,:,1), - zmax(2), ii, ij ) 
     112            CALL mpp_minloc( tsn(:,:,1,jp_sal),tmask(:,:,1), zzt, ii, ij ) 
    109113         ELSE 
    110114            ilocs = MINLOC( tsn(:,:,1,jp_sal), mask = tmask(:,:,1) == 1._wp ) 
Note: See TracChangeset for help on using the changeset viewer.