- Timestamp:
- 2012-07-11T13:22:58+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r3211 r3432 157 157 ! !------------------------------------------! 158 158 #if defined key_z_first 159 fsdept(:,:,1:jpkm1) = fsdept_n(:,:,1:jpkm1) ! now local depths stored in fsdep. arrays 160 fsdepw(:,:,1:jpkm1) = fsdepw_n(:,:,1:jpkm1) 161 fsde3w(:,:,1:jpkm1) = fsde3w_n(:,:,1:jpkm1) 162 ! 163 fse3t (:,:,1:jpkm1) = fse3t_n (:,:,1:jpkm1) ! vertical scale factors stored in fse3. arrays 164 fse3u (:,:,1:jpkm1) = fse3u_n (:,:,1:jpkm1) 165 fse3v (:,:,1:jpkm1) = fse3v_n (:,:,1:jpkm1) 166 fse3f (:,:,1:jpkm1) = fse3f_n (:,:,1:jpkm1) 167 fse3w (:,:,1:jpkm1) = fse3w_n (:,:,1:jpkm1) 168 fse3uw(:,:,1:jpkm1) = fse3uw_n(:,:,1:jpkm1) 169 fse3vw(:,:,1:jpkm1) = fse3vw_n(:,:,1:jpkm1) 159 ! DCSE_NEMO: can't use implicit loop over k here because the domzgr_substitute.h90 160 ! file causes the line below to be expanded to: 161 ! gdept_1(1:jpkm1,:,:) = (gdept(1:jpkm1,:,:)*(1.+sshn(:,:)*mut(1:jpkm1,:,:))) 162 ! which contains non-conforming array expressions. 163 DO jj=1,jpj,1 164 DO ji=1,jpi,1 165 DO jk=1,jpk,1 166 fsdept(ji,jj,jk) = fsdept_n(ji,jj,jk) ! now local depths stored in fsdep. arrays 167 END DO 168 END DO 169 END DO 170 DO jj=1,jpj,1 171 DO ji=1,jpi,1 172 DO jk=1,jpk,1 173 fsdepw(ji,jj,jk) = fsdepw_n(ji,jj,jk) 174 END DO 175 END DO 176 END DO 177 DO jj=1,jpj,1 178 DO ji=1,jpi,1 179 DO jk=1,jpk,1 180 fsde3w(ji,jj,jk) = fsde3w_n(ji,jj,jk) 181 END DO 182 END DO 183 END DO 184 ! 185 DO jj=1,jpj,1 186 DO ji=1,jpi,1 187 DO jk=1,jpk,1 188 fse3t (ji,jj,jk) = fse3t_n (ji,jj,jk) ! vertical scale factors stored in fse3. arrays 189 END DO 190 END DO 191 END DO 192 DO jj=1,jpj,1 193 DO ji=1,jpi,1 194 DO jk=1,jpk,1 195 fse3u (ji,jj,jk) = fse3u_n (ji,jj,jk) 196 END DO 197 END DO 198 END DO 199 DO jj=1,jpj,1 200 DO ji=1,jpi,1 201 DO jk=1,jpk,1 202 fse3v (ji,jj,jk) = fse3v_n (ji,jj,jk) 203 END DO 204 END DO 205 END DO 206 DO jj=1,jpj,1 207 DO ji=1,jpi,1 208 DO jk=1,jpk,1 209 fse3f (ji,jj,jk) = fse3f_n (ji,jj,jk) 210 END DO 211 END DO 212 END DO 213 DO jj=1,jpj,1 214 DO ji=1,jpi,1 215 DO jk=1,jpk,1 216 fse3w (ji,jj,jk) = fse3w_n (ji,jj,jk) 217 END DO 218 END DO 219 END DO 220 221 222 DO jj=1,jpj,1 223 DO ji=1,jpi,1 224 DO jk=1,jpk,1 225 fse3uw(ji,jj,jk) = fse3uw_n(ji,jj,jk) 226 END DO 227 END DO 228 END DO 229 230 DO jj=1,jpj,1 231 DO ji=1,jpi,1 232 DO jk=1,jpk,1 233 fse3vw(ji,jj,jk) = fse3vw_n(ji,jj,jk) 234 END DO 235 END DO 236 END DO 170 237 #else 171 238 DO jk = 1, jpkm1
Note: See TracChangeset
for help on using the changeset viewer.