Skip to content

Commit 25beed2

Browse files
fangjianfangjian
fangjian
authored and
fangjian
committed
minor correction on tgvsolver
1 parent 7c94c6e commit 25beed2

3 files changed

Lines changed: 114 additions & 14 deletions

File tree

‎miniapps/tgv_solver_3d/state.ref_128‎

Lines changed: 100 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,100 @@
1+
0 0.0000000000000000 0.12499999999996174 0.37499999999488698
2+
1 1.0000000000000000E-003 0.12499953071277797 0.37499863119528576
3+
2 2.0000000000000000E-003 0.12499906142723267 0.37499734052322387
4+
3 3.0000000000000001E-003 0.12499859214342805 0.37499612797711518
5+
4 4.0000000000000001E-003 0.12499812286167343 0.37499499355547478
6+
5 5.0000000000000001E-003 0.12499765358230532 0.37499393725797714
7+
6 6.0000000000000001E-003 0.12499718430591819 0.37499295908495961
8+
7 7.0000000000000001E-003 0.12499671503318036 0.37499205903747740
9+
8 8.0000000000000002E-003 0.12499624576491108 0.37499123711744031
10+
9 9.0000000000000011E-003 0.12499577650223034 0.37499049332710233
11+
10 1.0000000000000002E-002 0.12499530724606374 0.37498982766986294
12+
11 1.1000000000000003E-002 0.12499483799781376 0.37498924014984564
13+
12 1.2000000000000004E-002 0.12499436875874617 0.37498873077158523
14+
13 1.3000000000000005E-002 0.12499389953036803 0.37498829954026042
15+
14 1.4000000000000005E-002 0.12499343031427904 0.37498794646235062
16+
15 1.5000000000000006E-002 0.12499296111217625 0.37498767154424323
17+
16 1.6000000000000007E-002 0.12499249192570040 0.37498747479329597
18+
17 1.7000000000000008E-002 0.12499202275665258 0.37498735621785523
19+
18 1.8000000000000009E-002 0.12499155360707047 0.37498731582650030
20+
19 1.9000000000000010E-002 0.12499108447865999 0.37498735362835206
21+
20 2.0000000000000011E-002 0.12499061537345962 0.37498746963357643
22+
21 2.1000000000000012E-002 0.12499014629344375 0.37498766385234200
23+
22 2.2000000000000013E-002 0.12498967724050493 0.37498793629598659
24+
23 2.3000000000000013E-002 0.12498920821673662 0.37498828697601749
25+
24 2.4000000000000014E-002 0.12498873922393530 0.37498871590455912
26+
25 2.5000000000000015E-002 0.12498827026411298 0.37498922309456928
27+
26 2.6000000000000016E-002 0.12498780133905857 0.37498980855875019
28+
27 2.7000000000000017E-002 0.12498733245059222 0.37499047231097327
29+
28 2.8000000000000018E-002 0.12498686360054233 0.37499121436516575
30+
29 2.9000000000000019E-002 0.12498639479050826 0.37499203473599890
31+
30 3.0000000000000020E-002 0.12498592602204131 0.37499293343818763
32+
31 3.1000000000000021E-002 0.12498545729657723 0.37499391048724434
33+
32 3.2000000000000021E-002 0.12498498861558544 0.37499496589878706
34+
33 3.3000000000000022E-002 0.12498451998009355 0.37499609968853026
35+
34 3.4000000000000023E-002 0.12498405139132347 0.37499731187297747
36+
35 3.5000000000000024E-002 0.12498358285015605 0.37499860246879269
37+
36 3.6000000000000025E-002 0.12498311435737584 0.37499997149239361
38+
37 3.7000000000000026E-002 0.12498264591358588 0.37500141896122424
39+
38 3.8000000000000027E-002 0.12498217751921181 0.37500294489224623
40+
39 3.9000000000000028E-002 0.12498170917454515 0.37500454930303739
41+
40 4.0000000000000029E-002 0.12498124087967877 0.37500623221112228
42+
41 4.1000000000000029E-002 0.12498077263448849 0.37500799363408949
43+
42 4.2000000000000030E-002 0.12498030443871293 0.37500983358980522
44+
43 4.3000000000000031E-002 0.12497983629189796 0.37501175209612536
45+
44 4.4000000000000032E-002 0.12497936819325756 0.37501374917094515
46+
45 4.5000000000000033E-002 0.12497890014202934 0.37501582483195733
47+
46 4.6000000000000034E-002 0.12497843213711375 0.37501797909711709
48+
47 4.7000000000000035E-002 0.12497796417723200 0.37502021198425572
49+
48 4.8000000000000036E-002 0.12497749626088207 0.37502252351116949
50+
49 4.9000000000000037E-002 0.12497702838649806 0.37502491369522073
51+
50 5.0000000000000037E-002 0.12497656055211215 0.37502738255413676
52+
51 5.1000000000000038E-002 0.12497609275585499 0.37502993010507923
53+
52 5.2000000000000039E-002 0.12497562499546319 0.37503255636513283
54+
53 5.3000000000000040E-002 0.12497515726856288 0.37503526135123494
55+
54 5.4000000000000041E-002 0.12497468957257352 0.37503804508022442
56+
55 5.5000000000000042E-002 0.12497422190490641 0.37504090756811886
57+
56 5.6000000000000043E-002 0.12497375426266037 0.37504384883161807
58+
57 5.7000000000000044E-002 0.12497328664284370 0.37504686888572142
59+
58 5.8000000000000045E-002 0.12497281904248315 0.37504996774625476
60+
59 5.9000000000000045E-002 0.12497235145820414 0.37505314542789531
61+
60 6.0000000000000046E-002 0.12497188388682458 0.37505640194566764
62+
61 6.1000000000000047E-002 0.12497141632489267 0.37505973731371089
63+
62 6.2000000000000048E-002 0.12497094876895323 0.37506315154535913
64+
63 6.3000000000000042E-002 0.12497048121538416 0.37506664465424738
65+
64 6.4000000000000043E-002 0.12497001366064098 0.37507021665292339
66+
65 6.5000000000000044E-002 0.12496954610105623 0.37507386755378974
67+
66 6.6000000000000045E-002 0.12496907853303740 0.37507759736823437
68+
67 6.7000000000000046E-002 0.12496861095281558 0.37508140610769641
69+
68 6.8000000000000047E-002 0.12496814335681772 0.37508529378249567
70+
69 6.9000000000000047E-002 0.12496767574134479 0.37508926040258306
71+
70 7.0000000000000048E-002 0.12496720810280398 0.37509330597751839
72+
71 7.1000000000000049E-002 0.12496674043771912 0.37509743051591754
73+
72 7.2000000000000050E-002 0.12496627274248798 0.37510163402535074
74+
73 7.3000000000000051E-002 0.12496580501370057 0.37510591651384417
75+
74 7.4000000000000052E-002 0.12496533724816614 0.37511027798753083
76+
75 7.5000000000000053E-002 0.12496486944257798 0.37511471845273581
77+
76 7.6000000000000054E-002 0.12496440159384881 0.37511923791472412
78+
77 7.7000000000000055E-002 0.12496393369905512 0.37512383637771318
79+
78 7.8000000000000055E-002 0.12496346575543192 0.37512851384599749
80+
79 7.9000000000000056E-002 0.12496299776030585 0.37513327032208060
81+
80 8.0000000000000057E-002 0.12496252971119393 0.37513810580920420
82+
81 8.1000000000000058E-002 0.12496206160579457 0.37514302030795987
83+
82 8.2000000000000059E-002 0.12496159344210499 0.37514801381982082
84+
83 8.3000000000000060E-002 0.12496112521815415 0.37515308634492134
85+
84 8.4000000000000061E-002 0.12496065693232325 0.37515823788199149
86+
85 8.5000000000000062E-002 0.12496018858312687 0.37516346843011245
87+
86 8.6000000000000063E-002 0.12495972016943949 0.37516877798695109
88+
87 8.7000000000000063E-002 0.12495925169014827 0.37517416654956909
89+
88 8.8000000000000064E-002 0.12495878314462720 0.37517963411398431
90+
89 8.9000000000000065E-002 0.12495831453230338 0.37518518067612566
91+
90 9.0000000000000066E-002 0.12495784585296521 0.37519080623035578
92+
91 9.1000000000000067E-002 0.12495737710670637 0.37519651077100746
93+
92 9.2000000000000068E-002 0.12495690829374485 0.37520229429110447
94+
93 9.3000000000000069E-002 0.12495643941455443 0.37520815678348912
95+
94 9.4000000000000070E-002 0.12495597046999414 0.37521409823972929
96+
95 9.5000000000000070E-002 0.12495550146109116 0.37522011865123156
97+
96 9.6000000000000071E-002 0.12495503238925176 0.37522621800784006
98+
97 9.7000000000000072E-002 0.12495456325596352 0.37523239630033922
99+
98 9.8000000000000073E-002 0.12495409406305148 0.37523865351684493
100+
99 9.9000000000000074E-002 0.12495362481257724 0.37524498964658615

‎miniapps/tgv_solver_3d/tgvsolver.F90‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -163,7 +163,7 @@ subroutine program_init
163163

164164
alfa_filter=0.49d0
165165

166-
maxstep=10
166+
maxstep=100
167167

168168
write(*,'(3(A,I0))')' ** dimension set: ',im,' x ',jm,' x ',km
169169

@@ -223,10 +223,10 @@ subroutine flowfield_init
223223
enddo
224224
enddo
225225

226-
qrhs=0.d0
227-
228226
nstep=0
227+
229228
time =0.d0
229+
230230
deltat=1.d-3
231231

232232
print*,' ** flow field initilised'
@@ -669,7 +669,7 @@ subroutine diffusion
669669

670670
skk=num1d3*(s11+s22+s33)
671671

672-
miu=miucal(tmp(i,j,k))
672+
miu=miucal(tmp(i,j,k))/reynolds
673673
miu2=2.d0*miu
674674
hcc=(miu/prandtl)/const5
675675

@@ -776,10 +776,10 @@ subroutine diffusion
776776
df(:,n)=fds%central(fds_compact_k,f=ff(:,n),dim=km)
777777
enddo
778778

779-
qrhs(i,j,0:km,2)=qrhs(i,j,0:km,2)+df(0:km,1)/dz
780-
qrhs(i,j,0:km,3)=qrhs(i,j,0:km,3)+df(0:km,2)/dz
781-
qrhs(i,j,0:km,4)=qrhs(i,j,0:km,4)+df(0:km,3)/dz
782-
qrhs(i,j,0:km,5)=qrhs(i,j,0:km,5)+df(0:km,4)/dz
779+
qrhs(i,j,0:km,2)=qrhs(i,j,0:km,2)+df(0:km,2)/dz
780+
qrhs(i,j,0:km,3)=qrhs(i,j,0:km,3)+df(0:km,3)/dz
781+
qrhs(i,j,0:km,4)=qrhs(i,j,0:km,4)+df(0:km,4)/dz
782+
qrhs(i,j,0:km,5)=qrhs(i,j,0:km,5)+df(0:km,5)/dz
783783
!
784784
enddo
785785
enddo
@@ -973,6 +973,8 @@ subroutine rhscal(comptime)
973973
call cpu_time(tstart)
974974
endif
975975
!
976+
qrhs=0.d0
977+
976978
call gradcal
977979
!
978980
call convection
@@ -1033,8 +1035,6 @@ subroutine rk3(comptime)
10331035
endif
10341036
!
10351037
do rkstep=1,3
1036-
!
1037-
qrhs=0.d0
10381038

10391039
call bchomo
10401040

‎src/pp.F90‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -2003,10 +2003,10 @@ subroutine fieldview(flowfile,outputfile,viewmode,inputfile)
20032003
call h5_read2dfrom3d( p_2d,im,jm,km, 'p',flowfile,kslice=0)
20042004
call h5_read2dfrom3d( t_2d,im,jm,km, 't',flowfile,kslice=0)
20052005
!
2006-
call writeprvbin(outputfile,x_2d,'x',y_2d,'y',ro_2d,'ro',u_2d,'u', &
2007-
v_2d,'v',p_2d,'p',t_2d,'t',im,jm)
2008-
! call tecbin(outputfile,x_2d,'x',y_2d,'y',ro_2d,'ro',u_2d,'u', &
2009-
! v_2d,'v',p_2d,'p',t_2d,'t')
2006+
! call writeprvbin(outputfile,x_2d,'x',y_2d,'y',ro_2d,'ro',u_2d,'u', &
2007+
! v_2d,'v',p_2d,'p',t_2d,'t',im,jm)
2008+
call tecbin(outputfile,x_2d,'x',y_2d,'y',ro_2d,'ro',u_2d,'u', &
2009+
v_2d,'v',p_2d,'p',t_2d,'t')
20102010
elseif(viewmode=='3d') then
20112011
allocate(x(0:im,0:jm,0:km),y(0:im,0:jm,0:km),z(0:im,0:jm,0:km))
20122012
elseif(viewmode=='1d') then

0 commit comments

Comments
 (0)