Changeset 3921 for palm/trunk
- Timestamp:
- Apr 18, 2019 2:21:10 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3918 r3921 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Undo accidentally commented initialization 31 ! 32 ! 3918 2019-04-18 13:33:11Z suehring 30 33 ! Set green fraction to zero also at vertical surfaces 31 34 ! … … 4462 4465 surf_usm_h%lai(m) = building_pars_f%pars_xy(ind_lai_r,j,i) 4463 4466 4464 !IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= &4465 !building_pars_f%fill ) THEN4466 !surf_usm_h%rho_c_wall(nzb_wall,m) = &4467 !building_pars_f%pars_xy(ind_hc1,j,i)4468 !surf_usm_h%rho_c_wall(nzb_wall+1,m) = &4469 !building_pars_f%pars_xy(ind_hc1,j,i)4470 !ENDIF4471 !4472 !4473 !IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= &4474 !building_pars_f%fill ) &4475 !surf_usm_h%rho_c_wall(nzb_wall+2,m) = &4476 !building_pars_f%pars_xy(ind_hc2,j,i)4477 !4478 !IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= &4479 !building_pars_f%fill ) &4480 !surf_usm_h%rho_c_wall(nzb_wall+3,m) = &4481 !building_pars_f%pars_xy(ind_hc3,j,i)4467 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & 4468 building_pars_f%fill ) THEN 4469 surf_usm_h%rho_c_wall(nzb_wall,m) = & 4470 building_pars_f%pars_xy(ind_hc1,j,i) 4471 surf_usm_h%rho_c_wall(nzb_wall+1,m) = & 4472 building_pars_f%pars_xy(ind_hc1,j,i) 4473 ENDIF 4474 4475 4476 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= & 4477 building_pars_f%fill ) & 4478 surf_usm_h%rho_c_wall(nzb_wall+2,m) = & 4479 building_pars_f%pars_xy(ind_hc2,j,i) 4480 4481 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= & 4482 building_pars_f%fill ) & 4483 surf_usm_h%rho_c_wall(nzb_wall+3,m) = & 4484 building_pars_f%pars_xy(ind_hc3,j,i) 4482 4485 4483 4486 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & … … 4498 4501 building_pars_f%pars_xy(ind_hc3,j,i) 4499 4502 4500 !IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= &4501 !building_pars_f%fill ) THEN4502 !surf_usm_h%rho_c_window(nzb_wall,m) = &4503 !building_pars_f%pars_xy(ind_hc1,j,i)4504 !surf_usm_h%rho_c_window(nzb_wall+1,m) = &4505 !building_pars_f%pars_xy(ind_hc1,j,i)4506 !ENDIF4507 !IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= &4508 !building_pars_f%fill ) &4509 !surf_usm_h%rho_c_window(nzb_wall+2,m) = &4510 !building_pars_f%pars_xy(ind_hc2,j,i)4511 !4512 !IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= &4513 !building_pars_f%fill ) &4514 !surf_usm_h%rho_c_window(nzb_wall+3,m) = &4515 !building_pars_f%pars_xy(ind_hc3,j,i)4503 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & 4504 building_pars_f%fill ) THEN 4505 surf_usm_h%rho_c_window(nzb_wall,m) = & 4506 building_pars_f%pars_xy(ind_hc1,j,i) 4507 surf_usm_h%rho_c_window(nzb_wall+1,m) = & 4508 building_pars_f%pars_xy(ind_hc1,j,i) 4509 ENDIF 4510 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= & 4511 building_pars_f%fill ) & 4512 surf_usm_h%rho_c_window(nzb_wall+2,m) = & 4513 building_pars_f%pars_xy(ind_hc2,j,i) 4514 4515 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= & 4516 building_pars_f%fill ) & 4517 surf_usm_h%rho_c_window(nzb_wall+3,m) = & 4518 building_pars_f%pars_xy(ind_hc3,j,i) 4516 4519 4517 4520 IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= & … … 4754 4757 building_pars_f%pars_xy(ind_lai_w,j,i) 4755 4758 4756 !IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= &4757 !building_pars_f%fill ) THEN4758 !surf_usm_v(l)%rho_c_wall(nzb_wall,m) = &4759 !building_pars_f%pars_xy(ind_hc1,j,i)4760 !surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) = &4761 !building_pars_f%pars_xy(ind_hc1,j,i)4762 !ENDIF4763 !4764 !4765 !IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= &4766 !building_pars_f%fill ) &4767 !surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) = &4768 !building_pars_f%pars_xy(ind_hc2,j,i)4769 !4770 !IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= &4771 !building_pars_f%fill ) &4772 !surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) = &4773 !building_pars_f%pars_xy(ind_hc3,j,i)4759 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & 4760 building_pars_f%fill ) THEN 4761 surf_usm_v(l)%rho_c_wall(nzb_wall,m) = & 4762 building_pars_f%pars_xy(ind_hc1,j,i) 4763 surf_usm_v(l)%rho_c_wall(nzb_wall+1,m) = & 4764 building_pars_f%pars_xy(ind_hc1,j,i) 4765 ENDIF 4766 4767 4768 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= & 4769 building_pars_f%fill ) & 4770 surf_usm_v(l)%rho_c_wall(nzb_wall+2,m) = & 4771 building_pars_f%pars_xy(ind_hc2,j,i) 4772 4773 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= & 4774 building_pars_f%fill ) & 4775 surf_usm_v(l)%rho_c_wall(nzb_wall+3,m) = & 4776 building_pars_f%pars_xy(ind_hc3,j,i) 4774 4777 4775 4778 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & … … 4790 4793 building_pars_f%pars_xy(ind_hc3,j,i) 4791 4794 4792 !IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= &4793 !building_pars_f%fill ) THEN4794 !surf_usm_v(l)%rho_c_window(nzb_wall,m) = &4795 !building_pars_f%pars_xy(ind_hc1,j,i)4796 !surf_usm_v(l)%rho_c_window(nzb_wall+1,m) = &4797 !building_pars_f%pars_xy(ind_hc1,j,i)4798 !ENDIF4799 !IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= &4800 !building_pars_f%fill ) &4801 !surf_usm_v(l)%rho_c_window(nzb_wall+2,m) = &4802 !building_pars_f%pars_xy(ind_hc2,j,i)4803 !4804 !IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= &4805 !building_pars_f%fill ) &4806 !surf_usm_v(l)%rho_c_window(nzb_wall+3,m) = &4807 !building_pars_f%pars_xy(ind_hc3,j,i)4795 IF ( building_pars_f%pars_xy(ind_hc1,j,i) /= & 4796 building_pars_f%fill ) THEN 4797 surf_usm_v(l)%rho_c_window(nzb_wall,m) = & 4798 building_pars_f%pars_xy(ind_hc1,j,i) 4799 surf_usm_v(l)%rho_c_window(nzb_wall+1,m) = & 4800 building_pars_f%pars_xy(ind_hc1,j,i) 4801 ENDIF 4802 IF ( building_pars_f%pars_xy(ind_hc2,j,i) /= & 4803 building_pars_f%fill ) & 4804 surf_usm_v(l)%rho_c_window(nzb_wall+2,m) = & 4805 building_pars_f%pars_xy(ind_hc2,j,i) 4806 4807 IF ( building_pars_f%pars_xy(ind_hc3,j,i) /= & 4808 building_pars_f%fill ) & 4809 surf_usm_v(l)%rho_c_window(nzb_wall+3,m) = & 4810 building_pars_f%pars_xy(ind_hc3,j,i) 4808 4811 4809 4812 IF ( building_pars_f%pars_xy(ind_tc1,j,i) /= &
Note: See TracChangeset
for help on using the changeset viewer.