subroutine rdguerin(nopr,iout,&
   u10,mu,mc,mus,muc,c40,c04,c22,c12,c30,thetapmax)

! ******************************
   integer, parameter :: ndatguerin=25
   real :: u10(ndatguerin)
   real :: mu(ndatguerin),mc(ndatguerin),mus(ndatguerin),muc(ndatguerin)
   real :: c40(ndatguerin),c04(ndatguerin),c22(ndatguerin),c12(ndatguerin)
   real :: c30(ndatguerin),thetapmax(ndatguerin)

! ******************************
! See Guerin, Capelle, Jean-Michael Hartmann
! Revisiting the Cox and Munk wave-slope statistics using IASI observations
! of the sea surface
! October 11, 2022 preprint https://arxiv.org/abs/2210.05456

! Data is from thhe paper's Table 1

! *****
! Output ascii file f.out
   idat=20
  open(idat,form='formatted',file='guerin.dat',status='unknown')

! *****
  do i=1,ndatguerin
     read(idat,*) u10val,a1,a2,a3,a4,b1,b2,b3,b4,b5,c1
    u10(i)=u10val
    mu(i)=a1/100.0
    mc(i)=a2/100.0
    mus(i)=a3/100.0
    muc(i)=a4/100.0
    c40(i)=b1
    c04(i)=b2
    c22(i)=b3
    c12(i)=b4
    c30(i)=b5
    thetapmax(i)=c1
  enddo

  close (idat)

! ******************************
   if (nopr .eq. 1) then 
    write(iout,fmt=300)
    300 format(/,2x,"rdguerin:  i,u10(i),mu(i),mc(i),mus(i),muc(i),",/,&
    2x,"c40(i),c04(i),c22(i),c12(i),c30(i),thetapmax(i)")
    do i=1,ndatguerin
     write(iout,fmt=310) i,u10(i),mu(i),mc(i),mus(i),muc(i),&
     c40(i),c04(i),c22(i),c12(i),c30(i),thetapmax(i)
    enddo
    310 format(2x,i3,11(1x,f10.4))
   endif

! ******************************
  return
  end