Changeset 1001 for palm/trunk/SOURCE/diffusion_s.f90

Timestamp:

Sep 13, 2012 2:08:46 PM (12 years ago)

Author:

raasch

Message:

leapfrog timestep scheme and upstream-spline advection scheme completely removed from the code,
reading of dt_fixed from restart file removed

File:

• palm/trunk/SOURCE/diffusion_s.f90 (modified) (5 diffs)

palm/trunk/SOURCE/diffusion_s.f90

@@ -3 +3 @@
 !------------------------------------------------------------------------------!
 ! Current revisions:
-! -----------------
+! ------------------
+! some arrays comunicated by module instead of parameter list
 !
 ! Former revisions:
@@ -54 +55 @@
 ! Call for all grid points
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_s( ddzu, ddzw, kh, s, s_flux_b, s_flux_t, &
-                            wall_s_flux, tend )
-
+    SUBROUTINE diffusion_s( s, s_flux_b, s_flux_t, wall_s_flux )
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -65 +66 @@
        INTEGER ::  i, j, k
        REAL    ::  vertical_gridspace
-       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
-       REAL    ::  tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)
        REAL    ::  wall_s_flux(0:4)
-       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
-       REAL, DIMENSION(:,:,:), POINTER ::  kh, s
+       REAL, DIMENSION(nysg:nyng,nxlg:nxrg) ::  s_flux_b, s_flux_t
+       REAL, DIMENSION(:,:,:), POINTER ::  s
 
        DO  i = nxl, nxr
@@ -166 +165 @@
 ! Call for grid point i,j
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_s_ij( i, j, ddzu, ddzw, kh, s, s_flux_b, s_flux_t, &
-                               wall_s_flux, tend )
-
+    SUBROUTINE diffusion_s_ij( i, j, s, s_flux_b, s_flux_t, wall_s_flux )
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -177 +176 @@
        INTEGER ::  i, j, k
        REAL    ::  vertical_gridspace
-       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
-       REAL    ::  tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)
        REAL    ::  wall_s_flux(0:4)
-       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
-       REAL, DIMENSION(:,:,:), POINTER ::  kh, s
+       REAL, DIMENSION(nysg:nyng,nxlg:nxrg) ::  s_flux_b, s_flux_t
+       REAL, DIMENSION(:,:,:), POINTER ::  s
 
 !