Changeset 9397


Ignore:
Timestamp:
2018-03-14T14:27:22+01:00 (2 years ago)
Author:
clem
Message:

dev_CNRS_2017: change ice routines so that reproducibility issue is solved for the ice (but there is still repro issue between r8996 and r9012 after 268 time steps of orca2-lim3-pisces)

Location:
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/iceitd.F90

    r8967 r9397  
    119119            ! Compute thickness change in each ice category 
    120120            DO ji = 1, npti 
    121                zdhice(ji,jl) = h_i_2d(ji,jl) - h_ib_2d(ji,jl) 
     121               IF( a_i_2d(ji,jl) >  epsi10) zdhice(ji,jl) = h_i_2d(ji,jl) - h_ib_2d(ji,jl) 
    122122            END DO 
    123123         END DO 
     
    690690         zalpha = 0.05_wp 
    691691         zhmax  = 3._wp * rn_himean 
    692          hi_max(0) = 0._wp 
     692         hi_max(0) = 0._wp  
    693693         DO jl = 1, jpl 
    694694            znum = jpl * ( zhmax+1 )**zalpha 
  • branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_ent.F90

    r8882 r9397  
    103103      END DO 
    104104 
    105       zeh_cum1(1:npti,0) = 0._wp  
     105      zeh_cum1(1:npti,0:nlay_i) = 0._wp  
    106106      ! new cumulative q*h => linear interpolation 
    107107      DO jk0 = 1, nlay_i+1 
Note: See TracChangeset for help on using the changeset viewer.