Hi Neks,<br><br>I am trying to write a routine that averages over the y-z plane, yielding an array that varies in the x dimension.  I have seen the planar_average_z routine in navier5.f that averages over x-y, and I have a planar_average_y routine that Paul sent me.  So, I am trying to write planar_average_x.  I have been able to proceed by analogy, except for the two lines in which the weighting coefficients are computed.<br>

<br>I planar_average_z,<br><br>yy = (1.-zgm1(j,2))/2.  ! = 1 for j=1, = 0 for j=ny1<br>aa = yy*area(i,1,1,e) + (1-yy)*area(i,1,3,e)  ! wgtd jacobian, fc 1&3<br><br>In planar_aveage_y,<br><br>zz = (1.-zgm1(k,3))/2.  ! = 1 for k=1, = 0 for k=nz1<br>

aa = zz*area(i,1,5,e) + (1-zz)*area(i,1,6,e)  ! wgtd jacobian<br><br>Can someone tell me what the analogous lines would be for planar_average_x, as well as what the area tensor is?<br><br>Thanks very much,<br><br>David<br>

<br>