Changeset 9415


Ignore:
Timestamp:
2018-03-21T16:51:02+01:00 (2 years ago)
Author:
clem
Message:

debug ice diffusion: avoid too small ice thicknesses by including a threshold (rn_himin) after ice dynamics. This should prevent the ice-ocean conduction flux from being abnormally strong

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icecor.F90

    r9124 r9415  
    6666      ENDIF 
    6767      ! 
    68       IF( kn == 2 ) THEN 
    69          !                          !----------------------------------------------------- 
    70          !                          !  thickness of the smallest category above himin    ! 
    71          !                          !----------------------------------------------------- 
    72          WHERE( a_i(:,:,1) >= epsi20 )   ;   h_i(:,:,1) = v_i (:,:,1) / a_i(:,:,1) 
    73          ELSEWHERE                       ;   h_i(:,:,1) = 0._wp 
    74          END WHERE 
    75          WHERE( h_i(:,:,1) < rn_himin )      a_i(:,:,1) = a_i (:,:,1) * h_i(:,:,1) / rn_himin 
    76          ! 
    77       ENDIF 
     68      !                             !----------------------------------------------------- 
     69      !                             !  ice thickness must exceed himin (for ice diff)    ! 
     70      !                             !----------------------------------------------------- 
     71      WHERE( a_i(:,:,:) >= epsi20 )   ;   h_i(:,:,:) = v_i(:,:,:) / a_i(:,:,:) 
     72      ELSEWHERE                       ;   h_i(:,:,:) = 0._wp 
     73      END WHERE 
     74      WHERE( h_i(:,:,:) < rn_himin )      a_i(:,:,:) = a_i(:,:,:) * h_i(:,:,:) / rn_himin 
     75      ! 
    7876      !                             !----------------------------------------------------- 
    7977      !                             !  ice concentration should not exceed amax          ! 
Note: See TracChangeset for help on using the changeset viewer.