- Timestamp:
- 2019-07-09T12:05:00+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GO8_package/src/OCE/DIA/diaprod.F90
r11110 r11227 163 163 ENDIF 164 164 165 IF( iom_use("uv") ) THEN 166 z3d(:,:,:) = 0.e0 167 DO jk = 1, jpkm1 168 DO jj = 2, jpjm1 169 DO ji = fs_2, fs_jpim1 ! vector opt. 170 z3d(ji,jj,jk) = 0.25 * ( un(ji-1,jj,jk) + un(ji,jj,jk) ) * ( vn(ji,jj-1,jk) + vn(ji,jj,jk) ) 171 END DO 172 END DO 173 END DO 174 CALL iom_put( "uv", z3d ) ! product of zonal velocity and meridional velocity at T points 175 ENDIF 176 177 IF( iom_use("uw") ) THEN 178 z3d(:,:,:) = 0.e0 179 DO jj = 2, jpjm1 180 DO ji = fs_2, fs_jpim1 ! vector opt. 181 z3d(ji,jj,1) = 0.5 * ( wn(ji,jj,1) + wn(ji+1,jj,1) ) * un(ji,jj,1) 182 END DO 183 END DO 184 DO jk = 2, jpkm1 185 DO jj = 2, jpjm1 186 DO ji = fs_2, fs_jpim1 ! vector opt. 187 z3d(ji,jj,jk) = 0.25 * ( wn(ji,jj,jk) + wn(ji+1,jj,jk) ) * ( un(ji,jj,jk-1) + un(ji,jj,jk) ) 188 END DO 189 END DO 190 END DO 191 CALL iom_put( "uw", z3d ) ! product of zonal velocity and vertical velocity at UW points 192 ENDIF 193 194 IF( iom_use("vw") ) THEN 195 z3d(:,:,:) = 0.e0 196 DO jj = 2, jpjm1 197 DO ji = fs_2, fs_jpim1 ! vector opt. 198 z3d(ji,jj,1) = 0.5 * ( wn(ji,jj,1) + wn(ji,jj+1,1) ) * vn(ji,jj,1) 199 END DO 200 END DO 201 DO jk = 2, jpkm1 202 DO jj = 2, jpjm1 203 DO ji = fs_2, fs_jpim1 ! vector opt. 204 z3d(ji,jj,jk) = 0.25 * ( wn(ji,jj,jk) + wn(ji,jj+1,jk) ) * ( vn(ji,jj,jk-1) + vn(ji,jj,jk) ) 205 END DO 206 END DO 207 END DO 208 CALL iom_put( "vw", z3d ) ! product of meriodional velocity and vertical velocity at VW points 209 ENDIF 210 165 211 IF( iom_use("urhop") ) THEN 166 212 z3d(:,:,:) = 0.e0
Note: See TracChangeset
for help on using the changeset viewer.