Ignore:
Timestamp:
2017-12-22T18:20:43+01:00 (3 years ago)
Author:
clem
Message:

fix issues but agrif + lim3 is still not restartable because of interpolation when ice thermodynamics is activated

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/NST_SRC/agrif_lim3_update.F90

    r9160 r9167  
    9393            ptab(i1:i2,j1:j2,jm+2) = v_s (i1:i2,j1:j2,jl) 
    9494            ptab(i1:i2,j1:j2,jm+3) = sv_i(i1:i2,j1:j2,jl) 
    95             ptab(i1:i2,j1:j2,jm+4) = oa_i (i1:i2,j1:j2,jl) 
    96             jm = jm + 5 
     95            ptab(i1:i2,j1:j2,jm+4) = oa_i(i1:i2,j1:j2,jl) 
     96            ptab(i1:i2,j1:j2,jm+5) = a_ip(i1:i2,j1:j2,jl) 
     97            ptab(i1:i2,j1:j2,jm+6) = v_ip(i1:i2,j1:j2,jl) 
     98            ptab(i1:i2,j1:j2,jm+7) = t_su(i1:i2,j1:j2,jl) 
     99            jm = jm + 8 
    97100            DO jk = 1, nlay_s 
    98101               ptab(i1:i2,j1:j2,jm) = e_s(i1:i2,j1:j2,jk,jl)   ;   jm = jm + 1 
     
    120123                     sv_i(ji,jj,jl) = ptab(ji,jj,jm+3) * tmask(ji,jj,1) 
    121124                     oa_i(ji,jj,jl) = ptab(ji,jj,jm+4) * tmask(ji,jj,1) 
     125                     a_ip(ji,jj,jl) = ptab(ji,jj,jm+5) * tmask(ji,jj,1) 
     126                     v_ip(ji,jj,jl) = ptab(ji,jj,jm+6) * tmask(ji,jj,1) 
     127                     t_su(ji,jj,jl) = ptab(ji,jj,jm+7) * tmask(ji,jj,1) 
    122128                  ENDIF 
    123129               END DO 
    124130            END DO 
    125             jm = jm + 5 
     131            jm = jm + 8 
    126132            ! 
    127133            DO jk = 1, nlay_s 
     
    147153         et_s(i1:i2,j1:j2) = SUM( SUM( e_s(i1:i2,j1:j2,:,:), dim=4 ), dim=3 ) 
    148154         et_i(i1:i2,j1:j2) = SUM( SUM( e_i(i1:i2,j1:j2,:,:), dim=4 ), dim=3 ) 
     155 
     156         at_ip(i1:i2,j1:j2) = SUM( a_ip(i1:i2,j1:j2,:), dim=3 ) ! melt ponds 
     157         vt_ip(i1:i2,j1:j2) = SUM( v_ip(i1:i2,j1:j2,:), dim=3 ) 
     158         ! 
     159         ato_i(i1:i2,j1:j2) = 1._wp - at_i(i1:i2,j1:j2)         ! open water fraction   
    149160          
    150161      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.