C+ C NAME: C Mk_V2D C PURPOSE: C Make a density map from the velocity map using the constant momentum C flux assumption C CATEGORY: C Data processing C CALLING SEQUENCE: subroutine Mk_V2D(NN,VV,DD,D1AU,RR) C INPUTS: C VV(NN) real Velocity map C nLng integer # longitudes C nLat integer # latitudes C D1AU real density at 1 AU C RR real Velocity map reference height C OUTPUTS: C DD(NN) real Density map C CALLS: C BadR4 C PROCEDURE: C Bad values (indicated by BadR4()) are not processed C MODIFICATION HISTORY: C NOV, 1995 B. Jackson (STEL,UCSD) C- integer NN real VV(*) ! Velocity real DD(*) ! Density real D1AU real RR real V0 /400./ ! Standard velocity Bad = BadR4() D0 = V0/RR D0 = D1AU*D0*D0 do I=1,NN if (VV(I) .eq. Bad) then DD(I) = Bad else DD(I) = D0/(VV(I)*VV(I)) end if end do return end