! This file is generated by 'scripts/generate.py' using
!    /src/xbeachlibrary/variables.def
!    /src/xbeachlibrary/params.def
!    src/xbeachlibray/templates/set_var.mako
! Advice: do not edit this file, but above mentioned files.



 case("x")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%x))
   spointer%x(:,:) = x_2d_double_ptr

 case("y")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%y))
   spointer%y(:,:) = x_2d_double_ptr

 case("xz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%xz))
   spointer%xz(:,:) = x_2d_double_ptr

 case("yz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%yz))
   spointer%yz(:,:) = x_2d_double_ptr

 case("xu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%xu))
   spointer%xu(:,:) = x_2d_double_ptr

 case("yu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%yu))
   spointer%yu(:,:) = x_2d_double_ptr

 case("xv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%xv))
   spointer%xv(:,:) = x_2d_double_ptr

 case("yv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%yv))
   spointer%yv(:,:) = x_2d_double_ptr

 case("dsu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsu))
   spointer%dsu(:,:) = x_2d_double_ptr

 case("dsv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsv))
   spointer%dsv(:,:) = x_2d_double_ptr

 case("dsz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsz))
   spointer%dsz(:,:) = x_2d_double_ptr

 case("dsc")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsc))
   spointer%dsc(:,:) = x_2d_double_ptr

 case("dnu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dnu))
   spointer%dnu(:,:) = x_2d_double_ptr

 case("dnv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dnv))
   spointer%dnv(:,:) = x_2d_double_ptr

 case("dnz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dnz))
   spointer%dnz(:,:) = x_2d_double_ptr

 case("dnc")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dnc))
   spointer%dnc(:,:) = x_2d_double_ptr

 case("dsdnui")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsdnui))
   spointer%dsdnui(:,:) = x_2d_double_ptr

 case("dsdnvi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsdnvi))
   spointer%dsdnvi(:,:) = x_2d_double_ptr

 case("dsdnzi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dsdnzi))
   spointer%dsdnzi(:,:) = x_2d_double_ptr

 case("alfaz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%alfaz))
   spointer%alfaz(:,:) = x_2d_double_ptr

 case("alfau")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%alfau))
   spointer%alfau(:,:) = x_2d_double_ptr

 case("alfav")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%alfav))
   spointer%alfav(:,:) = x_2d_double_ptr

 case("sdist")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%sdist))
   spointer%sdist(:,:) = x_2d_double_ptr

 case("ndist")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ndist))
   spointer%ndist(:,:) = x_2d_double_ptr

 case("dx")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%dx = x_0d_double_ptr

 case("dy")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%dy = x_0d_double_ptr

 case("xori")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%xori = x_0d_double_ptr

 case("yori")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%yori = x_0d_double_ptr

 case("alfa")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%alfa = x_0d_double_ptr

 case("posdwn")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%posdwn = x_0d_double_ptr

 case("thetamin")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%thetamin = x_0d_double_ptr

 case("thetamax")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%thetamax = x_0d_double_ptr

 case("nx")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%nx = x_0d_int_ptr

 case("ny")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%ny = x_0d_int_ptr

 case("zs01")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%zs01 = x_0d_double_ptr

 case("zs02")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%zs02 = x_0d_double_ptr

 case("zs03")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%zs03 = x_0d_double_ptr

 case("zs04")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%zs04 = x_0d_double_ptr

 case("xyzs01")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%xyzs01))
   spointer%xyzs01(:) = x_1d_double_ptr

 case("xyzs02")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%xyzs02))
   spointer%xyzs02(:) = x_1d_double_ptr

 case("xyzs03")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%xyzs03))
   spointer%xyzs03(:) = x_1d_double_ptr

 case("xyzs04")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%xyzs04))
   spointer%xyzs04(:) = x_1d_double_ptr

 case("tidelen")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%tidelen = x_0d_int_ptr

 case("windlen")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%windlen = x_0d_int_ptr

 case("zb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zb))
   spointer%zb(:,:) = x_2d_double_ptr

 case("zb0")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zb0))
   spointer%zb0(:,:) = x_2d_double_ptr

 case("theta")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%theta))
   spointer%theta(:) = x_1d_double_ptr

 case("theta_s")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%theta_s))
   spointer%theta_s(:) = x_1d_double_ptr

 case("ntheta")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%ntheta = x_0d_int_ptr

 case("ntheta_s")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%ntheta_s = x_0d_int_ptr

 case("dtheta")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%dtheta = x_0d_double_ptr

 case("dtheta_s")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%dtheta_s = x_0d_double_ptr

 case("theta0")
   call c_f_pointer(xptr, x_0d_double_ptr) ! skip shape in case of rank 0
   spointer%theta0 = x_0d_double_ptr

 case("thetamean")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%thetamean))
   spointer%thetamean(:,:) = x_2d_double_ptr

 case("Fx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Fx))
   spointer%Fx(:,:) = x_2d_double_ptr

 case("Fy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Fy))
   spointer%Fy(:,:) = x_2d_double_ptr

 case("Sxy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Sxy))
   spointer%Sxy(:,:) = x_2d_double_ptr

 case("Syy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Syy))
   spointer%Syy(:,:) = x_2d_double_ptr

 case("Sxx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Sxx))
   spointer%Sxx(:,:) = x_2d_double_ptr

 case("n")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%n))
   spointer%n(:,:) = x_2d_double_ptr

 case("H")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%H))
   spointer%H(:,:) = x_2d_double_ptr

 case("cgx")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cgx))
   spointer%cgx(:,:,:) = x_3d_double_ptr

 case("cgy")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cgy))
   spointer%cgy(:,:,:) = x_3d_double_ptr

 case("cx")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cx))
   spointer%cx(:,:,:) = x_3d_double_ptr

 case("cy")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cy))
   spointer%cy(:,:,:) = x_3d_double_ptr

 case("ctheta")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ctheta))
   spointer%ctheta(:,:,:) = x_3d_double_ptr

 case("ee")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ee))
   spointer%ee(:,:,:) = x_3d_double_ptr

 case("thet")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%thet))
   spointer%thet(:,:,:) = x_3d_double_ptr

 case("costh")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%costh))
   spointer%costh(:,:,:) = x_3d_double_ptr

 case("sinth")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%sinth))
   spointer%sinth(:,:,:) = x_3d_double_ptr

 case("sigt")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%sigt))
   spointer%sigt(:,:,:) = x_3d_double_ptr

 case("rr")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%rr))
   spointer%rr(:,:,:) = x_3d_double_ptr

 case("cgx_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cgx_s))
   spointer%cgx_s(:,:,:) = x_3d_double_ptr

 case("cgy_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%cgy_s))
   spointer%cgy_s(:,:,:) = x_3d_double_ptr

 case("ctheta_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ctheta_s))
   spointer%ctheta_s(:,:,:) = x_3d_double_ptr

 case("ee_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ee_s))
   spointer%ee_s(:,:,:) = x_3d_double_ptr

 case("thet_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%thet_s))
   spointer%thet_s(:,:,:) = x_3d_double_ptr

 case("costh_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%costh_s))
   spointer%costh_s(:,:,:) = x_3d_double_ptr

 case("sinth_s")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%sinth_s))
   spointer%sinth_s(:,:,:) = x_3d_double_ptr

 case("k")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%k))
   spointer%k(:,:) = x_2d_double_ptr

 case("c")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%c))
   spointer%c(:,:) = x_2d_double_ptr

 case("cg")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cg))
   spointer%cg(:,:) = x_2d_double_ptr

 case("sigm")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%sigm))
   spointer%sigm(:,:) = x_2d_double_ptr

 case("wm")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%wm))
   spointer%wm(:,:) = x_2d_double_ptr

 case("hh")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hh))
   spointer%hh(:,:) = x_2d_double_ptr

 case("zs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zs))
   spointer%zs(:,:) = x_2d_double_ptr

 case("zs0")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zs0))
   spointer%zs0(:,:) = x_2d_double_ptr

 case("tideinpt")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%tideinpt))
   spointer%tideinpt(:) = x_1d_double_ptr

 case("tideinpz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%tideinpz))
   spointer%tideinpz(:,:) = x_2d_double_ptr

 case("windinpt")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%windinpt))
   spointer%windinpt(:) = x_1d_double_ptr

 case("windvelts")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%windvelts))
   spointer%windvelts(:) = x_1d_double_ptr

 case("winddirts")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%winddirts))
   spointer%winddirts(:) = x_1d_double_ptr

 case("windxts")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%windxts))
   spointer%windxts(:) = x_1d_double_ptr

 case("windyts")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%windyts))
   spointer%windyts(:) = x_1d_double_ptr

 case("windsu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%windsu))
   spointer%windsu(:,:) = x_2d_double_ptr

 case("windnv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%windnv))
   spointer%windnv(:,:) = x_2d_double_ptr

 case("dzsdt")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzsdt))
   spointer%dzsdt(:,:) = x_2d_double_ptr

 case("dzsdx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzsdx))
   spointer%dzsdx(:,:) = x_2d_double_ptr

 case("dzsdy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzsdy))
   spointer%dzsdy(:,:) = x_2d_double_ptr

 case("dzbdx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzbdx))
   spointer%dzbdx(:,:) = x_2d_double_ptr

 case("dzbdy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzbdy))
   spointer%dzbdy(:,:) = x_2d_double_ptr

 case("dzbdt")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzbdt))
   spointer%dzbdt(:,:) = x_2d_double_ptr

 case("dzbnow")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzbnow))
   spointer%dzbnow(:,:) = x_2d_double_ptr

 case("uu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%uu))
   spointer%uu(:,:) = x_2d_double_ptr

 case("vu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vu))
   spointer%vu(:,:) = x_2d_double_ptr

 case("uv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%uv))
   spointer%uv(:,:) = x_2d_double_ptr

 case("vv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vv))
   spointer%vv(:,:) = x_2d_double_ptr

 case("qx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%qx))
   spointer%qx(:,:) = x_2d_double_ptr

 case("qy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%qy))
   spointer%qy(:,:) = x_2d_double_ptr

 case("sedero")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%sedero))
   spointer%sedero(:,:) = x_2d_double_ptr

 case("dcbdx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dcbdx))
   spointer%dcbdx(:,:) = x_2d_double_ptr

 case("dcbdy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dcbdy))
   spointer%dcbdy(:,:) = x_2d_double_ptr

 case("dcsdx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dcsdx))
   spointer%dcsdx(:,:) = x_2d_double_ptr

 case("dcsdy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dcsdy))
   spointer%dcsdy(:,:) = x_2d_double_ptr

 case("ui")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ui))
   spointer%ui(:,:) = x_2d_double_ptr

 case("vi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vi))
   spointer%vi(:,:) = x_2d_double_ptr

 case("E")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%E))
   spointer%E(:,:) = x_2d_double_ptr

 case("R")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%R))
   spointer%R(:,:) = x_2d_double_ptr

 case("urms")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%urms))
   spointer%urms(:,:) = x_2d_double_ptr

 case("D")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%D))
   spointer%D(:,:) = x_2d_double_ptr

 case("Qb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Qb))
   spointer%Qb(:,:) = x_2d_double_ptr

 case("ust")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ust))
   spointer%ust(:,:) = x_2d_double_ptr

 case("ueu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ueu))
   spointer%ueu(:,:) = x_2d_double_ptr

 case("vev")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vev))
   spointer%vev(:,:) = x_2d_double_ptr

 case("vmagu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmagu))
   spointer%vmagu(:,:) = x_2d_double_ptr

 case("vmageu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmageu))
   spointer%vmageu(:,:) = x_2d_double_ptr

 case("vmagv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmagv))
   spointer%vmagv(:,:) = x_2d_double_ptr

 case("vmagev")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmagev))
   spointer%vmagev(:,:) = x_2d_double_ptr

 case("u")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%u))
   spointer%u(:,:) = x_2d_double_ptr

 case("v")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%v))
   spointer%v(:,:) = x_2d_double_ptr

 case("ue")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ue))
   spointer%ue(:,:) = x_2d_double_ptr

 case("ve")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ve))
   spointer%ve(:,:) = x_2d_double_ptr

 case("ue_sed")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ue_sed))
   spointer%ue_sed(:,:) = x_2d_double_ptr

 case("ve_sed")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ve_sed))
   spointer%ve_sed(:,:) = x_2d_double_ptr

 case("hold")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hold))
   spointer%hold(:,:) = x_2d_double_ptr

 case("wetu")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%wetu))
   spointer%wetu(:,:) = x_2d_int_ptr

 case("wetv")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%wetv))
   spointer%wetv(:,:) = x_2d_int_ptr

 case("wetz")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%wetz))
   spointer%wetz(:,:) = x_2d_int_ptr

 case("wete")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%wete))
   spointer%wete(:,:) = x_2d_int_ptr

 case("hu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hu))
   spointer%hu(:,:) = x_2d_double_ptr

 case("hv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hv))
   spointer%hv(:,:) = x_2d_double_ptr

 case("hum")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hum))
   spointer%hum(:,:) = x_2d_double_ptr

 case("hvm")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hvm))
   spointer%hvm(:,:) = x_2d_double_ptr

 case("vmag")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmag))
   spointer%vmag(:,:) = x_2d_double_ptr

 case("ccg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ccg))
   spointer%ccg(:,:,:) = x_3d_double_ptr

 case("uwf")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%uwf))
   spointer%uwf(:,:) = x_2d_double_ptr

 case("vwf")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vwf))
   spointer%vwf(:,:) = x_2d_double_ptr

 case("ustr")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ustr))
   spointer%ustr(:,:) = x_2d_double_ptr

 case("usd")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%usd))
   spointer%usd(:,:) = x_2d_double_ptr

 case("bi")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%bi))
   spointer%bi(:) = x_1d_double_ptr

 case("DR")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%DR))
   spointer%DR(:,:) = x_2d_double_ptr

 case("umean")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%umean))
   spointer%umean(:,:) = x_2d_double_ptr

 case("vmean")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmean))
   spointer%vmean(:,:) = x_2d_double_ptr

 case("ur")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ur))
   spointer%ur(:,:) = x_2d_double_ptr

 case("vardx")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%vardx = x_0d_int_ptr

 case("D15")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%D15))
   spointer%D15(:) = x_1d_double_ptr

 case("D50")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%D50))
   spointer%D50(:) = x_1d_double_ptr

 case("D90")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%D90))
   spointer%D90(:) = x_1d_double_ptr

 case("sedcal")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%sedcal))
   spointer%sedcal(:) = x_1d_double_ptr

 case("ucrcal")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%ucrcal))
   spointer%ucrcal(:) = x_1d_double_ptr

 case("Tsg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Tsg))
   spointer%Tsg(:,:,:) = x_3d_double_ptr

 case("Susg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Susg))
   spointer%Susg(:,:,:) = x_3d_double_ptr

 case("Svsg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Svsg))
   spointer%Svsg(:,:,:) = x_3d_double_ptr

 case("Subg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Subg))
   spointer%Subg(:,:,:) = x_3d_double_ptr

 case("Svbg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Svbg))
   spointer%Svbg(:,:,:) = x_3d_double_ptr

 case("ceqbg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ceqbg))
   spointer%ceqbg(:,:,:) = x_3d_double_ptr

 case("ceqsg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ceqsg))
   spointer%ceqsg(:,:,:) = x_3d_double_ptr

 case("ua")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ua))
   spointer%ua(:,:) = x_2d_double_ptr

 case("BR")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%BR))
   spointer%BR(:,:) = x_2d_double_ptr

 case("kb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%kb))
   spointer%kb(:,:) = x_2d_double_ptr

 case("Tbore")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Tbore))
   spointer%Tbore(:,:) = x_2d_double_ptr

 case("dzav")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzav))
   spointer%dzav(:,:) = x_2d_double_ptr

 case("maxzs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%maxzs))
   spointer%maxzs(:,:) = x_2d_double_ptr

 case("minzs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%minzs))
   spointer%minzs(:,:) = x_2d_double_ptr

 case("L1")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%L1))
   spointer%L1(:,:) = x_2d_double_ptr

 case("Sk")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Sk))
   spointer%Sk(:,:) = x_2d_double_ptr

 case("As")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%As))
   spointer%As(:,:) = x_2d_double_ptr

 case("gwhead")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwhead))
   spointer%gwhead(:,:) = x_2d_double_ptr

 case("gwheadb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwheadb))
   spointer%gwheadb(:,:) = x_2d_double_ptr

 case("gwlevel")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwlevel))
   spointer%gwlevel(:,:) = x_2d_double_ptr

 case("gwheight")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwheight))
   spointer%gwheight(:,:) = x_2d_double_ptr

 case("gwbottom")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwbottom))
   spointer%gwbottom(:,:) = x_2d_double_ptr

 case("gwu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwu))
   spointer%gwu(:,:) = x_2d_double_ptr

 case("gwv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwv))
   spointer%gwv(:,:) = x_2d_double_ptr

 case("gwqx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwqx))
   spointer%gwqx(:,:) = x_2d_double_ptr

 case("gwqy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwqy))
   spointer%gwqy(:,:) = x_2d_double_ptr

 case("gww")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gww))
   spointer%gww(:,:) = x_2d_double_ptr

 case("gwcurv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gwcurv))
   spointer%gwcurv(:,:) = x_2d_double_ptr

 case("dinfil")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dinfil))
   spointer%dinfil(:,:) = x_2d_double_ptr

 case("infil")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%infil))
   spointer%infil(:,:) = x_2d_double_ptr

 case("gw0back")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%gw0back))
   spointer%gw0back(:,:) = x_2d_double_ptr

 case("Kx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Kx))
   spointer%Kx(:,:) = x_2d_double_ptr

 case("Ky")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Ky))
   spointer%Ky(:,:) = x_2d_double_ptr

 case("Kz")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Kz))
   spointer%Kz(:,:) = x_2d_double_ptr

 case("Kzinf")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Kzinf))
   spointer%Kzinf(:,:) = x_2d_double_ptr

 case("kturb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%kturb))
   spointer%kturb(:,:) = x_2d_double_ptr

 case("ero")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ero))
   spointer%ero(:,:,:) = x_3d_double_ptr

 case("depo_im")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%depo_im))
   spointer%depo_im(:,:,:) = x_3d_double_ptr

 case("depo_ex")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%depo_ex))
   spointer%depo_ex(:,:,:) = x_3d_double_ptr

 case("nd")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%nd))
   spointer%nd(:,:) = x_2d_int_ptr

 case("pbbed")
   call c_f_pointer(xptr, x_4d_double_ptr, shape(spointer%pbbed))
   spointer%pbbed(:,:,:,:) = x_4d_double_ptr

 case("dzbed")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%dzbed))
   spointer%dzbed(:,:,:) = x_3d_double_ptr

 case("z0bed")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%z0bed))
   spointer%z0bed(:,:) = x_2d_double_ptr

 case("ureps")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ureps))
   spointer%ureps(:,:) = x_2d_double_ptr

 case("vreps")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vreps))
   spointer%vreps(:,:) = x_2d_double_ptr

 case("urepb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%urepb))
   spointer%urepb(:,:) = x_2d_double_ptr

 case("vrepb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vrepb))
   spointer%vrepb(:,:) = x_2d_double_ptr

 case("umwci")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%umwci))
   spointer%umwci(:,:) = x_2d_double_ptr

 case("vmwci")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vmwci))
   spointer%vmwci(:,:) = x_2d_double_ptr

 case("rolthick")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%rolthick))
   spointer%rolthick(:,:) = x_2d_double_ptr

 case("zswci")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zswci))
   spointer%zswci(:,:) = x_2d_double_ptr

 case("pres")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%pres))
   spointer%pres(:,:) = x_2d_double_ptr

 case("dU")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dU))
   spointer%dU(:,:) = x_2d_double_ptr

 case("dV")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dV))
   spointer%dV(:,:) = x_2d_double_ptr

 case("wb")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%wb))
   spointer%wb(:,:) = x_2d_double_ptr

 case("ws")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ws))
   spointer%ws(:,:) = x_2d_double_ptr

 case("wscrit")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%wscrit))
   spointer%wscrit(:,:) = x_2d_double_ptr

 case("bedfriccoef")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%bedfriccoef))
   spointer%bedfriccoef(:,:) = x_2d_double_ptr

 case("taubx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%taubx))
   spointer%taubx(:,:) = x_2d_double_ptr

 case("tauby")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%tauby))
   spointer%tauby(:,:) = x_2d_double_ptr

 case("Df")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Df))
   spointer%Df(:,:) = x_2d_double_ptr

 case("Dp")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Dp))
   spointer%Dp(:,:) = x_2d_double_ptr

 case("Sutot")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Sutot))
   spointer%Sutot(:,:) = x_2d_double_ptr

 case("Svtot")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Svtot))
   spointer%Svtot(:,:) = x_2d_double_ptr

 case("cctot")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cctot))
   spointer%cctot(:,:) = x_2d_double_ptr

 case("wi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%wi))
   spointer%wi(:,:) = x_2d_double_ptr

 case("dUi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dUi))
   spointer%dUi(:,:) = x_2d_double_ptr

 case("zi")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zi))
   spointer%zi(:,:) = x_2d_double_ptr

 case("nuh")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%nuh))
   spointer%nuh(:,:) = x_2d_double_ptr

 case("cf")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cf))
   spointer%cf(:,:) = x_2d_double_ptr

 case("cfu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cfu))
   spointer%cfu(:,:) = x_2d_double_ptr

 case("cfv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cfv))
   spointer%cfv(:,:) = x_2d_double_ptr

 case("D50top")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%D50top))
   spointer%D50top(:,:) = x_2d_double_ptr

 case("D90top")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%D90top))
   spointer%D90top(:,:) = x_2d_double_ptr

 case("structdepth")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%structdepth))
   spointer%structdepth(:,:) = x_2d_double_ptr

 case("zs0fac")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%zs0fac))
   spointer%zs0fac(:,:,:) = x_3d_double_ptr

 case("tdisch")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%tdisch))
   spointer%tdisch(:) = x_1d_double_ptr

 case("pdisch")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%pdisch))
   spointer%pdisch(:,:) = x_2d_double_ptr

 case("pntdisch")
   call c_f_pointer(xptr, x_1d_int_ptr, shape(spointer%pntdisch))
   spointer%pntdisch(:) = x_1d_int_ptr

 case("qdisch")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%qdisch))
   spointer%qdisch(:,:) = x_2d_double_ptr

 case("idrift")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%idrift))
   spointer%idrift(:) = x_1d_double_ptr

 case("jdrift")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%jdrift))
   spointer%jdrift(:) = x_1d_double_ptr

 case("tdriftb")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%tdriftb))
   spointer%tdriftb(:) = x_1d_double_ptr

 case("tdrifte")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%tdrifte))
   spointer%tdrifte(:) = x_1d_double_ptr

 case("runup")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%runup))
   spointer%runup(:) = x_1d_double_ptr

 case("Hrunup")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%Hrunup))
   spointer%Hrunup(:) = x_1d_double_ptr

 case("xHrunup")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%xHrunup))
   spointer%xHrunup(:) = x_1d_double_ptr

 case("istruct")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%istruct))
   spointer%istruct(:) = x_1d_double_ptr

 case("iwl")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%iwl))
   spointer%iwl(:) = x_1d_double_ptr

 case("strucslope")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%strucslope))
   spointer%strucslope(:) = x_1d_double_ptr

 case("Dc")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Dc))
   spointer%Dc(:,:) = x_2d_double_ptr

 case("ph")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ph))
   spointer%ph(:,:) = x_2d_double_ptr

 case("newstatbc")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%newstatbc = x_0d_int_ptr

 case("dobs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dobs))
   spointer%dobs(:,:) = x_2d_double_ptr

 case("sig2prior")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%sig2prior))
   spointer%sig2prior(:,:) = x_2d_double_ptr

 case("zbobs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zbobs))
   spointer%zbobs(:,:) = x_2d_double_ptr

 case("shobs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%shobs))
   spointer%shobs(:,:) = x_2d_double_ptr

 case("bwalpha")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%bwalpha))
   spointer%bwalpha(:,:) = x_2d_double_ptr

 case("dcmdo")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dcmdo))
   spointer%dcmdo(:,:) = x_2d_double_ptr

 case("dassim")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dassim))
   spointer%dassim(:,:) = x_2d_double_ptr

 case("cobs")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%cobs))
   spointer%cobs(:,:) = x_2d_double_ptr

 case("shipxCG")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipxCG))
   spointer%shipxCG(:) = x_1d_double_ptr

 case("shipyCG")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipyCG))
   spointer%shipyCG(:) = x_1d_double_ptr

 case("shipzCG")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipzCG))
   spointer%shipzCG(:) = x_1d_double_ptr

 case("shipFx")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipFx))
   spointer%shipFx(:) = x_1d_double_ptr

 case("shipFy")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipFy))
   spointer%shipFy(:) = x_1d_double_ptr

 case("shipFz")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipFz))
   spointer%shipFz(:) = x_1d_double_ptr

 case("shipMx")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipMx))
   spointer%shipMx(:) = x_1d_double_ptr

 case("shipMy")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipMy))
   spointer%shipMy(:) = x_1d_double_ptr

 case("shipMz")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipMz))
   spointer%shipMz(:) = x_1d_double_ptr

 case("shipphi")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipphi))
   spointer%shipphi(:) = x_1d_double_ptr

 case("shipchi")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shipchi))
   spointer%shipchi(:) = x_1d_double_ptr

 case("shippsi")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%shippsi))
   spointer%shippsi(:) = x_1d_double_ptr

 case("vegtype")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%vegtype))
   spointer%vegtype(:,:) = x_2d_int_ptr

 case("nsecveg")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%nsecveg))
   spointer%nsecveg(:,:) = x_2d_int_ptr

 case("nsecvegmax")
   call c_f_pointer(xptr, x_0d_int_ptr) ! skip shape in case of rank 0
   spointer%nsecvegmax = x_0d_int_ptr

 case("Cdveg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Cdveg))
   spointer%Cdveg(:,:,:) = x_3d_double_ptr

 case("ahveg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ahveg))
   spointer%ahveg(:,:,:) = x_3d_double_ptr

 case("bveg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%bveg))
   spointer%bveg(:,:,:) = x_3d_double_ptr

 case("Nveg")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%Nveg))
   spointer%Nveg(:,:,:) = x_3d_double_ptr

 case("Dveg")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Dveg))
   spointer%Dveg(:,:) = x_2d_double_ptr

 case("Fvegu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Fvegu))
   spointer%Fvegu(:,:) = x_2d_double_ptr

 case("Fvegv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%Fvegv))
   spointer%Fvegv(:,:) = x_2d_double_ptr

 case("ududx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ududx))
   spointer%ududx(:,:) = x_2d_double_ptr

 case("vdvdy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vdvdy))
   spointer%vdvdy(:,:) = x_2d_double_ptr

 case("udvdx")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%udvdx))
   spointer%udvdx(:,:) = x_2d_double_ptr

 case("vdudy")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vdudy))
   spointer%vdudy(:,:) = x_2d_double_ptr

 case("viscu")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%viscu))
   spointer%viscu(:,:) = x_2d_double_ptr

 case("viscv")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%viscv))
   spointer%viscv(:,:) = x_2d_double_ptr

 case("setbathy")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%setbathy))
   spointer%setbathy(:,:,:) = x_3d_double_ptr

 case("tsetbathy")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%tsetbathy))
   spointer%tsetbathy(:) = x_1d_double_ptr

 case("breaking")
   call c_f_pointer(xptr, x_2d_int_ptr, shape(spointer%breaking))
   spointer%breaking(:,:) = x_2d_int_ptr

 case("fw")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%fw))
   spointer%fw(:,:) = x_2d_double_ptr

 case("sigz")
   call c_f_pointer(xptr, x_1d_double_ptr, shape(spointer%sigz))
   spointer%sigz(:) = x_1d_double_ptr

 case("uz")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%uz))
   spointer%uz(:,:,:) = x_3d_double_ptr

 case("vz")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%vz))
   spointer%vz(:,:,:) = x_3d_double_ptr

 case("ustz")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ustz))
   spointer%ustz(:,:,:) = x_3d_double_ptr

 case("nutz")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%nutz))
   spointer%nutz(:,:,:) = x_3d_double_ptr

 case("dzs0dn")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%dzs0dn))
   spointer%dzs0dn(:,:) = x_2d_double_ptr

 case("ccz")
   call c_f_pointer(xptr, x_3d_double_ptr, shape(spointer%ccz))
   spointer%ccz(:,:,:) = x_3d_double_ptr

 case("refA")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%refA))
   spointer%refA(:,:) = x_2d_double_ptr

 case("ca")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ca))
   spointer%ca(:,:) = x_2d_double_ptr

 case("zs1")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%zs1))
   spointer%zs1(:,:) = x_2d_double_ptr

 case("taubx_add")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%taubx_add))
   spointer%taubx_add(:,:) = x_2d_double_ptr

 case("tauby_add")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%tauby_add))
   spointer%tauby_add(:,:) = x_2d_double_ptr

 case("hhw")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hhw))
   spointer%hhw(:,:) = x_2d_double_ptr

 case("hhws")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hhws))
   spointer%hhws(:,:) = x_2d_double_ptr

 case("uws")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%uws))
   spointer%uws(:,:) = x_2d_double_ptr

 case("vws")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vws))
   spointer%vws(:,:) = x_2d_double_ptr

 case("hhwcins")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%hhwcins))
   spointer%hhwcins(:,:) = x_2d_double_ptr

 case("uwcins")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%uwcins))
   spointer%uwcins(:,:) = x_2d_double_ptr

 case("vwcins")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vwcins))
   spointer%vwcins(:,:) = x_2d_double_ptr

 case("ucan")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%ucan))
   spointer%ucan(:,:) = x_2d_double_ptr

 case("vcan")
   call c_f_pointer(xptr, x_2d_double_ptr, shape(spointer%vcan))
   spointer%vcan(:,:) = x_2d_double_ptr



! vim: filetype=fortran