Changeset 4408 for palm/trunk/SOURCE/virtual_measurement_mod.f90

Timestamp:

Feb 14, 2020 10:04:39 AM (4 years ago)

Author:

gronemeier

Message:

write fill_value attribute in virtual-measurements module; enable character-array output in data-output module

File:

• palm/trunk/SOURCE/virtual_measurement_mod.f90 (modified) (14 diffs)

palm/trunk/SOURCE/virtual_measurement_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! write fill_value attribute
+!
+! 4406 2020-02-13 20:06:29Z knoop
 ! Bugix: removed oro_rel wrong loop bounds and removed unnecessary restart method
-! 
+!
 ! 4400 2020-02-10 20:32:41Z suehring
 ! Revision of the module:
@@ -268 +271 @@
        INTEGER(iwp) ::  start_coord_s = 0   !< start coordinate in NetCDF file for local soil observations
 
-       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  dim_t !< number observations individual for each trajectory 
+       INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  dim_t !< number observations individual for each trajectory
                                                          !< or station that are no _FillValues
 
@@ -287 +290 @@
        REAL(wp) ::  fill_nutm                            !< fill value for UTM coordinates in case of missing values
        REAL(wp) ::  fill_zar                             !< fill value for heigth coordinates in case of missing values
-       REAL(wp) ::  fillout = -9999.0                    !< fill value for output in case a observation is taken 
+       REAL(wp) ::  fillout = -9999.0                    !< fill value for output in case a observation is taken
                                                          !< e.g. from inside a building
        REAL(wp) ::  origin_x_obs                         !< origin of the observation in UTM coordiates in x-direction
@@ -306 +309 @@
 
     ! This need to be generalized
-!     CHARACTER(LEN=10) ::  char_fill = '_FillValue'
+    CHARACTER(LEN=10) ::  char_fill = '_FillValue'                 !< attribute name for fill value
     CHARACTER(LEN=9)  ::  char_long = 'long_name'                  !< attribute name for long_name
     CHARACTER(LEN=13) ::  char_standard = 'standard_name'          !< attribute name for standard_name
@@ -462 +465 @@
 
  END SUBROUTINE vm_check_parameters
- 
+
 !------------------------------------------------------------------------------!
 ! Description:
@@ -893 +896 @@
 !
 !-- Allocate flag array. This dummy array is used to identify grid points
-!-- where virtual measurements should be taken. Please note, in order to 
+!-- where virtual measurements should be taken. Please note, in order to
 !-- include also the surrounding grid points of the original coordinate
 !-- ghost points are required.
@@ -1000 +1003 @@
        ALLOCATE( vmea(l)%var_atts(1:vmea(l)%nmeas) )
 !
-!--    Store the variable names in a data structures, which assigns further 
-!--    attributes to this name. Further, for data output reasons, create a 
+!--    Store the variable names in a data structures, which assigns further
+!--    attributes to this name. Further, for data output reasons, create a
 !--    string of output variables, which will be written into the attribute
 !--    data_content.
@@ -1014 +1017 @@
 !--    define the grid-index space on each subdomain where virtual measurements
 !--    should be taken. Note, the entire coordinate array (on the entire model
-!--    domain) won't be stored as this would exceed memory requirements, 
+!--    domain) won't be stored as this would exceed memory requirements,
 !--    particularly for trajectories.
        IF ( vmea(l)%nmeas > 0 )  THEN
@@ -1344 +1347 @@
     DEALLOCATE( meas_flag )
 !
-!-- Close input file for virtual measurements. 
+!-- Close input file for virtual measurements.
     CALL close_input_file( pids_id )
 !
@@ -1398 +1401 @@
 !-- Note, start coordinates are initialized with zero for sake of simplicity
 !-- in summation. However, in NetCDF the start coordinates must be >= 1,
-!-- so that a one needs to be added at the end. 
+!-- so that a one needs to be added at the end.
     DO  l = 1, vmea_general%nvm
        DO  n  = 0, myid - 1
@@ -2237 +2240 @@
                                       value = TRIM( vmea(l)%var_atts(n)%coordinates ) )
 
-!           return_value = dom_def_att( vmea(l)%nc_filename,                     &
-!                                       variable_name = variable_name,           &
-!                                       attribute_name = char_fill,              &
-!                                       value = vmea(l)%var_atts(n)%fill_value )
+          return_value = dom_def_att( vmea(l)%nc_filename,                     &
+                                      variable_name = variable_name,           &
+                                      attribute_name = char_fill,              &
+                                      value = REAL( vmea(l)%var_atts(n)%fill_value, KIND=4 ) )
 
        ENDDO  ! loop over variables per site
@@ -2345 +2348 @@
           DEALLOCATE( output_values_1d_target )
 !
-!--       In case of sampled soil quantities, output also the respective 
+!--       In case of sampled soil quantities, output also the respective
 !--       coordinate arrays.
           IF ( vmea(l)%soil_sampling )  THEN
@@ -2412 +2415 @@
 !
 !--          Write the stations name
-             
+
           ENDIF
 
@@ -2910 +2913 @@
 !                     i = MERGE( vmea(l)%i(m), nxl, vmea(l)%i(m) > nxl )
 !                     i = MERGE( i           , nxr, i            < nxr )
-! 
+!
 !                     DO  mm = surf_lsm_h%start_index(j,i),                      &
 !                              surf_lsm_h%end_index(j,i)