Changeset 1001 for palm/trunk/SOURCE/diffusion_u.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_u.f90 (modified) (6 diffs)

palm/trunk/SOURCE/diffusion_u.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! arrays comunicated by module instead of parameter list
 !
 ! Former revisions:
@@ -68 +68 @@
 ! Call for all grid points
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_u( ddzu, ddzw, km, tend, u, usws, uswst, v, w )
-
+    SUBROUTINE diffusion_u
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -78 +79 @@
        INTEGER ::  i, j, k
        REAL    ::  kmym, kmyp, kmzm, kmzp
-       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
-       REAL    ::  tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)
-       REAL, DIMENSION(:,:),   POINTER ::  usws, uswst
-       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
+
        REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs
 
@@ -93 +91 @@
 
        DO  i = nxlu, nxr
-          DO  j = nys,nyn
+          DO  j = nys, nyn
 !
 !--          Compute horizontal diffusion
@@ -221 +219 @@
 ! Call for grid point i,j
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_u_ij( i, j, ddzu, ddzw, km, tend, u, usws, &
-                               uswst, v, w )
-
+    SUBROUTINE diffusion_u_ij( i, j )
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -232 +230 @@
        INTEGER ::  i, j, k
        REAL    ::  kmym, kmyp, kmzm, kmzp
-       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
-       REAL    ::  tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)
-       REAL, DIMENSION(nzb:nzt+1)      ::  usvs
-       REAL, DIMENSION(:,:),   POINTER ::  usws, uswst
-       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
+
+       REAL, DIMENSION(nzb:nzt+1) ::  usvs
 
 !