program hardcode_vcoord
  integer, parameter :: lev=26, ilev=lev+1
  real*8, parameter :: hyai(ilev)=(/ 0.00219406700000001, 0.00489520900000001, 0.009882418, 0.01805201, &
    0.02983724, 0.0446233400000002, 0.0616058700000002, 0.0785124300000004, &
    0.0773127100000002, 0.0759013100000003, 0.0742408600000002, &
    0.0722874400000002, 0.0699893299999998, 0.06728574, 0.06410509, &
    0.0603632200000002, 0.0559611100000001, 0.0507822500000001, &
    0.0446896000000001, 0.0375219099999999, 0.0290894900000001, 0.02084739, &
    0.01334443, 0.00708499000000001, 0.00252136, 0, 0 /)

  real*8, parameter :: hybi(ilev)=(/ 0, 0, 0, 0, 0, 0, 0, 0, 0.01505309, 0.03276228, 0.05359622, &
    0.0781062700000006, 0.1069411, 0.140863700000001, 0.180772, 0.227722, &
    0.282956200000001, 0.347936400000002, 0.4243822, 0.514316800000003, &
    0.620120200000002, 0.723535500000004, 0.817676800000001, &
    0.896215300000001, 0.953476100000003, 0.9851122, 1 /)
  real*8, parameter :: hyam(lev)= (/0.00354463800000002, 0.00738881300000002, 0.0139672100000001, &
    0.02394463, 0.0372302900000001, 0.0531146000000002, 0.0700591400000001, &
    0.0779125700000001, 0.0766070100000005, 0.0750710800000003, 0.07326415, &
    0.0711383900000002, 0.0686375400000006, 0.06569541, 0.06223416, &
    0.0581621700000001, 0.0533716900000001, 0.0477359300000002, &
    0.0411057500000001, 0.0333056999999999, 0.0249684400000001, &
    0.0170959100000001, 0.0102147100000001, 0.00480317500000002, 0.001260685, 0/)

  real*8, parameter :: hybm(lev)=(/ 0, 0, 0, 0, 0, 0, 0, 0.00752654999999999, 0.0239076800000001, &
    0.0431792499999999, 0.0658512300000003, 0.0925236599999999, 0.1239024, &
    0.1608178, 0.2042469, 0.2553391, 0.3154463, 0.386159300000001, &
    0.469349500000002, 0.567218400000002, 0.671827799999999, &
    0.770606100000002, 0.856946000000001, 0.9248458, 0.969294100000003, &
    0.992556100000005 /)

  character(len=80) :: fname
  integer :: k


  fname="cami-26.fbin"
  open(UNIT=11,FILE=fname,form='unformatted',status='unknown',iostat=ios)

  fname="camm-26.fbin"
  open(UNIT=12,FILE=fname,form='unformatted',status='unknown',iostat=ios)

  do k=1,ilev
     write(11) hyai(k)
     write(11) hybi(k)
  end do
  
  do k=1,lev
     write(12) hyam(k)
     write(12) hybm(k)
  end do

  close(11)
  close(12)

end program hardcode_vcoord
