Changeset 10787 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynnxt.F90
- Timestamp:
- 2019-03-21T09:14:33+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynnxt.F90
r10768 r10787 358 358 ! 359 359 IF(.NOT.ln_linssh ) THEN 360 hu _b(:,:) = e3u(:,:,1,ikt ) * umask(:,:,1)361 hv _b(:,:) = e3v(:,:,1,ikt ) * vmask(:,:,1)360 hu(:,:,ikt) = e3u(:,:,1,ikt ) * umask(:,:,1) 361 hv(:,:,ikt) = e3v(:,:,1,ikt ) * vmask(:,:,1) 362 362 DO jk = 2, jpkm1 363 hu _b(:,:) = hu_b(:,:) + e3u(:,:,jk,ikt ) * umask(:,:,jk)364 hv _b(:,:) = hv_b(:,:) + e3v(:,:,jk,ikt ) * vmask(:,:,jk)363 hu(:,:,ikt) = hu(:,:,ikt) + e3u(:,:,jk,ikt ) * umask(:,:,jk) 364 hv(:,:,ikt) = hv(:,:,ikt) + e3v(:,:,jk,ikt ) * vmask(:,:,jk) 365 365 END DO 366 r1_hu _b(:,:) = ssumask(:,:) / ( hu_b(:,:) + 1._wp - ssumask(:,:) )367 r1_hv _b(:,:) = ssvmask(:,:) / ( hv_b(:,:) + 1._wp - ssvmask(:,:) )366 r1_hu(:,:,ikt) = ssumask(:,:) / ( hu(:,:,ikt) + 1._wp - ssumask(:,:) ) 367 r1_hv(:,:,ikt) = ssvmask(:,:) / ( hv(:,:,ikt) + 1._wp - ssvmask(:,:) ) 368 368 ENDIF 369 369 ! … … 378 378 vb_b(:,:) = vb_b(:,:) + e3v(:,:,jk,ikt ) * vv(:,:,jk,ikt ) * vmask(:,:,jk) 379 379 END DO 380 un_b(:,:) = un_b(:,:) * r1_hu _a(:,:)381 vn_b(:,:) = vn_b(:,:) * r1_hv _a(:,:)382 ub_b(:,:) = ub_b(:,:) * r1_hu _b(:,:)383 vb_b(:,:) = vb_b(:,:) * r1_hv _b(:,:)380 un_b(:,:) = un_b(:,:) * r1_hu(:,:,kNp1) 381 vn_b(:,:) = vn_b(:,:) * r1_hv(:,:,kNp1) 382 ub_b(:,:) = ub_b(:,:) * r1_hu(:,:,ikt) 383 vb_b(:,:) = vb_b(:,:) * r1_hv(:,:,ikt) 384 384 ! 385 385 IF( .NOT.ln_dynspg_ts ) THEN ! output the barotropic currents
Note: See TracChangeset
for help on using the changeset viewer.