Ignore:
Timestamp:
2018-12-19T22:54:16+01:00 (2 years ago)
Author:
smasson
Message:

trunk: merge back dev_r10164_HPC09_ESIWACE_PREP_MERGE@10424 into the trunk

File:
1 edited

Legend:

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

    r10068 r10425  
    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 ) THEN 
    91          CALL mpp_max_multiple( zmax(:), 3 )    ! max over the global domain 
    92       ENDIF 
    93       ! 
    94       IF( MOD( kt, nwrite ) == 1 .AND. lwp ) THEN 
    95          WRITE(numout,*) ' ==>> time-step= ', kt, ' vt_i max: ',  zmax(1), ' |u_ice| max: ', zmax(2), ' tm_i min: ', -zmax(3) 
    96       ENDIF 
    97  
     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         CALL mpp_max( "stpctl", zmax )                                   ! max over the global domain 
     94      END IF 
    9895      !                                            !==  run statistics  ==!   ("run.stat" file) 
    99       IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 
    100       IF( lwm ) THEN 
     96      IF( ln_ctl .AND. lwm ) THEN 
     97         IF(lwp) WRITE(numrun,9500) kt, zmax(1), zmax(2), - zmax(3) 
    10198         istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) ) 
    10299         istatus = NF90_PUT_VAR( idrun,   idu, (/ zmax(2)/), (/kt/), (/1/) ) 
     
    106103      END IF 
    107104      ! 
    108 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) 
    109106      ! 
    110107   END SUBROUTINE stp_ctl 
Note: See TracChangeset for help on using the changeset viewer.