#include <bits/stdc++.h>
const int mod=1e9+7;
using namespace std;
int n,m,K,f[2][4096][73],p=0,o=1,C[150][150],fac[150],_k[150];
char a[144][144],b[144][144];
void qmo(int &x){
	x+=(x>>31)&mod;
}
void init(int N){
	for(int i=0;i<=N;i++) for(int j=C[i][0]=1;j<=i;j++) qmo(C[i][j]=C[i-1][j]+C[i-1][j-1]-mod);
	for(int i=fac[0]=1;i<=N;i++) fac[i]=1ll*i*fac[i-1]%mod;
	for(int i=_k[0]=1;i<=N;i++) qmo(_k[i]=_k[i-1]+_k[i-1]-mod);
}
void sol(){
	scanf("%d%d%d",&n,&m,&K);
	for(int i=0;i<n;i++) scanf("%s",a[i]);
	if(n<m){
		for(int i=0;i<n;i++) for(int j=0;j<m;j++) b[j][i]=a[i][j];
		swap(a,b),swap(n,m);
	}
	int cnt=0;
	for(int i=0;i<n;i++) for(int j=0;j<m;j++) cnt+=a[i][j]=='.';
	for(int S=0;S<(1<<m);S++) for(int j=0;j<=K;j++) f[0][S][j]=0;
	f[0][0][0]=1,p=0,o=1;
	for(int i=0;i<n;i++) for(int j=0;j<m;j++){
		for(int S=0;S<(1<<m);S++) for(int j=0;j<=K;j++) f[o][S][j]=0;
		for(int S=0;S<(1<<m);S++) if(a[i][j]=='.'||!((S>>j)&1)){
			if((S>>j)&1) for(int k=0;k<K;k++) qmo(f[o][S^(1<<j)][k+1]+=f[p][S][k]-mod);
			else{
				for(int k=0;k<=K;k++) qmo(f[o][S][k]+=f[p][S][k]-mod);
				if(a[i][j]=='.') for(int k=0;k<=K;k++) qmo(f[o][S^(1<<j)][k]+=f[p][S][k]-mod);
				if(a[i][j]=='.'&&j&&((S>>(j-1))&1)){
					for(int k=0;k<K;k++) qmo(f[o][S^(1<<(j-1))][k+1]+=f[p][S][k]-mod);
				}
			}
		}
		o^=1,p^=1;
	}
	int res=0;
	for(int j=0;j<=K;j++){
		int w=1ll*C[cnt-2*j][2*(K-j)]*fac[2*(K-j)]%mod*C[K][j]%mod*fac[j]%mod*f[p][0][j]*_k[j]%mod;
		if(j&1) qmo(res-=w);
		else qmo(res+=w-mod);
	}
	printf("%d\n",res);
}
int main(){
	init(144);
	int T;
	scanf("%d",&T);
	while(T--) sol();
}