Changeset 9866


Ignore:
Timestamp:
2018-07-02T18:19:04+02:00 (2 years ago)
Author:
clem
Message:

solve problem described in ticket #2107

Location:
NEMO/trunk/src/ICE
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icedyn_adv_umx.F90

    r9604 r9866  
    231231            ztra       = ztrd(ji,jj)  - (  zfu_ho(ji,jj) - zfu_ho(ji-1,jj  )   & 
    232232               &                         + zfv_ho(ji,jj) - zfv_ho(ji  ,jj-1) ) * r1_e1e2t(ji,jj)   
    233             ptc(ji,jj) = ptc(ji,jj) + pdt * ztra 
     233            ptc(ji,jj) = ( ptc(ji,jj) + pdt * ztra ) * tmask(ji,jj,1) 
    234234         END DO 
    235235      END DO 
  • NEMO/trunk/src/ICE/icedyn_rhg_evp.F90

    r9660 r9866  
    498498            CALL agrif_interp_ice( 'V' ) 
    499499#endif 
    500 !clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
    501 !            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     500            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
    502501            ! 
    503502            DO jj = 2, jpjm1 
     
    547546            CALL agrif_interp_ice( 'U' ) 
    548547#endif 
    549 !clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
    550 !            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     548            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
    551549            ! 
    552550         ELSE ! odd iterations 
     
    598596            CALL agrif_interp_ice( 'U' ) 
    599597#endif 
    600 !clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
    601 !            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
     598            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 
    602599            ! 
    603600            DO jj = 2, jpjm1 
     
    647644            CALL agrif_interp_ice( 'V' ) 
    648645#endif 
    649 !clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries 
    650 !            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
     646            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 
    651647            ! 
    652648         ENDIF 
  • NEMO/trunk/src/ICE/icestp.F90

    r9725 r9866  
    2121   !! 
    2222   !! SI3 has been made possible by a handful of persons who met as working group 
    23    !!      (in France, Belgium, UK and Italy) 
     23   !!      (from France, Belgium, UK and Italy) 
    2424   !!    * Clement Rousset, Martin Vancoppenolle & Gurvan Madec (LOCEAN, France) 
    2525   !!    * Matthieu Chevalier & David Salas (Meteo France, France) 
Note: See TracChangeset for help on using the changeset viewer.