Ignore:
Timestamp:
2018-11-25T15:24:21+01:00 (22 months ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: action 5b: by default, suppress global communication in stpctl, see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/stpctl.F90

    r10314 r10358  
    3333 
    3434   INTEGER  ::   idrun, idtime, idssh, idu, ids1, ids2, istatus 
     35   LOGICAL  ::   lsomeoce 
    3536   !!---------------------------------------------------------------------- 
    3637   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    6061      INTEGER, INTENT(inout) ::   kindic   ! error indicator 
    6162      !! 
    62       INTEGER  ::   ji, jj, jk             ! dummy loop indices 
    63       INTEGER  ::   ih(2)                  ! local integers 
    64       INTEGER  ::   iu(3)                  !   -       - 
    65       INTEGER  ::   is1(3)                 !   -       - 
    66       INTEGER  ::   is2(3)                 !   -       - 
    67       REAL(wp) ::   zzz                    ! local real  
    68       INTEGER , DIMENSION(3) ::   ilocu, ilocs1, ilocs2 
    69       INTEGER , DIMENSION(2) ::   iloch 
     63      INTEGER                ::   ji, jj, jk          ! dummy loop indices 
     64      INTEGER, DIMENSION(2)  ::   ih                  ! min/max loc indices 
     65      INTEGER, DIMENSION(3)  ::   iu, is1, is2        ! min/max loc indices 
     66      REAL(wp)               ::   zzz                 ! local real  
    7067      REAL(wp), DIMENSION(5) ::   zmax 
    7168      CHARACTER(len=20) :: clname 
     
    7774         WRITE(numout,*) '~~~~~~~' 
    7875         !                                ! open time.step file 
    79          CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     76         IF( lwm ) CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    8077         !                                ! open run.stat file 
    81          CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    82  
    83          IF( lwm ) THEN 
     78         IF( ln_ctl .AND. lwm ) THEN 
     79            CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    8480            clname = 'run.stat.nc' 
    8581            IF( .NOT. Agrif_Root() )   clname = TRIM(Agrif_CFixed())//"_"//TRIM(clname) 
     
    9288            istatus = NF90_ENDDEF(idrun) 
    9389         ENDIF 
    94           
    9590      ENDIF 
     91      IF( kt == nit000 )   lsomeoce = COUNT( ssmask(:,:) == 1._wp ) > 0 
    9692      ! 
    97       IF(lwp) THEN                        !==  current time step  ==!   ("time.step" file) 
     93      IF(lwm) THEN                        !==  current time step  ==!   ("time.step" file) 
    9894         WRITE ( numstp, '(1x, i8)' )   kt 
    9995         REWIND( numstp ) 
     
    111107      zmax(5) = REAL( nstop , wp )                                            ! stop indicator 
    112108      ! 
    113       IF( lk_mpp ) THEN 
    114          CALL mpp_max( "stpctl", zmax )    ! max over the global domain 
    115          ! 
     109      IF( lk_mpp .AND. ln_ctl ) THEN 
     110         CALL mpp_max( "stpctl", zmax )          ! max over the global domain 
    116111         nstop = NINT( zmax(5) )                 ! nstop indicator sheared among all local domains 
    117112      ENDIF 
    118       ! 
    119       IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 
    120          WRITE(numout,*) ' ==>> time-step= ', kt, ' |ssh| max: ',   zmax(1), ' |U| max: ', zmax(2),   & 
    121             &                                     ' S min: '    , - zmax(3), ' S max: ', zmax(4) 
    122       ENDIF 
    123       ! 
    124       IF (  zmax(1) >   15._wp .OR.   &                    ! too large sea surface height ( > 15 m ) 
    125          &  zmax(2) >   10._wp .OR.   &                    ! too large velocity ( > 10 m/s) 
    126          &  zmax(3) >=   0._wp .OR.   &                    ! negative or zero sea surface salinity 
    127          &  zmax(4) >= 100._wp .OR.   &                    ! too large sea surface salinity ( > 100 ) 
    128          &  zmax(4) <    0._wp .OR.   &                    ! too large sea surface salinity (keep this line for sea-ice) 
    129          &  ISNAN( zmax(1) + zmax(2) + zmax(3) )  ) THEN   ! NaN encounter in the tests 
    130          IF( lk_mpp ) THEN 
    131             CALL mpp_maxloc( 'stpctl', ABS(sshn)        , ssmask(:,:)  , zzz, ih  ) 
    132             CALL mpp_maxloc( 'stpctl', ABS(un)          , umask (:,:,:), zzz, iu  ) 
    133             CALL mpp_minloc( 'stpctl', tsn(:,:,:,jp_sal), tmask (:,:,:), zzz, is1 ) 
    134             CALL mpp_maxloc( 'stpctl', tsn(:,:,:,jp_sal), tmask (:,:,:), zzz, is2 ) 
    135          ELSE 
    136             iloch  = MINLOC( ABS( sshn(:,:)   )                               ) 
    137             ilocu  = MAXLOC( ABS( un  (:,:,:) )                               ) 
    138             ilocs1 = MINLOC( tsn(:,:,:,jp_sal) , mask = tmask(:,:,:) == 1._wp ) 
    139             ilocs2 = MAXLOC( tsn(:,:,:,jp_sal) , mask = tmask(:,:,:) == 1._wp ) 
    140             ih(1)  = iloch (1) + nimpp - 1   ;   ih(2)  = iloch (2) + njmpp - 1 
    141             iu(1)  = ilocu (1) + nimpp - 1   ;   iu(2)  = ilocu (2) + njmpp - 1   ;   iu(3)  = ilocu (3) 
    142             is1(1) = ilocs1(1) + nimpp - 1   ;   is1(2) = ilocs1(2) + njmpp - 1   ;   is1(3) = ilocs1(3) 
    143             is2(1) = ilocs2(1) + nimpp - 1   ;   is2(2) = ilocs2(2) + njmpp - 1   ;   is2(3) = ilocs2(3) 
    144          ENDIF 
    145          IF(lwp) THEN 
    146             WRITE(numout,cform_err) 
    147             WRITE(numout,*) ' stp_ctl: |ssh| > 10 m  or  |U| > 10 m/s  or  S <= 0  or  S >= 100  or  NaN encounter in the tests' 
    148             WRITE(numout,*) ' ======= ' 
    149             WRITE(numout,9100) kt,   zmax(1), ih(1) , ih(2) 
    150             WRITE(numout,9200) kt,   zmax(2), iu(1) , iu(2) , iu(3) 
    151             WRITE(numout,9300) kt, - zmax(3), is1(1), is1(2), is1(3) 
    152             WRITE(numout,9400) kt,   zmax(4), is2(1), is2(2), is2(3) 
    153             WRITE(numout,*) 
    154             WRITE(numout,*) '          output of last computed fields in output.abort.nc file' 
    155          ENDIF 
    156          kindic = -3 
    157          ! 
    158          nstop = nstop + 1                            ! increase nstop by 1 (on all local domains) 
    159          CALL dia_wri_state( 'output.abort', kt )     ! create an output.abort file 
    160          ! 
    161       ENDIF 
    162 9100  FORMAT (' kt=',i8,'   |ssh| max: ',1pg11.4,', at  i j  : ',2i5) 
    163 9200  FORMAT (' kt=',i8,'   |U|   max: ',1pg11.4,', at  i j k: ',3i5) 
    164 9300  FORMAT (' kt=',i8,'   S     min: ',1pg11.4,', at  i j k: ',3i5) 
    165 9400  FORMAT (' kt=',i8,'   S     max: ',1pg11.4,', at  i j k: ',3i5) 
    166       ! 
    167       !                                            !==  run statistics  ==!   ("run.stat" file) 
    168       IF(lwp) WRITE(numrun,9500) kt, zmax(1), zmax(2), -zmax(3), zmax(4) 
    169       IF( lwm ) THEN 
     113      !                                   !==  run statistics  ==!   ("run.stat" files) 
     114      IF( ln_ctl .AND. lwm ) THEN 
     115         WRITE(numrun,9500) kt, zmax(1), zmax(2), -zmax(3), zmax(4) 
    170116         istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) ) 
    171117         istatus = NF90_PUT_VAR( idrun,   idu, (/ zmax(2)/), (/kt/), (/1/) ) 
     
    175121         IF( kt == nitend         ) istatus = NF90_CLOSE(idrun) 
    176122      END IF 
     123      !                                   !==  error handling  ==! 
     124      IF( ( ln_ctl .OR. lsomeoce ) .AND. (   &             ! have use mpp_max (because ln_ctl=.T.) or contains some ocean points 
     125         &  zmax(1) >   15._wp .OR.   &                    ! too large sea surface height ( > 15 m ) 
     126         &  zmax(2) >   10._wp .OR.   &                    ! too large velocity ( > 10 m/s) 
     127         &  zmax(3) >=   0._wp .OR.   &                    ! negative or zero sea surface salinity 
     128         &  zmax(4) >= 100._wp .OR.   &                    ! too large sea surface salinity ( > 100 ) 
     129         &  zmax(4) <    0._wp .OR.   &                    ! too large sea surface salinity (keep this line for sea-ice) 
     130         &  ISNAN( zmax(1) + zmax(2) + zmax(3) ) ) ) THEN   ! NaN encounter in the tests 
     131         IF( lk_mpp .AND. ln_ctl ) THEN 
     132            CALL mpp_maxloc( 'stpctl', ABS(sshn)        , ssmask(:,:)  , zzz, ih  ) 
     133            CALL mpp_maxloc( 'stpctl', ABS(un)          , umask (:,:,:), zzz, iu  ) 
     134            CALL mpp_minloc( 'stpctl', tsn(:,:,:,jp_sal), tmask (:,:,:), zzz, is1 ) 
     135            CALL mpp_maxloc( 'stpctl', tsn(:,:,:,jp_sal), tmask (:,:,:), zzz, is2 ) 
     136         ELSE 
     137            ih(:)  = MAXLOC( ABS( sshn(:,:)   )                              ) + (/ nimpp - 1, njmpp - 1    /) 
     138            iu(:)  = MAXLOC( ABS( un  (:,:,:) )                              ) + (/ nimpp - 1, njmpp - 1, 0 /) 
     139            is1(:) = MINLOC( tsn(:,:,:,jp_sal), mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 
     140            is2(:) = MAXLOC( tsn(:,:,:,jp_sal), mask = tmask(:,:,:) == 1._wp ) + (/ nimpp - 1, njmpp - 1, 0 /) 
     141         ENDIF 
     142         IF( numout == 6 )   &   ! force to open ocean.output file 
     143            CALL ctl_opn( numout, 'ocean.output', 'APPEND', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. ) 
     144 
     145         WRITE(numout,cform_err) 
     146         WRITE(numout,*) ' stp_ctl: |ssh| > 10 m  or  |U| > 10 m/s  or  S <= 0  or  S >= 100  or  NaN encounter in the tests' 
     147         WRITE(numout,*) ' ======= ' 
     148         IF( lk_mpp .AND. .NOT. ln_ctl ) WRITE(numout,*) 'E R R O R message from sub-domain: ', narea 
     149         WRITE(numout,9100) kt,   zmax(1), ih(1) , ih(2) 
     150         WRITE(numout,9200) kt,   zmax(2), iu(1) , iu(2) , iu(3) 
     151         WRITE(numout,9300) kt, - zmax(3), is1(1), is1(2), is1(3) 
     152         WRITE(numout,9400) kt,   zmax(4), is2(1), is2(2), is2(3) 
     153         WRITE(numout,*) 
     154         WRITE(numout,*) '          output of last computed fields in output.abort.nc file' 
     155          
     156         CALL dia_wri_state( 'output.abort' )     ! create an output.abort file 
     157          
     158         IF( ln_ctl ) THEN 
     159            kindic = -3 
     160            nstop = nstop + 1                            ! increase nstop by 1 (on all local domains) 
     161         ELSE 
     162            CALL ctl_stop() 
     163            CALL mppstop(ld_force_abort = .true.) 
     164         ENDIF 
     165         ! 
     166      ENDIF 
    177167      ! 
     1689100  FORMAT (' kt=',i8,'   |ssh| max: ',1pg11.4,', at  i j  : ',2i5) 
     1699200  FORMAT (' kt=',i8,'   |U|   max: ',1pg11.4,', at  i j k: ',3i5) 
     1709300  FORMAT (' kt=',i8,'   S     min: ',1pg11.4,', at  i j k: ',3i5) 
     1719400  FORMAT (' kt=',i8,'   S     max: ',1pg11.4,', at  i j k: ',3i5) 
    1781729500  FORMAT(' it :', i8, '    |ssh|_max: ', D23.16, ' |U|_max: ', D23.16,' S_min: ', D23.16,' S_max: ', D23.16) 
    179173      ! 
Note: See TracChangeset for help on using the changeset viewer.