-
Notifications
You must be signed in to change notification settings - Fork 1
/
Intensitymax_parsing.f90
59 lines (47 loc) · 1.74 KB
/
Intensitymax_parsing.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
!12345----------------------------------------------
!>@file:Intensitymax.f90
!
!>@author:Changwan Sun/[email protected]
!>date: 2016. 06. 27. (월) 10:38:28 KST
!
!2345----------------------------------------------
Program Intensitymax
implicit none
Complex :: Pixel(1000)
Integer :: TI,i,j,Maxlocation
Real :: Al=0.009896
Real :: Lat, Long, Alt, H_offset,H_nadir, Start_range,Intensity(1000)
Open(10, File='W58.15Chandrayaan_N00_N40.dat',form='unformatted',status='old')
Open(11,File='W58.15_parsing_12_15.txt',status='replace')
!Open(12,File="Kaguya_direction.txt",status='replace')
j=0
100 Read(10,End=200) TI,Lat,Long,Alt,H_offset,H_nadir,Start_range,Pixel
! Write(11,*)'TI=',TI
! Write(11,*)'Lat=',Lat
! Write(11,*)'Long=',Long
! Write(11,*)'Alt=',Alt
! Write(11,*)'H_offset=',H_offset
! Write(11,*) 'H_nadir=',H_nadir
! Write(11,*)'Start_range=',Start_range
!!DTM Latittude!!
If((12.000244+Al) <= Lat.AND.Lat <=(15.000000-Al)) then
Do i=1,1000
Intensity(i)=20*log10(ABS(pixel(i)))
!Write(11,*)'intensity(',i,')',Intensity(i)
End do
!-----Max intensity and location in the intensity array-----!
! Write(11,*)'Max=',Maxval(Intensity)
! Write(11,*)'Maxlocation=',25*Maxloc(Intensity)
!-----------------------------------------------------------!
!----------Kaguya direction checkout------------------------!
!Write(12,*) Lat, Long, Alt, TI
!---------------------------------------------------------------!
Write(11,*)Lat,Long,Alt,25*Maxloc(Intensity),Start_range
!----------------------------------------------------------------!
End if
j=j+1
!Write(11,*)'j=',j,'times'
Goto 100
200 Close(10)
Stop
End Program Intensitymax