Changeset 10787 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DOM/domvvl.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/DOM/domvvl.F90
r10767 r10787 616 616 tilde_e3t_n(:,:,:) = tilde_e3t_a(:,:,:) 617 617 ENDIF 618 gdept_b(:,:,:) = gdept_n(:,:,:)619 gdepw_b(:,:,:) = gdepw_n(:,:,:)620 618 621 619 !! TO BE DELETED 620 !!$ gdept_b(:,:,:) = gdept_n(:,:,:) 621 !!$ gdepw_b(:,:,:) = gdepw_n(:,:,:) 622 622 623 !!$ e3t_n(:,:,:) = e3t_a(:,:,:) 623 624 !!$ e3u_n(:,:,:) = e3u_a(:,:,:) … … 643 644 644 645 ! t- and w- points depth (set the isf depth as it is in the initial step) 645 gdept _n(:,:,1) = 0.5_wp * e3w(:,:,1,kNnn_2lev)646 gdepw _n(:,:,1) = 0.0_wp647 gde3w _n(:,:,1) = gdept_n(:,:,1) - sshn(:,:)646 gdept(:,:,1,kNnn_2lev) = 0.5_wp * e3w(:,:,1,kNnn_2lev) 647 gdepw(:,:,1,kNnn_2lev) = 0.0_wp 648 gde3w(:,:,1) = gdept(:,:,1,kNnn_2lev) - sshn(:,:) 648 649 DO jk = 2, jpk 649 650 DO jj = 1,jpj … … 652 653 ! 1 for jk = mikt 653 654 zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk)) 654 gdepw _n(ji,jj,jk) = gdepw_n(ji,jj,jk-1) + e3t(ji,jj,jk-1,kNnn)655 gdept _n(ji,jj,jk) = zcoef * ( gdepw_n(ji,jj,jk) + 0.5 * e3w(ji,jj,jk,kNnn_2lev) ) &656 & + (1-zcoef) * ( gdept _n(ji,jj,jk-1) + e3w(ji,jj,jk,kNnn_2lev) )657 gde3w _n(ji,jj,jk) = gdept_n(ji,jj,jk) - sshn(ji,jj)655 gdepw(ji,jj,jk,kNnn_2lev) = gdepw(ji,jj,jk-1,kNnn_2lev) + e3t(ji,jj,jk-1,kNnn) 656 gdept(ji,jj,jk,kNnn_2lev) = zcoef * ( gdepw(ji,jj,jk,kNnn_2lev ) + 0.5 * e3w(ji,jj,jk,kNnn_2lev) ) & 657 & + (1-zcoef) * ( gdept(ji,jj,jk-1,kNnn_2lev) + e3w(ji,jj,jk,kNnn_2lev) ) 658 gde3w(ji,jj,jk) = gdept(ji,jj,jk,kNnn_2lev) - sshn(ji,jj) 658 659 END DO 659 660 END DO … … 662 663 ! Local depth and Inverse of the local depth of the water 663 664 ! ------------------------------------------------------- 664 hu_n(:,:) = hu_a(:,:) ; r1_hu_n(:,:) = r1_hu_a(:,:) 665 hv_n(:,:) = hv_a(:,:) ; r1_hv_n(:,:) = r1_hv_a(:,:) 666 ! 667 ht_n(:,:) = e3t(:,:,1,kNnn) * tmask(:,:,1) 665 !! TO BE DELETED 666 !!$ hu_n(:,:) = hu_a(:,:) ; r1_hu_n(:,:) = r1_hu_a(:,:) 667 !!$ hv_n(:,:) = hv_a(:,:) ; r1_hv_n(:,:) = r1_hv_a(:,:) 668 !! TO BE DELETED 669 ! 670 ht(:,:) = e3t(:,:,1,kNnn) * tmask(:,:,1) 668 671 DO jk = 2, jpkm1 669 ht _n(:,:) = ht_n(:,:) + e3t(:,:,jk,kNnn) * tmask(:,:,jk)672 ht(:,:) = ht(:,:) + e3t(:,:,jk,kNnn) * tmask(:,:,jk) 670 673 END DO 671 674
Note: See TracChangeset
for help on using the changeset viewer.