Changeset 4319 for palm/trunk


Ignore:
Timestamp:
Dec 3, 2019 6:14:41 PM (5 years ago)
Author:
suehring
Message:

Further revise check for surface fractions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/netcdf_data_input_mod.f90

    r4313 r4319  
    2525! -----------------
    2626! $Id$
     27! Further revise check for surface fractions
     28!
     29! 4313 2019-11-27 14:07:00Z suehring
    2730! Checks for surface fractions revised
    2831!
     
    37143717                ENDIF
    37153718!
    3716 !--             Sum of relative fractions must be 1. Note, due to type
    3717 !--             type conversions due to reading, the sum of surface fractions
     3719!--             Sum of relative fractions must be 1. Note, attributed to type
     3720!--             conversions due to reading, the sum of surface fractions
    37183721!--             might be not exactly 1. Hence, the sum is check with a
    37193722!--             tolerance. Later, in the land-surface model, the relative
    3720 !--             fractions are normalized to one.
    3721                 IF ( ANY( surface_fraction_f%frac(:,j,i) /=                    &
    3722                           surface_fraction_f%fill ) )  THEN
     3723!--             fractions are normalized to one. Actually, surface fractions
     3724!--             shall be _FillValue at building grid points, however, in order
     3725!--             to relax this requirement and allow that surface-fraction can
     3726!--             also be zero at these grid points, only perform this check
     3727!--             at locations where some vegetation, pavement or water is defined.
     3728                IF ( vegetation_type_f%var(j,i) /= vegetation_type_f%fill  .OR.&
     3729                     pavement_type_f%var(j,i)   /= pavement_type_f%fill    .OR.&
     3730                     water_type_f%var(j,i)      /= water_type_f%fill )  THEN
    37233731                   IF ( SUM ( surface_fraction_f%frac(0:2,j,i) ) >             &
    37243732                        1.0_wp + 1E-8_wp  .OR.                                 &
     
    37273735                      WRITE( message_string, * )                               &
    37283736                                    'The sum of all land-surface fractions ' //&
    3729                                     'must equal 1. ', i, j
     3737                                    'must equal 1. (i, j) = ', i, j
    37303738                      CALL message( 'netcdf_data_input_mod', 'PA0566',         &
    37313739                                     2, 2, myid, 6, 0 )
Note: See TracChangeset for help on using the changeset viewer.