The output (z vs t of DMS) was a constant field with initial value. So the calculation was not done. It turns out that the issue arises from the subroutine where you deal with the RHS of the equation:
In my "dms.F90":
subroutine do(self,_ARGUMENTS_DO_)
class (type_uvic_dms),intent(in) :: self
_DECLARE_ARGUMENTS_DO_
real(rk) :: dms
_LOOP_BEGIN_
_GET_(self%id_dms,dms)
dms=dms+2
_SET_ODE_(self%id_dms,self%q)
_LOOP_END_
end subroutine do
I used to name this subroutine "do_dms" as I thought it can be customizable. It seems, however, that it cannot be changed. If I changed to something else other than "do", the computation within the subroutine was not done. Not sure why. Anyway, I will just keep the subroutine name as "do" and not worry about it...
0 件のコメント:
コメントを投稿