C **********************************************************
C
C PLOTS X-Z CROSS-SECTION OF ZONAL AVR  STR. FUN.       
C TRANSPORT ACROSS ATALTIC IS CALCULATED ON ZIG-ZAG LINES    
C (May 93- improvements by glm for MED)
C ### (Nov 96- correct error in zigzag integration) #####
c ******************** NSATL MODEL **************************
      PROGRAM PLXZ
C
C                  files  rec/file  1st  
C                      5      36
      PARAMETER (NFILE=1,NRPF=36,NSTRT=1 ,IPLTS=1)
C
C
      PARAMETER (IM=128,JM=128,KB=16, IMJM=IM+JM)
C
C                  Y,Z   grid for interpolation  
C if LY=JM integrate on model grid instead of lat.
      PARAMETER ( LY=113, LZ=35, LYB=LY+4)   
C     PARAMETER ( LY= JM, LZ=35, LYB=LY+4)   
      PARAMETER (AWEST=-100.,AEAST=30.,ASOUTH=-33.,ANORTH=80.)
C
C PARAMETERS AND ARRAYS FOR FILLING AREA
      PARAMETER (LMX2=2000,LMX3=3000)
      DIMENSION DST(LMX2),NST(LMX3),XBOT(LYB),YBOT(LYB),KBOT(LYB)
C
      DIMENSION AA(LY,LZ),AX(LY,LZ),AZ(LY,LZ),EY(LY),ZS(LZ)
      DIMENSIO