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

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Enable output of qsurf and ssurf
+! 
+! 4442 2020-03-04 19:21:13Z suehring
 ! Change order of dimension in surface array %frac to allow for better 
 ! vectorization.
@@ -77 +80 @@
     USE averaging,                                                             &
         ONLY:  e_av, ghf_av, lpt_av, lwp_av, ol_av, p_av, pc_av, pr_av, pt_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,    &
+               q_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qsurf_av, qsws_av,     &
+               qv_av, r_a_av, s_av, shf_av, ssurf_av,                          &
+               ssws_av, ts_av, tsurf_av, u_av,                                 &
                us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
 
@@ -95 +99 @@
 
     USE indices,                                                               &
-        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt 
+        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,           &
+               topo_top_ind
 
     USE kinds
@@ -232 +237 @@
                 ql_vp_av = 0.0_wp
 
+             CASE ( 'qsurf*' )
+                IF ( .NOT. ALLOCATED( qsurf_av ) )  THEN
+                   ALLOCATE( qsurf_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                qsurf_av = 0.0_wp
+
              CASE ( 'qsws*' )
                 IF ( .NOT. ALLOCATED( qsws_av ) )  THEN
@@ -261 +272 @@
                 ENDIF
                 shf_av = 0.0_wp
-                
+
+             CASE ( 'ssurf*' )
+                IF ( .NOT. ALLOCATED( ssurf_av ) )  THEN
+                   ALLOCATE( ssurf_av(nysg:nyng,nxlg:nxrg) )
+                ENDIF
+                ssurf_av = 0.0_wp
+
              CASE ( 'ssws*' )
                 IF ( .NOT. ALLOCATED( ssws_av ) )  THEN
@@ -587 +604 @@
              ENDIF
 
+          CASE ( 'qsurf*' )
+             IF ( ALLOCATED( qsurf_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)
+                         qsurf_av(j,i) = qsurf_av(j,i) +                       &
+                                         surf_def_h(0)%q_surface(m)
+                      ELSEIF ( match_lsm  .AND.  .NOT. match_usm )  THEN
+                         m = surf_lsm_h%end_index(j,i)
+                         qsurf_av(j,i) = qsurf_av(j,i) +                       &
+                                         surf_lsm_h%q_surface(m)
+                      ELSEIF ( match_usm )  THEN
+                         m = surf_usm_h%end_index(j,i)
+                         qsurf_av(j,i) = qsurf_av(j,i) +                       &
+                                         surf_usm_h%q_surface(m)
+                      ENDIF
+                   ENDDO
+                ENDDO
+             ENDIF
+
           CASE ( 'qsws*' )
 !
@@ -703 +748 @@
              ENDIF
 
+          CASE ( 'ssurf*' )
+             IF ( ALLOCATED( ssurf_av ) ) THEN
+                DO  i = nxl, nxr
+                   DO  j = nys, nyn
+                      k = topo_top_ind(j,i,0)
+                      ssurf_av(j,i) = ssurf_av(j,i) + s(k,j,i)
+                   ENDDO
+                ENDDO
+             ENDIF
+
           CASE ( 'ssws*' )
              IF ( ALLOCATED( ssws_av ) ) THEN