Ignore:
Timestamp:
2019-03-20T19:58:40+01:00 (20 months ago)
Author:
clem
Message:

fix tickets #2256 and #2257

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/ICE/icethd.F90

    r10534 r10785  
    244244         ! 
    245245      END DO 
    246       ! update ice age (in case a_i changed, i.e. becomes 0 or lateral melting) 
    247       oa_i(:,:,:) = o_i(:,:,:) * a_i(:,:,:) 
    248246 
    249247      IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icethd', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) 
     
    418416         CALL tab_2d_1d( npti, nptidx(1:npti), sst_1d(1:npti), sst_m ) 
    419417         CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d(1:npti), sss_m ) 
    420  
     418         ! 
     419         ! to update ice age 
     420         CALL tab_2d_1d( npti, nptidx(1:npti), o_i_1d (1:npti), o_i (:,:,kl) ) 
     421         CALL tab_2d_1d( npti, nptidx(1:npti), oa_i_1d(1:npti), oa_i(:,:,kl) ) 
     422         ! 
    421423         ! --- Change units of e_i, e_s from J/m2 to J/m3 --- ! 
    422424         DO jk = 1, nlay_i 
     
    443445         sv_i_1d(1:npti) = s_i_1d (1:npti) * v_i_1d (1:npti) 
    444446         v_ip_1d(1:npti) = h_ip_1d(1:npti) * a_ip_1d(1:npti) 
     447         oa_i_1d(1:npti) = o_i_1d (1:npti) * a_i_1d (1:npti) 
    445448          
    446449         CALL tab_1d_2d( npti, nptidx(1:npti), at_i_1d(1:npti), at_i             ) 
     
    516519         CALL tab_1d_2d( npti, nptidx(1:npti), sv_i_1d(1:npti), sv_i(:,:,kl) ) 
    517520         CALL tab_1d_2d( npti, nptidx(1:npti), v_ip_1d(1:npti), v_ip(:,:,kl) ) 
     521         CALL tab_1d_2d( npti, nptidx(1:npti), oa_i_1d(1:npti), oa_i(:,:,kl) ) 
    518522         ! 
    519523      END SELECT 
Note: See TracChangeset for help on using the changeset viewer.