Main dependencies of ja_anglecheck_swath_graphics execution
(using Makefile_vrrms_redhat as a guide):

rt      = r-theta coordinate (relative to storm center)
first   = first interpolation, etc.
fourier = wavenumbers 0 and 1 
guess   = 1st guess wind field used to help dealiasing schemes

********************************************************************************

./wind_interpolate_rt_auto_first_defreckle g

(The following are NOT used, but listed in makefile:
 tv_map_sub.o x_to_tiff_subs.o tiff_subs.o diver3.o 
 miscellaneous_scan.o ctme.o )

>>>PROGRAM LEASTWIND <--> dnewbigr_abcdonly_auto_first.f

    >>CALL GETARG <--> FORTRAN intrinsic

    >>CALL OPENRAMFILE <--> readramfile.f

         IF(fileid .eq. 'ram1')
          CALL CTMEL <--> readramfile.f
          CALL CTMEL <--> readramfile.f

    >>CALL ABCD <--> abcd_scan_rt4.f

        >IF(IGRAPHICS.EQ.1) 
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c

        >IF(IUSEWIND.EQ.1) 
          CALL READWIND <--> readwind.f
              
             CALL READHEADERS <--> readfile.f
             CALL LDPLN <--> readfile.f

        >CALL SWEEPIN <--> newsweepread_auto_first_disc_defreckle.f

            IF(ICALLSWEEPIN.EQ.1)
             IF(IUSEWIND.EQ.1)
              CALL VRLOCATER <--> newsweepread_auto_first_disc_defreckle.f            
              CALL VRLOCATER <--> newsweepread_auto_first_disc_defreckle.f
              CALL VRLOCATER <--> newsweepread_auto_first_disc_defreckle.f

            IF(ICALLSWEEPIN.EQ.0) 

             IF(RAMFL .EQ. 'Y') 
              CALL OPENRAMFILE <--> readramfile.f   
       
             IF(IATTEN .EQ. 'Y') 
              CALL MAKEATTENTABLE1 <--> makeattentable1.f

             IF(ITAPEOPEN.EQ.1) 
              CALL TAPEOPEN <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f
             ELSE
              CALL OPENRADARFILE <--> radardiscread1.f

            IF(ICALLSWEEPIN.EQ.0) 

             CALL READTAPE <--> newsweepread_auto_first_disc_defreckle.f

                IF(ITAPEOPEN.EQ.0)
                 CALL RADARDISCREAD1 <--> radardiscread1.f                   
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f
                ELSE
                 CALL RADARTAPEREAD <--> tapedummy.f

                CALL ANGLECOR <--> paulsanglecor.f

                CALL CONVERTANGLES <--> convertangles.f

                IF(RAMFL .EQ. 'Y') 

                 CALL READRAMFILE <--> readramfile.f

                    IF(irec.lt.2 .or. irec.gt.maxrec .or. timein.lt.0.0)
                     CALL CTMEL <--> readramfile.f

                    IF(rlat .le. -999.)
                     CALL CTMEL <--> readramfile.f
                     CALL CTMEL <--> readramfile.f

                CALL RAYCOR1 <--> raycor1.f

                   IF (iattenflag.EQ.1)
                    CALL INTERVENATTEN1 <--> intervenatten1.f

                   IF (iunfoldflag.EQ.1)
                    CALL UNFLD1 <--> unfld1.f

             IF(IFLIGHTCHECK.EQ.0) 
              CALL FINDIDAYS <--> findidays.f
              CALL FINDIDAYS <--> findidays.f

             IF(TIMEMSTIME.LT.-72000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             IF(LUNR.NE.80)

              CALL TAPESEARCH <--> tapedummy.f

             ELSE

              CALL DISCSEARCH <--> radardiscread1.f

                 IF(key .le. itime(kindradar))
                  CALL CTMEL <--> readramfile.f
                  CALL CTMEL <--> readramfile.f 

                 CALL DISCENTRY <--> radardiscread1.f
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f

                 IF(ioldsweep(kindradar) .ne. numsweep(kindradar))
                  CALL CTMEL <--> readramfile.f

             CALL READTAPE <--> newsweepread_auto_first_disc_defreckle.f

             IF((TIME-STIME).LT.-60000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

            CALL READTAPE <--> newsweepread_auto_first_disc_defreckle.f

            IF(TIME.GT.STIME.AND.IFILECHECK.EQ.0.AND.ITAPEOPEN.EQ.1)
             CALL TAPECONTROL <--> tapedummy.f
             CALL TAPECONTROL <--> tapedummy.f

            IF(NUM_RAYS.GT.0.AND.IUSEWIND.EQ.1)
             CALL VRLOCATER <--> newsweepread_auto_first_disc_defreckle.f

            IF(NUM_RAYS.GT.0)
             CALL VRLOCATER <--> newsweepread_auto_first_disc_defreckle.f

            IF((ISWEEPSET.EQ.2.AND.SWEEPTEST.LT.0.).OR.NUM_RAYS.EQ.800)
             CALL CTME1 <--> readramfile.f

             IF(TRAKAZM.GT.91..AND.TRAKAZM.LT.269.)
              CALL FIND_GROUND <--> newsweepread_auto_first_disc_defreckle.f

            IF(TRAKAZM.GT.91..AND.TRAKAZM.LT.269.)
             CALL FIND_GROUND <--> newsweepread_auto_first_disc_defreckle.f

            IF(MOD(ITIME,60).EQ.0 .AND. itime.NE.OLDTIME)
             CALL CTME1 <--> readramfile.f

            IF(RDTAIIN(I).GT.RATEST)THEN

             IF(J.LE.NUM_RAYS.AND.AZ_EL(1,J).LT.180.)
              CALL GROUNDCHECK <--> newsweepread_auto_first_disc_defreckle.f

             IF(J.GE.1)
              IF(AZ_EL(1,J).GT.180.)
               CALL GROUNDCHECK <--> newsweepread_auto_first_disc_defreckle.f

            CALL DEFRECKLE <--> defreckle.f

            CALL UNFLDNEW <--> unfldnew.f

            CALL AZIUNFLD <--> aziunfld.f

        >IF(IGRAPHICS.EQ.1)

          CALL SOLO_CONVT <--> solo_non_convt.f
             CALL GET_RADII <--> solo_non_convt.f 
             FUNCTION AZDIF <--> azdif.f
             CALL WCLIP <--> wclip.f
             CALL WCLIP <--> wclip.f
             CALL WEDGINT <--> wedgint.f
             CALL WCLIP <--> wclip.f
          CALL SHOW_IZ <--> show_iz.f
             CALL RASX_DISPLAY_IMAGE <--> bigger_raster_x_functions.c
             CALL X_TVLINE <--> bigger_raster_x_functions.c
          CALL SOLO_CONVT <--> solo_non_convt.f
          CALL SHOW_IZ <--> show_iz.f

        >IF(IRAMSFILES.EQ.1)
          CALL READRAMFILE <--> readramfile.f

        >IF(IREADABCD.NE.1.AND.TIMES.GE.TSTART.AND.TIMES.LE.TEND)
          CALL VRLOCATE <--> abcd_scan_rt4.f

    >>CALL WRITEABCD <--> files_rho4.f
    >>CALL WRITEABCD <--> files_rho4.f

cp SUM*first* SUMFILES/

+++++++++++++
OUTCOME: Summations of products needed for least squares wind 
         solution (derived from composited observations) are 
         written to a file. Formulated in r-theta(wv#0,1)-z 
         relative to storm center. Aircraft motion removal, 
         Bargen-Brown unfolding, defreckling, and surface 
         removal are accomplished.
+++++++++++++

********************************************************************************

./wind_fourier_auto_first

(The following are NOT used, but listed in makefile: 
 miscellaneous_scan_fourier.o diver3_rt.o iiposi.o) 
	
>>>PROGRAM LEASTWIND <--> dnewbigr_fourier_auto_first.f

    >>CALL WZERO <--> dnewbigr_fourier_auto_first.f
        >CALL READABCD <--> files_rho.f

    >>IF(IDBZONLY.NE.1)

       CALL SPARSE <--> sparse_rt_fourier_nofilled_auto.f

         >CALL FUN <--> sparse_rt_fourier_nofilled_auto.f

         >CALL DFUN <--> sparse_rt_fourier_nofilled_auto.f

         >CALL SBINI <--> itpack.f
             CALL IVFILL <--> itpack.f
             CALL IVFILL <--> itpack.f
             CALL VFILL <--> itpack.f

         >CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
             CALL SORTER <--> sparse_rt_fourier_nofilled_auto.f
             IF(INDX(I).GE.INDX(L))
              IF(INDX(I).GT.0.AND.INDX(L).GT.0)
               IF(IN.GT.0.AND.INN.GT.0)
                CALL SBSIJ <--> itpack.f

         >CALL READABCD <--> files_rho.f

         >IF(U2(I,K).GT.10.1)
           CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f

         >IF(U2(I,K).GT.FLAG)
           IF(FACTZZ.EQ.0.)
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           ELSEIF(K.EQ.1)
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
           IF(I.LE.1)
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f

         >IF(INPOSITION.GT.0..AND.INDEX4(I,J,K).NE.0.AND.
             ABS(XWEIGHT(I,J,K1)).LT.XACC.AND.
             ABS(YWEIGHT(I,J,K1)).LT.YACC.AND.
             ABS(ZWEIGHT(I,J,K1)).LT.ZACC.AND.
             SUMTESTER.GT..1.AND.ABS(SUM12(I,J,K1)).LT.SINTOL)
           IF(NSETSET.EQ.1) 
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f
            IF(IP.LE.NCOEF(K))
             IF(NP.LE.NCOEF(K))
              CALL SBSIJ <--> itpack.f

         >IF(U2(I,K).GT.FLAG)
           IF(INDEXNN1.GT.0.OR.INDEXNN2.GT.0)
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f
            CALL MATRIXSUM <--> sparse_rt_fourier_nofilled_auto.f

         >CALL SBEND <--> itpack.f

       CALL WRITEHEADER <--> files_rho.f

       CALL WRPLN <--> files_rho.f

       CALL SOLVESPARSE <--> solvesparse_fourier_auto.f

         >CALL DFAULT <--> itpack.f

         >CALL SSORCG <--> itpack.f
             CALL ECHOUT <--> itpack.f
             CALL ECHALL <--> itpack.f
             CALL IVFILL <--> itpack.f
             CALL VFILL <--> itpack.f
             CALL SBELM <--> itpack.f
             CALL IVFILL <--> itpack.f
             CALL PRBNDX <--> itpack.f
             CALL PERMAT <--> itpack.f
             CALL PERVEC <--> itpack.f
             CALL PERVEC <--> itpack.f
             CALL VFILL <--> itpack.f
             CALL SCAL <--> itpack.f
             CALL VFILL <--> itpack.f
             CALL OMEG <--> itpack.f
             CALL DCOPY <--> itpack.f
             CALL DCOPY <--> itpack.f
             CALL PFSOR <--> itpack.f
             CALL VEVMW <--> itpack.f
             CALL ITSRCG <--> itpack.f
             CALL ITSRCG <--> itpack.f
             CALL UNSCAL <--> itpack.f
             CALL PERMAT <--> itpack.f
             CALL PERVEC <--> itpack.f
             CALL PERVEC <--> itpack.f
             CALL PERRORR <--> itpack.f

       CALL WRITEHEADER <--> files_rho.f

       CALL WRPLN <--> files_rho.f
       
    >>CALL WRITEHEADER <--> files_rho.f

    >>CALL WRPLN <--> files_rho.f

+++++++++++++
OUTCOME: Solve for wv#0,1 Fourier coefficients of ur,vt,w
+++++++++++++

********************************************************************************

./indexwind_rt_fourier_auto

>>>PROGRAM FILLOLD <--> indexwind_rt_fourier_auto.f

    >>CALL READHEADER <--> indexwind_rt_fourier_auto.f

    >>CALL LDPLN <--> indexwind_rt_fourier_auto.f

    >>CALL NEWWIND <--> indexwind_rt_fourier_auto.f
        >CALL FUN <--> indexwind_rt_fourier_auto.f

    >>CALL WRITEHEADER <--> indexwind_rt_fourier_auto.f

    >>CALL DIVER3 <--> indexwind_rt_fourier_auto.f

    >>CALL WRPLN <--> indexwind_rt_fourier_auto.f

+++++++++++++
OUTCOME: Fourier coefficients used to construct ur,vt,w (r,theta,z)
+++++++++++++

********************************************************************************

./convert_rt_xy_fourier_auto

>>>PROGRAM CONVERT_RT_XY <--> convert_rt_xy_fourier_auto.f

    >>CALL READHEADER <--> convert_rt_xy_fourier_auto.f

    >>CALL LDPLN <--> convert_rt_xy_fourier_auto.f

    >>IF(IVRVTMEAN.GT.0)
       CALL RTMEAN <--> convert_rt_xy_fourier_auto.f

    >>IF(IUVMEAN.EQ.1)
       CALL UVMEAN <--> convert_rt_xy_fourier_auto.f
      ELSEIF(IUVMEAN.EQ.2)
       CALL LUVMEAN <--> convert_rt_xy_fourier_auto.f

    >>IF(IUVMEAN.EQ.1.AND.IVRVTMEAN.NE.1)
       CALL NUMTEST <--> convert_rt_xy_fourier_auto.f

    >>IF(INUMTEST.EQ.1)
       CALL NUMTEST <--> convert_rt_xy_fourier_auto.f

    >>CALL CONVERT <--> convert_rt_xy_fourier_auto.f
        >CALL POLAR_CART <--> convert_rt_xy_fourier_auto.f
        >CALL POLAR_CART <--> convert_rt_xy_fourier_auto.f
        >CALL POLAR_CART <--> convert_rt_xy_fourier_auto.f
        >CALL POLAR_CART <--> convert_rt_xy_fourier_auto.f
        >CALL POLAR_CART <--> convert_rt_xy_fourier_auto.f

    >>CALL WRITEHEADER <--> convert_rt_xy_fourier_auto.f

    >>CALL WRPLN <--> convert_rt_xy_fourier_auto.f

+++++++++++++
OUTCOME: Convert from ur,vt,w (r,theta,z) to u,v,w (Cartesian)
+++++++++++++

********************************************************************************

./convert_wind_fourier_auto

>>>PROGRAM FILLOLD <--> convert_wind_fourier_auto.f

    >>CALL READHEADER <--> convert_wind_fourier_auto.f  

    >>CALL LDPLN <--> convert_wind_fourier_auto.f

    >>IF(IREMOVE.GT.0)
       CALL REMOVE_RZMEAN <--> convert_wind_fourier_auto.f

    >>IF(IREMOVE.EQ.1)
       IF(IREMOV.EQ.1)
        CALL REMOVE_CONSTANT <--> convert_wind_fourier_auto.f
       ELSE
        CALL REMOVE_CONSTANT2 <--> convert_wind_fourier_auto.f

    >>CALL WRITEHEADER <--> convert_wind_fourier_auto.f

    >>IF(IUPDN.NE.2)
       CALL DIVER <--> convert_wind_fourier_auto.f
       CALL VMOTN <--> convert_wind_fourier_auto.f

    >>CALL WRPLN <--> convert_wind_fourier_auto.f

+++++++++++++
OUTCOME: Convert 32-bit real winds to 16-bit integer winds for use
         in HRD packages.
+++++++++++++

********************************************************************************

./anglecheck

(The following are NOT used, but listed in makefile:
 diver3.o) 

>>>PROGRAM LEASTWIND <--> dnewbig_interpolate_auto_guess_xyrt.f

    >>CALL GETARG <--> FORTRAN intrinsic

    >>IF(IRAMSFILES.EQ.1)
       CALL OPENRAMFILE <--> readramfile.f

          IF(fileid .eq. 'ram1')
           CALL CTMEL <--> readramfile.f
           CALL CTMEL <--> readramfile.f

    >>CALL ABCD <--> abcd_shell.f

        >IF(IUSEWIND.EQ.1)
          CALL READWIND <--> readwind_auto_guess_fourier.f
             CALL READHEADERS <--> readfile.f
             CALL LDPLN <--> readfile.f

        >CALL SWEEPIN <--> sweepread_anglecheck.f

            IF(ICALLSWEEPIN.EQ.1)
             IF(IUSEWIND.EQ.1)
              CALL VRLOCATER <--> sweepread_anglecheck.f          
              CALL VRLOCATER <--> sweepread_anglecheck.f

            IF(ICALLSWEEPIN.EQ.0) 

             IF(RAMFL .EQ. 'Y') 
              CALL OPENRAMFILE <--> readramfile.f   

             IF(IATTEN .EQ. 'Y') 
              CALL MAKEATTENTABLE1 <--> makeattentable1.f

             IF(ITAPEOPEN.EQ.1) 
              CALL TAPEOPEN <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f
             ELSE
              CALL OPENRADARFILE <--> radardiscread1.f

            IF(ICALLSWEEPIN.EQ.0) 

             CALL READTAPE <--> sweepread_anglecheck.f

                IF(ITAPEOPEN.EQ.0)
                 CALL RADARDISCREAD1 <--> radardiscread1.f                   
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f
                ELSE
                 CALL RADARTAPEREAD <--> tapedummy.f

                CALL ANGLECOR <--> paulsanglecor.f

                CALL CONVERTANGLES <--> convertangles.f

                IF(RAMFL .EQ. 'Y') 

                 CALL READRAMFILE <--> readramfile.f

                    IF(irec.lt.2 .or. irec.gt.maxrec .or. timein.lt.0.0)
                     CALL CTMEL <--> readramfile.f

                    IF(rlat .le. -999.)
                     CALL CTMEL <--> readramfile.f
                     CALL CTMEL <--> readramfile.f

                CALL RAYCOR1 <--> raycor1.f

                   IF (iattenflag.EQ.1)
                    CALL INTERVENATTEN1 <--> intervenatten1.f

                   IF (iunfoldflag.EQ.1)
                    CALL UNFLD1 <--> unfld1.f

                CALL RAYCOR1 <--> raycor1.f

             IF(IFLIGHTCHECK.EQ.0) 
              CALL FINDIDAYS <--> findidays.f
              CALL FINDIDAYS <--> findidays.f

             IF(TIMEMSTIME.LT.-72000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             IF(LUNR.NE.80)

              CALL TAPESEARCH <--> tapedummy.f

             ELSE

              CALL DISCSEARCH <--> radardiscread1.f

                 IF(key .le. itime(kindradar))
                  CALL CTMEL <--> readramfile.f
                  CALL CTMEL <--> readramfile.f 

                 CALL DISCENTRY <--> radardiscread1.f
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f

                 IF(ioldsweep(kindradar) .ne. numsweep(kindradar))
                  CALL CTMEL <--> readramfile.f

             CALL READTAPE <--> sweepread_anglecheck.f

             IF((TIME-STIME).LT.-60000.)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

            CALL READTAPE <--> sweepread_anglecheck.f

            IF(TIME.GT.STIME.AND.IFILECHECK.EQ.0.AND.ITAPEOPEN.EQ.1)
             CALL TAPECONTROL <--> tapedummy.f
             CALL TAPECONTROL <--> tapedummy.f

            IF(NUM_RAYS.GT.0.AND.IUSEWIND.EQ.1)
             CALL VRLOCATER <--> sweepread_anglecheck.f

            IF(NUM_RAYS.GT.0)
             CALL VRLOCATER <--> sweepread_anglecheck.f

            IF((ISWEEPSET.EQ.2.AND.SWEEPTEST.LT.0.).OR.NUM_RAYS.EQ.800)
             CALL CTME1 <--> readramfile.f

             IF(TRAKAZM.GT.91..AND.TRAKAZM.LT.269.)
              CALL FIND_GROUND <--> sweepread_anglecheck.f

            IF(TRAKAZM.GT.91..AND.TRAKAZM.LT.269.)
             CALL FIND_GROUND <--> sweepread_anglecheck.f

            IF(MOD(ITIME,60).EQ.0 .AND. itime.NE.OLDTIME)
             CALL CTME1 <--> readramfile.f

            IF(RDTAIIN(I).GT.RATEST)THEN

             IF(J.LE.NUM_RAYS.AND.AZ_EL(1,J).LT.180.)
              CALL GROUNDCHECK <--> sweepread_anglecheck.f

             IF(J.GE.1)
              IF(AZ_EL(1,J).GT.180.)
               CALL GROUNDCHECK <--> sweepread_anglecheck.f

            CALL DEFRECKLE <--> defreckle.f

            CALL UNFLDNEW <--> unfldnew.f

            CALL AZIUNFLD <--> aziunfldguess.f

        >IF(IRAMSFILES.EQ.1)

          CALL READRAMFILE <--> readramfile.f

             IF(irec.lt.2 .or. irec.gt.maxrec .or. timein.lt.0.0)
              CALL CTMEL <--> readramfile.f

             IF(rlat .le. -999.)
              CALL CTMEL <--> readramfile.f
              CALL CTMEL <--> readramfile.f

        >IF(IREADABCD.NE.1.AND.TIMES.GE.TSTART.AND.TIMES.LE.TEND)

          CALL VRLOCATE <--> abcd_shell.f

          IF(VR(I).GT.-800..OR.DBZ(I).GT.DBZCHK.AND.
             ZR(I).GT.ZRCHECK.AND.ZR(I).GT..075.
             AND.RANGE(I).GT.RLOW.AND.RANGE(I).LT.RHIGH)
           IF(IGAMMATEST.EQ.1)
            CALL DIFFOBLONG <--> miscellaneous_scan.f
               CALL ROTATE <--> miscellaneous_scan.f
            CALL DIFFDISTXYRT <--> miscellaneous_scan.f

    >>IF(IWRITEABCD.EQ.1)
       CALL WRITEABCD <--> files_rho_13_real4.f

+++++++++++++
OUTCOME: Esentially requires no flow into or out of flight track
         (tilt correction) and level surface (azimuth correction).  
         Range delay is y-intercept of best fit line to sea surface.  
         Looks for unrealistic dBZ/Vr in first few bins and sets ibst.
+++++++++++++

********************************************************************************

./anglefix

>>>PROGRAM ANGLEFIX

+++++++++++++
OUTCOME: Create new job files with angle corrections. (No roll,pitch,
         drift corrections attempted)
+++++++++++++

********************************************************************************

./anglemover

<mv script>

********************************************************************************

./wind_interpolate_auto_guess_fourier_defreckle_passthrough g

(The following are NOT used, but listed in makefile:
solo_non_convt.o tv_map_sub.o wedgint.o x_to_tiff_subs.o 
tiff_subs.o wclip.o show_iz.o azdif.o diver3.o 
miscellaneous_scan.o unfldnew.o ctme.o) 

>>>PROGRAM LEASTWIND <--> dnewbig_interpolate_auto_guess_xyrt.f

    >>CALL GETARG <--> FORTRAN intrinsic

    >>IF(IRAMSFILES.EQ.1)
       CALL OPENRAMFILE <--> readramfile.f

          IF(fileid .eq. 'ram1')
           CALL CTMEL <--> readramfile.f
           CALL CTMEL <--> readramfile.f

    >>CALL ABCD <--> abcd_scan_interpolate4_xyrt_graphics.f

        >IF(IGRAPHICS.EQ.1) 
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c

        >IF(IUSEWIND.EQ.1) 
          CALL READWIND <--> readwind_auto_guess_fourier.f
              
             CALL READHEADERS <--> readfile.f
             CALL LDPLN <--> readfile.f

        >CALL SWEEPIN <--> sweepread_auto_guess_disc_defreckle_passthrough.f
            
            IF(ICALLSWEEPIN.EQ.1)
             IF(IUSEWIND.EQ.1)
              CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f
              CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f
              CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f

            IF(ICALLSWEEPIN.EQ.0)

             IF(RAMFL .EQ. 'Y') 
              CALL OPENRAMFILE <--> readramfile.f   
       
             IF(IATTEN .EQ. 'Y') 
              CALL MAKEATTENTABLE1 <--> makeattentable1.f

             IF(ITAPEOPEN.EQ.1) 
              CALL TAPEOPEN <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f
             ELSE
              CALL OPENRADARFILE <--> radardiscread1.f

            IF(ICALLSWEEPIN.EQ.0) 

             CALL READTAPE <--> sweepread_auto_guess_disc_defreckle_passthrough.f

                IF(ITAPEOPEN.EQ.0)
                 CALL RADARDISCREAD1 <--> radardiscread1.f                   
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f
                ELSE
                 CALL RADARTAPEREAD <--> tapedummy.f

                CALL ANGLECOR <--> paulsanglecor.f

                CALL RAYCOR1 <--> raycor1.f

                   IF (iattenflag.EQ.1)
                    CALL INTERVENATTEN1 <--> intervenatten1.f

                   IF (iunfoldflag.EQ.1)
                    CALL UNFLD1 <--> unfld1.f

                CALL RAYCOR1 

                CALL CONVERTANGLES <--> convertangles.f

                IF(RAMFL .EQ. 'Y') 

                 CALL READRAMFILE <--> readramfile.f

                    IF(irec.lt.2 .or. irec.gt.maxrec .or. timein.lt.0.0)
                     CALL CTMEL <--> readramfile.f

                    IF(rlat .le. -999.)
                     CALL CTMEL <--> readramfile.f
                     CALL CTMEL <--> readramfile.f

             IF(IFLIGHTCHECK.EQ.0) 
              CALL FINDIDAYS <--> findidays.f
              CALL FINDIDAYS <--> findidays.f

             IF(TIMEMSTIME.LT.-72000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             IF(LUNR.NE.80)

              CALL TAPESEARCH <--> tapedummy.f

             ELSE

              CALL DISCSEARCH <--> radardiscread1.f

                 IF(key .le. itime(kindradar))
                  CALL CTMEL <--> readramfile.f
                  CALL CTMEL <--> readramfile.f 

                 CALL DISCENTRY <--> radardiscread1.f
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f

                 IF(ioldsweep(kindradar) .ne. numsweep(kindradar))
                  CALL CTMEL <--> readramfile.f

             CALL READTAPE <--> sweepread_auto_guess_disc_defreckle_passthrough.f

             IF((TIME-STIME).LT.-60000.)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             CALL READTAPE <--> sweepread_anglecheck.f

             IF(TIME.GT.STIME.AND.IFILECHECK.EQ.0.AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f 

             IF(NUM_RAYS.GT.0.AND.IUSEWIND.EQ.1)
              CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f

             IF(NUM_RAYS.GT.0)
              CALL VRLOCATER <--> sweepread_auto_guess_disc_defreckle_passthrough.f

             IF((ISWEEPSET.EQ.2.AND.SWEEPTEST.LT.0.).OR.NUM_RAYS.EQ.800)
              CALL CTME1 <--> readramfile.f

              IF(TRAKAZM.GT.80..AND.TRAKAZM.LT.280.)
               CALL FIND_GROUND <--> sweepread_auto_guess_disc_defreckle_passthrough.f

             IF(TRAKAZM.GT.80..AND.TRAKAZM.LT.280.)
              CALL FIND_GROUND <--> sweepread_auto_guess_disc_defreckle_passthrough.f

             IF(MOD(ITIME,60).EQ.0 .AND. itime.NE.OLDTIME)
              CALL CTME1 <--> readramfile.f

             IF(RDTAIIN(I).GT.RATEST)THEN

              IF(J.LE.NUM_RAYS.AND.AZ_EL(1,J).LT.180.)
               CALL GROUNDCHECK <--> sweepread_auto_guess_disc_defreckle_passthrough.f

              IF(J.GE.1)
               IF(AZ_EL(1,J).GT.180.)
                CALL GROUNDCHECK <--> sweepread_auto_guess_disc_defreckle_passthrough.f

            CALL DEFRECKLE <--> defreckle.f

            CALL VARUNFLD <--> varunfldguess.f
               CALL SPARSEVAR <--> varunfldguess.f
                  CALL SBINI <--> ditp.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I+1,J).EQ.1.AND.
                     ABS(DVRADR(I,J)).LT..05)
                   CALL MATRIXSUM <--> varunfldguess.f
                      CALL SORTER <--> varunfldguess.f
                      IF(INDX(I).GE.INDX(L))
                       IF(INDX(I).GT.0.AND.INDX(L).GT.0)
                        IF(IN.GT.0.AND.INN.GT.0)
                         CALL SBSIJ <--> ditp.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I+1,J).EQ.1.AND.
                     INDEX4(I-1,J).EQ.1.AND.DR(I-1).EQ.DR(I))
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I,J+1).EQ.1.AND.
                     ABS(DVRADAZM(I,J)).LT.600.)
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I,J+1).EQ.1.AND.
                     INDEX4(I,J-1).EQ.1)
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1)THEN
                   CALL MATRIXSUM <--> varunfldguess.f
                  CALL SBEND <--> ditp.f
               CALL SOLVESPARSE <--> solvesparse_sweep.f
                  CALL DFAULT <--> ditp.f
                  CALL SSORCG <--> ditp.f

    >>IF(IWRITEABCD.EQ.1)
       CALL WRITEABCD <--> files_rho_13_real4.f

cp SUM*helped* SUMFILES/

cp jobfile_radar jobfile_radar_full

+++++++++++++
OUTCOME: QC done again. Bargen-Brown dealiasing again, but using
         Fourier wv#0,1 background to improve dealiasing.
+++++++++++++

********************************************************************************

./wind3_fill_auto_guess

(The following are NOT used, but listed in makefile:
 readramfile.o ctme.o)

>>>PROGRAM LEASTWIND <--> dnewbig_auto_guess.f

    >>IF(WTBACKU.GT.0..OR.WTBACKV.GT.0..OR.WTBACKW.GT.0.)
       IF(ISETZERO.NE.0)
        CALL READHEADER <--> files_rho_13_real.f
        CALL LDPLNW <--> files_rho_13_real.f
        CALL DIVER3 <--> diver3.f

    >>CALL DIVER <--> miscellaneous_scan_wind.f

    >>IF(IMULTI.EQ.0)
       CALL WCHECK <--> dnewbig_auto_guess.f
         >CALL READABCD <--> files_rho_13_real.f
      ELSE
       CALL WZERO <--> dnewbig_auto_guess.f
         >CALL READABCD <--> files_rho_13_real.f

    >>IF(IDBZONLY.NE.1)

       CALL SPARSE <--> sparsenew_fill.f

         >CALL SBINI <--> ditp.f
         >CALL BOUND5 <--> sparsenew_fill.f
         >CALL VALUESET <--> sparsenew_fill.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
             CALL SORTER <--> sparsenew_fill.f
             IF(INDX(I).GE.INDX(L))
              IF(INDX(I).GT.0.AND.INDX(L).GT.0)
               IF(IN.GT.0.AND.INN.GT.0)
                CALL SBSIJ <--> ditp.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(IABS(IBOUNDX).LE.1.AND.IABS(IBOUNDY).LE.1)
           CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(IABS(IBOUNDX).LE.1.AND.IABS(IBOUNDZ).LE.1)
           CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(IABS(IBOUNDY).LE.1.AND.IABS(IBOUNDZ).LE.1)
           CALL MATRIXSUM <--> sparsenew_fill.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(IBOUNDZ.EQ.-2000)
           CALL MATRIXSUM <--> sparsenew_fill.f
         >CALL MATRIXSUM <--> sparsenew_fill.f
         >CALL READABCD <--> files_rho_13_real.f
         >IF(U(I,J,K).GT.FLAG)
           IF(UFIT(I,J,K).GT.FLAG)
            CALL MATRIXSUM <--> sparsenew_fill.f
            CALL MATRIXSUM <--> sparsenew_fill.f
            CALL MATRIXSUM <--> sparsenew_fill.f
           ELSE
            CALL MATRIXSUM <--> sparsenew_fill.f
           IF(K.EQ.1.OR.(FACTZZ.LE.0..AND.FACTWP.LE.0.))
            IF(FACTZ.GT.0.)
             CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(K.EQ.KMAX.OR.K.EQ.KTTOP(I,J))
           IF(FACTZ.GT.0.)
            CALL MATRIXSUM <--> sparsenew_fill.f
         >IF(INPOSITIONU.GT.0.AND.INPOSITIONV.GT.0.AND.INPOSITIONW.GT.0)
           IF(BAND1(INPOSITIONU).GT.0..AND.BAND1(INPOSITIONV).GT.0.
              .AND.BAND1(INPOSITIONW).GT.0.)THEN
            CALL SBSIJ <--> ditp.f
            CALL SBSIJ <--> ditp.f
            CALL SBSIJ <--> ditp.f
            CALL SBSIJ <--> ditp.f
            CALL SBSIJ <--> ditp.f
            CALL SBSIJ <--> ditp.f
         >CALL SBEND <--> ditp.f

       CALL WRITEHEADER <--> files_rho_13_real.f

       CALL WRPLN <--> files_rho_13_real.f

       CALL SOLVESPARSE <--> solvesparse_itpack_auto.f
         >CALL DFAULT <--> ditp.f
         >CALL SSORCG <--> ditp.f

       CALL DIVER3 <--> diver3.f

       CALL WRITEHEADER <--> files_rho_13_real.f

       CALL DEFRECKLEXYZ <--> defrecklexyz.f

       CALL WRPLNEU <--> files_rho_13_real.f
         >IF(K.EQ.2.OR.K.EQ.3)
           CALL GETTIMES <--> files_rho_13_real.f
           CALL WRITETIME <--> files_rho_13_real.f
           CALL WRITETIME <--> files_rho_13_real.f
           CALL WRITETIME <--> files_rho_13_real.f

    >>IF(IDBZONLY.EQ.1)
       CALL WRITEHEADER <--> files_rho_13_real.f
       CALL MINMAX <--> solvesparse_itpack_auto.f
       CALL WRPLN <--> files_rho_13_real.f

+++++++++++++
OUTCOME: Synthesis with wv#0,1 as background.
+++++++++++++



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
MOST OF CODE BELOW IS USED FOR VERTICAL PROFILES
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$



********************************************************************************

./radazi_in

>>>PROGRAM RADAZI <--> radazi_in.f

+++++++++++++
OUTCOME: Determines azimuth of aircraft relative to storm center (in)
+++++++++++++

********************************************************************************

./radazi_out

>>>PROGRAM RADAZI <--> radazi_out.f

+++++++++++++
OUTCOME: Determines azimuth of aircraft relative to storm center (out)
+++++++++++++

********************************************************************************

cp aziradin aziradinout
cp inoutradazi.txt inoutradazi_save.txt
cp jobfile_radar_profile_in jobfile_radar_profile
cp jobfile_radar_in jobfile_radar

********************************************************************************

./wind_interpolate_rt_auto_profile_swath g

(in)

(The following are NOT used, but listed in makefile:
 tv_map_sub.o x_to_tiff_subs.o tiff_subs.o diver3.o 
 miscellaneous_scan.o unfld1.o unfldnew.o intervenatten1.o ctme.o) 

>>>PROGRAM LEASTWIND <--> dnewbigr_abcdonly_auto_guess_profile_swath.f

    >>CALL GETARG <--> FORTRAN intrinsic 

    >>IF(IGUESS.EQ.5)
       CALL READHEADER <--> files_rho4.f
       CALL LDPLNWD <--> files_rho4.f
       CALL WLOCATE <--> dnewbigr_abcdonly_auto_guess_profile_swath.f

    >>IF(IRAMSFILES.EQ.1)
       CALL OPENRAMFILE <--> readramfile.f

         >IF(fileid .eq. 'ram1')
           CALL CTMEL <--> readramfile.f
           CALL CTMEL <--> readramfile.f

    >>CALL ABCD <--> abcd_scan_rt4_512_swath.f

        >IF(IGRAPHICS.EQ.1) 
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c
          CALL RASX_OPEN_WINDOW <--> bigger_raster_x_functions.c

        >IF(IUSEWIND.EQ.1) 
          CALL READWIND <--> readwind_auto_guess_fourier.f
              
             CALL READHEADERS <--> readfile.f
             CALL LDPLN <--> readfile.f

        >CALL SWEEPIN <--> sweepread_vrrms.f

            IF(ICALLSWEEPIN.EQ.1)
             IF(IUSEWIND.EQ.1)
              CALL VRLOCATER <--> sweepread_vrrms.f             
              CALL VRLOCATER <--> sweepread_vrrms.f 
              CALL VRLOCATER <--> sweepread_vrrms.f

            IF(ICALLSWEEPIN.EQ.0) 

             IF(RAMFL .EQ. 'Y') 
              CALL OPENRAMFILE <--> readramfile.f   
       
             IF(IATTEN .EQ. 'Y') 
              CALL MAKEATTENTABLE1 <--> makeattentable1.f

             IF(ITAPEOPEN.EQ.1) 
              CALL TAPEOPEN <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f
             ELSE
              CALL OPENRADARFILE <--> radardiscread1.f

            IF(ICALLSWEEPIN.EQ.0) 

             CALL READTAPE <--> sweepread_vrrms.f 

                IF(ITAPEOPEN.EQ.0)
                 CALL RADARDISCREAD1 <--> radardiscread1.f                   
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f
                ELSE
                 CALL RADARTAPEREAD <--> tapedummy.f

                CALL ANGLECOR <--> paulsanglecor.f

                CALL CONVERTANGLES <--> convertangles.f

                IF(RAMFL .EQ. 'Y') 

                 CALL READRAMFILE <--> readramfile.f

                    IF(irec.lt.2 .or. irec.gt.maxrec .or. timein.lt.0.0)
                     CALL CTMEL <--> readramfile.f

                    IF(rlat .le. -999.)
                     CALL CTMEL <--> readramfile.f
                     CALL CTMEL <--> readramfile.f

                CALL RAYCOR1 <--> raycor1.f
                CALL RAYCOR1 <--> raycor1.f

             IF(IFLIGHTCHECK.EQ.0) 
              CALL FINDIDAYS <--> findidays.f
              CALL FINDIDAYS <--> findidays.f

             IF(TIMEMSTIME.LT.-72000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

             IF(LUNR.NE.80)

              CALL TAPESEARCH <--> tapedummy.f

             ELSE

              CALL DISCSEARCH <--> radardiscread1.f

                 IF(key .le. itime(kindradar))
                  CALL CTMEL <--> readramfile.f
                  CALL CTMEL <--> readramfile.f 

                 CALL DISCENTRY <--> radardiscread1.f
                    CALL DISKREAD1 <--> diskread1.f
                    CALL SAVESWEEP_DISC <--> radardiscread1.f
                    CALL CHKTM1 <--> radardiscread1.f
                    CALL DCOMPRESS1 <--> dcompress1.f
                       IF(CODE .NE. 0)
                        CALL UNPK81 <--> unpk81.f

                 IF(ioldsweep(kindradar) .ne. numsweep(kindradar))
                  CALL CTMEL <--> readramfile.f

             CALL READTAPE <--> sweepread_vrrms.f 

             IF((TIME-STIME).LT.-60000..AND.ITAPEOPEN.EQ.1)
              CALL TAPECONTROL <--> tapedummy.f
              CALL TAPECONTROL <--> tapedummy.f

            CALL READTAPE <--> sweepread_vrrms.f 

            IF(TIME.GT.STIME.AND.IFILECHECK.EQ.0.AND.ITAPEOPEN.EQ.1)
             CALL TAPECONTROL <--> tapedummy.f
             CALL TAPECONTROL <--> tapedummy.f

             IF(NUM_RAYS.GT.0.AND.IUSEWIND.EQ.1)            
              CALL VRLOCATER <--> sweepread_vrrms.f 

             IF(NUM_RAYS.GT.0.AND.IUSEWIND.EQ.1)
              CALL VRLOCATER <--> sweepread_vrrms.f 

             IF(NUM_RAYS.GT.0)
              CALL VRLOCATER <--> sweepread_vrrms.f

             IF((ISWEEPSET.EQ.2.AND.SWEEPTEST.LT.0.).OR.NUM_RAYS.EQ.800)
              CALL CTME1 <--> readramfile.f

              IF(TRAKAZM.GT.90..AND.TRAKAZM.LT.271.)
               CALL FIND_GROUND <--> sweepread_vrrms.f

             IF(TRAKAZM.GT.90..AND.TRAKAZM.LT.271.)
              CALL FIND_GROUND <--> sweepread_vrrms.f

             IF(MOD(ITIME,60).EQ.0 .AND. itime.NE.OLDTIME)
              CALL CTME1 <--> readramfile.f

             IF(RDTAIIN(I).GT.RATEST)THEN

              IF(J.LE.NUM_RAYS.AND.AZ_EL(1,J).LT.180.)
               CALL GROUNDCHECK <--> sweepread_vrrms.f

              IF(J.GE.1)
               IF(AZ_EL(1,J).GT.180.)
                CALL GROUNDCHECK <--> sweepread_vrrms.f

            CALL DEFRECKLE <--> defreckle.f

            CALL VARUNFLD <--> varunfldguess.f
               CALL SPARSEVAR <--> varunfldguess.f
                  CALL SBINI <--> ditp.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I+1,J).EQ.1.AND.
                     ABS(DVRADR(I,J)).LT..05)
                   CALL MATRIXSUM <--> varunfldguess.f
                      CALL SORTER <--> varunfldguess.f
                      IF(INDX(I).GE.INDX(L))
                       IF(INDX(I).GT.0.AND.INDX(L).GT.0)
                        IF(IN.GT.0.AND.INN.GT.0)
                         CALL SBSIJ <--> ditp.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I+1,J).EQ.1.AND.
                     INDEX4(I-1,J).EQ.1.AND.DR(I-1).EQ.DR(I))
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I,J+1).EQ.1.AND.
                     ABS(DVRADAZM(I,J)).LT.600.)
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1.AND.INDEX4(I,J+1).EQ.1.AND.
                     INDEX4(I,J-1).EQ.1)
                   CALL MATRIXSUM <--> varunfldguess.f
                  IF(INDEX4(I,J).EQ.1)THEN
                   CALL MATRIXSUM <--> varunfldguess.f
                  CALL SBEND <--> ditp.f
               CALL SOLVESPARSE <--> solvesparse_sweep.f
                  CALL DFAULT <--> ditp.f
                  CALL SSORCG <--> ditp.f

        >IF(IGRAPHICS.EQ.1)

          CALL SOLO_CONVT <--> solo_non_convt.f
             CALL GET_RADII <--> solo_non_convt.f 
             FUNCTION AZDIF <--> azdif.f
             CALL WCLIP <--> wclip.f
             CALL WCLIP <--> wclip.f
             CALL WEDGINT <--> wedgint.f
             CALL WCLIP <--> wclip.f
          CALL SHOW_IZ <--> show_iz.f
             CALL RASX_DISPLAY_IMAGE <--> bigger_raster_x_functions.c
             CALL X_TVLINE <--> bigger_raster_x_functions.c
          CALL SOLO_CONVT <--> solo_non_convt.f
          CALL SHOW_IZ <--> show_iz.f

        >IF(IRAMSFILES.EQ.1)
          CALL READRAMFILE <--> readramfile.f

        >IF(IREADABCD.NE.1.AND.TIMES.GE.TSTART.AND.TIMES.LE.TEND)
          CALL VRLOCATE <--> abcd_scan_rt4_512_swath.f

    >>IF(IWRITEABCD.EQ.1)
       CALL WRITEABCD <--> files_rho4.f
       CALL WRITEABCD <--> files_rho4.f

cp SUM*profile* SUMFILES/

+++++++++++++
OUTCOME: ??
+++++++++++++

********************************************************************************

./wind_rt_auto_guess_profile

(in)

(The following are NOT used, but listed in makefile:
 miscellaneous_scan_rt.o)

>>>PROGRAM LEASTWIND <--> dnewbigr_auto_guess_profile.f

    >>CALL GETARG <--> FORTRAN intrinsic

    >>IF(IGUESS.EQ.1)
       CALL READHEADER <--> files_rho.f
       CALL LDPLNW <--> files_rho.f
       CALL DIVER3 <--> diver3_rt.f

    >>CALL WZERO <--> dnewbigr_auto_guess_profile.f
        >CALL READABCD <--> files_rho.f
    >>CALL WRITEHEADER <--> files_rho.f
    >>CALL WRPLN <--> files_rho.f

    >>IF(IDBZONLY.NE.1)

       CALL SPARSE <--> sparsenew_rt.f
         >CALL SBINI <--> itpack.f
         >IF(IBOUNDZ.EQ.-2)
           CALL MATRIXSUM <--> sparsenew_rt.f
              CALL SORTER <--> sparsenew_rt.f
              IF((INDX(I).GE.INDX(L))
               IF(INDX(I).GT.0.AND.INDX(L).GT.0)
                IF(IN.GT.0.AND.INN.GT.0)
                 CALL SBSIJ <--> itpack.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >IF(IABS(IBOUNDX).LE.1.AND.IABS(IBOUNDY).LE.1)
           CALL MATRIXSUM <--> sparsenew_rt.f
         >IF(IABS(IBOUNDX).LE.1.AND.IABS(IBOUNDZ).LE.1)
           CALL MATRIXSUM <--> sparsenew_rt.f
         >IF(IABS(IBOUNDY).LE.1.AND.IABS(IBOUNDZ).LE.1)
           CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f
         >IF(IBOUNDZ.EQ.-2000)
           CALL MATRIXSUM <--> sparsenew_rt.f
         >CALL MATRIXSUM <--> sparsenew_rt.f        
         >CALL READABCD <--> files_rho.f
         >IF(K.EQ.1.AND.U(I,J,K).GT.FLAG.AND.
            (FACTWTOP.LE.0.01.OR.INDEX4(I,J,K).EQ.1))
           CALL MATRIXSUM <--> sparsenew_rt.f  
         >IF(BAND1(INPOSITIONU).GT.0..AND.BAND1(INPOSITIONV).GT.0.
            .AND.BAND1(INPOSITIONW).GT.0.)
           IF(IWZERO.EQ.1)
            CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
           CALL SBSIJ <--> itpack.f
         >CALL SBEND <--> itpack.f

       CALL WRITEHEADER <--> files_rho.f
       CALL WRPLN <--> files_rho.f
       CALL SOLVESPARSE <--> solvesparse_itpack.f
         >CALL DFAULT <--> itpack.f
         >CALL SSORCG <--> itpack.f
       CALL SETINDEX4 <--> dnewbigr_auto_guess_profile.f
       CALL DIVER3 <--> diver3_rt.f

    >>CALL FLAGBOTTOM <--> dnewbigr_auto_guess_profile.f
    >>CALL WRITEHEADER <--> files_rho.f
    >>CALL WRPLN <--> files_rho.f
    >>CALL WRITEHEADER <--> files_rho.f
    >>CALL DEFRECKLEXZ <--> defrecklexz.f
    >>CALL WCHANGE <--> dnewbigr_auto_guess_profile.f
        >CALL FALLDIF <--> dnewbigr_auto_guess_profile.f
    >>CALL INOROUT <--> dnewbigr_auto_guess_profile.f
        >CALL GETTIMES <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
        >CALL WRITETIME <--> dnewbigr_auto_guess_profile.f
    >>CALL WRPLN <--> files_rho.f

+++++++++++++
OUTCOME: ??
+++++++++++++

********************************************************************************

cp profile_full.w profile_full_in_save.w
cp jobfile_radar_profile_out jobfile_radar_profile
cp jobfile_radar_out jobfile_radar
cp aziradout aziradinout

********************************************************************************

./wind_interpolate_rt_auto_profile_swath g

(out)

See above (in) for code description.

cp SUM*profile* SUMFILES/

+++++++++++++
OUTCOME: ??
+++++++++++++

********************************************************************************

./wind_rt_auto_guess_profile

(out)

See above (in) for code description.

+++++++++++++
OUTCOME: ??
+++++++++++++

********************************************************************************

cp jobfile_radar_full jobfile_radar
cp *.w wfiles/
