##// END OF EJS Templates
schain
RSS
bfield2 supports general path to call bfmodel
rflores -
r1573:4add2f3d5a7c
parent child
Show More
Side by Side Unified
Context file:
r1573:4add2f3d5a7c Loading diff...:
Show file before | Show file after | Hide comments
@@ -17,7 +17,7 c usgs, ms 964, box 25046 federal center, denver, co 80225
17 c
17 c
18 c ===============================================================
18 c ===============================================================
19 character*8 filmod(17)
19 character*8 filmod(17)
20 character*9 fqual
20 c character*31 fqual
21 dimension gh1(220), gh2(220), gha(224), ext(3), dtemod(17)
21 dimension gh1(220), gh2(220), gha(224), ext(3), dtemod(17)
22 data ext /3*0./
22 data ext /3*0./
23 data filmod / 'dgrf45', 'dgrf50',
23 data filmod / 'dgrf45', 'dgrf50',
@@ -26,7 +26,7 c ===============================================================
26 3 'dgrf85', 'dgrf90', 'dgrf95',
26 3 'dgrf85', 'dgrf90', 'dgrf95',
27 4 'dgrf00', 'dgrf05', 'dgrf10',
27 4 'dgrf00', 'dgrf05', 'dgrf10',
28 5 'dgrf15', 'igrf20', 'igrf20s'/
28 5 'dgrf15', 'igrf20', 'igrf20s'/
29 data fqual/"/bfmodel/"/
29 c data fqual/"/usr/local/lib/faraday/bfmodel/"/
30 data dtemod / 1945., 1950., 1955., 1960.,
30 data dtemod / 1945., 1950., 1955., 1960.,
31 1 1965., 1970., 1975., 1980., 1985., 1990., 1995., 2000.,
31 1 1965., 1970., 1975., 1980., 1985., 1990., 1995., 2000.,
32 2 2005.,2010.,2015.,2020.,2025./
32 2 2005.,2010.,2015.,2020.,2025./
@@ -36,6 +36,13 c data a2/40680925./, b2/40408588./
36 c
36 c
37 data rtd/57.29577951/
37 data rtd/57.29577951/
38 data tmp/0./,lp/0/
38 data tmp/0./,lp/0/
39 character(1024) :: fqual_temp
40 character(:), allocatable :: fqual
41 call get_path(fqual_temp)
42 c write(*,*) "L_BEF: ", fqual_temp, "L_BEF_end"
43 fqual = TRIM(fqual_temp)
44 c write(*,*) "L: ", fqual, "L_end"
45
39 flat=90. - theta*rtd
46 flat=90. - theta*rtd
40 flon=rtd*phi
47 flon=rtd*phi
41 c*****if previous time is not equal to current time
48 c*****if previous time is not equal to current time
@@ -62,8 +69,9 c write(*,-)tm,l
62 if(l .ne. lp)then
69 if(l .ne. lp)then
63 c*********if previous epoch not the same, read in new coefs
70 c*********if previous epoch not the same, read in new coefs
64 c write(*,fmt='(" read coefs"))')
71 c write(*,fmt='(" read coefs"))')
65 c write(*,*) fqual//filmod(l)
72 c write(*,*) "filmod", fqual//filmod(l)
66 call getshc2 (iu, fqual//filmod(l), nmax1, erad, gh1, ier)
73 call getshc2 (iu, fqual//filmod(l), nmax1, erad, gh1, ier)
74 c write(*,*) "AA: ", ier, filmod(l)
67 if(ier .ne. 0)then
75 if(ier .ne. 0)then
68 write(*,fmt='(" geobfield: read error=",i2," on ",a)')
76 write(*,fmt='(" geobfield: read error=",i2," on ",a)')
69 1 ier,filmod(l)
77 1 ier,filmod(l)
Show file before | Show file after | Hide comments
@@ -87,6 +87,18 setup(
87 ],
87 ],
88 )
88 )
89
89
90 main_path = os.getcwd()
91 child_path = '/schainf/mkfact/bfmodel/'
92 #print("main_path: ", main_path)
93 my_str = ' '+"ppath = "+"\""+main_path+"\"\n"
94 my_str_2 = ' '+"cpath = "+"\""+child_path+"\""
95
96 #'''
97 with open('./schainf/mkfact/get_path_1.f', 'r') as path1: data1 = path1.read()
98 with open('./schainf/mkfact/get_path_2.f', 'r') as path2: data2 = path2.read()
99 with open('./schainf/mkfact/get_path.f', 'w') as final: final.write(data1+my_str+my_str_2+data2)
100 #'''
101
90 from numpy.distutils.core import Extension, setup
102 from numpy.distutils.core import Extension, setup
91
103
92 setup(name='schainpy',
104 setup(name='schainpy',
@@ -97,7 +109,8 setup(name='schainpy',
97 "schainf/mkfact/lmdif1.f",
109 "schainf/mkfact/lmdif1.f",
98 "schainf/mkfact/mkfact.f",
110 "schainf/mkfact/mkfact.f",
99 "schainf/mkfact/r1mach.f",
111 "schainf/mkfact/r1mach.f",
100 "schainf/mkfact/bfield2.f"])
112 "schainf/mkfact/bfield2.f",
113 "schainf/mkfact/get_path.f"])
101 ]
114 ]
102 )
115 )
103
116
General Comments 0
You need to be logged in to leave comments. Login now