Changeset 2031 for palm/trunk/SOURCE/sum_up_3d_data.f90

Timestamp:

Oct 21, 2016 3:11:58 PM (8 years ago)

Author:

knoop

Message:

Renamed variable rho to rho_ocean, rho_init to rho_ocean_init and rho_av to rho_ocean_av

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
 ! 
 ! Former revisions:
@@ -130 +130 @@
     USE arrays_3d,                                                             &
         ONLY:  dzw, e, nr, ol, p, pt, precipitation_rate, q, qc, ql, ql_c,     &
-               ql_v, qr, qsws, rho, s, sa, shf, ssws, ts, u, us, v, vpt, w, z0,&
+               ql_v, qr, qsws, rho_ocean, s, sa, shf, ssws, ts, u, us, v, vpt, w, z0,&
                z0h, z0q
 
@@ -136 +136 @@
         ONLY:  e_av, lpt_av, lwp_av, nr_av, ol_av, p_av, pc_av, pr_av, prr_av, &
                precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
-               ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av,  &
+               ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_ocean_av, s_av, sa_av,  &
                shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, &
                z0h_av, z0q_av
@@ -323 +323 @@
                 qv_av = 0.0_wp
 
-             CASE ( 'rho' )
-                IF ( .NOT. ALLOCATED( rho_av ) )  THEN
-                   ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-                ENDIF
-                rho_av = 0.0_wp
+             CASE ( 'rho_ocean' )
+                IF ( .NOT. ALLOCATED( rho_ocean_av ) )  THEN
+                   ALLOCATE( rho_ocean_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                rho_ocean_av = 0.0_wp
 
              CASE ( 's' )
@@ -657 +657 @@
              ENDDO
 
-          CASE ( 'rho' )
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   DO  k = nzb, nzt+1
-                      rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i)
+          CASE ( 'rho_ocean' )
+             DO  i = nxlg, nxrg
+                DO  j = nysg, nyng
+                   DO  k = nzb, nzt+1
+                      rho_ocean_av(k,j,i) = rho_ocean_av(k,j,i) + rho_ocean(k,j,i)
                    ENDDO
                 ENDDO