Changeset 261

Timestamp:

06/17/19 16:32:14 (5 years ago)

Author:

aquiquet

Message:

Mass conservation in double precision to deal with short timesteps

Location:

trunk/SOURCES

Files:

• conserv-mass-adv-diff_sept2009_mod.f90 (modified) (2 diffs)
• resol_adv_diff_2D-sept2009.f90 (modified) (3 diffs)

trunk/SOURCES/conserv-mass-adv-diff_sept2009_mod.f90

@@ -36 +36 @@
 integer :: itesti
 integer :: itour
-
+double precision, dimension(nx,ny) :: H_local
 
 contains
@@ -503 +503 @@
 
 
-call resol_adv_diff_2D_vect (Dmx,Dmy,advmx,advmy,Hp,i_Hp,bilmass,vieuxH,H)
-
+call resol_adv_diff_2D_vect (dble(Dmx),dble(Dmy),dble(advmx),dble(advmy),dble(Hp),i_Hp,dble(bilmass),dble(vieuxH),H_local)
+
+H(:,:) = real(H_local(:,:))
 !$OMP PARALLEL
 !$OMP DO

trunk/SOURCES/resol_adv_diff_2D-sept2009.f90

@@ -18 +18 @@
 
 implicit none
-real :: omega   !< parametre schema temporel de la resolution partie diffusion
-real :: mu_adv  !< parametre schema temporel de la resolution advection
-real :: upwind  !< schema spatial pour l'advection
+double precision :: omega   !< parametre schema temporel de la resolution partie diffusion
+double precision :: mu_adv  !< parametre schema temporel de la resolution advection
+double precision :: upwind  !< schema spatial pour l'advection
 
 ! tableaux de travail. resolution M H = Frelax
-real,dimension(nx,ny) :: crelax      !< diagnonale de M
-real,dimension(nx,ny) :: arelax      !< sous diagonale selon x
-real,dimension(nx,ny) :: brelax      !< sur diagonale selon x
-real,dimension(nx,ny) :: drelax      !< sous diagonale selon y
-real,dimension(nx,ny) :: erelax      !< sur diagonale selon y
-real,dimension(nx,ny) :: frelax      !< vecteur
-real,dimension(nx,ny) :: c_west      !< sur demi mailles Ux
-real,dimension(nx,ny) :: c_east      !< sur demi mailles Ux
-real,dimension(nx,ny) :: c_north     !< sur demi mailles Uy
-real,dimension(nx,ny) :: c_south     !< sur demi mailles Uy
-
-real,dimension(nx,ny) :: bdx         !< pente socle
-real,dimension(nx,ny) :: bdy         !< pente socle
-
-real,dimension(nx,ny) :: hdx         !< pente epaisseur
-real,dimension(nx,ny) :: hdy         !< pente epaisseur
+double precision,dimension(nx,ny) :: crelax      !< diagnonale de M
+double precision,dimension(nx,ny) :: arelax      !< sous diagonale selon x
+double precision,dimension(nx,ny) :: brelax      !< sur diagonale selon x
+double precision,dimension(nx,ny) :: drelax      !< sous diagonale selon y
+double precision,dimension(nx,ny) :: erelax      !< sur diagonale selon y
+double precision,dimension(nx,ny) :: frelax      !< vecteur
+double precision,dimension(nx,ny) :: c_west      !< sur demi mailles Ux
+double precision,dimension(nx,ny) :: c_east      !< sur demi mailles Ux
+double precision,dimension(nx,ny) :: c_north     !< sur demi mailles Uy
+double precision,dimension(nx,ny) :: c_south     !< sur demi mailles Uy
+
+double precision,dimension(nx,ny) :: bdx         !< pente socle
+double precision,dimension(nx,ny) :: bdy         !< pente socle
+
+double precision,dimension(nx,ny) :: hdx         !< pente epaisseur
+double precision,dimension(nx,ny) :: hdy         !< pente epaisseur
 
 
@@ -95 +95 @@
 implicit none
 
-real,dimension(nx,ny), intent(in) :: Dfx          !< terme diffusif selon x
-real,dimension(nx,ny), intent(in) :: Dfy          !< terme diffusif selon y
-real,dimension(nx,ny), intent(in) :: Advx         !< terme advectif selon x
-real,dimension(nx,ny), intent(in) :: Advy         !< terme advectif selon y
-real,dimension(nx,ny), intent(in) :: vieuxH       !< ancienne valeur de H
-real,dimension(nx,ny), intent(out):: newH         !< nouvelle valeur de H
-
-real,dimension(nx,ny), intent(in) :: bil          !< bilan de masse pour la colonne
-
-real,dimension(nx,ny), intent(in) :: H_presc      !< H value if prescribed
+double precision,dimension(nx,ny), intent(in) :: Dfx          !< terme diffusif selon x
+double precision,dimension(nx,ny), intent(in) :: Dfy          !< terme diffusif selon y
+double precision,dimension(nx,ny), intent(in) :: Advx         !< terme advectif selon x
+double precision,dimension(nx,ny), intent(in) :: Advy         !< terme advectif selon y
+double precision,dimension(nx,ny), intent(in) :: vieuxH       !< ancienne valeur de H
+double precision,dimension(nx,ny), intent(out):: newH         !< nouvelle valeur de H
+
+double precision,dimension(nx,ny), intent(in) :: bil          !< bilan de masse pour la colonne
+
+double precision,dimension(nx,ny), intent(in) :: H_presc      !< H value if prescribed
 integer,dimension(nx,ny), intent(in) :: i_Hpresc  !< 1 if H is prescribed on this node, else 0
 
@@ -126 +126 @@
 !!$real,dimension(nx,ny) :: hdy         !< pente epaisseur
 
-real :: frdx,frdy                    !< pour calcul frelax : termes diffusion
-real :: fraxw,fraxe,frays,frayn      !< termes advection
-
-real,dimension(nx,ny) :: deltah      ! dans calcul relax
-real :: delh                         ! dans calcul relax
-real :: testh                        ! dans calcul relax
+double precision :: frdx,frdy                    !< pour calcul frelax : termes diffusion
+double precision :: fraxw,fraxe,frays,frayn      !< termes advection
+
+double precision,dimension(nx,ny) :: deltah      ! dans calcul relax
+double precision :: delh                         ! dans calcul relax
+double precision :: testh                        ! dans calcul relax
 
 logical :: stopp
 integer :: ntour
-real :: reste
+double precision :: reste