deall_telemac3d.f

Go to the documentation of this file.

!                   **************************
                    SUBROUTINE deall_telemac3d
!                   **************************
!
!
!***********************************************************************
! TELEMAC3D   V7P2                                   02/08/2017
!***********************************************************************
!
!brief    Memory allocation of structures, aliases, blocks...
!
!history  J-M HERVOUET (LNHE)
!+        02/04/2009
!+
!+   T3D_FILES(T3DGEO)%LU REPLACES NGEO
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE tpxo, ONLY : dealloc_tpxo
      USE tel4del, ONLY : dealloc_tel4del
      USE meteo_telemac, ONLY : dealloc_meteo
      USE gotm_coupling, ONLY : dealloc_gotm
      USE oilspill
      USE declarations_telemac3d
      USE declarations_gaia, ONLY : num_isusp_icla,fludp,fludpt,fluer
      USE declarations_waqtel, ONLY : waqprocess
#if defined HAVE_AED2
      USE t3d_aed2, ONLY : clean_aed2_models
#endif
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER :: I
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!

!
!-----------------------------------------------------------------------
!

      CALL bief_deallobj(varsor )
      CALL bief_deallobj(varcl )
      CALL bief_deallobj(varso3)
      CALL bief_deallobj(bl_fc )
      CALL bief_deallobj(bl_fn )
      CALL bief_deallobj(bl_s0f)
      CALL bief_deallobj(bl_fsc)
      CALL bief_deallobj(bl_bol)
      CALL bief_deallobj(bl_borl)
      CALL bief_deallobj(fn3d  )
      CALL bief_deallobj(fc3d  )
      CALL bief_deallobj(prive  )
      CALL bief_deallobj(prive2d)
      CALL bief_deallobj(viscnu)

!
!     Deallocate the mesh structure
      CALL dealmesh(mesh2d)
      CALL dealmesh(mesh3d)
      NULLIFY(ikle2)
      NULLIFY(ikle3)
      NULLIFY(x)
      NULLIFY(y)
      NULLIFY(z)
      NULLIFY(nelem2)
      NULLIFY(nelem3)
      NULLIFY(nelmax2)
      NULLIFY(nelmax3)
      NULLIFY(nptfr2)
      NULLIFY(nptfr3)
      NULLIFY(nptfrx2)
      NULLIFY(nptfrx3)
      NULLIFY(typelm2)
      NULLIFY(typelm3)
      NULLIFY(npoin2)
      NULLIFY(npoin3)
      NULLIFY(npmax2)
      NULLIFY(npmax3)
      NULLIFY(mxptvs2)
      NULLIFY(mxptvs3)
      NULLIFY(lv)


      CALL bief_deallobj(zprop )
      CALL bief_deallobj(zt    )
      CALL bief_deallobj(zn    )

      CALL bief_deallobj(gradzf)
      CALL bief_deallobj(gradzs)
      CALL bief_deallobj(gradzn)
      CALL bief_deallobj(dssudt)

      CALL bief_deallobj(zstar )
      CALL bief_deallobj(zplane)
      CALL bief_deallobj(zchar )
      CALL bief_deallobj(transf_plane)

      CALL bief_deallobj(lihbor)
      CALL bief_deallobj(numliq)
      CALL bief_deallobj(numliq_elm)
      CALL bief_deallobj(boundary_colour)
      CALL bief_deallobj(limpro)

      CALL bief_deallobj(nliqbed)
      CALL bief_deallobj(bedflu)

      CALL bief_deallobj(mask  )
      CALL bief_deallobj(mask_3d)
      CALL bief_deallobj(masktr)
      CALL bief_deallobj(maskel)
      CALL bief_deallobj(maskpt)
      CALL bief_deallobj(maskbr)

      CALL bief_deallobj(un)
      CALL bief_deallobj(vn)
      CALL bief_deallobj(uc)
      CALL bief_deallobj(vc)
      CALL bief_deallobj(ud)
      CALL bief_deallobj(vd)
      CALL bief_deallobj(u)
      CALL bief_deallobj(v)
      CALL bief_deallobj(s0u)
      CALL bief_deallobj(s0v)
      CALL bief_deallobj(s1u)
      CALL bief_deallobj(s1v)
      CALL bief_deallobj(smu)
      CALL bief_deallobj(smv)
      CALL bief_deallobj(uborf   )
      CALL bief_deallobj(vborf   )
      CALL bief_deallobj(uborl   )
      CALL bief_deallobj(vborl   )
      CALL bief_deallobj(ubors   )
      CALL bief_deallobj(vbors   )
      CALL bief_deallobj(auborf  )
      CALL bief_deallobj(avborf  )
      CALL bief_deallobj(auborl  )
      CALL bief_deallobj(uetcal  )
      CALL bief_deallobj(avborl  )
      CALL bief_deallobj(aubors  )
      CALL bief_deallobj(avbors  )
      CALL bief_deallobj(buborf  )
      CALL bief_deallobj(buborl  )
      CALL bief_deallobj(bubors  )
      CALL bief_deallobj(bvborf  )
      CALL bief_deallobj(bvborl  )
      CALL bief_deallobj(bvbors  )
      CALL bief_deallobj(uborsave)
      CALL bief_deallobj(vborsave)
      CALL bief_deallobj(wborsave)

      CALL bief_deallobj(hbtide )
      CALL bief_deallobj(ubtide )
      CALL bief_deallobj(vbtide )
      CALL bief_deallobj(numtide)

      CALL bief_deallobj(uetcar  )
      CALL bief_deallobj(ipbot )

      CALL bief_deallobj(liubof)
      CALL bief_deallobj(liubol)
      CALL bief_deallobj(liubos)
      CALL bief_deallobj(livbof)
      CALL bief_deallobj(livbol)
      CALL bief_deallobj(livbos)

      IF(thomfr) THEN
        CALL bief_deallobj(lihbor_user)
        CALL bief_deallobj(liubol_user)
        CALL bief_deallobj(livbol_user)
      ENDIF

      CALL bief_deallobj(wn)
      CALL bief_deallobj(wc)
      CALL bief_deallobj(wd)
      CALL bief_deallobj(wconv )
      CALL bief_deallobj(w)
      CALL bief_deallobj(s0w)
      CALL bief_deallobj(s1w)
      CALL bief_deallobj(wborf   )
      CALL bief_deallobj(wborl   )
      CALL bief_deallobj(wbors   )
      CALL bief_deallobj(awborf  )
      CALL bief_deallobj(awborl  )
      CALL bief_deallobj(awbors  )
      CALL bief_deallobj(bwborf  )
      CALL bief_deallobj(bwborl  )
      CALL bief_deallobj(bwbors  )
      CALL bief_deallobj(liwbof)
      CALL bief_deallobj(liwbol)
      CALL bief_deallobj(liwbos)
      CALL bief_deallobj(ws    )
      CALL bief_deallobj(wel   )

      CALL bief_deallobj(dp    )
      CALL bief_deallobj(dpn   )
      CALL bief_deallobj(ph    )
      CALL bief_deallobj(pborf )
      CALL bief_deallobj(pborl )
      CALL bief_deallobj(pbors )
      CALL bief_deallobj(lipbof)
      CALL bief_deallobj(lipbol)
      CALL bief_deallobj(lipbos)
!     block
      CALL bief_deallobj(trn   )
      CALL bief_deallobj(tac   )
      CALL bief_deallobj(ta    )
      CALL bief_deallobj(s0ta  )
      CALL bief_deallobj(s1ta  )
      CALL bief_deallobj(taborf)
      CALL bief_deallobj(taborl)
      CALL bief_deallobj(tabors)
      CALL bief_deallobj(atabof)
      CALL bief_deallobj(atabol)
      CALL bief_deallobj(atabos)
      CALL bief_deallobj(btabof)
      CALL bief_deallobj(btabol)
      CALL bief_deallobj(btabos)
      CALL bief_deallobj(litabf)
      CALL bief_deallobj(litabl)
      CALL bief_deallobj(litabs)
      CALL bief_deallobj(trborsave)
      CALL bief_deallobj(ta_sce )

      CALL bief_deallobj(akn   )
      CALL bief_deallobj(akc   )
      CALL bief_deallobj(ak    )
      CALL bief_deallobj(s0ak  )
      CALL bief_deallobj(s1ak  )
      CALL bief_deallobj(epn   )
      CALL bief_deallobj(epc   )
      CALL bief_deallobj(ep    )
      CALL bief_deallobj(s0ep  )
      CALL bief_deallobj(s1ep  )
      CALL bief_deallobj(kborf )
      CALL bief_deallobj(kborl )
      CALL bief_deallobj(kborsave)
      CALL bief_deallobj(kbors )
      CALL bief_deallobj(eborf )
      CALL bief_deallobj(eborl )
      CALL bief_deallobj(eborsave)
      CALL bief_deallobj(ebors )
      CALL bief_deallobj(akborf)
      CALL bief_deallobj(akborl)
      CALL bief_deallobj(akbors)
      CALL bief_deallobj(aeborf)
      CALL bief_deallobj(aeborl)
      CALL bief_deallobj(aebors)
      CALL bief_deallobj(bkborf)
      CALL bief_deallobj(bkborl)
      CALL bief_deallobj(bkbors)
      CALL bief_deallobj(beborf)
      CALL bief_deallobj(beborl)
      CALL bief_deallobj(bebors)
      CALL bief_deallobj(likbof)
      CALL bief_deallobj(likbol)
      CALL bief_deallobj(likbos)
      CALL bief_deallobj(liebof)
      CALL bief_deallobj(liebol)
      CALL bief_deallobj(liebos)
!     SPALART
      CALL bief_deallobj(nun   )
      CALL bief_deallobj(nuc   )
      CALL bief_deallobj(nu    )
      CALL bief_deallobj(s0nu  )
      CALL bief_deallobj(s1nu  )
      CALL bief_deallobj(nuborf)
      CALL bief_deallobj(nuborl)
      CALL bief_deallobj(nubors)
      CALL bief_deallobj(nuborsave)
      CALL bief_deallobj(anuborf)
      CALL bief_deallobj(anuborl)
      CALL bief_deallobj(anubors)
      CALL bief_deallobj(bnuborf)
      CALL bief_deallobj(bnuborl)
      CALL bief_deallobj(bnubors)
      CALL bief_deallobj(linubof)
      CALL bief_deallobj(linubol)
      CALL bief_deallobj(linubos)
      CALL bief_deallobj(wdist)
      CALL bief_deallobj(strain)
      CALL bief_deallobj(rotan)
!
      CALL bief_deallobj(u2d   )
      CALL bief_deallobj(un2d   )
      CALL bief_deallobj(fu    )
      CALL bief_deallobj(v2d   )
      CALL bief_deallobj(vn2d   )
      CALL bief_deallobj(fv    )
      CALL bief_deallobj(ubor2d)
      CALL bief_deallobj(vbor2d)
      CALL bief_deallobj(flbor   )
      CALL bief_deallobj(flblim  )
      CALL bief_deallobj(uconv )
      CALL bief_deallobj(vconv )
      CALL bief_deallobj(uconvc)
      CALL bief_deallobj(vconvc)
      CALL bief_deallobj(wsconv)
      CALL bief_deallobj(dm1)

      CALL bief_deallobj(h)
      CALL bief_deallobj(hn)
      CALL bief_deallobj(hprop )
      CALL bief_deallobj(nuwave)
      CALL bief_deallobj(hbor)
      CALL bief_deallobj(dh)
      CALL bief_deallobj(smh   )

      CALL bief_deallobj(pluie  )
      IF(ncsize.LE.1.AND.rain) THEN
        ALLOCATE(parapluie%R(1))
      ENDIF
      CALL bief_deallobj(parapluie)
      CALL bief_deallobj(maxz )
      CALL bief_deallobj(tmaxz)

      CALL bief_deallobj(fxh   )
      CALL bief_deallobj(fyh   )

      CALL bief_deallobj(dirmoy )
      CALL bief_deallobj(hm0    )
      CALL bief_deallobj(tpr5   )
      CALL bief_deallobj(orbvel )

      CALL bief_deallobj(volu   )
      CALL bief_deallobj(volun  )
      CALL bief_deallobj(volu3d )
      IF(ncsize.GT.1) THEN
        CALL bief_deallobj(volupar)
        CALL bief_deallobj(volunpar)
        CALL bief_deallobj(volu3dpar)
      ELSE
        NULLIFY(volupar%R)
        ALLOCATE(volupar%R(1))
        CALL bief_deallobj(volupar)
        NULLIFY(volunpar%R)
        ALLOCATE(volunpar%R(1))
        CALL bief_deallobj(volunpar)
        NULLIFY(volu3dpar%R)
        ALLOCATE(volu3dpar%R(1))
        CALL bief_deallobj(volu3dpar)
      ENDIF
      CALL bief_deallobj(volut  )
      CALL bief_deallobj(volu2d)
      CALL bief_deallobj(v2dpar)
      CALL bief_deallobj(unsv2d)
      CALL bief_deallobj(unsv3d)
      CALL bief_deallobj(fluint )
      CALL bief_deallobj(fluext )
      IF(ncsize.GT.1) THEN
        CALL bief_deallobj(fluextpar)
      ELSE
        NULLIFY(fluextpar%R)
        ALLOCATE(fluextpar%R(1))
        CALL bief_deallobj(fluextpar)
      ENDIF
      CALL bief_deallobj(flint2)

      CALL bief_deallobj(viscvi)
      CALL bief_deallobj(viscta)
      CALL bief_deallobj(rotat   )
      CALL bief_deallobj(deltar  )
      CALL bief_deallobj(rho)
      CALL bief_deallobj(rhopot)
      CALL bief_deallobj(ri      )
      CALL bief_deallobj(rugof )
      CALL bief_deallobj(rugol   )
      CALL bief_deallobj(cf    )
      CALL bief_deallobj(windspd)
      CALL bief_deallobj(rho4aed2)
      CALL bief_deallobj(temp4aed2)
      CALL bief_deallobj(sali4aed2)
      CALL bief_deallobj(thick4aed2)
      CALL bief_deallobj(benth4aed2)

      CALL dealloc_meteo()

      CALL bief_deallobj(masini )
      CALL bief_deallobj(masse  )
      CALL bief_deallobj(massen )
      CALL bief_deallobj(flucum )
      CALL bief_deallobj(flux   )
      CALL bief_deallobj(fluxb  )
      CALL bief_deallobj(zf    )
      CALL bief_deallobj(zfe   )
      CALL bief_deallobj(dhhn  )

      CALL bief_deallobj(xflot )
      CALL bief_deallobj(yflot )
      CALL bief_deallobj(zflot )
      CALL bief_deallobj(shpflo)
      CALL bief_deallobj(shzflo)
      CALL bief_deallobj(tagflo)
      CALL bief_deallobj(clsflo)
      CALL bief_deallobj(eltflo)
      CALL bief_deallobj(etaflo)

      CALL bief_deallobj(u_sce  )
      CALL bief_deallobj(v_sce  )
      CALL bief_deallobj(w_sce  )
      CALL bief_deallobj(ak_sce )
      CALL bief_deallobj(ep_sce )
      CALL bief_deallobj(nu_sce )
!     block
      CALL bief_deallobj(itrav3)

      CALL bief_deallobj(msupg)
      CALL bief_deallobj(mmurd)
      CALL bief_deallobj(murd_tf)
      CALL bief_deallobj(mdiff)
      CALL bief_deallobj(mtra1)
      CALL bief_deallobj(mtra2)
      CALL bief_deallobj(mat2d )
      CALL bief_deallobj(mbor2d)
      CALL bief_deallobj(matr2h)

      CALL bief_deallobj(trav3 )
      CALL bief_deallobj(tb2   )
      CALL bief_deallobj(sem3d )
      CALL bief_deallobj(sources)
      CALL bief_deallobj(trav2 )
      CALL bief_deallobj(sem2d )

      CALL bief_deallobj(te1)
      CALL bief_deallobj(te2)
      CALL bief_deallobj(te3)
      CALL bief_deallobj(zflats)
      CALL bief_deallobj(zconv  )
      CALL bief_deallobj(grazco )
      CALL bief_deallobj(flodel )
      CALL bief_deallobj(flopar )
      CALL bief_deallobj(flulim )
      CALL bief_deallobj(svide  )
      CALL bief_deallobj(wchu   )
      IF(s3d_sedi) THEN
        CALL bief_deallobj(trav1  )
        DEALLOCATE(s3d_conc)
        DEALLOCATE(s3d_epai)
        DEALLOCATE(s3d_toce)
        DO i=1,s3d_ncouch
          ALLOCATE(s3d_laytoce%ADR(i)%P%R(1))
          ALLOCATE(s3d_laythi%ADR(i)%P%R(1))
        ENDDO
        CALL bief_deallobj(s3d_laytoce)
        CALL bief_deallobj(s3d_laythi)
        CALL bief_deallobj(s3d_fludp)
        CALL bief_deallobj(s3d_fludpt)
        CALL bief_deallobj(s3d_zf_s)
        CALL bief_deallobj(s3d_esomt)
      ENDIF
      CALL bief_deallobj(s3d_wchu)
      CALL bief_deallobj(s3d_dmoy)
      CALL bief_deallobj(s3d_cref)
      CALL bief_deallobj(zref)
      CALL bief_deallobj(s3d_hdep)
      CALL bief_deallobj(s3d_fluer)
      CALL bief_deallobj(s3d_pdepo)
      CALL bief_deallobj(zr)
      CALL bief_deallobj(npf)

      CALL bief_deallobj(s3d_ivide)
      CALL bief_deallobj(s3d_temp)

      CALL bief_deallobj(s3d_epaico)
      CALL bief_deallobj(s3d_epainco)
      CALL bief_deallobj(s3d_pvsco)
      CALL bief_deallobj(s3d_pvsnco)
      CALL bief_deallobj(s3d_fluerc)
      CALL bief_deallobj(s3d_fluernc)
      CALL bief_deallobj(s3d_fludptc)
      CALL bief_deallobj(s3d_fludptnc)
      CALL bief_deallobj(s3d_fludpc)
      CALL bief_deallobj(s3d_fludpnc)
      CALL bief_deallobj(s3d_wcs)

      CALL bief_deallobj(altbus)
      CALL bief_deallobj(csbus )
      CALL bief_deallobj(cebus )
      CALL bief_deallobj(angbus)
      CALL bief_deallobj(lbus  )
      CALL bief_deallobj(dbus  )
      CALL bief_deallobj(lrgbus)
      CALL bief_deallobj(haubus)
      CALL bief_deallobj(secbus)
      CALL bief_deallobj(clpbus)
      CALL bief_deallobj(cv    )
      CALL bief_deallobj(c56   )
      CALL bief_deallobj(cv5   )
      CALL bief_deallobj(c5    )
      CALL bief_deallobj(ctrash)
      CALL bief_deallobj(fricbus)
      CALL bief_deallobj(longbus)
      CALL bief_deallobj(circ  )
      CALL bief_deallobj(delbus)
      CALL bief_deallobj(aabus )
      CALL bief_deallobj(tbus  )
!
      CALL bief_deallobj(wip    )
      CALL bief_deallobj(wipdx  )
      CALL bief_deallobj(wipdy  )
      CALL bief_deallobj(ustokes)
      CALL bief_deallobj(vstokes)
      CALL bief_deallobj(wstokes)
      CALL bief_deallobj(wst1   )
      CALL bief_deallobj(wst2   )
      CALL bief_deallobj(us2d   )
      CALL bief_deallobj(vs2d   )
      CALL bief_deallobj(dus2d  )
      CALL bief_deallobj(dvs2d  )
      CALL bief_deallobj(dus2dx )
      CALL bief_deallobj(dvs2dy )
      CALL bief_deallobj(ustx   )
      CALL bief_deallobj(vsty   )
      CALL bief_deallobj(fdx    )
      CALL bief_deallobj(fdy    )
      CALL bief_deallobj(fbx    )
      CALL bief_deallobj(fby    )
      CALL bief_deallobj(cfwc   )
      CALL bief_deallobj(fdk    )
      CALL bief_deallobj(fwx    )
      CALL bief_deallobj(fwy    )
      CALL bief_deallobj(hborcli)
      CALL bief_deallobj(uborcli)
      CALL bief_deallobj(vborcli)
!
!     Deallocate some global variables allocated
!     at some point in the code
!
        IF(ALLOCATED(calcflu)) DEALLOCATE(calcflu)
        IF(ALLOCATED(calcrain)) DEALLOCATE(calcrain)
!
      IF(ALLOCATED(particules)) THEN
        DO i=1,nflot_max
!!          DEALLOCATE(PARTICULES(I)%COMPO) !????
!!          DEALLOCATE(PARTICULES(I)%HAP)   !????
        ENDDO
        DEALLOCATE(particules)
      ENDIF
!
!     TEL4DEL
!
      CALL dealloc_tel4del()
!
!     TPXO
!
      CALL dealloc_tpxo()
!
!     GOTM
!
      IF(iturbv.EQ.6) THEN
        CALL dealloc_gotm()
      ENDIF
!
!     Breach
!
!
!     Deallocating arrays allocated in lecdon_TELEMAC3D
!
      DEALLOCATE(frtype)
      DEALLOCATE(debimp)
      DEALLOCATE(cotimp)
      DEALLOCATE(vitimp)
      DEALLOCATE(flux_boundaries)
      DEALLOCATE(tracer)
      DEALLOCATE(pts_curves)
      DEALLOCATE(sta_dis_curves)
      DEALLOCATE(profvel)
      DEALLOCATE(verprovel)
      DEALLOCATE(verprotra)
      DEALLOCATE(dirflu)
      DEALLOCATE(okq3)
      DEALLOCATE(oksl3)
      DEALLOCATE(oktr3)
      DEALLOCATE(okvit3)
      DEALLOCATE(bnd_tide)
      DEALLOCATE(isce)
      DEALLOCATE(ksce)
      DEALLOCATE(xsce)
      DEALLOCATE(ysce)
      DEALLOCATE(zsce)
      DEALLOCATE(qsce)
      DEALLOCATE(qsce2)
      DEALLOCATE(tasce)
      DEALLOCATE(usce)
      DEALLOCATE(vsce)
      DEALLOCATE(wsce)
      DEALLOCATE(okdebsce)
      DEALLOCATE(oktrsce)
      DEALLOCATE(schcta)
      DEALLOCATE(train)
      DEALLOCATE(trac0)
      DEALLOCATE(t0ac)
      DEALLOCATE(betac)
      DEALLOCATE(dnutah)
      DEALLOCATE(dnutav)
      DEALLOCATE(s_adv)
      DEALLOCATE(nom_adv)
      DEALLOCATE(slvdta)
      DEALLOCATE(nametrac)
      DEALLOCATE(bedflo)
      DEALLOCATE(bedqarea)
      DEALLOCATE(optadv_tr)

      IF(deja_murd3d_pos) THEN
        DEALLOCATE(indic_murd3d_pos)
      ENDIF

      IF(ALLOCATED(qz)) DEALLOCATE(qz)
      ! Resetting initial values of variables
      deja_murd3d_pos = .false.
      dejalu_trisou = .false.

      ! Temporary
      IF(.NOT.inclus(coupling,'GAIA  ')) THEN
        CALL bief_deallobj(fludp)
        CALL bief_deallobj(fludpt)
        CALL bief_deallobj(fluer)
        DEALLOCATE(num_isusp_icla)
      ENDIF

      IF( 13*int(waqprocess/13).EQ.waqprocess ) THEN
        IF( inclus(coupling,'WAQTEL')) THEN
#if defined HAVE_AED2
          CALL clean_aed2_models()
#endif
        ENDIF
      ENDIF
!
!=======================================================================
!
! WRITES OUT TO LISTING :
!
      IF(listin) THEN
        WRITE(lu,23)
      ENDIF
23    FORMAT(1x,///,21x,'*************************************',/,
     &21x,              '*    END OF MEMORY ORGANIZATION:    *',/,
     &21x,              '*************************************',/)
!
!-----------------------------------------------------------------------
!
      RETURN
!
!
      END SUBROUTINE deall_telemac3d