#define MA_U(j)       mean1d(j,1,1)
#define MA_V(j)       mean1d(j,1,2)
#define MA_W(j)       mean1d(j,1,3)
#define MA_P(j)       mean1d(j,1,4)
#define MA_R(j)       mean1d(j,1,5)
#define MA_UU(j)      mean1d(j,1,6)
#define MA_VV(j)      mean1d(j,1,7)
#define MA_WW(j)      mean1d(j,1,8)
#define MA_PP(j)      mean1d(j,1,9)
#define MA_RR(j)      mean1d(j,1,10)
#define MA_UV(j)      mean1d(j,1,11)
#define MA_UW(j)      mean1d(j,1,12)
#define MA_VW(j)      mean1d(j,1,13)
#define MA_Ux(j)      mean1d(j,1,14)
#define MA_Uy(j)      mean1d(j,1,15)
#define MA_Uz(j)      mean1d(j,1,16)
#define MA_Vx(j)      mean1d(j,1,17)
#define MA_Vy(j)      mean1d(j,1,18)
#define MA_Vz(j)      mean1d(j,1,19)
#define MA_Wx(j)      mean1d(j,1,20) 
#define MA_Wy(j)      mean1d(j,1,21) 
#define MA_Wz(j)      mean1d(j,1,22) 
#define MA_Ux2(j)     mean1d(j,1,23)
#define MA_Uy2(j)     mean1d(j,1,24)
#define MA_Uz2(j)     mean1d(j,1,25)
#define MA_Vx2(j)     mean1d(j,1,26) 
#define MA_Vy2(j)     mean1d(j,1,27) 
#define MA_Vz2(j)     mean1d(j,1,28) 
#define MA_Wx2(j)     mean1d(j,1,29) 
#define MA_Wy2(j)     mean1d(j,1,30) 
#define MA_Wz2(j)     mean1d(j,1,31)
#define MA_VxUy(j)    mean1d(j,1,32) 
#define MA_WxUz(j)    mean1d(j,1,33) 
#define MA_WyVz(j)    mean1d(j,1,34) 
#define MA_UXVY(j)    mean1d(j,1,35)
#define MA_UxWz(j)    mean1d(j,1,36) 
#define MA_VyWz(j)    mean1d(j,1,37)
#define MA_T(j)       mean1d(j,1,38)
#define MA_TT(j)      mean1d(j,1,39)
#define MA_RU(j)      mean1d(j,1,40)
#define MA_RV(j)      mean1d(j,1,41)
#define MA_RW(j)      mean1d(j,1,42)
#define MA_RP(j)      mean1d(j,1,43)
#define MA_RT(j)      mean1d(j,1,44)
#define MA_RUU(j)     mean1d(j,1,45)
#define MA_RVV(j)     mean1d(j,1,46)
#define MA_RWW(j)     mean1d(j,1,47)
#define MA_RUV(j)     mean1d(j,1,48)
#define MA_RUW(j)     mean1d(j,1,49)
#define MA_RVW(j)     mean1d(j,1,50)
#define MA_PU(j)      mean1d(j,1,51)
#define MA_PV(j)      mean1d(j,1,52)
#define MA_PW(j)      mean1d(j,1,53)
#define MA_TU(j)      mean1d(j,1,54)
#define MA_TV(j)      mean1d(j,1,55)
#define MA_TW(j)      mean1d(j,1,56)
#define MA_VIS(j)     mean1d(j,1,57)
#define MA_VIS2(j)    mean1d(j,1,58)
#define MA_Rx(j)      mean1d(j,1,59) 
#define MA_Ry(j)      mean1d(j,1,60) 
#define MA_Rz(j)      mean1d(j,1,61) 
#define MA_Px(j)      mean1d(j,1,62)
#define MA_Py(j)      mean1d(j,1,63)
#define MA_Pz(j)      mean1d(j,1,64)
#define MA_Tx(j)      mean1d(j,1,65)
#define MA_Ty(j)      mean1d(j,1,66)
#define MA_Tz(j)      mean1d(j,1,67)
#define MA_Tkk(j)     mean1d(j,1,68)
#define MA_UkPk(j)    mean1d(j,1,69)
#define MA_RTT(j)     mean1d(j,1,70)
#define MA_URx(j)     mean1d(j,1,71) 
#define MA_URy(j)     mean1d(j,1,72) 
#define MA_URz(j)     mean1d(j,1,73) 
#define MA_VRx(j)     mean1d(j,1,74)
#define MA_VRy(j)     mean1d(j,1,75)
#define MA_VRz(j)     mean1d(j,1,76)
#define MA_WRx(j)     mean1d(j,1,77)
#define MA_WRy(j)     mean1d(j,1,78)
#define MA_WRz(j)     mean1d(j,1,79)
#define MA_TAUXkk(j)  mean1d(j,1,80)
#define MA_TAUYkk(j)  mean1d(j,1,81)
#define MA_TAUZkk(j)  mean1d(j,1,82)

#define MA_RUUx(j)    mean1d(j,1,84)
#define MA_RUUy(j)    mean1d(j,1,85)
#define MA_RUUz(j)    mean1d(j,1,86)
#define MA_RVVx(j)    mean1d(j,1,87)
#define MA_RVVy(j)    mean1d(j,1,88)
#define MA_RVVz(j)    mean1d(j,1,89)
#define MA_RWWx(j)    mean1d(j,1,90)
#define MA_RWWy(j)    mean1d(j,1,91)
#define MA_RWWz(j)    mean1d(j,1,92)
#define MA_RUVx(j)    mean1d(j,1,93)
#define MA_RUVy(j)    mean1d(j,1,94)
#define MA_RUVz(j)    mean1d(j,1,95)
#define MA_RUWx(j)    mean1d(j,1,96)
#define MA_RUWy(j)    mean1d(j,1,97)
#define MA_RUWz(j)    mean1d(j,1,98)
#define MA_RVWx(j)    mean1d(j,1,99)
#define MA_RVWy(j)    mean1d(j,1,100)
#define MA_RVWz(j)    mean1d(j,1,101)
#define MA_PUx(j)     mean1d(j,1,102)
#define MA_PUy(j)     mean1d(j,1,103)
#define MA_PUz(j)     mean1d(j,1,104)
#define MA_PVX(j)     mean1d(j,1,105)
#define MA_PVY(j)     mean1d(j,1,106)
#define MA_PVZ(j)     mean1d(j,1,107)
#define MA_PWx(j)     mean1d(j,1,108)
#define MA_PWY(j)     mean1d(j,1,109)
#define MA_PWz(j)     mean1d(j,1,110)
#define MA_UPx(j)     mean1d(j,1,111)
#define MA_UPy(j)     mean1d(j,1,112)
#define MA_UPz(j)     mean1d(j,1,113)
#define MA_VPx(j)     mean1d(j,1,114)
#define MA_VPy(j)     mean1d(j,1,115)
#define MA_VPz(j)     mean1d(j,1,116)
#define MA_WPx(j)     mean1d(j,1,117)
#define MA_WPy(j)     mean1d(j,1,118)
#define MA_WPz(j)     mean1d(j,1,119)
#define MA_TAUXkUk(j) mean1d(j,1,120)
#define MA_TAUYkUk(j) mean1d(j,1,121)
#define MA_TAUZkUk(j) mean1d(j,1,122)
#define MA_TAUXkVk(j) mean1d(j,1,123)
#define MA_TAUYkVk(j) mean1d(j,1,124)
#define MA_TAUZkVk(j) mean1d(j,1,125)
#define MA_TAUXkWk(j) mean1d(j,1,126)
#define MA_TAUYkWk(j) mean1d(j,1,127)
#define MA_TAUZkWk(j) mean1d(j,1,128)
#define MA_UTAUXkk(j) mean1d(j,1,129)
#define MA_VTAUXkk(j) mean1d(j,1,130)
#define MA_WTAUXkk(j) mean1d(j,1,131)
#define MA_UTAUYkk(j) mean1d(j,1,132)
#define MA_VTAUYkk(j) mean1d(j,1,133)
#define MA_WTAUYkk(j) mean1d(j,1,134)
#define MA_UTAUZkk(j) mean1d(j,1,135)
#define MA_VTAUZkk(j) mean1d(j,1,136)
#define MA_WTAUZkk(j) mean1d(j,1,137)
#define MA_RUUUkk(j)  mean1d(j,1,138)
#define MA_RVVUkk(j)  mean1d(j,1,139)
#define MA_RWWUkk(j)  mean1d(j,1,140)
#define MA_RUVUkk(j)  mean1d(j,1,141)
#define MA_RUWUkk(j)  mean1d(j,1,142)
#define MA_RVWUkk(j)  mean1d(j,1,143)

#define MA_ZPx(j)     mean1d(j,1,153)
#define MA_ZPy(j)     mean1d(j,1,154)
#define MA_ZPz(j)     mean1d(j,1,155)
#define MA_RUx(j)     mean1d(j,1,156)
#define MA_RUy(j)     mean1d(j,1,157)
#define MA_RUz(j)     mean1d(j,1,158)
#define MA_RVx(j)     mean1d(j,1,159)
#define MA_RVy(j)     mean1d(j,1,160)
#define MA_RVz(j)     mean1d(j,1,161) 
#define MA_RWx(j)     mean1d(j,1,162) 
#define MA_RWy(j)     mean1d(j,1,163) 
#define MA_RWz(j)     mean1d(j,1,164)
#define MA_ZTAUXkk(j) mean1d(j,1,165)
#define MA_ZTAUYkk(j) mean1d(j,1,166)
#define MA_ZTAUZkk(j) mean1d(j,1,167)

#define MA_PHI1(j)    mean1d(j,1,171)
#define MA_PHI4(j)    mean1d(j,1,172)
#define MA_PHI6(j)    mean1d(j,1,173)
#define MA_PTx(j)     mean1d(j,1,174)
#define MA_PTy(j)     mean1d(j,1,175)
#define MA_PTz(j)     mean1d(j,1,176)
#define MA_TPx(j)     mean1d(j,1,177)
#define MA_TPy(j)     mean1d(j,1,178)
#define MA_TPz(j)     mean1d(j,1,179)
#define MA_UTkk(j)    mean1d(j,1,180)
#define MA_VTkk(j)    mean1d(j,1,181)
#define MA_WTkk(j)    mean1d(j,1,182)
#define MA_TTAUXkk(j) mean1d(j,1,183)
#define MA_TTAUYkk(j) mean1d(j,1,184)
#define MA_TTAUZkk(j) mean1d(j,1,185)
#define MA_PHI2(j)    mean1d(j,1,186)
#define MA_PHI3(j)    mean1d(j,1,187)
#define MA_PHI5(j)    mean1d(j,1,188)
#define MA_UPHI(j)    mean1d(j,1,189)
#define MA_VPHI(j)    mean1d(j,1,190)
#define MA_WPHI(j)    mean1d(j,1,191)
#define MA_UUX(j)     mean1d(j,1,192)
#define MA_UUY(j)     mean1d(j,1,193)
#define MA_UUZ(j)     mean1d(j,1,194)
#define MA_UVX(j)     mean1d(j,1,195)
#define MA_UVY(j)     mean1d(j,1,196)
#define MA_UVZ(j)     mean1d(j,1,197)
#define MA_UWx(j)     mean1d(j,1,198) 
#define MA_UWY(j)     mean1d(j,1,199)
#define MA_UWz(j)     mean1d(j,1,200) 
#define MA_VUX(j)     mean1d(j,1,201)
#define MA_VUY(j)     mean1d(j,1,202)
#define MA_VUZ(j)     mean1d(j,1,203)
#define MA_VVx(j)     mean1d(j,1,204) 
#define MA_VVy(j)     mean1d(j,1,205) 
#define MA_VVZ(j)     mean1d(j,1,206)
#define MA_VWx(j)     mean1d(j,1,207) 
#define MA_VWY(j)     mean1d(j,1,208)
#define MA_VWz(j)     mean1d(j,1,209) 
#define MA_WUX(j)     mean1d(j,1,210)
#define MA_WUY(j)     mean1d(j,1,211)
#define MA_WUZ(j)     mean1d(j,1,212)
#define MA_WVx(j)     mean1d(j,1,213) 
#define MA_WVy(j)     mean1d(j,1,214)
#define MA_WVZ(j)     mean1d(j,1,215)
#define MA_WWx(j)     mean1d(j,1,216) 
#define MA_WWY(j)     mean1d(j,1,217)
#define MA_WWz(j)     mean1d(j,1,218) 
#define MA_TAUxx(j)   mean1d(j,1,219)
#define MA_TAUyy(j)   mean1d(j,1,220)
#define MA_TAUzz(j)   mean1d(j,1,221)
#define MA_TAUxy(j)   mean1d(j,1,222)
#define MA_TAUxz(j)   mean1d(j,1,223)
#define MA_TAUyz(j)   mean1d(j,1,224)

#define MA_FLT_RU(j)  mean1d(j,1,225)
#define MA_FLT_RV(j)  mean1d(j,1,226)
#define MA_FLT_RW(j)  mean1d(j,1,227)
#define MA_FLT_RUU(j) mean1d(j,1,228)
#define MA_FLT_RVV(j) mean1d(j,1,229)
#define MA_FLT_RWW(j) mean1d(j,1,230)

#define MA_R3(j)      mean1d(j,1,231)
#define MA_R4(j)      mean1d(j,1,232)
#define MA_U3(j)      mean1d(j,1,233)
#define MA_U4(j)      mean1d(j,1,234)
#define MA_V3(j)      mean1d(j,1,235)
#define MA_V4(j)      mean1d(j,1,236)
#define MA_W3(j)      mean1d(j,1,237)
#define MA_W4(j)      mean1d(j,1,238)
#define MA_P3(j)      mean1d(j,1,239)
#define MA_P4(j)      mean1d(j,1,240)
#define MA_T3(j)      mean1d(j,1,241)
#define MA_T4(j)      mean1d(j,1,242)

#define MA_RUTTx(j)   mean1d(j,1,243)
#define MA_RVTTy(j)   mean1d(j,1,244)
#define MA_RWTTz(j)   mean1d(j,1,245)
#define MA_PTUx(j)    mean1d(j,1,246)
#define MA_PTVy(j)    mean1d(j,1,247)
#define MA_PTWz(j)    mean1d(j,1,248)
#define MA_TUx(j)     mean1d(j,1,249)
#define MA_TVy(j)     mean1d(j,1,250)
#define MA_TWz(j)     mean1d(j,1,251)

#define MA_MOMENTUM_SIZE 251

! Scalar

#define MS_S(j)       mean1d_sc(j,1,1)
#define MS_S2(j)      mean1d_sc(j,1,2)
#define MS_SU(j)      mean1d_sc(j,1,3)
#define MS_SV(j)      mean1d_sc(j,1,4)
#define MS_SW(j)      mean1d_sc(j,1,5)
#define MS_RS(j)      mean1d_sc(j,1,6)

#define MS_RSU(j)     mean1d_sc(j,1,7)
#define MS_RSV(j)     mean1d_sc(j,1,8)
#define MS_RSW(j)     mean1d_sc(j,1,9)

#define MS_RSS(j)     mean1d_sc(j,1,10)

#define MS_RSx(j)     mean1d_sc(j,1,11)
#define MS_RSy(j)     mean1d_sc(j,1,12)
#define MS_SRx(j)     mean1d_sc(j,1,13)
#define MS_SRy(j)     mean1d_sc(j,1,14)
#define MS_RSSx(j)    mean1d_sc(j,1,15)
#define MS_RSSy(j)    mean1d_sc(j,1,16)

#define MS_RUSS(j)    mean1d_sc(j,1,17)
#define MS_RVSS(j)    mean1d_sc(j,1,18)
#define MS_RWSS(j)    mean1d_sc(j,1,19)

#define MS_RSUx(j)    mean1d_sc(j,1,20)
#define MS_RSVy(j)    mean1d_sc(j,1,21)
#define MS_RSSUx(j)   mean1d_sc(j,1,22)
#define MS_RSSVy(j)   mean1d_sc(j,1,23)

#define MS_Sx(j)      mean1d_sc(j,1,24)
#define MS_Sy(j)      mean1d_sc(j,1,25)

#define MS_Fx(j)      mean1d_sc(j,1,26)
#define MS_Fy(j)      mean1d_sc(j,1,27)
#define MS_Fz(j)      mean1d_sc(j,1,28)
#define MS_SEPS(j)    mean1d_sc(j,1,29)

#define MS_SFx(j)     mean1d_sc(j,1,30)
#define MS_SFy(j)     mean1d_sc(j,1,31)
#define MS_Fxx(j)     mean1d_sc(j,1,32)
#define MS_Fyy(j)     mean1d_sc(j,1,33)

#define MS_FkdkS(j)    mean1d_sc(j,1,34)

#define MS_RSUy(j)    mean1d_sc(j,1,35)
#define MS_RSVx(j)    mean1d_sc(j,1,36)
#define MS_RSWx(j)    mean1d_sc(j,1,37)
#define MS_RSWy(j)    mean1d_sc(j,1,38)

#define MS_FkUk(j)    mean1d_sc(j,1,39)
#define MS_FkVk(j)    mean1d_sc(j,1,40)
#define MS_FkWk(j)    mean1d_sc(j,1,41)

#define MS_TAUxkSk(j) mean1d_sc(j,1,42)
#define MS_TAUykSk(j) mean1d_sc(j,1,43)
#define MS_TAUzkSk(j) mean1d_sc(j,1,44)

#define MS_PSx(j)     mean1d_sc(j,1,45)
#define MS_PSy(j)     mean1d_sc(j,1,46)
#define MS_PSz(j)     mean1d_sc(j,1,47)

#define MS_RUUS(j)    mean1d_sc(j,1,48)
#define MS_RVVS(j)    mean1d_sc(j,1,49)
#define MS_RWWS(j)    mean1d_sc(j,1,50)
#define MS_RUVS(j)    mean1d_sc(j,1,51)
#define MS_RUWS(j)    mean1d_sc(j,1,52)
#define MS_RVWS(j)    mean1d_sc(j,1,53)

#define MS_RUUSx(j)    mean1d_sc(j,1,54)
#define MS_RUVSx(j)    mean1d_sc(j,1,55)
#define MS_RUWSx(j)    mean1d_sc(j,1,56)
#define MS_RUVSy(j)    mean1d_sc(j,1,57)
#define MS_RVVSy(j)    mean1d_sc(j,1,58)
#define MS_RVWSy(j)    mean1d_sc(j,1,59)

#define MS_FkdkU(j)    mean1d_sc(j,1,60)
#define MS_FkdkV(j)    mean1d_sc(j,1,61)
#define MS_FkdkW(j)    mean1d_sc(j,1,62)

#define MS_TAUxxx(j)   mean1d_sc(j,1,63)
#define MS_TAUxyy(j)   mean1d_sc(j,1,64)
#define MS_TAUxyx(j)   mean1d_sc(j,1,65)
#define MS_TAUyyy(j)   mean1d_sc(j,1,66)
#define MS_TAUxzx(j)   mean1d_sc(j,1,67)
#define MS_TAUyzy(j)   mean1d_sc(j,1,68)
#define MS_STAUxkk(j)  mean1d_sc(j,1,69)
#define MS_STAUykk(j)  mean1d_sc(j,1,70)
#define MS_STAUzkk(j)  mean1d_sc(j,1,71)
#define MS_SPx(j)      mean1d_sc(j,1,72)
#define MS_SPy(j)      mean1d_sc(j,1,73)
#define MS_SPz(j)      mean1d_sc(j,1,74)

#define MS_GAMMA(j)    mean1d_sc(j,1,75)
#define MS_S3(j)       mean1d_sc(j,1,76)
#define MS_S4(j)       mean1d_sc(j,1,77)

#define MS_SCALAR_SIZE 77

