2013年11月14日木曜日

subroutine "do" cannot be renamed.

Since yesterday, I started adding parameters and equations in my code. The compilation was successful, but the output of the model simulation was not promising:

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 件のコメント:

コメントを投稿