!----- AGPL --------------------------------------------------------------------
!                                                                               
!  Copyright (C)  Stichting Deltares, 2015-2022.                                
!                                                                               
!  This file is part of Delft3D (D-Flow Flexible Mesh component).               
!                                                                               
!  Delft3D is free software: you can redistribute it and/or modify              
!  it under the terms of the GNU Affero General Public License as               
!  published by the Free Software Foundation version 3.                         
!                                                                               
!  Delft3D  is distributed in the hope that it will be useful,                  
!  but WITHOUT ANY WARRANTY; without even the implied warranty of               
!  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the                
!  GNU Affero General Public License for more details.                          
!                                                                               
!  You should have received a copy of the GNU Affero General Public License     
!  along with Delft3D.  If not, see <http://www.gnu.org/licenses/>.             
!                                                                               
!  contact: delft3d.support@deltares.nl                                         
!  Stichting Deltares                                                           
!  P.O. Box 177                                                                 
!  2600 MH Delft, The Netherlands                                               
!                                                                               
!  All indications and logos of, and references to, "Delft3D",                  
!  "D-Flow Flexible Mesh" and "Deltares" are registered trademarks of Stichting 
!  Deltares, and remain the property of Stichting Deltares. All rights reserved.
!                                                                               
!-------------------------------------------------------------------------------

  select case(var_name)
  case("DFM_COMM_DFMWORLD")
     call c_f_pointer(xptr, x_0d_int_ptr)
     DFM_COMM_DFMWORLD = x_0d_int_ptr
  case("iglobal_s")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(iglobal_s))
     iglobal_s(:) = x_1d_int_ptr
  case("hwav")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(hwav))
     hwav(:) = x_1d_double_ptr
  case("twav")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(twav))
     twav(:) = x_1d_double_ptr
  case("Uorb")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(Uorb))
     Uorb(:) = x_1d_double_ptr
  case("bodsed")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(stmpar%morlyr%state%bodsed))
     stmpar%morlyr%state%bodsed(:,:) = x_2d_double_ptr
  case("dpsed")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stmpar%morlyr%state%dpsed))
     stmpar%morlyr%state%dpsed(:) = x_1d_double_ptr
  case("msed")
     call c_f_pointer(xptr, x_3d_double_ptr, shape(stmpar%morlyr%state%msed))
     stmpar%morlyr%state%msed(:,:,:) = x_3d_double_ptr
  case("thlyr")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(stmpar%morlyr%state%thlyr))
     stmpar%morlyr%state%thlyr(:,:) = x_2d_double_ptr
  case("rsedeq")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%rsedeq))
     sedtra%rsedeq(:,:) = x_2d_double_ptr
  case("sbcx")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sbcx))
     sedtra%sbcx(:,:) = x_2d_double_ptr
  case("sbcy")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sbcy))
     sedtra%sbcy(:,:) = x_2d_double_ptr
  case("sbwx")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sbwx))
     sedtra%sbwx(:,:) = x_2d_double_ptr
  case("sbwy")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sbwy))
     sedtra%sbwy(:,:) = x_2d_double_ptr
  case("sscx")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sscx))
     sedtra%sscx(:,:) = x_2d_double_ptr
  case("sscy")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sscy))
     sedtra%sscy(:,:) = x_2d_double_ptr
  case("sswx")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sswx))
     sedtra%sswx(:,:) = x_2d_double_ptr
  case("sswy")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(sedtra%sswy))
     sedtra%sswy(:,:) = x_2d_double_ptr
  case("taucr")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stmpar%sedpar%taucr))
     stmpar%sedpar%taucr(:) = x_1d_double_ptr
  case("tetacr")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stmpar%sedpar%tetacr))
     stmpar%sedpar%tetacr(:) = x_1d_double_ptr
  case("infilt")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(infilt))
     infilt(:) = x_1d_double_ptr
  case("infiltcap")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(infiltcap))
     infiltcap(:) = x_1d_double_ptr
  case("shx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shx))
     shx(:) = x_1d_double_ptr
  case("shy")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shy))
     shy(:) = x_1d_double_ptr
  case("shi")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shi))
     shi(:) = x_1d_double_ptr
  case("zsp")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zsp))
     zsp(:) = x_1d_double_ptr
  case("zsp0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zsp0))
     zsp0(:) = x_1d_double_ptr
  case("zspc")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zspc))
     zspc(:) = x_1d_double_ptr
  case("zspc0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zspc0))
     zspc0(:) = x_1d_double_ptr
  case("v0ship")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(v0ship))
     v0ship(:) = x_1d_double_ptr
  case("v1ship")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(v1ship))
     v1ship(:) = x_1d_double_ptr
  case("qinship")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qinship))
     qinship(:) = x_1d_double_ptr
  case("vicushp")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vicushp))
     vicushp(:) = x_1d_double_ptr
  case("shL")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shL))
     shL(:) = x_1d_double_ptr
  case("shB")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shB))
     shB(:) = x_1d_double_ptr
  case("shd")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(shd))
     shd(:) = x_1d_double_ptr
  case("stuw")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stuw))
     stuw(:) = x_1d_double_ptr
  case("fstuw")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(fstuw))
     fstuw(:) = x_1d_double_ptr
  case("stuwmx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stuwmx))
     stuwmx(:) = x_1d_double_ptr
  case("roer")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(roer))
     roer(:) = x_1d_double_ptr
  case("froer")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(froer))
     froer(:) = x_1d_double_ptr
  case("roermx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(roermx))
     roermx(:) = x_1d_double_ptr
  case("wx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(wx))
     wx(:) = x_1d_double_ptr
  case("wy")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(wy))
     wy(:) = x_1d_double_ptr
  case("rain")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(rain))
     rain(:) = x_1d_double_ptr
  case("evap")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(evap))
     evap(:) = x_1d_double_ptr
  case("numlatsg")
     call c_f_pointer(xptr, x_0d_int_ptr)
     numlatsg = x_0d_int_ptr
  case("qplat")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qplat))
     qplat(:) = x_1d_double_ptr
  case("qqlat")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qqlat))
     qqlat(:) = x_1d_double_ptr
  case("balat")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(balat))
     balat(:) = x_1d_double_ptr
  case("n1latsg")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(n1latsg))
     n1latsg(:) = x_1d_int_ptr
  case("n2latsg")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(n2latsg))
     n2latsg(:) = x_1d_int_ptr
  case("nnlat")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(nnlat))
     nnlat(:) = x_1d_int_ptr
  case("kclat")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(kclat))
     kclat(:) = x_1d_int_ptr
  case("qext")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qext))
     qext(:) = x_1d_double_ptr
  case("qextreal")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qextreal))
     qextreal(:) = x_1d_double_ptr
  case("vextcum")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vextcum))
     vextcum(:) = x_1d_double_ptr
  case("rainuni")
     call c_f_pointer(xptr, x_0d_double_ptr)
     rainuni = x_0d_double_ptr
  case("Qsunmap")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(Qsunmap))
     Qsunmap(:) = x_1d_double_ptr
  case("zbndz")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zbndz))
     zbndz(:) = x_1d_double_ptr
  case("zbndu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zbndu))
     zbndu(:) = x_1d_double_ptr
  case("zbndq")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zbndq))
     zbndq(:) = x_1d_double_ptr
  case("turkin1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(turkin1))
     turkin1(:) = x_1d_double_ptr
  case("vicwws")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vicwws))
     vicwws(:) = x_1d_double_ptr
  case("rnveg")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(rnveg))
     rnveg(:) = x_1d_double_ptr
  case("diaveg")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(diaveg))
     diaveg(:) = x_1d_double_ptr
  case("cfuveg")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(cfuveg))
     cfuveg(:) = x_1d_double_ptr
  case("alfaveg")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(alfaveg))
     alfaveg(:) = x_1d_double_ptr
  case("stemdens")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stemdens))
     stemdens(:) = x_1d_double_ptr
  case("stemdiam")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stemdiam))
     stemdiam(:) = x_1d_double_ptr
  case("stemheight")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(stemheight))
     stemheight(:) = x_1d_double_ptr
  case("Cdvegsp")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(Cdvegsp))
     Cdvegsp(:) = x_1d_double_ptr
  case("zws")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zws))
     zws(:) = x_1d_double_ptr
  case("kbot")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(kbot))
     kbot(:) = x_1d_int_ptr
  case("ktop")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(ktop))
     ktop(:) = x_1d_int_ptr
  case("Lbot")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(Lbot))
     Lbot(:) = x_1d_int_ptr
  case("Ltop")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(Ltop))
     Ltop(:) = x_1d_int_ptr
  case("s0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(s0))
     s0(:) = x_1d_double_ptr
  case("s1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(s1))
     s1(:) = x_1d_double_ptr
  case("a0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(a0))
     a0(:) = x_1d_double_ptr
  case("a1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(a1))
     a1(:) = x_1d_double_ptr
  case("vol0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vol0))
     vol0(:) = x_1d_double_ptr
  case("vol1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vol1))
     vol1(:) = x_1d_double_ptr
  case("vol1_f")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vol1_f))
     vol1_f(:) = x_1d_double_ptr
  case("hs")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(hs))
     hs(:) = x_1d_double_ptr
  case("ucx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ucx))
     ucx(:) = x_1d_double_ptr
  case("ucy")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ucy))
     ucy(:) = x_1d_double_ptr
  case("ucz")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ucz))
     ucz(:) = x_1d_double_ptr
  case("ucmag")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ucmag))
     ucmag(:) = x_1d_double_ptr
  case("sa0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(sa0))
     sa0(:) = x_1d_double_ptr
  case("sa1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(sa1))
     sa1(:) = x_1d_double_ptr
  case("satop")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(satop))
     satop(:) = x_1d_double_ptr
  case("sabot")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(sabot))
     sabot(:) = x_1d_double_ptr
  case("tem0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(tem0))
     tem0(:) = x_1d_double_ptr
  case("tem1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(tem1))
     tem1(:) = x_1d_double_ptr
  case("u1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(u1))
     u1(:) = x_1d_double_ptr
  case("u_to_umain")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(u_to_umain))
     u_to_umain(:) = x_1d_double_ptr
  case("q1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(q1))
     q1(:) = x_1d_double_ptr
  case("q1_main")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(q1_main))
     q1_main(:) = x_1d_double_ptr
  case("au")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(au))
     au(:) = x_1d_double_ptr
  case("hu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(hu))
     hu(:) = x_1d_double_ptr
  case("frcu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(frcu))
     frcu(:) = x_1d_double_ptr
  case("viusp")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(viusp))
     viusp(:) = x_1d_double_ptr
  case("diusp")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(diusp))
     diusp(:) = x_1d_double_ptr
  case("taus")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(taus))
     taus(:) = x_1d_double_ptr
  case("ndxi")
     call c_f_pointer(xptr, x_0d_int_ptr)
     ndxi = x_0d_int_ptr
  case("ndx1db")
     call c_f_pointer(xptr, x_0d_int_ptr)
     ndx1db = x_0d_int_ptr
  case("kfs")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(kfs))
     kfs(:) = x_1d_int_ptr
  case("bare")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(bare))
     bare(:) = x_1d_double_ptr
  case("ba_mor")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ba_mor))
     ba_mor(:) = x_1d_double_ptr
  case("bl")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(bl))
     bl(:) = x_1d_double_ptr
  case("lnx1D")
     call c_f_pointer(xptr, x_0d_int_ptr)
     lnx1D = x_0d_int_ptr
  case("lnxi")
     call c_f_pointer(xptr, x_0d_int_ptr)
     lnxi = x_0d_int_ptr
  case("lnx1Db")
     call c_f_pointer(xptr, x_0d_int_ptr)
     lnx1Db = x_0d_int_ptr
  case("lnx")
     call c_f_pointer(xptr, x_0d_int_ptr)
     lnx = x_0d_int_ptr
  case("ln")
     call c_f_pointer(xptr, x_2d_int_ptr, shape(ln))
     ln(:,:) = x_2d_int_ptr
  case("LLkkk")
     call c_f_pointer(xptr, x_2d_int_ptr, shape(LLkkk))
     LLkkk(:,:) = x_2d_int_ptr
  case("lncn")
     call c_f_pointer(xptr, x_2d_int_ptr, shape(lncn))
     lncn(:,:) = x_2d_int_ptr
  case("iadv")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(iadv))
     iadv(:) = x_1d_int_ptr
  case("dx")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(dx))
     dx(:) = x_1d_double_ptr
  case("wu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(wu))
     wu(:) = x_1d_double_ptr
  case("wu_mor")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(wu_mor))
     wu_mor(:) = x_1d_double_ptr
  case("wu1D2D")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(wu1D2D))
     wu1D2D(:) = x_1d_double_ptr
  case("hh1D2D")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(hh1D2D))
     hh1D2D(:) = x_1d_double_ptr
  case("bob")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(bob))
     bob(:,:) = x_2d_double_ptr
  case("bob0")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(bob0))
     bob0(:,:) = x_2d_double_ptr
  case("xu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(xu))
     xu(:) = x_1d_double_ptr
  case("yu")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(yu))
     yu(:) = x_1d_double_ptr
  case("vort")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(vort))
     vort(:) = x_1d_double_ptr
  case("xzw")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(xzw))
     xzw(:) = x_1d_double_ptr
  case("yzw")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(yzw))
     yzw(:) = x_1d_double_ptr
  case("xk")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(xk))
     xk(:) = x_1d_double_ptr
  case("yk")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(yk))
     yk(:) = x_1d_double_ptr
  case("zk")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zk))
     zk(:) = x_1d_double_ptr
  case("kn")
     call c_f_pointer(xptr, x_2d_int_ptr, shape(kn))
     kn(:,:) = x_2d_int_ptr
  case("numk")
     call c_f_pointer(xptr, x_0d_int_ptr)
     numk = x_0d_int_ptr
  case("zbnd1d2d1")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zbnd1d2d1))
     zbnd1d2d1(:) = x_1d_double_ptr
  case("zbnd1d2d0")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zbnd1d2d0))
     zbnd1d2d0(:) = x_1d_double_ptr
  case("zcrest1d2d")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(zcrest1d2d))
     zcrest1d2d(:) = x_1d_double_ptr
  case("edgenumbers1d2d")
     call c_f_pointer(xptr, x_1d_int_ptr, shape(edgenumbers1d2d))
     edgenumbers1d2d(:) = x_1d_int_ptr
  case("kbnd1d2d")
     call c_f_pointer(xptr, x_2d_int_ptr, shape(kbnd1d2d))
     kbnd1d2d(:,:) = x_2d_int_ptr
  case("width_1d")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(width_1d))
     width_1d(:) = x_1d_double_ptr
  case("qzeta_1d2d")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qzeta_1d2d))
     qzeta_1d2d(:) = x_1d_double_ptr
  case("qlat_1d2d")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qlat_1d2d))
     qlat_1d2d(:) = x_1d_double_ptr
  case("qtotal_1d2d")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(qtotal_1d2d))
     qtotal_1d2d(:) = x_1d_double_ptr
  case("is_sumvalsnd")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(is_sumvalsnd))
     is_sumvalsnd(:,:) = x_2d_double_ptr
  case("is_maxvalsnd")
     call c_f_pointer(xptr, x_2d_double_ptr, shape(is_maxvalsnd))
     is_maxvalsnd(:,:) = x_2d_double_ptr
  case("is_dtint")
     call c_f_pointer(xptr, x_0d_double_ptr)
     is_dtint = x_0d_double_ptr
  case("ndx2d")
     call c_f_pointer(xptr, x_0d_int_ptr)
     ndx2d = x_0d_int_ptr
  case("ndx")
     call c_f_pointer(xptr, x_0d_int_ptr)
     ndx = x_0d_int_ptr
  case("xz")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(xz))
     xz(:) = x_1d_double_ptr
  case("yz")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(yz))
     yz(:) = x_1d_double_ptr
  case("ba")
     call c_f_pointer(xptr, x_1d_double_ptr, shape(ba))
     ba(:) = x_1d_double_ptr
  end select