gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

USER_PROBLEM: Correct a few comments.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2024-04-21 11:10:19 -03:00
parent a8be8f1710
commit 9e45431c89
1 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
sources/user_problem.F90
View File

@ -23,7 +23,7 @@
!! !!
!! module: USER_PROBLEM !! module: USER_PROBLEM
!! !!
!! This module provides subroutines to setup custom problem. !! This module provides subroutines to set up a custom problem.
!! !!
!!***************************************************************************** !!*****************************************************************************
! !
@ -91,14 +91,14 @@ module user_problem
! subroutine INITIALIZE_USER_PROBLEM: ! subroutine INITIALIZE_USER_PROBLEM:
! ---------------------------------- ! ----------------------------------
! !
! Subroutine initializes user problem. It could read problem parameters which ! Subroutine initializes the user problem. It can read problem parameters
! are used in all subroutines defining this specific problem. ! that are used in all subroutines defining this specific problem.
! !
! Arguments: ! Arguments:
! !
! problem - the problem name ! problem - the problem name
! rcount - the run count for restarted jobs ! rcount - the run count for restarted jobs
! verbose - a logical flag turning the information printing; ! verbose - a logical flag turning on information printing;
! status - an integer flag for error return value; ! status - an integer flag for error return value;
! !
!=============================================================================== !===============================================================================
@ -847,7 +847,7 @@ module user_problem
end do ! k = 1, nn end do ! k = 1, nn
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! calculate magnetic field components from vector potential ! calculate velocity components from stream function
! !
call curl(dh(1:3), pot(:,:,:,:), qpert(ivx:ivz,:,:,:)) call curl(dh(1:3), pot(:,:,:,:), qpert(ivx:ivz,:,:,:))
@ -1885,7 +1885,7 @@ module user_problem
if (resistive) then if (resistive) then
! diffusion of Bx through ! diffusion of Bx
! !
rterms(7) = rterms(7) & rterms(7) = rterms(7) &
#if NDIMS == 3 #if NDIMS == 3
@ -1896,7 +1896,7 @@ module user_problem
pdata%q(ibx,nb:ne,nm:ni, : ))) * dxz pdata%q(ibx,nb:ne,nm:ni, : ))) * dxz
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! diffusion of magnetic energy through the lower Y boundary ! diffusion of magnetic energy
! !
rterms(21) = rterms(21) & rterms(21) = rterms(21) &
#if NDIMS == 3 #if NDIMS == 3
@ -2139,10 +2139,10 @@ module user_problem
! ------------------- ! -------------------
! !
! Function calculates the logarithm of the hyperbolic cosine, which is ! Function calculates the logarithm of the hyperbolic cosine, which is
! the result of the integration of tanh(x). Direct calculation using ! the result of integrating tanh(x). Direct calculation using Fortran
! Fortran intrinsic subroutines fails for large values of x, therefore ! intrinsic subroutines fails for large values of x; therefore,
! the logarithm of cosh is approximated as |x| + log(1/2) for ! the logarithm of cosh is approximated as |x| + log(1/2) for |x| greater
! |x| > threshold. ! than a certain threshold.
! !
! Arguments: ! Arguments:
! !