deall_sisyphe.f

Go to the documentation of this file.

!                   **************************
                    SUBROUTINE deall_sisyphe
!                   **************************
!
!
!***********************************************************************
! SISYPHE   V7P1                                   19/05/2016
!***********************************************************************
!
!brief    Memory deallocation of structures, aliases, blocks...
!
!Author  R-S MOURADI (LNHE)
!
!history  P TASSI, F CORDIER, S PAVAN
!+        18/09/2018
!+        V8P0
!     + Deallocates memory for SANFRA
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      USE declarations_sisyphe
!
      IMPLICIT NONE

!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER :: I, K
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!
      CALL outbief(mesh)
!     Deallocate the mesh structure
      CALL dealmesh(mesh)

      NULLIFY(ikle)
      NULLIFY(x)
      NULLIFY(y)
      NULLIFY(nelem)
      NULLIFY(nelmax)
      NULLIFY(nptfr)
      NULLIFY(nptfrx)
      NULLIFY(typelm)
      NULLIFY(npoin)
      NULLIFY(npmax)
      NULLIFY(mxptvs)
      NULLIFY(lv)
      CALL bief_deallobj(varsor)
!
!     Deallocating all the blocks first
      CALL bief_deallobj(s     )
      CALL bief_deallobj(e     )
      CALL bief_deallobj(z     )
      CALL bief_deallobj(del_z )
      CALL bief_deallobj(zf_c  )
      CALL bief_deallobj(zf_s  )
      CALL bief_deallobj(esomt )
      CALL bief_deallobj(emax  )
      CALL bief_deallobj(q     )
      CALL bief_deallobj(qu    )
      CALL bief_deallobj(qv    )
      CALL bief_deallobj(del_qu)
      CALL bief_deallobj(del_qv)
      CALL bief_deallobj(del_uw)
      CALL bief_deallobj(u2d   )
      CALL bief_deallobj(v2d   )
      CALL bief_deallobj(qs    )
      CALL bief_deallobj(qsx   )
      CALL bief_deallobj(qsy   )
      CALL bief_deallobj(qs_c  )
      CALL bief_deallobj(qsxc  )
      CALL bief_deallobj(qsyc  )
      CALL bief_deallobj(qs_s  )
      CALL bief_deallobj(qsxs  )
      CALL bief_deallobj(qsys  )
      CALL bief_deallobj(hiding)
      CALL bief_deallobj(zf    )
      CALL bief_deallobj(zr    )
      CALL bief_deallobj(radsec)
      CALL bief_deallobj(zref  )
      CALL bief_deallobj(chestr)
      CALL bief_deallobj(coefpn)
      CALL bief_deallobj(cf    )
      CALL bief_deallobj(tob   )
      CALL bief_deallobj(tobw  )
      CALL bief_deallobj(mu    )
      CALL bief_deallobj(ksp   )
      CALL bief_deallobj(ks    )
      CALL bief_deallobj(ksr   )
      CALL bief_deallobj(thetaw)
      CALL bief_deallobj(fw    )
      CALL bief_deallobj(uw    )
      CALL bief_deallobj(hw    )
      CALL bief_deallobj(tw    )
      CALL bief_deallobj(dzf_gf)
      CALL bief_deallobj(acladm)
      CALL bief_deallobj(unladm)
      CALL bief_deallobj(hcpl  )
      CALL bief_deallobj(ecpl  )
      CALL bief_deallobj(elay  )
      CALL bief_deallobj(estrat)
      CALL bief_deallobj(kx    )
      CALL bief_deallobj(ky    )
      CALL bief_deallobj(kz    )
      CALL bief_deallobj(uconv )
      CALL bief_deallobj(vconv )
      CALL bief_deallobj(unorm )
      CALL bief_deallobj(disp  )
      CALL bief_deallobj(disp_c)
      CALL bief_deallobj(maskb )
      CALL bief_deallobj(mask  )
      CALL bief_deallobj(afbor )
      CALL bief_deallobj(bfbor )
      CALL bief_deallobj(flbor )
      CALL bief_deallobj(q2bor )
!     BOUNDARY FLUX FOR CALL TO CVDFTR
      CALL bief_deallobj(flbor_sis )
      CALL bief_deallobj(flbortra  )
      CALL bief_deallobj(cstaeq)
!     MAK ADDITION
      CALL bief_deallobj(csratio)
      CALL bief_deallobj(hn    )
      CALL bief_deallobj(hclip )
      CALL bief_deallobj(hprop )
      CALL bief_deallobj(volu2d)
      CALL bief_deallobj(v2dpar)
      CALL bief_deallobj(unsv2d)
      CALL bief_deallobj(mpm_aray)
      CALL bief_deallobj(flulim  )

      ! MSK
      CALL bief_deallobj(maskel)
      CALL bief_deallobj(msktmp)
      CALL bief_deallobj(maskpt)
!     FOR MIXED SEDIMENTS
      CALL bief_deallobj(fluer_vase )
      CALL bief_deallobj(toce_mixte )
      CALL bief_deallobj(ms_sable   )
      CALL bief_deallobj(ms_vase    )
      CALL bief_deallobj( liebor)
      CALL bief_deallobj( liqbor)
      CALL bief_deallobj( limtec)
      CALL bief_deallobj( numliq)
      CALL bief_deallobj( clt   )
      CALL bief_deallobj( clu   )
      CALL bief_deallobj( clv   )
      CALL bief_deallobj( limdif)
      CALL bief_deallobj( licbor)
      CALL bief_deallobj( lihbor)
      CALL bief_deallobj( boundary_colour)
      CALL bief_deallobj( limpro)
      CALL bief_deallobj( indic )
      CALL bief_deallobj( it1   )
      CALL bief_deallobj( it2   )
      CALL bief_deallobj( it3   )
      CALL bief_deallobj( it4   )
! NUMBER OF LAYERS
      CALL bief_deallobj( nlayer)
      CALL bief_deallobj(breach)
      CALL bief_deallobj(ifamas)

      DEALLOCATE(avail)
      DEALLOCATE(es)
      DEALLOCATE(es_sable)
      DEALLOCATE(es_vase)
      DEALLOCATE(conc)
      DEALLOCATE(ivide)
      DEALLOCATE(sanfra)

      CALL bief_deallobj(masktr)
      CALL bief_deallobj(ebor  )
      CALL bief_deallobj(qbor  )
      DO i=1,nsicla
        DO k=1,nomblay
          ALLOCATE(avai%ADR(k+(i-1)*nomblay)%P%R(1))
        ENDDO
      ENDDO
      CALL bief_deallobj(avai  )
      DO k=1,nomblay
        ALLOCATE(laythi%ADR(k)%P%R(1))
      ENDDO
      CALL bief_deallobj(laythi)
      DO k=1,nomblay
        ALLOCATE(layconc%ADR(k)%P%R(1))
      ENDDO
      CALL bief_deallobj(layconc)
!
      CALL bief_deallobj(qscl  )
      CALL bief_deallobj(qscl_c)
      CALL bief_deallobj(qsclxc)
      CALL bief_deallobj(qsclyc)
      CALL bief_deallobj(zfcl  )
      CALL bief_deallobj(zfcl_c)
!
      CALL bief_deallobj(cbor  )
      CALL bief_deallobj(qscl_s)
      CALL bief_deallobj(qsclxs)
      CALL bief_deallobj(qsclys)
      CALL bief_deallobj(zfcl_s)
      CALL bief_deallobj(fludp )
      CALL bief_deallobj(fludpt)
      CALL bief_deallobj(fluer )
      CALL bief_deallobj(fluert)
      CALL bief_deallobj(cs    )
      CALL bief_deallobj(ctild )
      CALL bief_deallobj(cst   )
!
      CALL bief_deallobj(zfcl_ms)
      CALL bief_deallobj(calfa_cl)
      CALL bief_deallobj(salfa_cl)
      CALL bief_deallobj(flbcla)

      CALL bief_deallobj(am1_s)
      CALL bief_deallobj(am2_s)
      CALL bief_deallobj(mbor )

      CALL bief_deallobj( w1 )
      CALL bief_deallobj( te1)
      CALL bief_deallobj( te2)
      CALL bief_deallobj( te3)

      CALL bief_deallobj(varcl)
      CALL bief_deallobj(prive)
      CALL bief_deallobj(tb   )
      CALL bief_deallobj(tb2  )

      IF(nadvar.GT.0) THEN
        CALL bief_deallobj(advar )
      ENDIF

      ! Resetting variable
      init_fluxpr = .true.
      IF(deja_rfc) THEN
        deja_rfc = .false.
        DEALLOCATE(infic_rfc)
        DEALLOCATE(time_rfc)
      ENDIF
      IF(deja_flusec) THEN
        deja_flusec = .false.
        old_method_flusec = .false.
        DEALLOCATE(nseg)
        DEALLOCATE(liste)
        DEALLOCATE(volnegs)
        DEALLOCATE(volposs)
        DEALLOCATE(volnegc)
        DEALLOCATE(volposc)
        DEALLOCATE(flx)
        DEALLOCATE(volneg)
        DEALLOCATE(volpos)
        DEALLOCATE(flxs)
        DEALLOCATE(flxc)
      ENDIF
      IF(deja_flusec2) THEN
        old_method_flusec = .false.
        deja_flusec2 = .false.
        DO i = 1,numberoflines_flusec2
          DEALLOCATE(fluxlinedata_flusec2(i)%SECTIONIDS)
          DEALLOCATE(fluxlinedata_flusec2(i)%DIRECTION)
        ENDDO
        DEALLOCATE(fluxlinedata_flusec2)
        DEALLOCATE(flux_flusec2)
        DEALLOCATE(volflux_flusec2)
      ENDIF

      ! Lecdon
      DEALLOCATE(soldis)
      DEALLOCATE(okcgl)
      DEALLOCATE(okqgl)
      DEALLOCATE(cbor_classe)
      DEALLOCATE(ctrlsc)
!
!-----------------------------------------------------------------------
!
!=======================================================================
!
! 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_sisyphe