Changeset 10512


Ignore:
Timestamp:
2019-01-14T19:53:17+01:00 (19 months ago)
Author:
clem
Message:

change the default options for ice advection

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

Legend:

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

    r10425 r10512  
    146146         END DO 
    147147         ! --- 
    148          CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice + correction on ice thickness 
     148         CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice 
    149149 
    150150         ! diagnostics: divergence at T points  
     
    167167         !CALL RANDOM_NUMBER(v_ice(:,:)) ; v_ice(:,:) = v_ice(:,:) * 0.1 + rn_vice * 0.9 * vmask(:,:,1) 
    168168         ! --- 
    169          CALL ice_dyn_adv   ( kt )   ;   CALL Hbig( zhi_max, zhs_max, zhip_max )   ! -- advection of ice + correction on ice thickness 
     169         CALL ice_dyn_adv   ( kt )                                       ! -- advection of ice 
    170170 
    171171         ! diagnostics: divergence at T points  
  • NEMO/trunk/src/ICE/icedyn_adv_umx.F90

    r10475 r10512  
    3939 
    4040   ! if T interpolated at u/v points is negative, then interpolate T at u/v points using the upstream scheme 
    41    LOGICAL ::   ll_neg = .TRUE. 
     41   !   clem: if set to true, the 2D test case "diagonal advection" does not work (I do not understand why)  
     42   LOGICAL ::   ll_neg = .FALSE. 
    4243    
    4344   ! alternate directions for upstream 
     
    4849    
    4950   ! prelimiter: use it to avoid overshoot in H 
    50    LOGICAL ::   ll_prelimiter_zalesak = .TRUE.  ! from: Zalesak(1979) eq. 14 => better for 1D. Not well defined in 2D 
     51   !   clem: if set to true, the 2D test case "diagnoal advection" does not work (I do not understand why) 
     52   LOGICAL ::   ll_prelimiter_zalesak = .FALSE.  ! from: Zalesak(1979) eq. 14 => better for 1D. Not well defined in 2D 
    5153 
    5254 
Note: See TracChangeset for help on using the changeset viewer.