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

Timestamp:

Feb 8, 2018 1:24:35 PM (6 years ago)

Author:

suehring

Message:

Output of ground-heat flux at natural- and urban-type surfaces in one output variable; enable restart data of _av variables that belong to both land- and urban-surface model

File:

• palm/trunk/SOURCE/data_output_2d.f90 (modified) (2 diffs)

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable output of ground-heat flux also at urban surfaces.
+! 
+! 2743 2018-01-12 16:03:39Z suehring
 ! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
 ! 
@@ -545 +548 @@
                 ENDIF
                 IF ( mode == 'xy' )  level_z = zu
+
+             CASE ( 'ghf*_xy' )        ! 2d-array
+                IF ( av == 0 )  THEN
+                   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%ghf(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%frac(0,m)     *          &
+                                            surf_usm_h%wghf_eb(m)        +      &
+                                            surf_usm_h%frac(1,m)     *          &
+                                            surf_usm_h%wghf_eb_green(m)  +      &
+                                            surf_usm_h%frac(2,m)     *          &
+                                            surf_usm_h%wghf_eb_window(m)
+                   ENDDO
+                ELSE
+                   DO  i = nxl, nxr
+                      DO  j = nys, nyn
+                         local_pf(i,j,nzb+1) = ghf_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+
+                resorted = .TRUE.
+                two_d = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
 
              CASE ( 'ol*_xy' )        ! 2d-array