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