[580] | 1 | !-------------------------------------------------------------------------- |
---|
| 2 | !---------------------------- caldyn_slow_hydro ---------------------------------- |
---|
| 3 | DO l = ll_begin, ll_end |
---|
| 4 | !DIR$ SIMD |
---|
| 5 | DO ij=ij_begin_ext, ij_end_ext |
---|
| 6 | uu = .5*(rhodz(ij,l)+rhodz(ij+t_right,l))*u(ij+u_right,l) |
---|
| 7 | hflux(ij+u_right,l) = uu*le_de(ij+u_right) |
---|
| 8 | uu = .5*(rhodz(ij,l)+rhodz(ij+t_lup,l))*u(ij+u_lup,l) |
---|
| 9 | hflux(ij+u_lup,l) = uu*le_de(ij+u_lup) |
---|
| 10 | uu = .5*(rhodz(ij,l)+rhodz(ij+t_ldown,l))*u(ij+u_ldown,l) |
---|
| 11 | hflux(ij+u_ldown,l) = uu*le_de(ij+u_ldown) |
---|
| 12 | END DO |
---|
| 13 | END DO |
---|
| 14 | DO l = ll_begin, ll_end |
---|
| 15 | !DIR$ SIMD |
---|
| 16 | DO ij=ij_begin, ij_end |
---|
| 17 | ke=0.d0 |
---|
| 18 | ke = ke + le_de(ij+u_rup)*u(ij+u_rup,l)**2 |
---|
| 19 | ke = ke + le_de(ij+u_lup)*u(ij+u_lup,l)**2 |
---|
| 20 | ke = ke + le_de(ij+u_left)*u(ij+u_left,l)**2 |
---|
| 21 | ke = ke + le_de(ij+u_ldown)*u(ij+u_ldown,l)**2 |
---|
| 22 | ke = ke + le_de(ij+u_rdown)*u(ij+u_rdown,l)**2 |
---|
| 23 | ke = ke + le_de(ij+u_right)*u(ij+u_right,l)**2 |
---|
| 24 | BERNI(ij,l)=ke*(.25/Ai(ij)) |
---|
| 25 | END DO |
---|
| 26 | END DO |
---|
| 27 | IF(zero) THEN |
---|
| 28 | DO l = ll_begin, ll_end |
---|
| 29 | !DIR$ SIMD |
---|
| 30 | DO ij=ij_begin, ij_end |
---|
| 31 | du(ij+u_right,l) = ne_right*(berni(ij,l)-berni(ij+t_right,l)) ! minus gradient |
---|
| 32 | du(ij+u_lup,l) = ne_lup*(berni(ij,l)-berni(ij+t_lup,l)) ! minus gradient |
---|
| 33 | du(ij+u_ldown,l) = ne_ldown*(berni(ij,l)-berni(ij+t_ldown,l)) ! minus gradient |
---|
| 34 | END DO |
---|
| 35 | END DO |
---|
| 36 | ELSE |
---|
| 37 | DO l = ll_begin, ll_end |
---|
| 38 | !DIR$ SIMD |
---|
| 39 | DO ij=ij_begin, ij_end |
---|
| 40 | du(ij+u_right,l) = du(ij+u_right,l) + ne_right*(berni(ij,l)-berni(ij+t_right,l)) ! minus gradient |
---|
| 41 | du(ij+u_lup,l) = du(ij+u_lup,l) + ne_lup*(berni(ij,l)-berni(ij+t_lup,l)) ! minus gradient |
---|
| 42 | du(ij+u_ldown,l) = du(ij+u_ldown,l) + ne_ldown*(berni(ij,l)-berni(ij+t_ldown,l)) ! minus gradient |
---|
| 43 | END DO |
---|
| 44 | END DO |
---|
| 45 | END IF |
---|
| 46 | !---------------------------- caldyn_slow_hydro ---------------------------------- |
---|
| 47 | !-------------------------------------------------------------------------- |
---|