b'Magnetics in the mountainsFeatureCOSI = COS(AI)SIND = SIN(AD)COSD = COS(AD)COSA = COS(AZIMR)SINA = SIN(AZIMR)! initialize operation factors!IFLAG = 1CFAC = CMPLX(1.0_8,0.0_8,8)!! process each ROWDO IROW = 1,NYFMODE=1CALL GRDROW(LuHFREQ,MODE,NXF,IROW,ITYPEF,C1)IF(MOD(IROW,50) .EQ. 1) THENCALL GSTATuS(IROW,NYF)ENDIF! set operation factors for this ROWIX = IROW - 1IF (IROW NYQX) IX = IX - NXPWX = IX * DELWXWXSQ = WX * WXWXCD = WX * COSDXLANC = 1.0XARG = IX * XLFACIF (XARG .NE. 0.0) XLANC = SIN(XARG) / XARGXHANN = 0.5 * (1.0 + COS(WX))!! process each COLuMNDO ICOL = 1,NYP!! branch if first COLuMN OF FIRST ROW - DC componentIF (IFLAG.EQ.1) THENIFLAG = 2ELSE!! set operation factors for each COLuMNIY = ICOL - 1IF (ICOL .GT. NYQY) IY = IY - NYPWY = IY * DELWYWYSQ = WY * WYWNuM = SQRT(WXSQ + WYSQ)CFAC = CMPLX(1.0_8,0.0_8,8) !! LanczosIF (LOP(1)) THENYLANC = 1.0YARG = IY * YLFACIF (YARG .NE. 0.0) YLANC = SIN(YARG) / YARGRLANC = XLANC * YLANCCFAC = CFAC * RLANCEND IF!! HanningIF (LOP(2)) THENYHANN = 0.5 * (1.0 + COS(WY))RHANN = XHANN * YHANNCFAC = CFAC * RHANNEND IF!! PSEuDO-GRAVITY TRANSFORMIF (LOP(10)) THENWYSD = WY * SIND! NB COMMENTED OuT LINES DO SAME AS ACTIVE LINES! COMMENTED LINES ARE IN uSGS FORMAT! RPOLER = WNuM*SINI**2-((COSI * (WXCD + WYSD))**2)/WNuM! RPOLEI = 2*SINI*COSI * (WXCD + WYSD)! RDENOM = 1.0/(RPOLER**2+RPOLEI**2) 43 PREVIEW JuNE 2020'