Changeset 2963 for palm/trunk/SOURCE/netcdf_data_input_mod.f90

Timestamp:

Apr 12, 2018 2:47:44 PM (6 years ago)

Author:

suehring

Message:

Minor revision of static input file checks, bugfix in initialization of surface-fractions in LSM; minor bugfix in initialization of albedo at window-surfaces; for clearer access of albedo and emissivity introduce index for vegetation/wall, pavement/green-wall and water/window surfaces

File:

• palm/trunk/SOURCE/netcdf_data_input_mod.f90 (modified) (7 diffs)

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - Revise checks for static input variables.
+! - Introduce index for vegetation/wall, pavement/green-wall and water/window 
+!   surfaces, for clearer access of surface fraction, albedo, emissivity, etc. .
+! 
+! 2958 2018-04-11 15:38:13Z suehring
 ! Synchronize longitude and latitude between nested model domains, values are
 ! taken from the root model.
@@ -126 +131 @@
 
     USE pegrid
+
+    USE surface_mod,                                                           &
+        ONLY:  ind_pav_green, ind_veg_wall, ind_wat_win
 !
 !-- Define type for dimensions. 
@@ -2980 +2988 @@
 !--    More advanced methods, e.g. a nearest neighbor algorithm as used in GIS
 !--    systems might be implemented later. 
+!--    Please note, if no terrain height is provided, it is set to 0. 
        IF ( ANY( terrain_height_f%var == terrain_height_f%fill ) )  THEN
           message_string = 'NetCDF variable orography_2D is not ' //           &
@@ -3012 +3021 @@
 !--    Check for minimum requirement of surface-classification data in case
 !--    static input file is used.
-       IF ( .NOT. vegetation_type_f%from_file  .OR.                            &
-            .NOT. pavement_type_f%from_file    .OR.                            &
-            .NOT. building_type_f%from_file    .OR.                            &
-            .NOT. water_type_f%from_file       .OR.                            &
-            .NOT. soil_type_f%from_file             )  THEN
+       IF ( ( .NOT. vegetation_type_f%from_file  .OR.                          &
+              .NOT. pavement_type_f%from_file    .OR.                          &
+              .NOT. water_type_f%from_file       .OR.                          &
+              .NOT. soil_type_f%from_file             ) .OR.                   &
+             ( urban_surface  .AND.  .NOT. building_type_f%from_file ) )  THEN
           message_string = 'Minimum requirement for surface classification ' //&
                            'is not fulfilled. At least ' //                    &
                            'vegetation_type, pavement_type, ' //               &
-                           'building_type, soil_type and water_type are '//    &
-                           'required.'
+                           'soil_type and water_type are '//                   &
+                           'required. If urban-surface model is applied, ' //  &
+                           'also building_type ist required' 
           CALL message( 'netcdf_data_input_mod', 'NDI016', 1, 2, 0, 6, 0 )
        ENDIF
@@ -3129 +3139 @@
                   building_type_f%var(j,i)   == building_type_f%fill    .AND.  &
                   water_type_f%var(j,i)      == water_type_f%fill )  THEN
-                message_string = 'At least one of the paramters '       //     &
+                WRITE( message_string, * ) 'At least one of the parameters '// &
                                  'vegetation_type, pavement_type, '     //     &
                                  'building_type, or water_type must be set '// &
-                                 'to a non-missing value'
+                                 'to a non-missing value. Grid point: ', j, i
                 CALL message( 'netcdf_data_input_mod', 'NDI025', 2, 2, 0, 6, 0 )
              ENDIF
@@ -3188 +3198 @@
 !--          surface vegetation fraction is zero. 
              IF ( ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill  .AND.&
-                 ( surface_fraction_f%frac(0,j,i) == 0.0_wp .OR.               &
-                   surface_fraction_f%frac(0,j,i) == surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_veg_wall,j,i) == 0.0_wp .OR.    &
+                   surface_fraction_f%frac(ind_veg_wall,j,i) ==                &
+                                                     surface_fraction_f%fill ) &
                   )  .OR.                                                      &
                   ( pavement_type_f%var(j,i) /= pavement_type_f%fill     .AND. &
-                 ( surface_fraction_f%frac(1,j,i) == 0.0_wp .OR.               &
-                   surface_fraction_f%frac(1,j,i) == surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_pav_green,j,i) == 0.0_wp .OR.   &
+                   surface_fraction_f%frac(ind_pav_green,j,i) ==               &
+                                                     surface_fraction_f%fill ) &
                   )  .OR.                                                      &
                   ( water_type_f%var(j,i) /= water_type_f%fill           .AND. &
-                 ( surface_fraction_f%frac(2,j,i) == 0.0_wp .OR.               &
-                   surface_fraction_f%frac(2,j,i) == surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_wat_win,j,i) == 0.0_wp .OR.     &
+                   surface_fraction_f%frac(ind_wat_win,j,i) ==                 &
+                                                     surface_fraction_f%fill ) &
                   ) )  THEN 
                 WRITE( message_string, * ) 'Mismatch in setting of '     //    &
@@ -3210 +3223 @@
 !--          surface vegetation fraction is non-zero. 
              IF ( ( vegetation_type_f%var(j,i) == vegetation_type_f%fill  .AND.&
-                 ( surface_fraction_f%frac(0,j,i) /= 0.0_wp .AND.              &
-                   surface_fraction_f%frac(0,j,i) /= surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_veg_wall,j,i) /= 0.0_wp .AND.   &
+                   surface_fraction_f%frac(ind_veg_wall,j,i) /=                &
+                                                     surface_fraction_f%fill ) &
                   )  .OR.                                                      &
                   ( pavement_type_f%var(j,i) == pavement_type_f%fill     .AND. &
-                 ( surface_fraction_f%frac(1,j,i) /= 0.0_wp .AND.              &
-                   surface_fraction_f%frac(1,j,i) /= surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_pav_green,j,i) /= 0.0_wp .AND.  &
+                   surface_fraction_f%frac(ind_pav_green,j,i) /=               &
+                                                     surface_fraction_f%fill ) &
                   )  .OR.                                                      &
                   ( water_type_f%var(j,i) == water_type_f%fill           .AND. &
-                 ( surface_fraction_f%frac(2,j,i) /= 0.0_wp .AND.              &
-                   surface_fraction_f%frac(2,j,i) /= surface_fraction_f%fill ) &
+                 ( surface_fraction_f%frac(ind_wat_win,j,i) /= 0.0_wp .AND.    &
+                   surface_fraction_f%frac(ind_wat_win,j,i) /=                 &
+                                                     surface_fraction_f%fill ) &
                   ) )  THEN 
                 WRITE( message_string, * ) 'Mismatch in setting of '     //    &