How did convert this code to the python

void srcn(short int** pd,int* nse,int il,int jl)
{
	int i,j,k,l;
	double std,avg,temp;
	for(i=rn;i<il-rn;i++){
		for(j=rn;j<jl-rn;j++){
			avg=0;
			temp=0;
			std=0;
			for(k=-rn;k<=rn;k++){
				for(l=-rn;l<=rn;l++){
					avg+=(double)pd[i+k][j+l];
				}
			}
			avg=avg/((double)sn*sn);
			for(k=-rn;k<=rn;k++){
				for(l=-rn;l<=rn;l++){
				temp+=((double)pd[i+k][j+l]-avg)*((double)pd[i+k][j+l]-avg);
				}
			}
			std=pow(temp/((double)sn*sn),0.5);
			if(std<c*drt){
				nse[0]++;
			}
			else{
				nse[1]++;
			}
		}
	}
	
}