数値粘性項を評価する.
粘性係数は arareprm 中で, \nu = \frac{\alpha \DeltaX}{\Delta t} で決める. \alpha < 1/8 である. CReSS マニュアル参照.
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax use differentiate, only: ss_dx_fs, fs_dx_ss, ss_dz_sf, sf_dz_ss use arareprm, only: Cm, NuH, NuV
real(8), intent(in) :: ss_PotTemp_n(DimXMin:DimXMax, DimZMin:DimZMax) !温位
real(8), intent(out) :: ss_NumDiffScalar_n(DimXMin:DimXMax, DimZMin:DimZMax) !数値粘性項