Changeset 12294 for NEMO/trunk
- Timestamp:
- 2019-12-30T19:46:59+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/SBC/sbcisf.F90
r11536 r12294 738 738 CHARACTER(len=1), INTENT(in ) :: cd_ptin ! point of variable in/out 739 739 ! 740 INTEGER :: ji, jj, jk ! loop index741 INTEGER :: ikt, ikb ! top and bottom index of the tbl740 INTEGER :: ji, jj, jk ! loop index 741 INTEGER :: ikt, ikb ! top and bottom index of the tbl 742 742 REAL(wp) :: ze3, zhk 743 743 REAL(wp), DIMENSION(jpi,jpj) :: zhisf_tbl ! thickness of the tbl 744 REAL(wp), DIMENSION(jpi,jpj) :: zvarout 744 745 !!---------------------------------------------------------------------- 745 746 … … 749 750 SELECT CASE ( cd_ptin ) 750 751 CASE ( 'U' ) ! compute U in the top boundary layer at T- point 752 ! 753 zvarout(:,:)=0._wp 754 ! 751 755 DO jj = 1,jpj 752 756 DO ji = 1,jpi … … 764 768 DO jk = ikt, ikb - 1 765 769 ze3 = e3u_n(ji,jj,jk) 766 pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3770 zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 767 771 END DO 768 772 769 773 ! level partially include in ice shelf boundary layer 770 774 zhk = SUM( e3u_n(ji, jj, ikt:ikb - 1)) / zhisf_tbl(ji,jj) 771 pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,ikb) * (1._wp - zhk)775 zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,ikb) * (1._wp - zhk) 772 776 END DO 773 777 END DO … … 775 779 DO ji = 2, jpi 776 780 !!gm a wet-point only average should be used here !!! 777 pvarout(ji,jj) = 0.5_wp * ( pvarout(ji,jj) + pvarout(ji-1,jj))781 pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji-1,jj)) 778 782 END DO 779 783 END DO … … 781 785 782 786 CASE ( 'V' ) ! compute V in the top boundary layer at T- point 787 ! 788 zvarout(:,:)=0._wp 789 ! 783 790 DO jj = 1,jpj 784 791 DO ji = 1,jpi … … 796 803 DO jk = ikt, ikb - 1 797 804 ze3 = e3v_n(ji,jj,jk) 798 pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3805 zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 799 806 END DO 800 807 801 808 ! level partially include in ice shelf boundary layer 802 809 zhk = SUM( e3v_n(ji, jj, ikt:ikb - 1)) / zhisf_tbl(ji,jj) 803 pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,ikb) * (1._wp - zhk)810 zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,ikb) * (1._wp - zhk) 804 811 END DO 805 812 END DO … … 807 814 DO ji = 2, jpi 808 815 !!gm a wet-point only average should be used here !!! 809 pvarout(ji,jj) = 0.5_wp * ( pvarout(ji,jj) + pvarout(ji,jj-1))816 pvarout(ji,jj) = 0.5_wp * (zvarout(ji,jj) + zvarout(ji,jj-1)) 810 817 END DO 811 818 END DO
Note: See TracChangeset
for help on using the changeset viewer.