b'Magnetics in the mountainsFeature! RPOLEI = COSI * (WXCD + WYSD) / WNuM**2RPOLER = SINI*WNuMRPOLEI = COSI * (WXCD + WYSD)CPOLE=CMPLX(RPOLER,RPOLEI)! CPOLE=CMPLX(RPOLER,RPOLEI)/RDENOMCPOLE = 1.0 / (CPOLE * CPOLE)CFAC = CFAC * CPOLE*CMTOGEND IF!! PSEuDO-MAGNETICS TRANSFORMIF (LOP(11)) THENWYSD = WY * SIND! NB COMMENTED OuT LINES DO SAME AS ACTIVE LINES! COMMENTED LINES ARE IN uSGS FORMAT! RPOLER = SINI*WNuM! RPOLEI = COSI * (WXCD + WYSD)RPOLER = WNuM*SINI**2-((COSI * (WXCD + WYSD))**2)/WNuMRPOLEI = 2*SINI*COSI * (WXCD + WYSD)CPOLE=CMPLX(RPOLER,RPOLEI)!CPOLE = 1.0 / (CPOLE * CPOLE)CFAC = CFAC * CPOLE/CMTOGEND IF!! multiply coefficients by operation factorEND IFC1(ICOL) = C1(ICOL)* CFACEND DO!! write updated ROW C1 to output grid! YOuR CODEEND DODEALLOCATE(C1)RETuRNEND SuBROuTINE INVF01JuNE 2020 PREVIEW 44'