- Timestamp:
- 2018-03-27T15:57:02+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/CANAL/MY_SRC/stpctl.F90
r9403 r9441 26 26 USE wet_dry, ONLY : ll_wd, ssh_ref ! reference depth for negative bathy 27 27 28 USE netcdf ! NetCDF library 28 29 IMPLICIT NONE 29 30 PRIVATE 30 31 31 32 PUBLIC stp_ctl ! routine called by step.F90 33 34 INTEGER :: idrun, idtime, idssh, idu, ids, istatus 32 35 !!---------------------------------------------------------------------- 33 36 !! NEMO/OPA 4.0 , NEMO Consortium (2017) … … 75 78 ! ! open run.stat file 76 79 CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 80 81 IF( lwm ) THEN 82 istatus = NF90_CREATE( 'run.stat.nc', NF90_CLOBBER, idrun ) 83 istatus = NF90_DEF_DIM( idrun, 'time', NF90_UNLIMITED, idtime ) 84 istatus = NF90_DEF_VAR( idrun, 'abs_ssh_max', NF90_DOUBLE, (/ idtime /), idssh ) 85 istatus = NF90_DEF_VAR( idrun, 'abs_u_max', NF90_DOUBLE, (/ idtime /), idu ) 86 istatus = NF90_ENDDEF(idrun) 87 ENDIF 88 77 89 ENDIF 78 90 ! … … 89 101 ENDIF 90 102 zmax(2) = MAXVAL( ABS( un(:,:,:) ) ) ! velocity max (zonal only) 91 zmax(3) = MAXVAL( -tsn(:,:,:,jp_sal) , mask = tmask(:,:,:) == 1._wp ) ! minus salinity max 103 ! zmax(3) = MAXVAL( -tsn(:,:,:,jp_sal) , mask = tmask(:,:,:) == 1._wp ) ! minus salinity max 104 zmax(3) = 0.0_wp 92 105 zmax(4) = REAL( nstop , wp ) ! stop indicator 93 106 ! … … 109 122 CALL mpp_maxloc( ABS(sshn) , ssmask(:,:) , zzz, iih, ijh ) 110 123 CALL mpp_maxloc( ABS(un) , umask (:,:,:), zzz, iiu, iju, iku ) 111 124 ! CALL mpp_minloc( tsn(:,:,:,jp_sal), tmask (:,:,:), zzz, iis, ijs, iks ) 112 125 ELSE 113 126 iloch = MINLOC( ABS( sshn(:,:) ) ) 114 127 ilocu = MAXLOC( ABS( un (:,:,:) ) ) 115 128 ! ilocs = MINLOC( tsn(:,:,:,jp_sal) , mask = tmask(:,:,:) == 1._wp ) 116 129 iih = iloch(1) + nimpp - 1 ; ijh = iloch(2) + njmpp - 1 117 130 iiu = ilocu(1) + nimpp - 1 ; iju = ilocu(2) + njmpp - 1 ; iku = ilocu(3) 118 131 ! iis = ilocs(1) + nimpp - 1 ; ijs = ilocs(2) + njmpp - 1 ; iks = ilocu(3) 119 132 ENDIF 120 133 IF(lwp) THEN … … 139 152 ! 140 153 ! !== run statistics ==! ("run.stat" file) 141 IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 154 ! IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 155 IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2) 156 IF( lwm ) THEN 157 istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) ) 158 istatus = NF90_PUT_VAR( idrun, idu, (/ zmax(2)/), (/kt/), (/1/) ) 159 ! istatus = NF90_PUT_VAR( idrun, ids, (/-zmax(3)/), (/kt/), (/1/) ) 160 IF( MOD( kt , 100 ) == 0 ) istatus = NF90_SYNC(idrun) 161 IF( kt == nitend ) istatus = NF90_CLOSE(idrun) 162 END IF 142 163 ! 143 9400 FORMAT(' it :', i8, ' |ssh|_max: ', e16.10, ' |U|_max: ',e16.10,' S_min: ',e16.10) 164 !9400 FORMAT(' it :', i8, ' |ssh|_max: ', D23.16, ' |U|_max: ', D23.16,' S_min: ', D23.16) 165 9400 FORMAT(' it :', i8, ' |ssh|_max: ', D23.16, ' |U|_max: ', D23.16) 144 166 ! 145 167 END SUBROUTINE stp_ctl
Note: See TracChangeset
for help on using the changeset viewer.