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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Add possibility to initialize surface sensible and latent heat fluxes via
+! a static driver.
+! 
+! 4493 2020-04-10 09:49:43Z pavelkrc
 ! Overwrite u_init, v_init, pt_init, q_init and s_init with hom for all
 ! cyclic_fill-cases, not only for turbulent_inflow = .TRUE.
@@ -1599 +1603 @@
 !--       Initialize roughness length. Note, z0 will be only initialized at    
 !--       default-type surfaces. At natural or urban z0 is implicitly          
-!--       initialized bythe respective parameter lists.                        
+!--       initialized by the respective parameter lists.                        
 !--       Initialize horizontal surface elements.                              
           CALL init_single_surface_properties( surf_def_h(0)%z0,               &
@@ -1624 +1628 @@
        ENDIF
 !
-!--    Additional variables, e.g. shf, qsws, etc, can be initialized the 
+!--    Input surface sensible heat flux. 
+       IF ( check_existence( vars_pids, 'shf' ) )  THEN
+!
+!--       Read _FillValue attribute
+          CALL get_attribute( pids_id, char_fill, tmp_2d%fill,                 &
+                              .FALSE., 'shf' )
+!
+!--       Read variable
+          CALL get_variable( pids_id, 'shf', tmp_2d%var,                       &
+                             nxl, nxr, nys, nyn )
+!
+!--       Initialize heat flux. Note, shf will be only initialized at
+!--       default-type surfaces. At natural or urban shf is implicitly
+!--       initialized by the respective parameter lists.
+!--       Initialize horizontal surface elements.
+          CALL init_single_surface_properties( surf_def_h(0)%shf,              &
+                                               tmp_2d%var,                     &
+                                               surf_def_h(0)%ns,               &
+                                               tmp_2d%fill,                    &
+                                               surf_def_h(0)%i,                &
+                                               surf_def_h(0)%j )
+!
+!--       Initialize heat flux also at vertical surface elements.
+!--       Note, the actual 2D input arrays are only defined on the
+!--       subdomain. Therefore, pass the index arrays with their respective
+!--       offset values.
+          DO  l = 0, 3
+             CALL init_single_surface_properties(                              &
+                                         surf_def_v(l)%shf,                    &
+                                         tmp_2d%var,                           &
+                                         surf_def_v(l)%ns,                     &
+                                         tmp_2d%fill,                          &
+                                         surf_def_v(l)%i + surf_def_v(l)%ioff, &
+                                         surf_def_v(l)%j + surf_def_v(l)%joff )
+          ENDDO
+
+       ENDIF
+!
+!--    Input surface sensible heat flux. 
+       IF ( check_existence( vars_pids, 'qsws' ) )  THEN
+!
+!--       Read _FillValue attribute
+          CALL get_attribute( pids_id, char_fill, tmp_2d%fill,                 &
+                              .FALSE., 'qsws' )
+!
+!--       Read variable
+          CALL get_variable( pids_id, 'qsws', tmp_2d%var,                      &
+                             nxl, nxr, nys, nyn )
+!
+!--       Initialize latent heat flux. Note, qsws will be only initialized at
+!--       default-type surfaces. At natural or urban qsws is implicitly
+!--       initialized by the respective parameter lists.
+!--       Initialize horizontal surface elements.
+          CALL init_single_surface_properties( surf_def_h(0)%qsws,             &
+                                               tmp_2d%var,                     &
+                                               surf_def_h(0)%ns,               &
+                                               tmp_2d%fill,                    &
+                                               surf_def_h(0)%i,                &
+                                               surf_def_h(0)%j )
+!
+!--       Initialize latent heat flux also at vertical surface elements.
+!--       Note, the actual 2D input arrays are only defined on the
+!--       subdomain. Therefore, pass the index arrays with their respective
+!--       offset values.
+          DO  l = 0, 3
+             CALL init_single_surface_properties(                              &
+                                         surf_def_v(l)%qsws,                   &
+                                         tmp_2d%var,                           &
+                                         surf_def_v(l)%ns,                     &
+                                         tmp_2d%fill,                          &
+                                         surf_def_v(l)%i + surf_def_v(l)%ioff, &
+                                         surf_def_v(l)%j + surf_def_v(l)%joff )
+          ENDDO
+
+       ENDIF
+!
+!--    Additional variables, can be initialized the 
 !--    same way.
-                           
+
 !
 !--    Finally, close the input file and deallocate temporary arrays