b'Magnetics in the mountainsFeature! THE NEGATIVE CONTRASTS. JAM KARATCOMPLEX :: CPOLECOMPLEX(SELECTED_REAL_KIND(15)) :: CFACCOMPLEX(SELECTED_REAL_KIND(15)), ALLOCATABLE, DIMENSION(:) :: C1INTEGER :: NXP, NYPINTEGER :: NYQX,NYQYINTEGER :: IFLAG,ICOL,IY,IROW,IX! VALuES FOR LOP(I)! 1 = LANCOZ! 2 = HANNING! 3 = REDuCTION TO THE POLE! 4 = WITH LOW LATITuDE FILTER! 5 = DIRECTIONAL FILTER! 6 = uPWARD/DOWNWARD CONTINuATION! 7 = VERTICAL DERIVATIVE/INTEGRAL! 8 = BAND PASS FILTER! 9 = NOTCH FILTER! 10 = PSEuDO-GRAVITY TRANSFORM! 11 = PSEuDO-MAGNETIC TRANSFORM! 12 = Horizontal X derivative! 13 = Horizontal Y derivative! 14 = Horizontal slope! 15 = Hilbert Transform! 16 = REDuCTION TO THE EQuATOR! set up FT parameters! OPEN THE GRID CONTAINING THE FOuRIER COEFFICIENTS! YOuR CODENXP = NYF ! NuMBER OF ROWS IN GRIDGRID ARRAY STORED BY COLuMNS THEN ROWSNYP = NXF/2! NuMBER OF COMPLEX ELEMENTS IN A ROW OF THE GRID!ALLOCATE MEMORYALLOCATE(C1(NYP),STAT=IASTAT)IF (IASTAT.NE.0) THENIDuM=(2*NYP)*4WRITE (OuTMESS,(PROBLEM ALLOCATING ,I15, BYTES OF MEMORY TO COMPLEX &&ARRAY))IDuMCALL WERRORBOX(OuTMESS,OuTSuB)END IF!! CHECK THAT IT IS A FREQuENCY DOMAIN GRID - 2nd VALuE IN FIRST COLuMN IS ALWAYS! ZERO !!!!!!! READ A ROW FROM THE GRID INTO C1! YOuR CODEIF (ABS(AIMAG(C1(1))) EPSILON(AIMAG(C1(1)))) THENWRITE(OuTMESS,(A,G15.7))TRIM(FILEFREQ)// DOES NOT LOOK LIKE A FREQuENCY DOMAIN && GRID. CHECK YOuR INPuTS DC IMAGINARY IS ,AIMAG(C1(1))CALL WERRORBOX(OuTMESS,OuTSuB)END IFNYQX = NXP/2 + 1NYQY = NYP/2 + 1DELWX = TWOPI/NXPDELWY = TWOPI/NYPXLFAC = PI/ NYQXYLFAC = PI/ NYQY!! set operation constantsAI=0.0AD=0.0IF (LOP(10)) THENAI = PINCL * DEGRAAD = (-PDECL+90.0) * DEGRAELSE IF (LOP(11)) THENAI=PMINCL*DEGRAAD=(-PMDECL+90.0) * DEGRAEND IFSINI = SIN(AI)JuNE 2020 PREVIEW 42'