Changeset 6140 for trunk/NEMOGCM/NEMO/OPA_SRC/stpctl.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/stpctl.F90
r5930 r6140 16 16 USE oce ! ocean dynamics and tracers variables 17 17 USE dom_oce ! ocean space and time domain variables 18 USE c1d ! 1D vertical configuration 19 ! 18 20 USE in_out_manager ! I/O manager 19 21 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 20 22 USE lib_mpp ! distributed memory computing 21 USE c1d ! 1D vertical configuration22 23 23 24 IMPLICIT NONE … … 30 31 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 31 32 !!---------------------------------------------------------------------- 32 33 33 CONTAINS 34 34 … … 46 46 !! 47 47 !!---------------------------------------------------------------------- 48 INTEGER, INTENT( in ) :: kt! ocean time-step index49 INTEGER, INTENT( inout ) :: kindic! error indicator48 INTEGER, INTENT(in ) :: kt ! ocean time-step index 49 INTEGER, INTENT(inout) :: kindic ! error indicator 50 50 !! 51 INTEGER :: ji, jj, jk 52 INTEGER :: ii, ij, ik ! temporaryintegers53 REAL(wp) :: zumax, zsmin, zssh2 ! temporaryscalars54 INTEGER, DIMENSION(3) :: ilocu 55 INTEGER, DIMENSION(2) :: ilocs 51 INTEGER :: ji, jj, jk ! dummy loop indices 52 INTEGER :: ii, ij, ik ! local integers 53 REAL(wp) :: zumax, zsmin, zssh2 ! local scalars 54 INTEGER, DIMENSION(3) :: ilocu ! 55 INTEGER, DIMENSION(2) :: ilocs ! 56 56 !!---------------------------------------------------------------------- 57 57 ! 58 58 IF( kt == nit000 .AND. lwp ) THEN 59 59 WRITE(numout,*) … … 63 63 CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 64 64 ENDIF 65 65 ! 66 66 IF(lwp) WRITE ( numstp, '(1x, i8)' ) kt !* save the current time step in numstp 67 67 IF(lwp) REWIND( numstp ) ! -------------------------- 68 68 ! 69 69 ! !* Test maximum of velocity (zonal only) 70 70 ! ! ------------------------ … … 102 102 ENDIF 103 103 9400 FORMAT (' kt=',i6,' max abs(U): ',1pg11.4,', i j k: ',3i5) 104 104 ! 105 105 ! !* Test minimum of salinity 106 106 ! ! ------------------------ 107 107 !! zsmin = MINVAL( tsn(:,:,1,jp_sal), mask = tmask(:,:,1) == 1.e0 ) slower than the following loop on NEC SX5 108 zsmin = 100. e0108 zsmin = 100._wp 109 109 DO jj = 2, jpjm1 110 110 DO ji = 1, jpi … … 136 136 ENDIF 137 137 9500 FORMAT (' kt=',i6,' min SSS: ',1pg11.4,', i j: ',2i5) 138 139 138 ! 139 ! 140 140 IF( lk_c1d ) RETURN ! No log file in case of 1D vertical configuration 141 141 … … 151 151 IF(lwp) WRITE(numsol,9300) kt, zssh2, zumax, zsmin ! ssh statistics 152 152 ! 153 154 153 9200 FORMAT('it:', i8, ' iter:', i4, ' r: ',e16.10, ' b: ',e16.10 ) 155 154 9300 FORMAT(' it :', i8, ' ssh2: ', e16.10, ' Umax: ',e16.10,' Smin: ',e16.10)
Note: See TracChangeset
for help on using the changeset viewer.