Modelo 4/3 de deformacion atmosferica

PRO model43, hrad, angles, recorregut_feix, recorregut_Terra, alcades, TERRA=terra

; *********************************************************************
; DESCRIPTION: Programa per calcularalçades del feix del radar i
; distancies recorregudes amb un model de 4/3 de Terra
; equivalent. Mirar Doviak i Zrnic per formules
;*********************************************************************
; AUTHOR: XAVI LLORT, DAVID FORCADELL
; E-MAIL incidencies (24h): david(AT)grahi.upc.edu
; DATE: FEBRER 2005
; *********************************************************************; *********************************************************************
; VARIABLES
; hrad: Altura radar (km) - ENTRADA
; angles: Angles de les elevacions (en radians) - ENTRADA
;recorregut_feix: Distancia recorreguda en el feix de radar (en km) - ENTRADA (no necessaria si /TERRA)
; recorregut_Terra: Recorregut sobre la superfície de la Terra (en km) - ENTRADA/SORTIDA
; alcades:Alçades dels punts (en km) - SORTIDA

; *********************************************************************
; KEYWORDS
; TERRA indica que li passem recorregut a terra per calcular alcades.
; Si noindiquem res, calculem alcades i recorregut Terra a partir de Recorregut Feix.

; *********************************************************************
; CONSTANTS

r = 6371.1 ; Radi de laTerra (en km)
r_equi = 4.*r/3. ; Radi equivalent

; *********************************************************************
; CALCUL

num_elev = N_ELEMENTS(angles)
IF (num_elev EQ 1) THEN angles =[angles]
cos_angles = COS (angles) ; Per a optimitzacio
sin_angles = SIN (angles) ; Per a optimitzacio

IF NOT(KEYWORD_SET(terra)) THEN BEGIN ; Cas de les distancies en el feix

;************* Reform de les dades a un vector per angle *****************
dims = SIZE(recorregut_feix)
recorregut_feix = REFORM( recorregut_feix, FIX(N_ELEMENTS(recorregut_feix)/num_elev), num_elev)...