- Timestamp:
- 2019-12-10T12:57:49+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/SAS/diawri.F90
r10425 r12143 125 125 !! define all the NETCDF files and fields 126 126 !! At each time step call histdef to compute the mean if ncessary 127 !! Each n write time step, output the instantaneous or mean fields127 !! Each nn_write time step, output the instantaneous or mean fields 128 128 !!---------------------------------------------------------------------- 129 129 !! … … 138 138 REAL(wp) :: zsto, zout, zmax, zjulian ! local scalars 139 139 !!---------------------------------------------------------------------- 140 !141 IF( ln_timing ) CALL timing_start('dia_wri')142 140 ! 143 141 ! Output the initial state and forcings … … 147 145 ENDIF 148 146 ! 147 IF( nn_write == -1 ) RETURN ! we will never do any output 148 ! 149 IF( ln_timing ) CALL timing_start('dia_wri') 150 ! 149 151 ! 0. Initialisation 150 152 ! ----------------- … … 159 161 ENDIF 160 162 #if defined key_diainstant 161 zsto = n write * rdt163 zsto = nn_write * rdt 162 164 clop = "inst("//TRIM(clop)//")" 163 165 #else … … 165 167 clop = "ave("//TRIM(clop)//")" 166 168 #endif 167 zout = n write * rdt169 zout = nn_write * rdt 168 170 zmax = ( nitend - nit000 + 1 ) * rdt 169 171 … … 196 198 ! WRITE root name in date.file for use by postpro 197 199 IF(lwp) THEN 198 CALL dia_nam( clhstnam, n write,' ' )200 CALL dia_nam( clhstnam, nn_write,' ' ) 199 201 CALL ctl_opn( inum, 'date.file', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 200 202 WRITE(inum,*) clhstnam … … 204 206 ! Define the T grid FILE ( nid_T ) 205 207 206 CALL dia_nam( clhstnam, n write, 'grid_T' )208 CALL dia_nam( clhstnam, nn_write, 'grid_T' ) 207 209 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 208 210 CALL histbeg( clhstnam, jpi, glamt, jpj, gphit, & ! Horizontal grid: glamt and gphit … … 216 218 ! Define the U grid FILE ( nid_U ) 217 219 218 CALL dia_nam( clhstnam, n write, 'grid_U' )220 CALL dia_nam( clhstnam, nn_write, 'grid_U' ) 219 221 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam ! filename 220 222 CALL histbeg( clhstnam, jpi, glamu, jpj, gphiu, & ! Horizontal grid: glamu and gphiu … … 228 230 ! Define the V grid FILE ( nid_V ) 229 231 230 CALL dia_nam( clhstnam, n write, 'grid_V' ) ! filename232 CALL dia_nam( clhstnam, nn_write, 'grid_V' ) ! filename 231 233 IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam 232 234 CALL histbeg( clhstnam, jpi, glamv, jpj, gphiv, & ! Horizontal grid: glamv and gphiv … … 291 293 ! donne le nombre d'elements, et ndex la liste des indices a sortir 292 294 293 IF( lwp .AND. MOD( itmod, n write ) == 0 ) THEN295 IF( lwp .AND. MOD( itmod, nn_write ) == 0 ) THEN 294 296 WRITE(numout,*) 'dia_wri : write model outputs in NetCDF files at ', kt, 'time-step' 295 297 WRITE(numout,*) '~~~~~~ '
Note: See TracChangeset
for help on using the changeset viewer.