Here is the file test.inp, for 274 POPC lipids.(I made it following the example above, with a saved trajectory file for lipids without water and a psf file written in charm format.The number of atoms in trajectory file match those in psf ):

! the axial avg order parameter for an input C no. @C
! modified for C=C case; first file @FF last file @LF
! subdir @S


! Read initial data files
open unit 1 read form name top_all36_lipid.rtf
read rtf card unit 1
close unit 1
open unit 1 read form name par_all36_lipid.prm
read param card unit 1
close unit 1

! Read parameter file in command line

open read card unit 1 name popc.psf
read psf card unit 1
open read card unit 1 name popc.crd
read coor card unit 1

! number of lipids
set nl = 274
! number of files
calc n = 1

calc mxa @NL * 6 * 4
calc mxs 10 + @NL*2 * 4
calc mxt 2000 * @N

! unsaturation check; beta chain
set ub = 0
define unsat sele atom A 1 C2@C .and. chem CEL1 end
if ub lt ?NSEL set ub ?NSEL
! unsaturation check; gamma chain
set ug = 0
define unsat sele atom A 1 C3@C .and. chem CEL1 end
if ug lt ?NSEL set ug ?NSEL

correl maxtim @MXT maxa @MXA maxs @MXS
enter sr zero
if ub eq 0 enter ss dupl sr
enter sx dupl sr
if ug eq 0 enter sy dupl sr
set k 1
label elp
if ub eq 0 then
enter w@K vect z A @K H@{C}R A @K C2@C
enter a@K vect r A @K H@{C}R A @K C2@C
enter x@K vect z A @K H@{C}S A @K C2@C
enter b@K vect r A @K H@{C}S A @K C2@C
else
enter w@K vect z A @K H@{C}1 A @K C2@C
enter a@K vect r A @K H@{C}1 A @K C2@C
endif
if c gt 15 goto eskip
if ug eq 0 then
enter y@K vect z A @K H@{C}X A @K C3@C
enter c@K vect r A @K H@{C}X A @K C3@C
enter z@K vect z A @K H@{C}Y A @K C3@C
enter d@K vect r A @K H@{C}Y A @K C3@C
else
enter y@K vect z A @K H@{C}1 A @K C3@C
enter c@K vect r A @K H@{C}1 A @K C3@C
endif
label eskip
incr k by 1
if k le @NL goto elp

! OPEN FILES, FILL SERIES VIA TRAJ

open read unit 1 file name popclipidsnowater.dcd

! LOOP OVER LIPIDS; NORM, CALC SCD, ACCUM
set k 1
label clp
mantim w@K ratio a@K
mantim w@K squa
mantim w@K mult 1.5
mantim w@K shift -0.5
mantim sr add w@K
if ub eq 0 then
mantim x@K ratio b@K
mantim x@K squa
mantim x@K mult 1.5
mantim x@K shift -0.5
mantim ss add x@K
endif
if c gt 15 goto mskip
mantim y@K ratio c@K
mantim y@K squa
mantim y@K mult 1.5
mantim y@K shift -0.5
mantim sx add y@K
if ug eq 0 then
mantim z@K ratio d@K
mantim z@K squa
mantim z@K mult 1.5
mantim z@K shift -0.5
mantim sy add z@K
endif
label mskip
incr k by 1
if k le @NL goto clp

! use final mantime to set AVER, FLUC for output
open write unit 1 card name cd@{C}.txt
echu 1
mantim sr divi @NL
echo sr@C ?AVER ?FLUC
if ub eq 0 then
mantim ss divi @NL
echo ss@C ?AVER ?FLUC
endif
if c gt 15 goto oskip
mantim sx divi @NL
echo sx@C ?AVER ?FLUC
if ug eq 0 then
mantim sy divi @NL
echo sy@C ?AVER ?FLUC
endif
label oskip
close unit 1