Changeset 12859 for NEMO/releases/r4.0/r4.0-HEAD/src/SAS/stpctl.F90
- Timestamp:
- 2020-05-03T11:33:32+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/SAS/stpctl.F90
r10603 r12859 32 32 33 33 INTEGER :: idrun, idtime, idssh, idu, ids, istatus 34 LOGICAL :: lsomeoce35 34 !!---------------------------------------------------------------------- 36 35 !! NEMO/SAS 4.0 , NEMO Consortium (2018) … … 62 61 !!---------------------------------------------------------------------- 63 62 ! 64 ll_wrtstp = ( MOD( kt, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 65 ll_colruns = ll_wrtstp .AND. ( ln_ctl .OR. sn_cfctl%l_runstat ) 66 ll_wrtruns = ll_colruns .AND. lwm 67 IF( kt == nit000 .AND. lwp ) THEN 68 WRITE(numout,*) 69 WRITE(numout,*) 'stp_ctl : time-stepping control' 70 WRITE(numout,*) '~~~~~~~' 63 ll_wrtstp = ( MOD( kt-nit000, sn_cfctl%ptimincr ) == 0 ) .OR. ( kt == nitend ) 64 ll_colruns = ll_wrtstp .AND. sn_cfctl%l_runstat .AND. jpnij > 1 65 ll_wrtruns = ( ll_colruns .OR. jpnij == 1 ) .AND. lwm 66 ! 67 IF( kt == nit000 ) THEN 68 ! 69 IF( lwp ) THEN 70 WRITE(numout,*) 71 WRITE(numout,*) 'stp_ctl : time-stepping control' 72 WRITE(numout,*) '~~~~~~~' 73 ENDIF 71 74 ! ! open time.step file 72 75 IF( lwm ) CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 73 76 ! ! open run.stat file(s) at start whatever 74 77 ! ! the value of sn_cfctl%ptimincr 75 IF( l wm .AND. ( ln_ctl .OR. sn_cfctl%l_runstat )) THEN78 IF( ll_wrtruns ) THEN 76 79 CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 77 80 clname = 'run.stat.nc' … … 85 88 ENDIF 86 89 ENDIF 87 IF( kt == nit000 ) lsomeoce = COUNT( ssmask(:,:) == 1._wp ) > 088 90 ! 89 91 IF(lwm .AND. ll_wrtstp) THEN !== current time step ==! ("time.step" file) … … 92 94 ENDIF 93 95 ! !== test of extrema ==! 94 IF( ll_colruns ) THEN96 IF( ll_colruns .OR. jpnij == 1 ) THEN 95 97 zmax(1) = MAXVAL( vt_i (:,:) ) ! max ice thickness 96 98 zmax(2) = MAXVAL( ABS( u_ice(:,:) ) ) ! max ice velocity (zonal only) 97 99 zmax(3) = MAXVAL( -tm_i (:,:)+273.15_wp , mask = ssmask(:,:) == 1._wp ) ! min ice temperature 98 CALL mpp_max( "stpctl", zmax )! max over the global domain100 IF( ll_colruns ) CALL mpp_max( "stpctl", zmax ) ! max over the global domain 99 101 END IF 100 102 ! !== run statistics ==! ("run.stat" file)
Note: See TracChangeset
for help on using the changeset viewer.