Changeset 261 for trunk/SOURCES/resol_adv_diff_2D-sept2009.f90

Timestamp:

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

Author:

aquiquet

Message:

Mass conservation in double precision to deal with short timesteps

File:

• trunk/SOURCES/resol_adv_diff_2D-sept2009.f90 (modified) (3 diffs)

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