Changeset 1822 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Apr 7, 2016 7:49:42 AM (8 years ago)

Author:

hoffmann

Message:

changes in LPM and bulk cloud microphysics

File:

• palm/trunk/SOURCE/data_output_3d.f90 (modified) (7 diffs)

palm/trunk/SOURCE/data_output_3d.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-!
+! prr vertical dimensions set to nzb_do to nzt_do. Unused variables deleted.
 !
 ! Former revisions:
@@ -122 +122 @@
         
     USE control_parameters,                                                    &
-        ONLY:  avs_data_file, cloud_physics, do3d, do3d_avs_n,                 &
-               do3d_no, do3d_time_count, io_blocks, io_group,                  &
-               message_string, ntdim_3d,                                       &
-               nz_do3d, plot_3d_precision, psolver, simulated_time,            &
-               simulated_time_chr, skip_do_avs, time_since_reference_point
+        ONLY:  cloud_physics, do3d, do3d_no, do3d_time_count, io_blocks,       &
+               io_group, message_string, ntdim_3d, nz_do3d, psolver,           &
+               simulated_time, time_since_reference_point
         
     USE cpulog,                                                                &
@@ -132 +130 @@
         
     USE indices,                                                               &
-        ONLY:  nbgp, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg, nzt,  &
-               nzb
+        ONLY:  nbgp, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg, nzb
         
     USE kinds
@@ -163 +160 @@
     IMPLICIT NONE
 
-    CHARACTER (LEN=9) ::  simulated_time_mod  !< 
-
     INTEGER(iwp) ::  av        !< 
     INTEGER(iwp) ::  i         !< 
@@ -173 +168 @@
     INTEGER(iwp) ::  nzb_do    !< vertical lower limit for data output
     INTEGER(iwp) ::  nzt_do    !< vertical upper limit for data output
-    INTEGER(iwp) ::  pos       !< 
-    INTEGER(iwp) ::  prec      !< 
-    INTEGER(iwp) ::  psi       !< 
 
     LOGICAL      ::  found     !< 
@@ -374 +366 @@
                 DO  i = nxlg, nxrg
                    DO  j = nysg, nyng
-                      DO  k = nzb, nzt+1
+                      DO  k = nzb_do, nzt_do
                          local_pf(i,j,k) = prr(k,j,i)
                       ENDDO
@@ -383 +375 @@
                 DO  i = nxlg, nxrg
                    DO  j = nysg, nyng
-                      DO  k = nzb, nzt+1
+                      DO  k = nzb_do, nzt_do
                          local_pf(i,j,k) = prr_av(k,j,i)
                       ENDDO