Changeset 10512

Timestamp:

2019-01-14T19:53:17+01:00 (4 years ago)

Author:

clem

Message:

change the default options for ice advection

Location:

NEMO/trunk/src/ICE

Files:

• icedyn.F90 (modified) (2 diffs)
• icedyn_adv_umx.F90 (modified) (2 diffs)

NEMO/trunk/src/ICE/icedyn.F90

@@ -146 +146 @@
          END DO
          ! ---
-         CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice + correction on ice thickness
+         CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice
 
          ! diagnostics: divergence at T points 
@@ -167 +167 @@
          !CALL RANDOM_NUMBER(v_ice(:,:)) ; v_ice(:,:) = v_ice(:,:) * 0.1 + rn_vice * 0.9 * vmask(:,:,1)
          ! ---
-         CALL ice_dyn_adv   ( kt )   ;   CALL Hbig( zhi_max, zhs_max, zhip_max )   ! -- advection of ice + correction on ice thickness
+         CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice
 
          ! diagnostics: divergence at T points 

NEMO/trunk/src/ICE/icedyn_adv_umx.F90

@@ -39 +39 @@
 
    ! if T interpolated at u/v points is negative, then interpolate T at u/v points using the upstream scheme
-   LOGICAL ::   ll_neg = .TRUE.
+   !   clem: if set to true, the 2D test case "diagonal advection" does not work (I do not understand why) 
+   LOGICAL ::   ll_neg = .FALSE.
    
    ! alternate directions for upstream
@@ -48 +49 @@
    
    ! prelimiter: use it to avoid overshoot in H
-   LOGICAL ::   ll_prelimiter_zalesak = .TRUE.  ! from: Zalesak(1979) eq. 14 => better for 1D. Not well defined in 2D
+   !   clem: if set to true, the 2D test case "diagnoal advection" does not work (I do not understand why)
+   LOGICAL ::   ll_prelimiter_zalesak = .FALSE.  ! from: Zalesak(1979) eq. 14 => better for 1D. Not well defined in 2D