point.f
51 lines
| 1.6 KiB
| text/x-fortran
|
FortranFixedLexer
| r0 | C $Id: point.f 3304 2011-01-17 15:25:59Z brideout $ | |||
| C | ||||
| SUBROUTINE POINT(SR,SLAT,SLON,AZ,EL,RANGE,PR,GLAT,GLON) | ||||
| C | ||||
| C jmh - 1/80 ans fortran 66 | ||||
| C | ||||
| C POINT calculates the position of a point defined by the radar | ||||
| C line-of sight vector to that point. | ||||
| C | ||||
| C Input: | ||||
| C SR - distance of station from center of earth (km) | ||||
| C SLAT - geocentric latitude of station (deg) | ||||
| C SLON - longitude of station (deg) | ||||
| C AZ - radar azimuth (deg) | ||||
| C EL - radar elevation (deg) | ||||
| C RANGE - radar range (km) | ||||
| C | ||||
| C Output: | ||||
| C PR - distance from center of earth of observation point (km) | ||||
| C GLAT - observation point geocentric latitude (deg) | ||||
| C GLON - observation point longitude (deg) | ||||
| C | ||||
| C ...calculate "line-of-sight" station centered cartesian coords... | ||||
| C .. Scalar Arguments .. | ||||
| DOUBLE PRECISION AZ,EL,GLAT,GLON,PR,RANGE,SLAT,SLON,SR | ||||
| C .. | ||||
| C .. Local Scalars .. | ||||
| DOUBLE PRECISION RP,RR,RT,T | ||||
| C .. | ||||
| C .. Local Arrays .. | ||||
| DOUBLE PRECISION P(3),R(3),S(3) | ||||
| C .. | ||||
| C .. External Subroutines .. | ||||
| EXTERNAL CSCONV,VADD,VCTCNV | ||||
| C .. | ||||
| CALL CSCONV(RT,RP,RR,RANGE,90.0D0-EL,180.0D0-AZ,2) | ||||
| C | ||||
| C ...calculate "line-of-sight" earth centered cartesian coords | ||||
| C and "station" earth centered cartesian coords... | ||||
| CALL VCTCNV(R(1),R(2),R(3),S(1),S(2),S(3),RR,RT,RP,SR,90.0D0-SLAT, | ||||
| * SLON,2) | ||||
| C | ||||
| C ...calculate "observation-point" earth centered cartesian coords.. | ||||
| CALL VADD(S,R,P) | ||||
| C | ||||
| C ...calculate "observation-point" earth centered spherical coords.. | ||||
| CALL CSCONV(P(1),P(2),P(3),PR,T,GLON,1) | ||||
| GLAT = 90.0D0 - T | ||||
| RETURN | ||||
| C | ||||
| END | ||||
