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/SAS/stpctl.F90

    r10314 r10358  
    3232 
    3333   INTEGER  ::   idrun, idtime, idssh, idu, ids, istatus 
     34   LOGICAL  ::   lsomeoce 
    3435   !!---------------------------------------------------------------------- 
    3536   !! NEMO/SAS 4.0 , NEMO Consortium (2018) 
     
    5758      !! 
    5859      REAL(wp), DIMENSION(3) ::   zmax 
     60      CHARACTER(len=20) :: clname 
    5961      !!---------------------------------------------------------------------- 
    6062 
     
    6466         WRITE(numout,*) '~~~~~~~' 
    6567         !                                ! open time.step file 
    66          CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     68         IF( lwm ) CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    6769         !                                ! open run.stat file 
    68          CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    69  
    70          IF( lwm ) THEN 
     70         IF( ln_ctl .AND. lwm ) THEN 
     71            CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     72            clname = 'run.stat.nc' 
     73            IF( .NOT. Agrif_Root() )   clname = TRIM(Agrif_CFixed())//"_"//TRIM(clname) 
    7174            istatus = NF90_CREATE( 'run.stat.nc', NF90_CLOBBER, idrun ) 
    7275            istatus = NF90_DEF_DIM( idrun, 'time'     , NF90_UNLIMITED, idtime ) 
     
    7679            istatus = NF90_ENDDEF(idrun) 
    7780         ENDIF 
    78           
    7981      ENDIF 
     82      IF( kt == nit000 )   lsomeoce = COUNT( ssmask(:,:) == 1._wp ) > 0 
    8083      ! 
    81       IF(lwp) THEN                        !==  current time step  ==!   ("time.step" file) 
     84      IF(lwm) THEN                        !==  current time step  ==!   ("time.step" file) 
    8285         WRITE ( numstp, '(1x, i8)' )   kt 
    8386         REWIND( numstp ) 
    8487      ENDIF 
    8588      !                                   !==  test of extrema  ==! 
    86       zmax(1) = MAXVAL(      vt_i (:,:) )                                           ! max ice thickness 
    87       zmax(2) = MAXVAL( ABS( u_ice(:,:) ) )                                         ! max ice velocity (zonal only) 
    88       zmax(3) = MAXVAL(     -tm_i (:,:)+273.15_wp , mask = ssmask(:,:) == 1._wp )   ! min ice temperature 
    89       ! 
    90       IF( lk_mpp ) CALL mpp_max( "stpctl", zmax )    ! max over the global domain 
    91       ! 
    92       IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 
    93          WRITE(numout,*) ' ==>> time-step= ', kt, ' vt_i max: ',  zmax(1), ' |u_ice| max: ', zmax(2), ' tm_i min: ', -zmax(3) 
    94       ENDIF 
    95  
     89      IF( ln_ctl ) THEN   ! must be done by all processes because of the mpp_max 
     90         zmax(1) = MAXVAL(      vt_i (:,:) )                                           ! max ice thickness 
     91         zmax(2) = MAXVAL( ABS( u_ice(:,:) ) )                                         ! max ice velocity (zonal only) 
     92         zmax(3) = MAXVAL(     -tm_i (:,:)+273.15_wp , mask = ssmask(:,:) == 1._wp )   ! min ice temperature 
     93         IF( lk_mpp ) CALL mpp_max( "stpctl", zmax )                                   ! max over the global domain 
     94      END IF 
    9695      !                                            !==  run statistics  ==!   ("run.stat" file) 
    97       IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 
    98       IF( lwm ) THEN 
     96      IF( ln_ctl .AND. lwm ) THEN 
     97         IF(lwp) WRITE(numrun,9500) kt, zmax(1), zmax(2), - zmax(3) 
    9998         istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) ) 
    10099         istatus = NF90_PUT_VAR( idrun,   idu, (/ zmax(2)/), (/kt/), (/1/) ) 
     
    104103      END IF 
    105104      ! 
    106 9400  FORMAT(' it :', i8, '    vt_i_max: ', D23.16, ' |u|_max: ', D23.16,' tm_i_min: ', D23.16) 
     1059500  FORMAT(' it :', i8, '    vt_i_max: ', D23.16, ' |u|_max: ', D23.16,' tm_i_min: ', D23.16) 
    107106      ! 
    108107   END SUBROUTINE stp_ctl 
Note: See TracChangeset for help on using the changeset viewer.