Changeset 1001 for palm/trunk/SOURCE/diffusion_e.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_e.f90 (modified) (5 diffs)

palm/trunk/SOURCE/diffusion_e.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! most arrays comunicated by module instead of parameter list
 !
 ! Former revisions:
@@ -64 +64 @@
 ! Call for all grid points
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, var, &
-                            var_reference, rif, tend, zu, zw )
-
+    SUBROUTINE diffusion_e( var, var_reference )
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -75 +75 @@
 
        INTEGER ::  i, j, k
-       REAL            ::  dvar_dz, l_stable, phi_m, var_reference
-       REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt+1), &
-                           l_grid(1:nzt), zu(nzb:nzt+1), zw(nzb:nzt+1)
-       REAL, DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: diss, tend
-       REAL, DIMENSION(:,:), POINTER   ::  rif
-       REAL, DIMENSION(:,:,:), POINTER ::  e, km, var
+       REAL    ::  dvar_dz, l_stable, phi_m, var_reference
+
+       REAL, DIMENSION(:,:,:), POINTER ::  var
        REAL, DIMENSION(nzb+1:nzt,nys:nyn) ::  dissipation, l, ll
  
@@ -285 +282 @@
 ! Call for grid point i,j
 !------------------------------------------------------------------------------!
-    SUBROUTINE diffusion_e_ij( i, j, ddzu, dd2zu, ddzw, diss, e, km, l_grid, &
-                               var, var_reference, rif, tend, zu, zw )
-
+    SUBROUTINE diffusion_e_ij( i, j, var, var_reference )
+
+       USE arrays_3d
        USE control_parameters
        USE grid_variables
@@ -295 +292 @@
        IMPLICIT NONE
 
-       INTEGER         ::  i, j, k
-       REAL            ::  dvar_dz, l_stable, phi_m, var_reference
-       REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt+1), &
-                           l_grid(1:nzt), zu(nzb:nzt+1), zw(nzb:nzt+1)
-       REAL, DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: diss, tend
-       REAL, DIMENSION(:,:), POINTER   ::  rif
-       REAL, DIMENSION(:,:,:), POINTER ::  e, km, var
-       REAL, DIMENSION(nzb+1:nzt)    ::  dissipation, l, ll
+       INTEGER ::  i, j, k
+       REAL    ::  dvar_dz, l_stable, phi_m, var_reference
+
+       REAL, DIMENSION(:,:,:), POINTER ::  var
+       REAL, DIMENSION(nzb+1:nzt) ::  dissipation, l, ll