C1061: 编译器限制 : 块嵌套太深...怎么解决？？？

[zysg1314]

编译程序时出现了这个错误，希望大家给出回复。。。

#include "StdAfx.h"
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <vector>
#include <iomanip>
#include <iostream>
using namespace std;

double *TDMAsolver2(double*,double *,double *,double*,double *,double *,double*,int,int);
double PF (double a ,int b)//求平方
{
double jc=1;
for (int i=0;i<=b;i++)//是否有错？？
{
jc *=a;
}
return jc;
}


#define Do 1.8*PF(10,10)
#define q 0.001
#define g 9.8

int k,i,j,I,J;
int M=26,N=441;
int m=26,n=441;




double max (double e,double f)//比较大小
{
if(e>=f)
return e;
else
return f;
}

double yJ(int j,double yy) // 用J求y,在求冷却速度R时需要调用
{ double y;
y=(0.8-yy/2)+j*yy;
return y;
}


int main()
{

vector<vector<double> > gs(M+1,vector<double>(N+1));
vector<vector<double> > P(M+1,vector<double>(N+1));
for(int I=0;I<=M;I++)
for(int J=0;J<=N;J++)
P[I][J]=1;
vector<vector<double> > r(M+1,vector<double>(N+1));//r代表ρ
vector<vector<double> > Ue(M+1,vector<double>(N+1));//Ue代表Ueff，Ueff=Ul
for(int I=0;I<=M;I++)
for(int J=0;J<=N;J++)
Ue[I][J]=0.0062;
vector<vector<double> > rL(M+1,vector<double>(N+1));//rL代表ρL
for(int I=0;I<=M;I++)
for(int J=0;J<=N;J++)
rL[I][J]=6800;
vector<vector<double> > Vv(M+1,vector<double>(n+1));
vector<vector<double> > V(M+1,vector<double>(n+1));
for (int I=0,j=1;j<=n;j++)
V[I][j] =-0.02;
for (int I=1,j=1;j<=n;j++)
V[I][j] =-0.02;
for (int I=2,j=1;j<=n;j++)
V[I][j] =-0.02;
for (int I=3,j=1;j<=n;j++)
V[I][j] =-0.0326;
for (int I=4,j=1;j<=n;j++)
V[I][j] =-0.0334;
for (int I=5,j=1;j<=n;j++)
V[I][j] =-0.0342;
for (int I=6,j=1;j<=n;j++)
V[I][j] =-0.03505;
for (int I=7,j=1;j<=n;j++)
V[I][j] =-0.0359;
for (int I=8,j=1;j<=n;j++)
V[I][j] =-0.0367;
for (int I=9,j=1;j<=n;j++)
V[I][j] =-0.0375;
for (int I=10,j=1;j<=n;j++)
V[I][j] =-0.03825;
for (int I=11,j=1;j<=n;j++)
V[I][j] =-0.039;
for (int I=12,j=1;j<=n;j++)
V[I][j] =-0.03975;
for (int I=13,j=1;j<=n;j++)
V[I][j] =-0.0404;
for (int I=14,j=1;j<=n;j++)
V[I][j] =-0.041;
for (int I=15,j=1;j<=n;j++)
V[I][j] =-0.0416;
for (int I=16,j=1;j<=n;j++)
V[I][j] =-0.04215;
for (int I=17,j=1;j<=n;j++)
V[I][j] =-0.0426;
for (int I=18,j=1;j<=n;j++)
V[I][j] =-0.043;
for (int I=19,j=1;j<=n;j++)
V[I][j] =-0.0434;
for (int I=20,j=1;j<=n;j++)
V[I][j] =-0.04375;
for (int I=21,j=1;j<=n;j++)
V[I][j] =-0.044;
for (int I=22,j=1;j<=n;j++)
V[I][j] =-0.0442;
for (int I=23,j=1;j<=n;j++)
V[I][j] =-0.0444;
for (int I=24,j=1;j<=n;j++)
V[I][j] =-0.0445;
for (int I=25,j=1;j<=n;j++)
V[I][j] =-0.04455;
for (int I=26,j=1;j<=n;j++)
V[I][j] =-0.04455;


vector<vector<double> > Uu(m+1,vector<double>(N+1));
vector<vector<double> > U(m+1,vector<double>(N+1));/*U代表u,x方向的速度*/
for (int i=1,J=0;J<=N;J++) // 测试用给Uin初值
U[i][J] =0;
for (int i=2,J=0;J<=N;J++)
U[i][J] =0;
for (int i=3,J=0;J<=N;J++)
U[i][J] =0.00237;
for (int i=4,J=0;J<=N;J++)
U[i][J] =0.00232;
for (int i=5,J=0;J<=N;J++)
U[i][J] =0.00229;
for (int i=6,J=0;J<=N;J++)
U[i][J] =0.00222;
for (int i=7,J=0;J<=N;J++)
U[i][J] =0.00214;
for (int i=8,J=0;J<=N;J++)
U[i][J] =0.00204;
for (int i=9,J=0;J<=N;J++)
U[i][J] =0.00194;
for (int i=10,J=0;J<=N;J++)
U[i][J] =0.00183;
for (int i=11,J=0;J<=N;J++)
U[i][J] =0.00171;
for (int i=12,J=0;J<=N;J++)
U[i][J] =0.0016;
for (int i=13,J=0;J<=N;J++)
U[i][J] =0.00148;
for (int i=14,J=0;J<=N;J++)
U[i][J] =0.00136;
for (int i=15,J=0;J<=N;J++)
U[i][J] =0.00124;
for (int i=16,J=0;J<=N;J++)
U[i][J] =0.00112;
for (int i=17,J=0;J<=N;J++)
U[i][J] =0.00101;
for (int i=18,J=0;J<=N;J++)
U[i][J] =0.000893;
for (int i=19,J=0;J<=N;J++)
U[i][J] =0.000779;
for (int i=20,J=0;J<=N;J++)
U[i][J] =0.000666;
for (int i=21,J=0;J<=N;J++)
U[i][J] =0.000555;
for (int i=22,J=0;J<=N;J++)
U[i][J] =0.000444;
for (int i=23,J=0;J<=N;J++)
U[i][J] =0.000335;
for (int i=24,J=0;J<=N;J++)
U[i][J] =0.000227;
for (int i=25,J=0;J<=N;J++)
U[i][J] =0.000119 ;
for (int i=26,J=0;J<=N;J++)
U[i][J] =1.17e-005;
vector<vector<double> > T(M+1,vector<double>(N+1));
for (int I=0,J=0;J<=N;J++) // 测试用给T初值
T[I][J] =865.97;
for (int I=1,J=0;J<=N;J++)
T[I][J] =1038;
for (int I=2,J=0;J<=N;J++)
T[I][J] =1198.4;
for (int I=3,J=0;J<=N;J++)
T[I][J] =1332.1;
for (int I=4,J=0;J<=N;J++)
T[I][J] =1436.9;
for (int I=5,J=0;J<=N;J++)
T[I][J] =1493.8;
for (int I=6,J=0;J<=N;J++)
T[I][J] =1508;
for (int I=7,J=0;J<=N;J++)
T[I][J] =1512.8;
for (int I=8,J=0;J<=N;J++)
T[I][J] =1516.6;
for (int I=9,J=0;J<=N;J++)
T[I][J] =1519.6;
for (int I=10,J=0;J<=N;J++)
T[I][J] =1522.2;
for (int I=11,J=0;J<=N;J++)
T[I][J] =1524.4;
for (int I=12,J=0;J<=N;J++)
T[I][J] =1526.2;
for (int I=13,J=0;J<=N;J++)
T[I][J] =1527.8;
for (int I=14,J=0;J<=N;J++)
T[I][J] =1529.1;
for (int I=15,J=0;J<=N;J++)
T[I][J] =1530.2;
for (int I=16,J=0;J<=N;J++)
T[I][J] =1531.1;
for (int I=17,J=0;J<=N;J++)
T[I][J] =1531.9;
for (int I=18,J=0;J<=N;J++)
T[I][J] =1532.5;
for (int I=19,J=0;J<=N;J++)
T[I][J] =1533;
for (int I=20,J=0;J<=N;J++)
T[I][J] =1533.4;
for (int I=21,J=0;J<=N;J++)
T[I][J] =1533.7;
for (int I=22,J=0;J<=N;J++)
T[I][J] =1534;
for (int I=23,J=0;J<=N;J++)
T[I][J] =1534.2;
for (int I=24,J=0;J<=N;J++)
T[I][J] =1534.3;
for (int I=25,J=0;J<=N;J++)
T[I][J] =1534.3;
for (int I=26,J=0;J<=N;J++)
T[I][J] =1534.4;
初始场赋值的时候就出现了这个问题。。。。。。无语中，小弟新手