C+
C NAME:
C	MkVObstd
C PURPOSE:
C	Make pseudo observed point P velocities from a velocity map.  This program has not
C	been used so far in the time-dependent mode to check structures and their reconstruction
C	and it is not called in the main program
.
C CATEGORY:
C	Data processing
C CALLING SEQUENCE:
C	call MkVObstd(XCbeg,XCend,VMAP,nLng,nLat,XC,XL,NL,NLV,NLOSP1,
C		XLS,XCE,XP,XD,XE,IYRS,DOYS,VOBS,XLSV,XCEV,XPV,XDV,XEV,IYRSV,DOYSV)  
C INPUTS:
C       XCbeg           	real    	Beginning Carrington Rotation value
C       XCend           	real    	Ending Carrington Rotation value
C	VMAP(nLng,nLat)		real		Velocity map
C	nLng			integer		Number of longitudes
C	nLat			integer		Number of latitudes
C	XC(NL,NLOSP1)		integer		Projected Carr. rot. value of point on L.O.S.
C	XL(NL,NLOSP1)		integer		Heliographic lat. (in deg.) point on L.O.S.
C	NL			integer		Number of G-level data points
C	NLOSP1			integer		Number of L.O.S. distance segments + 1
C	XLS(NL)			real		Ecliptic longitude of Sun
C	XCE(NL)			real		Carrington variable of Earth
C	XP(NL)			real		Longitude diff. (Earth-Sun)
C	XD(NL)			real		Declination of Source
C	XE(NL)			real		Elongation of source
C	IYRS(NL)		integer		Year of source observation
C	DOYS(NL)		real		Day of year (and fraction) of source obs.
C OUTPUTS:
C	NLV			integer		Number of velocity data points
C	VOBS(NL)		real		Pseudo observed velocities for a given source
C	XLSV(NL)		real		Ecliptic longitude of Sun of p.v. source obs.
C	XCEV(NL)		real		Carrington variable of Earth of p.v. source obs.
C	XPV(NL)			real		Longitude diff. (Earth-Sun) of p.v. source obs. 
C	XDV(NL)			real		Declination of p.v. source
C	XEV(NL)			real		Elongation of p.v. source
C	IYRSV(NL)		integer		Year of p.v. source observation
C	DOYSV(NL)		real		Day of year (and fraction) of p.v. source obs.
C FUNCTIONS/SUBROUTINES:
C	Get2DVal
C PROCEDURE:
C MODIFICATION HISTORY:
C	NOV, 1995 B. Jackson (STEL,UCSD)
C-
	subroutine MkVObstd(XCbeg,Xcend,VMAP,nLng,nLat,XC,XL,NL,NLV,NLOSP1,
     &		XLS,XCE,XP,XD,XE,IYRS,DOYS,VOBS,XLSV,XCEV,XPV,XDV,XEV,IYRSV,DOYSV)
	real		VMAP(nLng,nLat)		! Velocity map
	integer		IYRS(NL),IYRSV(NL)	! Year of source observation
	real		DOYS(NL),DOYSV(NL),	! Day of year (and fraction) of source obs.
     &	        	XLS (NL),XLSV (NL),	! Ecliptic longitude of Sun
     &			XCE (NL),XCEV (NL),	! Carrington variable of Earth
     &			XP  (NL),XPV  (NL),	! Longitude diff. (Earth-Sun)
     &			XD  (NL),XDV  (NL),	! Declination of source
     &			XE  (NL),XEV  (NL),	! Elongation of source
     &			VOBS(NL)		! Pseudo observed velocities for a given source
	real		XC  (NLOSP1,NL),	! Projected Carr. rot. value of pt on L.O.S.  
     &			XL  (NLOSP1,NL)		! Heliographic lat. (in deg.) point on L.O.S.
	do I=1,NL
	    call Get2Dval(XCbeg,XCend,nLng,nLat,VMAP,XC(NLOSP1,I),XL(NLOSP1,I),VOBS(I))
	    IYRSV(I) = IYRS(I)
	    DOYSV(I) = DOYS(I)
	    XLSV(I)  = XLS(I)
	    XCEV(I)  = XCE(I)
	    XPV (I)  = XP (I)
	    XDV (I)  = XD (I)
	    XEV (I)  = XE (I)
	end do
	NLV = NL	
	return
	end