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

Timestamp:

Dec 10, 2019 6:25:02 PM (4 years ago)

Author:

suehring

Message:

New diagnostic output for 10-m wind speed; Diagnostic output of 2-m potential temperature moved to diagnostic output

File:

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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Move 2-m potential temperature output to diagnostic_output_quantities
+! 
+! 4182 2019-08-22 15:20:23Z scharf
 ! Corrected "Former revisions" section
 ! 
@@ -64 +67 @@
 
     USE arrays_3d,                                                             &
-        ONLY:  dzw, d_exner, e, heatflux_output_conversion, p,    &
-               pt, q, ql, ql_c, ql_v, s, u, v, vpt, w,                 &
+        ONLY:  dzw, d_exner, e, heatflux_output_conversion, p,                 &
+               pt, q, ql, ql_c, ql_v, s, u, v, vpt, w,                         &
                waterflux_output_conversion
 
     USE averaging,                                                             &
         ONLY:  e_av, ghf_av, lpt_av, lwp_av, ol_av, p_av, pc_av, pr_av, pt_av, &
-               pt_2m_av, q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsws_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
@@ -296 +299 @@
                 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' )
@@ -730 +727 @@
              ENDIF
 
-          CASE ( 'theta_2m*' )
-             IF ( ALLOCATED( pt_2m_av ) ) THEN    
-                DO  i = nxl, nxr
-                   DO  j = nys, nyn
-                      match_def = surf_def_h(0)%start_index(j,i) <=            &
-                                  surf_def_h(0)%end_index(j,i)
-                      match_lsm = surf_lsm_h%start_index(j,i) <=               &
-                                  surf_lsm_h%end_index(j,i)
-                      match_usm = surf_usm_h%start_index(j,i) <=               &
-                                  surf_usm_h%end_index(j,i)
-
-                      IF ( match_def )  THEN
-                         m = surf_def_h(0)%end_index(j,i)
-                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
-                                         surf_def_h(0)%pt_2m(m)
-                      ELSEIF ( match_lsm  .AND.  .NOT. match_usm )  THEN
-                         m = surf_lsm_h%end_index(j,i)
-                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
-                                         surf_lsm_h%pt_2m(m)
-                      ELSEIF ( match_usm )  THEN
-                         m = surf_usm_h%end_index(j,i)
-                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
-                                         surf_usm_h%pt_2m(m)
-                      ENDIF
-                   ENDDO
-                ENDDO
-             ENDIF
-             
-             
           CASE ( 't*' )
              IF ( ALLOCATED( ts_av ) ) THEN