Ignore:
Timestamp:
2020-07-30T15:40:36+02:00 (2 months ago)
Author:
mathiot
Message:

ticket #1900: add subroutine icb_utl_zavg, fix issue in icb_utl_getkb, add subroutine test_icb_utl_getkb to test awkward cases (will be removed after the review), simplify icbdyn.F90 icbthm.F90 according to the new routine created.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/tickets_icb_1900/src/OCE/ICB/icbdyn.F90

    r13357 r13359  
    332332          
    333333         ! compute mean velocity  
    334          zuo = 0.0 ; zvo = 0.0; zdep = 0.0 
    335          DO jk = 1, ikb-1 
    336             zuo = zuo + zuoce(jk) * ze3t(jk) 
    337             zvo = zvo + zvoce(jk) * ze3t(jk) 
    338             zdep = zdep + ze3t(jk) 
    339          END DO 
    340          zuo = (zuo + zuoce(ikb) * (zD - zdep)) / zD 
    341          zvo = (zvo + zvoce(ikb) * (zD - zdep)) / zD 
     334         CALL icb_utl_zavg(zuo, zuoce, ze3t, zD, ikb) 
     335         CALL icb_utl_zavg(zvo, zvoce, ze3t, zD, ikb) 
    342336      ELSE 
    343337         zuo = zssu 
Note: See TracChangeset for help on using the changeset viewer.