Changeset 12294


Ignore:
Timestamp:
2019-12-30T19:46:59+01:00 (4 weeks ago)
Author:
mathiot
Message:

fix #2359

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/sbcisf.F90

    r11536 r12294  
    738738      CHARACTER(len=1),           INTENT(in   ) :: cd_ptin ! point of variable in/out 
    739739      ! 
    740       INTEGER ::   ji, jj, jk                ! loop index 
    741       INTEGER ::   ikt, ikb                    ! top and bottom index of the tbl 
     740      INTEGER ::   ji, jj, jk                   ! loop index 
     741      INTEGER ::   ikt, ikb                     ! top and bottom index of the tbl 
    742742      REAL(wp) ::   ze3, zhk 
    743743      REAL(wp), DIMENSION(jpi,jpj) :: zhisf_tbl ! thickness of the tbl 
     744      REAL(wp), DIMENSION(jpi,jpj) :: zvarout 
    744745      !!---------------------------------------------------------------------- 
    745746       
     
    749750      SELECT CASE ( cd_ptin ) 
    750751      CASE ( 'U' ) ! compute U in the top boundary layer at T- point  
     752         ! 
     753         zvarout(:,:)=0._wp 
     754         ! 
    751755         DO jj = 1,jpj 
    752756            DO ji = 1,jpi 
     
    764768               DO jk = ikt, ikb - 1 
    765769                  ze3 = e3u_n(ji,jj,jk) 
    766                   pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 
     770                  zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 
    767771               END DO 
    768772 
    769773               ! level partially include in ice shelf boundary layer  
    770774               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) 
    772776            END DO 
    773777         END DO 
     
    775779            DO ji = 2, jpi 
    776780!!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)) 
    778782            END DO 
    779783         END DO 
     
    781785       
    782786      CASE ( 'V' ) ! compute V in the top boundary layer at T- point  
     787         ! 
     788         zvarout(:,:)=0._wp 
     789         ! 
    783790         DO jj = 1,jpj 
    784791            DO ji = 1,jpi 
     
    796803               DO jk = ikt, ikb - 1 
    797804                  ze3 = e3v_n(ji,jj,jk) 
    798                   pvarout(ji,jj) = pvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 
     805                  zvarout(ji,jj) = zvarout(ji,jj) + pvarin(ji,jj,jk) / zhisf_tbl(ji,jj) * ze3 
    799806               END DO 
    800807 
    801808               ! level partially include in ice shelf boundary layer  
    802809               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) 
    804811            END DO 
    805812         END DO 
     
    807814            DO ji = 2, jpi 
    808815!!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)) 
    810817            END DO 
    811818         END DO 
Note: See TracChangeset for help on using the changeset viewer.