Changeset 4258 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Oct 7, 2019 1:29:08 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r4245 r4258 28 28 ! ----------------- 29 29 ! $Id$ 30 ! - Add checks to ensure that relative fractions of walls, windowns and green 31 ! surfaces sum-up to one. 32 ! - Revise message calls dealing with local checks. 33 ! 34 ! 4245 2019-09-30 08:40:37Z pavelkrc 30 35 ! Initialize explicit per-surface parameters from building_surface_pars 31 36 ! … … 2355 2360 THEN 2356 2361 WRITE( message_string, * ) 'building_type = is out of ' // & 2357 'the valid range at (j,i) = ', j, i2358 CALL message( 'usm_check_parameters', 'PA0529', 2, 2, 0, 6, 0 )2362 'the valid range at (j,i) = ', j, i 2363 CALL message( 'usm_check_parameters', 'PA0529', 2, 2, myid, 6, 0 ) 2359 2364 ENDIF 2360 2365 ENDDO … … 4434 4439 building_pars_f%pars_xy(ind_lambda_surf,j,i) 4435 4440 4441 write(9,*) m, SUM( surf_usm_h%frac(:,m) ), "indiv", surf_usm_h%frac(0,m), surf_usm_h%frac(1,m), surf_usm_h%frac(2,m) 4436 4442 ENDDO 4443 flush(9) 4437 4444 4438 4445 … … 5028 5035 ENDDO 5029 5036 ENDIF 5037 ! 5038 !-- Run further checks to ensure that the respecitve material fractions are 5039 !-- prescribed properly. 5040 DO m = 1, surf_usm_h%ns 5041 IF ( SUM( surf_usm_h%frac(:,m) ) /= 1.0_wp ) THEN 5042 WRITE(message_string,*) 'The relative material fractions do ' // & 5043 'not sum-up to one at horizotal ' // & 5044 'surface. (i,j) = ', & 5045 surf_usm_h%i(m), surf_usm_h%j(m) 5046 CALL message( 'urban_surface_model_mod', 'PA0686', 2, 2, myid, 6, 0 ) 5047 ENDIF 5048 ENDDO 5049 5050 DO l = 0, 3 5051 DO m = 1, surf_usm_v(l)%ns 5052 IF ( SUM( surf_usm_v(l)%frac(:,m) ) /= 1.0_wp ) THEN 5053 WRITE(message_string,*) & 5054 'The relative material fractions do ' // & 5055 'not sum-up to one at vertical ' // & 5056 'surface. (i,j) = ', & 5057 surf_usm_v(l)%i(m), surf_usm_v(l)%j(m) 5058 CALL message( 'urban_surface_model_mod', 'PA0686', 2, 2, myid, 6, 0 ) 5059 ENDIF 5060 ENDDO 5061 ENDDO 5030 5062 ! 5031 5063 !-- Read the surface_types array. … … 5059 5091 surf_usm_h%i(m), surf_usm_h%j(m) 5060 5092 CALL message( 'urban_surface_model_mod', 'PA0503', & 5061 0, 0, 0, 6, 0 )5093 0, 0, myid, 6, 0 ) 5062 5094 ENDIF 5063 5095 IF ( surf_usm_h%z0h(m) >= surf_usm_h%z_mo(m) ) THEN … … 5071 5103 surf_usm_h%i(m), surf_usm_h%j(m) 5072 5104 CALL message( 'urban_surface_model_mod', 'PA0507', & 5073 0, 0, 0, 6, 0 )5105 0, 0, myid, 6, 0 ) 5074 5106 ENDIF 5075 5107 ENDDO … … 5087 5119 surf_usm_v(l)%j(m)+surf_usm_v(l)%joff 5088 5120 CALL message( 'urban_surface_model_mod', 'PA0503', & 5089 0, 0, 0, 6, 0 )5121 0, 0, myid, 6, 0 ) 5090 5122 ENDIF 5091 5123 IF ( surf_usm_v(l)%z0h(m) >= surf_usm_v(l)%z_mo(m) ) THEN … … 5100 5132 surf_usm_v(l)%j(m)+surf_usm_v(l)%joff 5101 5133 CALL message( 'urban_surface_model_mod', 'PA0507', & 5102 0, 0, 0, 6, 0 )5134 0, 0, myid, 6, 0 ) 5103 5135 ENDIF 5104 5136 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.