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_interp.F90

    r9160 r9167  
    165165            ptab(i1:i2,j1:j2,jm+3) = sv_i(i1:i2,j1:j2,jl) 
    166166            ptab(i1:i2,j1:j2,jm+4) = oa_i(i1:i2,j1:j2,jl) 
    167             jm = jm + 5 
     167            ptab(i1:i2,j1:j2,jm+5) = a_ip(i1:i2,j1:j2,jl) 
     168            ptab(i1:i2,j1:j2,jm+6) = v_ip(i1:i2,j1:j2,jl) 
     169            ptab(i1:i2,j1:j2,jm+7) = t_su(i1:i2,j1:j2,jl) 
     170            jm = jm + 8 
    168171            DO jk = 1, nlay_s 
    169172               ptab(i1:i2,j1:j2,jm) = e_s(i1:i2,j1:j2,jk,jl)   ;   jm = jm + 1 
     
    192195                     sv_i(ji,jj,jl) = ptab(ji,jj,jm+3) * tmask(ji,jj,1) 
    193196                     oa_i(ji,jj,jl) = ptab(ji,jj,jm+4) * tmask(ji,jj,1) 
     197                     a_ip(ji,jj,jl) = ptab(ji,jj,jm+5) * tmask(ji,jj,1) 
     198                     v_ip(ji,jj,jl) = ptab(ji,jj,jm+6) * tmask(ji,jj,1) 
     199                     t_su(ji,jj,jl) = ptab(ji,jj,jm+7) * tmask(ji,jj,1) 
    194200                  END DO 
    195201               END DO 
    196                jm = jm + 5 
     202               jm = jm + 8 
    197203               ! 
    198204               DO jk = 1, nlay_s 
     
    217223            jm = 1 
    218224            DO jl = 1, jpl 
    219                ztab(:,:,jm  ) = a_i  (:,:,jl) 
    220                ztab(:,:,jm+1) = v_i  (:,:,jl) 
    221                ztab(:,:,jm+2) = v_s  (:,:,jl) 
     225               ztab(:,:,jm  ) = a_i (:,:,jl) 
     226               ztab(:,:,jm+1) = v_i (:,:,jl) 
     227               ztab(:,:,jm+2) = v_s (:,:,jl) 
    222228               ztab(:,:,jm+3) = sv_i(:,:,jl) 
    223229               ztab(:,:,jm+4) = oa_i(:,:,jl) 
    224                jm = jm + 5 
     230               ztab(:,:,jm+5) = a_ip(:,:,jl) 
     231               ztab(:,:,jm+6) = v_ip(:,:,jl) 
     232               ztab(:,:,jm+7) = t_su(:,:,jl) 
     233               jm = jm + 8 
    225234               DO jk = 1, nlay_s 
    226235                  ztab(:,:,jm) = e_s(:,:,jk,jl) 
     
    324333                     v_s (ji,jj,jl) = ztab(ji,jj,jm+2) * tmask(ji,jj,1) 
    325334                     sv_i(ji,jj,jl) = ztab(ji,jj,jm+3) * tmask(ji,jj,1) 
    326                      oa_i (ji,jj,jl) = ztab(ji,jj,jm+4) * tmask(ji,jj,1) 
     335                     oa_i(ji,jj,jl) = ztab(ji,jj,jm+4) * tmask(ji,jj,1) 
     336                     a_ip(ji,jj,jl) = ztab(ji,jj,jm+5) * tmask(ji,jj,1) 
     337                     v_ip(ji,jj,jl) = ztab(ji,jj,jm+6) * tmask(ji,jj,1) 
     338                     t_su(ji,jj,jl) = ztab(ji,jj,jm+7) * tmask(ji,jj,1) 
    327339                  END DO 
    328340               END DO 
    329                jm = jm + 5 
     341               jm = jm + 8 
    330342               ! 
    331343               DO jk = 1, nlay_s 
     
    344356          
    345357         ! integrated values 
    346          vt_i (i1:i2,j1:j2) = SUM( v_i(i1:i2,j1:j2,:), dim=3 ) 
    347          vt_s (i1:i2,j1:j2) = SUM( v_s(i1:i2,j1:j2,:), dim=3 ) 
    348          at_i (i1:i2,j1:j2) = SUM( a_i(i1:i2,j1:j2,:), dim=3 ) 
     358         vt_i (i1:i2,j1:j2) = SUM(      v_i(i1:i2,j1:j2,:),            dim=3 ) 
     359         vt_s (i1:i2,j1:j2) = SUM(      v_s(i1:i2,j1:j2,:),            dim=3 ) 
     360         at_i (i1:i2,j1:j2) = SUM(      a_i(i1:i2,j1:j2,:),            dim=3 ) 
    349361         et_s(i1:i2,j1:j2)  = SUM( SUM( e_s(i1:i2,j1:j2,:,:), dim=4 ), dim=3 ) 
    350362         et_i(i1:i2,j1:j2)  = SUM( SUM( e_i(i1:i2,j1:j2,:,:), dim=4 ), dim=3 ) 
Note: See TracChangeset for help on using the changeset viewer.