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

Timestamp:

Dec 4, 2018 8:40:18 AM (5 years ago)

Author:

maronga

Message:

revised calculation of near surface air potential temperature

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added output of theta_2m
+! 
+! 3589 2018-11-30 15:09:51Z suehring
 ! Move the control parameter "salsa" from salsa_mod to control_parameters 
 ! (M. Kurppa)
@@ -249 +252 @@
     USE averaging,                                                             &
         ONLY:  e_av, ghf_av, lpt_av, lwp_av, ol_av, p_av, pc_av, pr_av, pt_av, &
-               q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsws_av, qv_av,        &
-               r_a_av, s_av, shf_av, ssws_av, ts_av, tsurf_av, u_av, us_av,    &
-               v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
+               pt_2m_av, q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsws_av,     &
+               qv_av, r_a_av, s_av, shf_av, ssws_av, ts_av, tsurf_av, u_av,    &
+               us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
 
     USE basic_constants_and_equations_mod,                                     &
@@ -497 +500 @@
                 ENDIF
                 vpt_av = 0.0_wp
+
+             CASE ( 'theta_2m*' )
+                IF ( .NOT. ALLOCATED( pt_2m_av ) )  THEN
+                   ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                pt_2m_av = 0.0_wp
 
              CASE ( 'w' )