#
# QSUB -eo
# QSUB -s /bin/csh
# QSUB -r ds218 -q prem -lT 40


set file = Y04345

echo start "$user"$$
set x = $TMPDIR
set y = /usr/tmp/$USER
mkdir $y
set local = $y/$file
set remote = /DSS/$file
cd $x

cat << EOF > tsrc.f
      program ds2181
      dimension data(144,73)
C-        READ (1)IV,IM,(DATA(I,J),I=1,144),J=1,73)
C-        AS I=1,144, LONG =0.0 TO 357.5 E, EASTWARD BY 2.5
C-        AS J=1,73,  LAT =90N TO 90S BY 2.5
C-
C-
C-      WHERE IM IS MONTH, IV IS A VARIABLE/LEVEL CODE
C-          IV 1-8 IS U AT 8 LEVELS (1000,850,700,500,300,250,200,100)
C-            9-16 IS V AT 8 LEVELS (M/S)
C-           17-24 IS T AT 8 LEVELS (DEG K)
C-           25-32 IS Z AT 8 LEVELS (M)
C-              DATA IS ORDERED 12 MONTHS FOR IV=1, THEN 12 FOR IV=2, ETC.
C-
      nr=0
    5 read (11,end=9)iv,im,((data(i,j),i=1,144),j=1,73)
      nr=nr+1
C select 1000 mb Z
      if(iv.ne.25) go to 5
      print 1001,nr,iv,im,data(1,1),data(1,37),
     2 data(144,72),data(144,73)
 1001 format(1x,i8,2i3,4f10.2)
      if(nr.gt.0) go to 5
      do 20 ie=12,144,12
      is=ie-11
      print 1003,(il,il=is,ie)
 1003 format(5x,12i10)
      do 18 j=73,1,-1
      print 1002,j,(data(i,j),i=is,ie)
 1002 format(1x,i3,1x,12f10.2)
   18 continue
   20 continue
      go to 5
    9 continue
      print 1004,nr
 1004 format(' stop recs =',i8)
      end
EOF

cft77 -e h tsrc.f

echo begin loading
segldr tsrc.o -L /lib,/usr/lib,/usr/local/lib -l ncarg,ncarm,ncaro

set opt = "OPTIONS=RETRY,NOREPLACE"
LREAD FORMAT=TR $opt REMOTE=$remote LOCAL=$local

assign -a$local fort.11

a.out

ls -al