Changeset 7753 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/depth_e3.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/depth_e3.F90
r7698 r7753 150 150 REAL(wp), DIMENSION(:,:,:), INTENT( out) :: pdept_3d, pdepw_3d ! depth = SUM( e3 ) [m] 151 151 ! 152 INTEGER :: jk , jj, ji! dummy loop indices152 INTEGER :: jk ! dummy loop indices 153 153 !!---------------------------------------------------------------------- 154 154 ! 155 !$OMP PARALLEL 156 !$OMP DO schedule(static) private(jj,ji) 157 DO jj = 1, jpj 158 DO ji = 1, jpi 159 pdepw_3d(ji,jj,1) = 0.0_wp 160 pdept_3d(ji,jj,1) = 0.5_wp * pe3w_3d(ji,jj,1) 161 END DO 155 pdepw_3d(:,:,1) = 0.0_wp 156 pdept_3d(:,:,1) = 0.5_wp * pe3w_3d(:,:,1) 157 DO jk = 2, jpk 158 pdepw_3d(:,:,jk) = pdepw_3d(:,:,jk-1) + pe3t_3d(:,:,jk-1) 159 pdept_3d(:,:,jk) = pdept_3d(:,:,jk-1) + pe3w_3d(:,:,jk ) 162 160 END DO 163 DO jk = 2, jpk164 !$OMP DO schedule(static) private(jj,ji)165 DO jj = 1, jpj166 DO ji = 1, jpi167 pdepw_3d(ji,jj,jk) = pdepw_3d(ji,jj,jk-1) + pe3t_3d(ji,jj,jk-1)168 pdept_3d(ji,jj,jk) = pdept_3d(ji,jj,jk-1) + pe3w_3d(ji,jj,jk )169 END DO170 END DO171 END DO172 !$OMP END PARALLEL173 161 ! 174 162 END SUBROUTINE e3_to_depth_3d
Note: See TracChangeset
for help on using the changeset viewer.