 Timestamp:
 Dec 3, 2019 6:14:41 PM (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/netcdf_data_input_mod.f90
r4313 r4319 25 25 !  26 26 ! $Id$ 27 ! Further revise check for surface fractions 28 ! 29 ! 4313 20191127 14:07:00Z suehring 27 30 ! Checks for surface fractions revised 28 31 ! … … 3714 3717 ENDIF 3715 3718 ! 3716 ! Sum of relative fractions must be 1. Note, dueto type3717 ! typeconversions due to reading, the sum of surface fractions3719 ! Sum of relative fractions must be 1. Note, attributed to type 3720 ! conversions due to reading, the sum of surface fractions 3718 3721 ! might be not exactly 1. Hence, the sum is check with a 3719 3722 ! tolerance. Later, in the landsurface 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 surfacefraction 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 3723 3731 IF ( SUM ( surface_fraction_f%frac(0:2,j,i) ) > & 3724 3732 1.0_wp + 1E8_wp .OR. & … … 3727 3735 WRITE( message_string, * ) & 3728 3736 'The sum of all landsurface fractions ' //& 3729 'must equal 1. ', i, j3737 'must equal 1. (i, j) = ', i, j 3730 3738 CALL message( 'netcdf_data_input_mod', 'PA0566', & 3731 3739 2, 2, myid, 6, 0 )
Note: See TracChangeset
for help on using the changeset viewer.