求助！LBM计算线性算子矩阵出现错误

GGGjjjjj

正在调试正压模式，执行到计算线性算子矩阵computing linear operator matrix时，执行csh文件时输出不了dat文件，为什么？
我的csh文件内容如下：
setenv LNHOME   /g6/zuojq/LBM/ln_solver.wm    # ROOT of LBM
setenv LBMDIR   $LNHOME/model                          # ROOT of model
setenv SYSTEM   pi                                     # execute system
setenv RUN      $LBMDIR/bin/$SYSTEM/lbm2.t42ml1cbaro   # Excutable file
setenv DIR      $LNHOME/data                           # Output Dir
setenv BSFILE   $LNHOME/bs/gt3/ncep.bs.t42l1           # Atm. BS File
setenv MATDIR   $DIR/mat                               # matrix directory
setenv MATNAM   MATBR.t42l1.dat                        # output matrix file name
ls -l $LNHOME/bs/gt3/ncep.bs.t42l1
#
#
setenv NLEVS     1  # no. of vertical levels
setenv NWSTR     0   # first zonal wave number
setenv NWEND    42   # truncate wave number
setenv KORD      4
setenv KTEF      1
#
setenv STM       1
setenv STD       1
#
if (! -e $DIR) mkdir -p $DIR
if (! -e $DIR/out) mkdir -p $DIR/out
if (! -e $MATDIR) mkdir -p $MATDIR
#
#/bin/rm -f $DIR/out/SYSOUT* $DIR/out/ERROUT
#/bin/rm -f $MATDIR/$MATNAM
#
echo ' '
echo ------ job started at `date`
echo ------ Error message  > $DIR/out/ERROUT
#
#       resolution      
#
@ maxw = 41 #$NWEND
@ izwv = $NWSTR
while ( $izwv <= $maxw )
#
#       loop for a unity (NMDIM) of vor.
#
echo VOR zonal wave M is $izwv
#
#      setup
#
cd $DIR
echo job started at `date` > $DIR/out/SYSOUT$izwv
/bin/rm -f $DIR/out/SYSIN
#
#      parameters
#
cat << END_OF_DATA >>! $DIR/out/SYSIN
&nmrun  run='make matrix'                                     &end

&nmtime start=0,$STM,$STD,0,0,0, end=0,$STM,$STD,0,30,0       &end

&nmdelt delt=40, tunit='MIN', inistp=0                        &end

&nmhdif order=$KORD, tefold=$KTEF, tunit='DAY'                &end

&nmdamp ddragv=10,
         tunit='DAY'                                           &end

&nmzmfct zmfctv=5.                                           &end

&nminit file='$BSFILE',
         DTBFR=1., DTAFTR=1., TUNIT='YEAR'                     &end

&nmchck ocheck=t, ockall=f                                    &end

&nmamat amatf='$MATDIR/$MATNAM', kmatw=$izwv                  &end
END_OF_DATA
echo "Processing wave $izwv, output to $MATDIR/$MATNAM"
#
#  run
#
( $RUN < $DIR/out/SYSIN >> $DIR/out/SYSOUT$izwv ) >& $DIR/out/ERROUT
#$RUN < $DIR/out/SYSIN >> $DIR/out/SYSOUT$izwv
echo job end at `date` >> $DIR/out/SYSOUT$izwv
#
   @ izwv++
end
#
echo ------ redist end at `date`
echo ' '