Changeset 3597 for palm/trunk/SOURCE/data_output_2d.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/data_output_2d.f90 (modified) (2 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added theta_2m
+! 
+! 3589 2018-11-30 15:09:51Z suehring
 ! Move the control parameter "salsa" from salsa_mod to control_parameters 
 ! (M. Kurppa)
@@ -1212 +1215 @@
                 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' )
                 flag_nr = 3