integer function c2i_convert( cval )
! 文字型を実数型に変換する
  implicit none
  character(*), intent(in) :: cval  ! 変換する文字

  read(cval,*) c2i_convert

  return
end function

real function c2r_convert( cval )
! 文字型を実数型に変換する
  implicit none
  character(*), intent(in) :: cval  ! 変換する文字

  read(cval,*) c2r_convert

  return
end function

integer function counter_day( stime, etime )
! 開始日から終了日までの日数をカウントする.
  implicit none
  type(dtime), intent(in) :: stime  ! 開始日
  type(dtime), intent(in) :: etime  ! 終了日
  integer, parameter, dimension(13) :: month=(/31,28,31,30,31,30, 31,31,30,31,30,31,29/)
  integer :: nt, nm, days, i, year_tmp, year_fact
  integer :: nsy, nsm, nsd, ney, nem, ned

  nt=etime%year_d-stime%year_d+1
  nm=etime%month_d-stime%month_d+1
  days=0
  nsy=stime%year_d
  nsm=stime%month_d
  nsd=stime%day_d
  ney=etime%year_d
  nem=etime%month_d
  ned=etime%day_d

!-- 一年分ある場合, 閏年で場合分けする.
  if(nt>2)then
     do i=2,nt-1
        year_tmp=nsy+i-1
        if(mod(year_tmp,4)==0)then
           year_fact=366
        else
           year_fact=365
        end if
        days=days+year_fact
     end do
  end if

!-- 上の処理を行っているので, あと足し合わせるのは開始年と終了年の日数のみ.
!-- まず, 1 ケ月まるまる存在する月を足し合わせる.
!-- ただし, 開始日が 12 月, 終了日が 1 月の場合は, 日を足す過程で足されるので,
!-- ここではカウントしない.
  if(nt>1)then  ! 年がまたがる場合
     if(nsm<12)then  ! 開始日が 12 月ではない場合.
        do i=nsm+1,12
           if(i==2.and.mod(nsy,4)==0)then
              days=days+month(13)
           else
              days=days+month(i)
           end if
        end do
     end if

     if(nem>1)then  ! 終了日が 1 月ではない場合.
        do i=1,nem-1
           if(i==2.and.mod(ney,4)==0)then
              days=days+month(13)
           else
              days=days+month(i)
           end if
        end do
     end if
  else
     if(nm>2)then  ! 年がまたがらず, 月がまたがり, かつ間に 1 ヶ月分
                   ! まるまる入っている場合.
        do i=nsm+1,nem-1
           if(i==2.and.mod(ney,4)==0)then
              days=days+month(13)
           else
              days=days+month(i)
           end if
        end do
     end if
  end if

!-- 後は, 開始月と終了月の日数を足し合わせる.
  if(nt/=1.or.nm/=1)then  ! 年も月もかぶっていない場合
     if(mod(nsy,4)==0.and.nsm==2)then
        days=days+month(13)-nsd
     else
        days=days+month(nsm)-nsd
     end if

     if(mod(ney,4)==0.and.nem==2)then
        days=days+ned
     else
        days=days+ned
     end if
  else
     days=days+ned-nsd
  end if

  days=days+1
  counter_day=days

  return
end function

character(100) function i2c_convert( ival, forma )
! 実数型を文字型に変換する
  implicit none
  integer, intent(in) :: ival  ! 変換する整数
  character(*), intent(in), optional :: forma  ! 指定するフォーマット
  character(100) :: tmp

  if(present(forma))then
     write(tmp,trim(forma)) ival
  else
     write(tmp,*) ival
  end if

  i2c_convert=tmp

  return
end function

character(100) function r2c_convert( rval, forma )
! 実数型を文字型に変換する
  implicit none
  real, intent(in) :: rval  ! 変換する実数
  character(*), intent(in), optional :: forma  ! 指定するフォーマット
  character(100) :: tmp

  if(present(forma))then
     write(tmp,trim(forma)) rval
  else
     write(tmp,*) rval
  end if

  r2c_convert=tmp

  return
end function

subroutine rand_make(L,output)
  ! 任意の桁数で擬似乱数を生成するサブルーチン
  ! 混合合同法という乱数生成アルゴリズムを用いて擬似乱数を生成
  ! $x_{n+1}=a\times x_{n}+b (mod \; L)$
  implicit none
  integer, intent(in) :: L  ! 出力する最大桁数 + 1 の数値
  integer, intent(inout) :: output  ! 出力される乱数
  integer :: a, b, x0, i, input
  integer, external :: time

  input=time()
  input=mod(input,L)
  write(*,*) input
  a=11
  b=12
  x0=input

  do i=1,10
     x0=a*x0+b
     x0=mod(x0,L)
  end do

  output=x0

end subroutine