- Timestamp:
- 2018-03-27T15:57:02+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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
Note: See TracChangeset
for help on using the changeset viewer.