Ignore:
Timestamp:
2018-07-05T17:43:44+02:00 (3 years ago)
Author:
clem
Message:

sea-ice: remove negative values mainly created by the advection schemes (UMx and Prather) + adapt BDY to ensure that all the fields are coherent with the rest of the code (max concentration etc)

File:
1 edited

Legend:

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

    r9604 r9880  
    9595      END SELECT 
    9696 
     97      !---------------------------- 
     98      ! Debug the advection schemes 
     99      !---------------------------- 
     100      ! clem: The 2 advection schemes above are not strictly positive. 
     101      !       In Prather, advected fields are bounded by 0 in the routine with a MAX(0,field) ==> likely conservation issues 
     102      !       In UMx    , advected fields are not bounded and negative values can appear. 
     103      !                   These values are usually very small but in some occasions they can also be non-negligible 
     104      !                   Therefore one needs to bound the advected fields by 0 (though this is not a clean fix) 
     105      ! ==> 1) remove negative ice areas and volumes (conservation is ensure) 
     106      CALL ice_var_zapsmall  
     107      ! ==> 2) remove remaining negative advected fields (conservation is not preserved) 
     108      WHERE( v_s (:,:,:)   < 0._wp )   v_s (:,:,:)   = 0._wp 
     109      WHERE( sv_i(:,:,:)   < 0._wp )   sv_i(:,:,:)   = 0._wp 
     110      WHERE( e_i (:,:,:,:) < 0._wp )   e_i (:,:,:,:) = 0._wp 
     111      WHERE( e_s (:,:,:,:) < 0._wp )   e_s (:,:,:,:) = 0._wp 
     112 
    97113      !------------ 
    98114      ! diagnostics 
Note: See TracChangeset for help on using the changeset viewer.