#include<bits/stdc++.h>

using namespace std;

#define il inline

#define gc getchar()

#define t(i) edge[i].to

#define n(i) edge[i].nxt

#define w(i) edge[i].wei

#define ri register int

#define rb register bool

#define rc register char

#define rp(i,x,y) for(ri i=x;i<=y;++i)

#define e(i,x) for(ri i=head[x];~i;i=n(i))

const int N=+,inf=1e9;

int S,T,ed_cnt=-,head[N],as,cur[N],dep[N],p,q,r,d;

int mov_x[]={,,,-},mov_y[]={,-,,};

struct ed{int to,nxt,wei;}edge[N<<];

il int read()

{

    rc ch=gc;ri x=;rb y=;

    while(ch!='-' && (ch>'' || ch<''))ch=gc;

    if(ch=='-')ch=gc,y=;

    while(ch>='' && ch<='')x=(x<<)+(x<<)+(ch^''),ch=gc;

    return y?x:-x;

}

il int nam(ri x,ri y,ri z){return ((x-)*q+y-)*r+z;}

il void print(ri dat)

{

    if(!dat)return void(printf("0 , 0 , 0"));

    if(dat==p*q*r+)return void(printf("0 , 0 , 0"));

    ri z=dat%r;if(!z)z+=r;dat-=z;dat/=r;dat+=;

    ri y=dat%q;if(!y)y+=q;dat-=y;dat/=q;dat+=;

    ri x=dat;

    printf("%d , %d , %d",x,y,z);

}

il void ad(ri x,ri y,ri z)

{

    //printf("(");print(y);printf(") -> (");print(x);printf(") : %d\n",z);

    //printf("%d -> %d : %d\n",y,x,z);

    edge[++ed_cnt]=(ed){x,head[y],z};head[y]=ed_cnt;edge[++ed_cnt]=(ed){y,head[x],};head[x]=ed_cnt;}

il bool bfs()

{

    queue<int>Q;Q.push(S);memset(dep,,sizeof(dep));dep[S]=;

    while(!Q.empty()){ri nw=Q.front();Q.pop();e(i,nw)if(w(i) && !dep[t(i)])dep[t(i)]=dep[nw]+,Q.push(t(i));}

    return dep[T];

}

il int dfs(ri nw,ri flow)

{

    if(nw==T || !flow)return flow;ri ret=;

    for(ri &i=cur[nw];~i;i=n(i))

        if(w(i) && dep[t(i)]==dep[nw]+)

        {ri tmp=dfs(t(i),min(flow,w(i)));ret+=tmp,w(i)-=tmp;w(i^)+=tmp,flow-=tmp;}

    return ret;

}

il int dinic(){ri ret=;while(bfs()){rp(i,S,T)cur[i]=head[i];while(int d=dfs(S,inf))ret+=d;}return ret;}

int main()

{

    //freopen("3227.in","r",stdin);freopen("3227.out","w",stdout);

    memset(head,-,sizeof(head));p=read();q=read();r=read();d=read();S=;T=p*q*r+;

    rp(i,,p)rp(j,,q){ad(nam(i,j,),S,read());}

    rp(k,,r){rp(i,,p)rp(j,,q)ad(nam(i,j,k),nam(i,j,k-),read());}

    rp(i,,p)rp(j,,q)ad(T,nam(i,j,r),inf);

    rp(i,,p)

        rp(j,,q)

            rp(t,,)

            {

                ri to_x=i+mov_x[t],to_y=j+mov_y[t];

                if(!to_x || !to_y || to_x>p || to_y>q)continue;

                rp(k,,r)

                {

                    ri tmp=k-d;if(tmp>){ad(nam(to_x,to_y,tmp),nam(i,j,k),inf);}

                    tmp=k+d+;if(tmp<=r){ad(nam(i,j,k+),nam(to_x,to_y,tmp),inf);}

                }

            }

    printf("%d\n",dinic());

    return ;

}