Changeset 9866

Timestamp:

2018-07-02T18:19:04+02:00 (5 years ago)

Author:

clem

Message:

solve problem described in ticket #2107

Location:

NEMO/trunk/src/ICE

Files:

• icedyn_adv_umx.F90 (modified) (1 diff)
• icedyn_rhg_evp.F90 (modified) (4 diffs)
• icestp.F90 (modified) (1 diff)

NEMO/trunk/src/ICE/icedyn_adv_umx.F90

@@ -231 +231 @@
             ztra       = ztrd(ji,jj)  - (  zfu_ho(ji,jj) - zfu_ho(ji-1,jj  )   &
                &                         + zfv_ho(ji,jj) - zfv_ho(ji  ,jj-1) ) * r1_e1e2t(ji,jj)  
-            ptc(ji,jj) = ptc(ji,jj) + pdt * ztra
+            ptc(ji,jj) = ( ptc(ji,jj) + pdt * ztra ) * tmask(ji,jj,1)
          END DO
       END DO

NEMO/trunk/src/ICE/icedyn_rhg_evp.F90

@@ -498 +498 @@
             CALL agrif_interp_ice( 'V' )
 #endif
-!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries
-!            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' )
+            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' )
             !
             DO jj = 2, jpjm1
@@ -547 +546 @@
             CALL agrif_interp_ice( 'U' )
 #endif
-!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries
-!            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' )
+            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' )
             !
          ELSE ! odd iterations
@@ -598 +596 @@
             CALL agrif_interp_ice( 'U' )
 #endif
-!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries
-!            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' )
+            IF( ln_bdy ) CALL bdy_ice_dyn( 'U' )
             !
             DO jj = 2, jpjm1
@@ -647 +644 @@
             CALL agrif_interp_ice( 'V' )
 #endif
-!clem: this call makes rheology numerically unstable, thus I removed it but one has to check the effect on the boundaries
-!            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' )
+            IF( ln_bdy ) CALL bdy_ice_dyn( 'V' )
             !
          ENDIF

NEMO/trunk/src/ICE/icestp.F90

@@ -21 +21 @@
    !!
    !! SI3 has been made possible by a handful of persons who met as working group
-   !!      (in France, Belgium, UK and Italy)
+   !!      (from France, Belgium, UK and Italy)
    !!    * Clement Rousset, Martin Vancoppenolle & Gurvan Madec (LOCEAN, France)
    !!    * Matthieu Chevalier & David Salas (Meteo France, France)