New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14866 for NEMO/branches/UKMO – NEMO

Changeset 14866 for NEMO/branches/UKMO


Ignore:
Timestamp:
2021-05-14T17:21:36+02:00 (3 years ago)
Author:
edblockley
Message:

Fixing reproducibility issue in melt-ponds scheme

Location:
NEMO/branches/UKMO/NEMO_4.0.4_GO8_package
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/cfgs/SHARED/namelist_ice_ref

    r14302 r14866  
    199199         rn_apnd_min =   0.15         !     minimum meltwater fraction contributing to pond growth (TOPO and LEV) 
    200200         rn_apnd_max =   0.85         !     maximum meltwater fraction contributing to pond growth (TOPO and LEV) 
    201          rn_pnd_flush=   0.01         !     pond flushing efficiency (tuning parameter) (LEV) 
     201         rn_pnd_flush=   0.1          !     pond flushing efficiency (tuning parameter) (LEV) 
    202202      ln_pnd_CST     = .false.        !  constant  melt ponds 
    203203         rn_apnd     =   0.2          !     prescribed pond fraction, at Tsu=0 degC 
  • NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/ICE/icestp.F90

    r14335 r14866  
    269269      ENDIF 
    270270      CALL ice_var_glo2eqv 
    271       CALL ice_var_agg(1) 
     271      CALL ice_var_agg(2) 
    272272      ! 
    273273      CALL ice_dyn_init                ! set ice dynamics parameters 
  • NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/ICE/icethd_pnd.F90

    r14727 r14866  
    582582       
    583583      zvolp(:,:) = 0._wp 
     584      zvolp_res(:,:) = 0._wp 
    584585 
    585586      DO jl = 1, jpl 
     
    647648                  zTavg = zTavg + t_su(ji,jj,jl)*a_i(ji,jj,jl) 
    648649               END DO 
    649 !! EWB: bug here - should be divided by at_i(ji,jj) but zTavg is not used anywhere 
    650                zTavg = zTavg / a_i(ji,jj,jl) !!! could get a division by zero here 
     650               zTavg = zTavg / at_i(ji,jj) !!! could get a division by zero here 
    651651          
    652652               DO jl = 1, jpl-1 
Note: See TracChangeset for help on using the changeset viewer.