fixed calculation of CurlRMS when grid(1)==1
This commit is contained in:
parent
a0f304021f
commit
ad32aa3ced
|
@ -691,6 +691,7 @@ real(pReal) function utilities_curlRMS()
|
|||
2.0_pReal*sum(real(curl_fourier)**2.0_pReal + aimag(curl_fourier)**2.0_pReal)
|
||||
enddo; enddo
|
||||
utilities_curlRMS = sqrt(utilities_curlRMS) * wgt
|
||||
if(grid(1) == 1_pInt) utilities_curlRMS = utilities_curlRMS * 0.5_pReal ! counted twice in case of grid(1) == 1
|
||||
|
||||
end function utilities_curlRMS
|
||||
|
||||
|
|
Loading…
Reference in New Issue