C+ C NAME: C MkD2V C PURPOSE: C Make a velocity map from a density map using the constant C momentum flux assumption C CATEGORY: C Data processing C CALLING SEQUENCE: subroutine MkD2V(N,DD,VV,V0) C INPUTS: C N integer # data points C DD(N) real nr^2 C V0 real Standard velocity C OUTPUTS: C VV(N) real Velocity C CALLS: C iArrR4ValuePresent, ArrR4Total, BadR4 C PROCEDURE: C Bad values (indicated by BadR4()) are not processed C MODIFICATION HISTORY: C MAY-1999, Paul Hick (UCSD/CASS; pphick@ucsd.edu) C- integer N real DD(N) ! nr^2 map real VV(N) ! Velocity map real V0 Bad = BadR4() F0 = ArrR4Total(-N,DD,I) F0 = F0/I ! Average normalized density F0 = F0*V0*V0 ! Average momentum flux do I=1,N ! Make velocity map if (DD(I) .eq. Bad) then VV(I) = Bad else VV(I) = sqrt(F0/DD(I)) end if end do return end