#include<deque>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<set>
#include<cmath>
#include<tuple>
#include<string>
#include<chrono>
#include<functional>
#include<iterator>
#include<random>
#include<unordered_set>
#include<unordered_map>
#include<array>
#include<map>
#include<bitset>
#include<iomanip>
#include<list>
#include <numeric>
using namespace std;
typedef unsigned long long int ulint;
typedef long long int llint;
typedef long double lldo;
#define mp make_pair
#define mt make_tuple
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define fir first
#define sec second
#define res resize
#define ins insert
#define era erase
#define RE return 0
//ios::sync_with_stdio(false);
//std::cin.tie(0);
//<< setprecision(20)
const int mod=(int)1000000007;
const llint big=(llint)(2.19e15)+1;
const long double pai=3.141592653589793238462643383279502884197;
const long double ena=2.71828182845904523536;
const long double eps=1e-7;
template <class T,class U>bool mineq(T& a,U b){if(a>b){a=b;return true;}return false;}
template <class T,class U>bool maxeq(T& a,U b){if(a<b){a=b;return true;}return false;}
template <class T> void soun(T& ar)
{sort(ar.begin(),ar.end());ar.erase(unique(ar.begin(),ar.end()),ar.end());}
llint gcd(llint a,llint b){if(a%b==0){return b;}else{return gcd(b,a%b);}}
llint lcm(llint a,llint b){return a/gcd(a,b) *b;}
template<class T,class U> auto LB(T& ve,U in){return lower_bound(ve.begin(),ve.end(),in);}
template<class T,class U> auto UB(T& ve,U in){return upper_bound(ve.begin(),ve.end(),in);}
template<class T,class U> auto LBI(T& ve,U in){return LB(ve,in)-ve.begin();}
template<class T,class U> auto UBI(T& ve,U in){return UB(ve,in)-ve.begin();}
template<class T> void SO(T& ve){sort(ve.begin(),ve.end());}
template<class T> void REV(T& ve){reverse(ve.begin(),ve.end());}
using pli=pair<llint,int>;
using daic=priority_queue<pli,vector<pli>,greater<pli>>;
int main(void){
	llint H,W,K,i,j;cin>>H>>W>>K;
	vector<string>grd(H);
	for(i=0;i<H;i++){cin>>grd[i];}
	bool ta=false,yo=false;
	for(j=0;j<W;j++){if(grd[0][j]=='#'&&grd[H-1][j]=='#'){ta=true;break;}}
	for(i=0;i<H;i++){if(grd[i][0]=='#'&&grd[i][W-1]=='#'){yo=true;break;}}
	if(ta&&yo){cout<<"1\n";RE;}
	llint kaz=0,ans=1;
	for(i=0;i<H;i++){
		for(j=0;j<W;j++){
			if(grd[i][j]=='#'){kaz++;}
		}
	}
	if((!ta)&&(!yo)){
		for(int h=0;h<60;h++){
			if((1LL<<h)&(K-1)){ans*=kaz;ans%=mod;}
			kaz*=kaz;kaz%=mod;
		}
		cout<<ans<<endl;RE;
	}
	if(ta){
		vector<string> iti(W,string(H,'.'));
		for(i=0;i<H;i++){
			for(j=0;j<W;j++){
				iti[j][i]=grd[i][j];
			}
		}
		swap(H,W);
		grd=iti;
	}
	llint ren=0;
	for(i=0;i<H;i++){if(grd[i][0]=='#'&&grd[i][W-1]=='#'){ren++;}}
	//行列累乗
	llint fe=0;
	for(i=0;i<H;i++){
		for(j=0;j<W;j++){
			if(grd[i][j]=='#'&&(j+1==W||grd[i][j+1]=='.')){fe++;}
		}
	}
	//2*2の行列累乗を書く
	array<array<llint,2>,2> gyo;
	gyo[0][0]=kaz;gyo[0][1]=-1;
	gyo[1][0]=0; gyo[1][1]=ren;
	//このK乗...
	array<array<llint,2>,2> re={};
	re[0][0]=1;re[1][1]=1;int e;
	//cerr<<"K="<<K<<" "<<kaz<<" "<<ren<<endl;
	for(int h=0;h<60;h++){
		if((1LL<<h)&(K-1)){
			//cerr<<"de";
			array<array<llint,2>,2> iti={{0,0,0,0}};
			for(i=0;i<2;i++){
				for(j=0;j<2;j++){
					for(e=0;e<2;e++){
						iti[i][j]+=re[i][e]*gyo[e][j];
						iti[i][j]%=mod;
					}
				}
			}
			re=iti;
		}
		array<array<llint,2>,2> iati={{0,0,0,0}};
		for(i=0;i<2;i++){
			for(j=0;j<2;j++){
				for(e=0;e<2;e++){
					iati[i][j]+=gyo[i][e]*gyo[e][j];
					iati[i][j]%=mod;
				}
			}
		}
		gyo=iati;
	}
	//cerr<<"re"<<re[0][0]<<endl;
	llint tna=0;
	for(i=0;i<H;i++){
		for(j=0;j<W-1;j++){
			if(grd[i][j]=='#'&&grd[i][j+1]=='#'){tna++;}
		}
	}
	ans=(re[0][0]+re[0][1]*tna)%mod;
	if(ans<0){ans+=mod;}
	cout<<ans<<endl;
	RE;
}

Measured memory_bytes: 3328