double G = lgord( cosc[i], nlat); double GM = lgord( cosc[i], nlat-1); double GP = lgord( cosc[i], nlat+1); double GT = (cosc[i]*cosc[i]-1.0) / (A*GP-B*GM); double D = lgord( cosc[i], nlat-1 ); D = D * D * FI * FI; gaussw[i] = C *( FI-0.5 ) / D; cosc[next] = 0.0; double C = 2.0; double D = lgord(cosc[next], nlat-1); D = D*D*FI*FI; gaussw[next] = C *( FI-0.5 ) / D;
double G = lgord( cosc[i], nlat); double GM = lgord( cosc[i], nlat-1); double GP = lgord( cosc[i], nlat+1); double GT = (cosc[i]*cosc[i]-1.0) / (A*GP-B*GM); double D = lgord( cosc[i], nlat-1 ); D = D * D * FI * FI; gaussw[i] = C *( FI-0.5 ) / D; cosc[next] = 0.0; double C = 2.0; double D = lgord(cosc[next], nlat-1); D = D*D*FI*FI; gaussw[next] = C *( FI-0.5 ) / D;
double G = lgord( cosc[i], nlat); double GM = lgord( cosc[i], nlat-1); double GP = lgord( cosc[i], nlat+1); double GT = (cosc[i]*cosc[i]-1.0) / (A*GP-B*GM); double D = lgord( cosc[i], nlat-1 ); D = D * D * FI * FI; gaussw[i] = C *( FI-0.5 ) / D; cosc[next] = 0.0; double C = 2.0; double D = lgord(cosc[next], nlat-1); D = D*D*FI*FI; gaussw[next] = C *( FI-0.5 ) / D;
double G = lgord( cosc[i], nlat); double GM = lgord( cosc[i], nlat-1); double GP = lgord( cosc[i], nlat+1); double GT = (cosc[i]*cosc[i]-1.0) / (A*GP-B*GM); double D = lgord( cosc[i], nlat-1 ); D = D * D * FI * FI; gaussw[i] = C *( FI-0.5 ) / D; cosc[next] = 0.0; double C = 2.0; double D = lgord(cosc[next], nlat-1); D = D*D*FI*FI; gaussw[next] = C *( FI-0.5 ) / D;