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

Timestamp:

Dec 10, 2019 6:25:02 PM (5 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/data_output_2d.f90 (modified) (2 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Move 2-m potential temperature output to diagnostic_output_quantities
+! 
+! 4329 2019-12-10 15:46:36Z motisi
 ! Renamed wall_flags_0 to wall_flags_static_0
 ! 
@@ -1001 +1004 @@
                 ENDIF
                 IF ( mode == 'xy' )  level_z = zu
-
-             CASE ( 'theta_2m*_xy' )        ! 2d-array
-                IF ( av == 0 )  THEN
-                   DO  m = 1, surf_def_h(0)%ns
-                      i = surf_def_h(0)%i(m)
-                      j = surf_def_h(0)%j(m)
-                      local_pf(i,j,nzb+1) = surf_def_h(0)%pt_2m(m)
-                   ENDDO
-                   DO  m = 1, surf_lsm_h%ns
-                      i = surf_lsm_h%i(m)
-                      j = surf_lsm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_lsm_h%pt_2m(m)
-                   ENDDO
-                   DO  m = 1, surf_usm_h%ns
-                      i = surf_usm_h%i(m)
-                      j = surf_usm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%pt_2m(m)
-                   ENDDO
-                ELSE
-                   IF ( .NOT. ALLOCATED( pt_2m_av ) ) THEN
-                      ALLOCATE( pt_2m_av(nysg:nyng,nxlg:nxrg) )
-                      pt_2m_av = REAL( fill_value, KIND = wp )
-                   ENDIF
-                   DO  i = nxl, nxr
-                      DO  j = nys, nyn
-                         local_pf(i,j,nzb+1) = pt_2m_av(j,i)
-                      ENDDO
-                   ENDDO
-                ENDIF
-                resorted = .TRUE.
-                two_d = .TRUE.
-                level_z(nzb+1) = zu(nzb+1)
 
              CASE ( 'w_xy', 'w_xz', 'w_yz' )