133 |
& globfld3d(i,bi,ip,j,bj,jp,k) * |
& globfld3d(i,bi,ip,j,bj,jp,k) * |
134 |
& sqrt(weightfld(k,bi,bj)) |
& sqrt(weightfld(k,bi,bj)) |
135 |
else |
else |
|
cph( |
|
|
print *, 'ph-nondim bef. ', k, j, i, |
|
|
& globfld3d(i,bi,ip,j,bj,jp,k), |
|
|
& weightfld(k,bi,bj) |
|
|
cph) |
|
136 |
cbuff(cbuffindex) = |
cbuff(cbuffindex) = |
137 |
& globfld3d(i,bi,ip,j,bj,jp,k) / |
& globfld3d(i,bi,ip,j,bj,jp,k) / |
138 |
& sqrt(weightfld(k,bi,bj)) |
& sqrt(weightfld(k,bi,bj)) |
|
cph( |
|
|
write(6,'(A,4I5,F10.2)'), 'ph-nondim aft. ', |
|
|
& k, j, i, cbuffindex, |
|
|
& cbuff(cbuffindex) |
|
|
cph) |
|
139 |
endif |
endif |
140 |
#else |
#else |
141 |
cbuff(cbuffindex) = globfld3d(i,bi,ip,j,bj,jp,k) |
cbuff(cbuffindex) = globfld3d(i,bi,ip,j,bj,jp,k) |