32 |
_RL gam1, gam2 |
_RL gam1, gam2 |
33 |
|
|
34 |
C magic parameters( M. Satoh, p. 314) |
C magic parameters( M. Satoh, p. 314) |
35 |
gam1 = 1./3.-2.*A1/B1 |
gam1 = 1. _d 0/3. _d 0 - 2. _d 0 *A1/B1 |
36 |
gam2 = (B2+6.*A1)/B1 |
gam2 = (B2+6. _d 0*A1)/B1 |
37 |
alpha1 = 3.*A2*gam1 |
alpha1 = 3. _d 0 *A2*gam1 |
38 |
alpha2 = 3.*A2*(gam1+gam2) |
alpha2 = 3. _d 0 *A2*(gam1+gam2) |
39 |
beta1 = A1*B1*(gam1-C1) |
beta1 = A1*B1*(gam1-C1) |
40 |
beta2 = A1*(B1*(gam1-C1) + 6.*A1 + 3.*A2) |
beta2 = A1*(B1*(gam1-C1) + 6. _d 0 *A1 + 3. _d 0 *A2) |
41 |
beta3 = A2*B1*gam1 |
beta3 = A2*B1*gam1 |
42 |
beta4 = A2*(B1*(gam1+gam2) - 3.*A1) |
beta4 = A2*(B1*(gam1+gam2) - 3. _d 0 *A1) |
43 |
|
|
44 |
C----------------------------------------------------------------------- |
C----------------------------------------------------------------------- |
45 |
C Initialize MY82 variables MYviscAr, MYdiffKrT |
C Initialize MY82 variables MYviscAr, MYdiffKrT |