Changeset 11245


Ignore:
Timestamp:
2019-07-10T17:27:04+02:00 (13 months ago)
Author:
girrmann
Message:

dev_r10984_HPC-13 : bugfix following [11241], see #2285

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10984_HPC-13_IRRMANN_BDY_optimization/src/OCE/DYN/dynspg_ts.F90

    r11241 r11245  
    567567               END DO 
    568568            END DO 
     569            !                    !* Update ocean depth (variable volume case only) 
     570            hu_e (2:jpim1,2:jpjm1) = hu_0(2:jpim1,2:jpjm1) + zsshu_a(2:jpim1,2:jpjm1)          
     571            hv_e (2:jpim1,2:jpjm1) = hv_0(2:jpim1,2:jpjm1) + zsshv_a(2:jpim1,2:jpjm1) 
     572            hur_e(2:jpim1,2:jpjm1) = ssumask(2:jpim1,2:jpjm1) / ( hu_e(2:jpim1,2:jpjm1) + 1._wp - ssumask(2:jpim1,2:jpjm1) ) 
     573            hvr_e(2:jpim1,2:jpjm1) = ssvmask(2:jpim1,2:jpjm1) / ( hv_e(2:jpim1,2:jpjm1) + 1._wp - ssvmask(2:jpim1,2:jpjm1) ) 
    569574         ENDIF    
    570575         !          
     
    679684            END DO 
    680685         ENDIF 
    681  
    682686        
    683          IF( .NOT. ln_linssh ) THEN 
    684             !                                             !* domain lateral boundary 
    685             CALL lbc_lnk_multi( 'dynspg_ts', ua_e   , 'U', -1._wp, va_e   , 'V', -1._wp & 
    686                  &                         , zsshu_a, 'U',  1._wp, zsshv_a, 'V',  1._wp ) 
     687         IF( .NOT. ln_linssh ) THEN                      !* domain lateral boundary 
     688            CALL lbc_lnk_multi( 'dynspg_ts', ua_e , 'U', -1._wp, va_e , 'V', -1._wp  & 
     689                 &                         , hu_e , 'U',  1._wp, hv_e , 'V',  1._wp & 
     690                 &                         , hur_e, 'U',  1._wp, hvr_e, 'V',  1._wp ) 
    687691            ! 
    688             hu_e (:,:) = hu_0(:,:) + zsshu_a(:,:)         !* Update ocean depth (variable volume case only) 
    689             hv_e (:,:) = hv_0(:,:) + zsshv_a(:,:) 
    690             hur_e(:,:) = ssumask(:,:) / ( hu_e(:,:) + 1._wp - ssumask(:,:) ) 
    691             hvr_e(:,:) = ssvmask(:,:) / ( hv_e(:,:) + 1._wp - ssvmask(:,:) ) 
    692             ! 
    693          ELSE 
    694             !                                             !* domain lateral boundary 
    695             CALL lbc_lnk_multi( 'dynspg_ts', ua_e   , 'U', -1._wp, va_e   , 'V', -1._wp ) 
     692         ELSE                                            !* domain lateral boundary 
     693            CALL lbc_lnk_multi( 'dynspg_ts', ua_e , 'U', -1._wp, va_e , 'V', -1._wp  ) 
    696694            ! 
    697695         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.