Changeset 9441
- Timestamp:
- 2018-03-27T15:57:02+02:00 (6 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM
- Files:
-
- 3 edited
-
CONFIG/TEST_CASES/CANAL/MY_SRC/stpctl.F90 (modified) (5 diffs)
-
NEMO/OPA_SRC/stpctl.F90 (modified) (3 diffs)
-
NEMO/SAS_SRC/stpctl.F90 (modified) (3 diffs)
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 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/stpctl.F90
r9210 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_DEF_VAR( idrun, 's_min', NF90_DOUBLE, (/ idtime /), ids ) 87 istatus = NF90_ENDDEF(idrun) 88 ENDIF 89 77 90 ENDIF 78 91 ! … … 141 154 ! !== run statistics ==! ("run.stat" file) 142 155 IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 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 143 163 ! 144 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) 145 165 ! 146 166 END SUBROUTINE stp_ctl -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/SAS_SRC/stpctl.F90
r9267 r9441 25 25 USE lib_mpp ! distributed memory computing 26 26 27 USE netcdf ! NetCDF library 27 28 IMPLICIT NONE 28 29 PRIVATE 29 30 30 31 PUBLIC stp_ctl ! routine called by step.F90 32 33 INTEGER :: idrun, idtime, idssh, idu, ids, istatus 31 34 !!---------------------------------------------------------------------- 32 35 !! NEMO/OPA 4.0 , NEMO Consortium (2017) … … 64 67 ! ! open run.stat file 65 68 CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 69 70 IF( lwm ) THEN 71 istatus = NF90_CREATE( 'run.stat.nc', NF90_CLOBBER, idrun ) 72 istatus = NF90_DEF_DIM( idrun, 'time', NF90_UNLIMITED, idtime ) 73 istatus = NF90_DEF_VAR( idrun, 'abs_ssh_max', NF90_DOUBLE, (/ idtime /), idssh ) 74 istatus = NF90_DEF_VAR( idrun, 'abs_u_max', NF90_DOUBLE, (/ idtime /), idu ) 75 istatus = NF90_DEF_VAR( idrun, 's_min', NF90_DOUBLE, (/ idtime /), ids ) 76 istatus = NF90_ENDDEF(idrun) 77 ENDIF 78 66 79 ENDIF 67 80 ! … … 84 97 85 98 ! !== run statistics ==! ("run.stat" file) 86 IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), -zmax(3) 99 IF(lwp) WRITE(numrun,9400) kt, zmax(1), zmax(2), - zmax(3) 100 IF( lwm ) THEN 101 istatus = NF90_PUT_VAR( idrun, idssh, (/ zmax(1)/), (/kt/), (/1/) ) 102 istatus = NF90_PUT_VAR( idrun, idu, (/ zmax(2)/), (/kt/), (/1/) ) 103 istatus = NF90_PUT_VAR( idrun, ids, (/-zmax(3)/), (/kt/), (/1/) ) 104 IF( MOD( kt , 100 ) == 0 ) istatus = NF90_SYNC(idrun) 105 IF( kt == nitend ) istatus = NF90_CLOSE(idrun) 106 END IF 87 107 ! 88 9400 FORMAT(' it :', i8, ' vt_i_max: ', e16.10, ' |u_ice|_max: ',e16.10,' tm_i_min: ',e16.10)108 9400 FORMAT(' it :', i8, ' |ssh|_max: ', D23.16, ' |U|_max: ', D23.16,' S_min: ', D23.16) 89 109 ! 90 110 END SUBROUTINE stp_ctl
Note: See TracChangeset
for help on using the changeset viewer.
