Changeset 3614 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Dec 10, 2018 7:05:46 AM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r3607 r3614 28 28 ! ----------------- 29 29 ! $Id$ 30 ! unused variables removed 31 ! 32 ! 3607 2018-12-07 11:56:58Z suehring 30 33 ! Output of radiation-related quantities migrated to radiation_model_mod. 31 34 ! … … 960 963 961 964 962 !-- arrays for time averages963 !-- Attention: the variable rad_net_av is also used in the 3d field variable in radiation_model_mod.f90. It may be better to rename it964 REAL(wp), DIMENSION(:), ALLOCATABLE :: wghf_eb_av !< average of wghf_eb965 REAL(wp), DIMENSION(:), ALLOCATABLE :: wshf_eb_av !< average of wshf_eb966 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: t_wall_av !< Average of t_wall967 REAL(wp), DIMENSION(:), ALLOCATABLE :: wghf_eb_green_av !< average of wghf_eb_green968 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: t_green_av !< Average of t_green969 REAL(wp), DIMENSION(:), ALLOCATABLE :: wghf_eb_window_av !< average of wghf_eb_window970 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: t_window_av !< Average of t_window971 REAL(wp), DIMENSION(:), ALLOCATABLE :: qsws_eb_av !< average of qsws_eb972 REAL(wp), DIMENSION(:), ALLOCATABLE :: qsws_veg_av !< average of qsws_veg973 REAL(wp), DIMENSION(:), ALLOCATABLE :: qsws_liq_av !< average of qsws_liq974 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: swc_av !< Average of swc975 976 977 965 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 978 966 !-- anthropogenic heat sources … … 1061 1049 1062 1050 #endif 1063 REAL(wp), DIMENSION(:), ALLOCATABLE, TARGET :: t_surf_wall_av !< average of wall surface temperature (K)1064 REAL(wp), DIMENSION(:), ALLOCATABLE, TARGET :: t_surf_window_av !< average of window surface temperature (K)1065 REAL(wp), DIMENSION(:), ALLOCATABLE, TARGET :: t_surf_green_av !< average of green wall surface temperature (K)1066 1067 !-- Temporal tendencies for time stepping1068 REAL(wp), DIMENSION(:), ALLOCATABLE :: tt_surface_wall_m !< surface temperature tendency of wall (K)1069 REAL(wp), DIMENSION(:), ALLOCATABLE :: tt_surface_window_m !< surface temperature tendency of window (K)1070 REAL(wp), DIMENSION(:), ALLOCATABLE :: tt_surface_green_m !< surface temperature tendency of green wall (K)1071 1051 1072 1052 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 1077 1057 #if defined( __nopointer ) 1078 1058 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_wall_h !< Wall temperature (K) 1079 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_wall_h_av !< Average of t_wall1080 1059 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_wall_h_p !< Prog. wall temperature (K) 1081 1060 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_window_h !< Window temperature (K) 1082 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_window_h_av !< Average of t_window1083 1061 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_window_h_p !< Prog. window temperature (K) 1084 1062 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_green_h !< Green temperature (K) 1085 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_green_h_av !< Average of t_green1086 1063 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_green_h_p !< Prog. green temperature (K) 1087 1064 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: swc_h !< soil water content green building layer … … 1096 1073 1097 1074 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_wall_v !< Wall temperature (K) 1098 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_wall_v_av !< Average of t_wall1099 1075 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_wall_v_p !< Prog. wall temperature (K) 1100 1076 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_window_v !< Window temperature (K) 1101 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_window_v_av !< Average of t_window1102 1077 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_window_v_p !< Prog. window temperature (K) 1103 1078 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_green_v !< Green temperature (K) 1104 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_green_v_av !< Average of t_green1105 1079 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_green_v_p !< Prog. green temperature (K) 1106 1080 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: swc_v !< Wall swc 1107 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: swc_v_av !< Average of swc1108 1081 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: swc_v_p !< Prog. swc 1109 1082 1110 1083 #else 1111 1084 REAL(wp), DIMENSION(:,:), POINTER :: t_wall_h, t_wall_h_p 1112 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_wall_h_ av, t_wall_h_1, t_wall_h_21085 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_wall_h_1, t_wall_h_2 1113 1086 REAL(wp), DIMENSION(:,:), POINTER :: t_window_h, t_window_h_p 1114 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_window_h_ av, t_window_h_1, t_window_h_21087 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_window_h_1, t_window_h_2 1115 1088 REAL(wp), DIMENSION(:,:), POINTER :: t_green_h, t_green_h_p 1116 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_green_h_ av, t_green_h_1, t_green_h_21089 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: t_green_h_1, t_green_h_2 1117 1090 REAL(wp), DIMENSION(:,:), POINTER :: swc_h, rootfr_h, wilt_h, fc_h, swc_sat_h, swc_h_p, swc_res_h 1118 1091 REAL(wp), DIMENSION(:,:), ALLOCATABLE, TARGET :: swc_h_1, rootfr_h_1, & … … 1121 1094 1122 1095 TYPE(t_wall_vertical), DIMENSION(:), POINTER :: t_wall_v, t_wall_v_p 1123 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_wall_v_ av, t_wall_v_1, t_wall_v_21096 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_wall_v_1, t_wall_v_2 1124 1097 TYPE(t_wall_vertical), DIMENSION(:), POINTER :: t_window_v, t_window_v_p 1125 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_window_v_ av, t_window_v_1, t_window_v_21098 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_window_v_1, t_window_v_2 1126 1099 TYPE(t_wall_vertical), DIMENSION(:), POINTER :: t_green_v, t_green_v_p 1127 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_green_v_ av, t_green_v_1, t_green_v_21100 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: t_green_v_1, t_green_v_2 1128 1101 TYPE(t_wall_vertical), DIMENSION(:), POINTER :: swc_v, swc_v_p 1129 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: swc_v_ av, swc_v_1, swc_v_21102 TYPE(t_wall_vertical), DIMENSION(0:3), TARGET :: swc_v_1, swc_v_2 1130 1103 #endif 1131 1132 !-- Wall temporal tendencies for time stepping1133 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tt_wall_m !< t_wall prognostic array1134 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tt_window_m !< t_window prognostic array1135 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tt_green_m !< t_green prognostic array1136 1104 1137 1105 !-- Surface and material parameters classes (surface_type) … … 2850 2818 REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: temp_pf !< temp array for urban surface output procedure 2851 2819 2852 CHARACTER (len=varnamelength) :: var , surfid2820 CHARACTER (len=varnamelength) :: var 2853 2821 INTEGER(iwp), PARAMETER :: nd = 5 2854 2822 CHARACTER(len=6), DIMENSION(0:nd-1), PARAMETER :: dirname = (/ '_roof ', '_south', '_north', '_west ', '_east ' /) … … 2856 2824 INTEGER(iwp), DIMENSION(0:nd-1), PARAMETER :: diridx = (/ -1, 1, 0, 3, 2 /) 2857 2825 !< index for surf_*_v: 0:3 = (North, South, East, West) 2858 INTEGER(iwp) :: ids,idsint,idsidx,is urf,isvf,isurfs,isurflt,ipcgb2859 INTEGER(iwp) :: i s,js,ks,i,j,k,iwl,istat, l, m2826 INTEGER(iwp) :: ids,idsint,idsidx,isvf 2827 INTEGER(iwp) :: i,j,k,iwl,istat, l, m 2860 2828 2861 2829 found = .TRUE. … … 7499 7467 INTEGER(iwp), DIMENSION(0:17, nysg:nyng, nxlg:nxrg) :: usm_par 7500 7468 REAL(wp), DIMENSION(1:14, nysg:nyng, nxlg:nxrg) :: usm_val 7501 INTEGER(iwp) :: k, l, d,iw, jw, kw, it, ip, ii, ij, m7469 INTEGER(iwp) :: k, l, iw, jw, kw, it, ip, ii, ij, m 7502 7470 INTEGER(iwp) :: i, j 7503 7471 INTEGER(iwp) :: nz, roof, dirwe, dirsn … … 8145 8113 IMPLICIT NONE 8146 8114 8147 INTEGER(iwp) :: i, j, k, l, d,m !< running indices8115 INTEGER(iwp) :: i, j, k, l, m !< running indices 8148 8116 8149 8117 INTEGER(iwp) :: i_off !< offset to determine index of surface element, seen from atmospheric grid point, for x … … 8153 8121 LOGICAL :: spinup !true during spinup 8154 8122 8155 REAL(wp) :: u1,v1,w1 !< near wall u,v,w8156 8123 REAL(wp) :: stend_wall !< surface tendency 8157 8124
Note: See TracChangeset
for help on using the changeset viewer.