- Timestamp:
- 2021-07-26T16:46:48+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_generic_obs/src/OCE/OBS/obs_oper.F90
r15089 r15144 178 178 ! Initialize daily mean for first timestep of the day 179 179 IF ( idayend == 1 .OR. kt == 0 ) THEN 180 DO_3D( 1, 1, 1, 1, 1, jpk ) 181 prodatqc%vdmean(ji,jj,jk,kvar) = 0.0 182 END_3D 183 ENDIF 184 185 DO_3D( 1, 1, 1, 1, 1, jpk ) 186 ! Increment field 1 for computing daily mean 187 prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 188 & + pvar(ji,jj,jk) 189 END_3D 180 DO jk = 1, jpk 181 DO jj = 1, jpj 182 DO ji = 1, jpi 183 prodatqc%vdmean(ji,jj,jk,kvar) = 0.0 184 END DO 185 END DO 186 END DO 187 ENDIF 188 189 DO jk = 1, jpk 190 DO jj = 1, jpj 191 DO ji = 1, jpi 192 ! Increment field 1 for computing daily mean 193 prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 194 & + pvar(ji,jj,jk) 195 END DO 196 END DO 197 END DO 190 198 191 199 ! Compute the daily mean at the end of day … … 194 202 IF (lwp) WRITE(numout,*) 'Calculating prodatqc%vdmean on time-step: ',kt 195 203 CALL FLUSH(numout) 196 DO_3D( 1, 1, 1, 1, 1, jpk ) 197 prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 198 & * zdaystp 199 END_3D 204 DO jk = 1, jpk 205 DO jj = 1, jpj 206 DO ji = 1, jpi 207 prodatqc%vdmean(ji,jj,jk,kvar) = prodatqc%vdmean(ji,jj,jk,kvar) & 208 & * zdaystp 209 END DO 210 END DO 211 END DO 200 212 ENDIF 201 213
Note: See TracChangeset
for help on using the changeset viewer.