How to calculate order parameter in coarse-gained simulation of a complex membrane

GROMACS version: 2022
GROMACS modification: No
dear GROMACS users
I am working on a project using Martini 3 force field. I have simulated a membrane containing DOPC, DOPE, DOPS, POSM and CHOL, and I want to calculate the lipid tail order parameter for the inner and outer layers. However, I have encountered some difficulties in this process.
Could any one please share with me the methods and scripts that used to calculate this parameter in coarse-gained simulation?

Check MARTINI website tutorials.