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

Timestamp:

Apr 30, 2020 4:29:59 PM (4 years ago)

Author:

suehring

Message:

Bugfix in plant-canopy model for output of averaged transpiration rate after a restart; Revise check for output for plant heating rate and rename error message number; Surface-data output: enable output of mixing ratio and passive scalar concentration at the surface; Surface-data input: Add possibility to prescribe surface sensible and latent heat fluxes via static input file

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable output of qsurf and ssurf
+! 
+! 4500 2020-04-17 10:12:45Z suehring
 ! Unify output conversion of sensible and latent heat flux
 ! 
@@ -123 +126 @@
     USE indices,                                                               &
         ONLY:  nbgp, nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg,       &
-               nzb, nzt, wall_flags_total_0
+               nzb, nzt,                                                       &
+               topo_top_ind,                                                   &
+               wall_flags_total_0
 
     USE kinds
@@ -701 +706 @@
                 IF ( mode == 'xy' )  level_z = zu
 
+             CASE ( 'qsurf*_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)%q_surface(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%q_surface(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%q_surface(m)
+                   ENDDO
+
+                ELSE
+                   IF ( .NOT. ALLOCATED( qsurf_av ) ) THEN
+                      ALLOCATE( qsurf_av(nysg:nyng,nxlg:nxrg) )
+                      qsurf_av = REAL( fill_value, KIND = wp )
+                   ENDIF
+                   DO  i = nxl, nxr
+                      DO  j = nys, nyn
+                         local_pf(i,j,nzb+1) = qsurf_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted       = .TRUE.
+                two_d          = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
              CASE ( 'qsws*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN
@@ -845 +885 @@
                 two_d = .TRUE.
                 level_z(nzb+1) = zu(nzb+1)
-                
+
+             CASE ( 'ssurf*_xy' )        ! 2d-array
+                IF ( av == 0 )  THEN
+                   DO  i = nxl, nxr
+                      DO  j = nys, nyn
+                         k = topo_top_ind(j,i,0)
+                         local_pf(i,j,nzb+1) = s(k+k,j,i)
+                      ENDDO
+                   ENDDO
+                ELSE
+                   IF ( .NOT. ALLOCATED( ssurf_av ) ) THEN
+                      ALLOCATE( ssurf_av(nysg:nyng,nxlg:nxrg) )
+                      ssurf_av = REAL( fill_value, KIND = wp )
+                   ENDIF
+                   DO  i = nxl, nxr
+                      DO  j = nys, nyn
+                         local_pf(i,j,nzb+1) = ssurf_av(j,i)
+                      ENDDO
+                   ENDDO
+                ENDIF
+                resorted       = .TRUE.
+                two_d          = .TRUE.
+                level_z(nzb+1) = zu(nzb+1)
+
              CASE ( 'ssws*_xy' )        ! 2d-array
                 IF ( av == 0 ) THEN