Changeset 2512 for palm/trunk/SOURCE/radiation_model_mod.f90
- Timestamp:
- Oct 4, 2017 8:26:59 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/radiation_model_mod.f90
r2504 r2512 25 25 ! ----------------- 26 26 ! $Id$ 27 ! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny 28 ! no output of ghost layer data 29 ! 30 ! 2504 2017-09-27 10:36:13Z maronga 27 31 ! Updates pavement types and albedo parameters 28 32 ! … … 2810 2814 LOGICAL :: two_d !< flag parameter that indicates 2D variables (horizontal cross sections) 2811 2815 2812 REAL(wp), DIMENSION(nxl g:nxrg,nysg:nyng,nzb:nzt+1) :: local_pf !<2816 REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) :: local_pf !< 2813 2817 2814 2818 found = .TRUE. … … 2818 2822 CASE ( 'rad_net*_xy' ) ! 2d-array 2819 2823 IF ( av == 0 ) THEN 2820 DO i = nxl g, nxrg2821 DO j = nys g, nyng2824 DO i = nxl, nxr 2825 DO j = nys, nyn 2822 2826 local_pf(i,j,nzb+1) = rad_net(j,i) 2823 2827 ENDDO 2824 2828 ENDDO 2825 2829 ELSE 2826 DO i = nxl g, nxrg2827 DO j = nys g, nyng2830 DO i = nxl, nxr 2831 DO j = nys, nyn 2828 2832 local_pf(i,j,nzb+1) = rad_net_av(j,i) 2829 2833 ENDDO … … 2836 2840 CASE ( 'rad_lw_in_xy', 'rad_lw_in_xz', 'rad_lw_in_yz' ) 2837 2841 IF ( av == 0 ) THEN 2838 DO i = nxl g, nxrg2839 DO j = nys g, nyng2842 DO i = nxl, nxr 2843 DO j = nys, nyn 2840 2844 DO k = nzb, nzt+1 2841 2845 local_pf(i,j,k) = rad_lw_in(k,j,i) … … 2844 2848 ENDDO 2845 2849 ELSE 2846 DO i = nxl g, nxrg2847 DO j = nys g, nyng2850 DO i = nxl, nxr 2851 DO j = nys, nyn 2848 2852 DO k = nzb, nzt+1 2849 2853 local_pf(i,j,k) = rad_lw_in_av(k,j,i) … … 2856 2860 CASE ( 'rad_lw_out_xy', 'rad_lw_out_xz', 'rad_lw_out_yz' ) 2857 2861 IF ( av == 0 ) THEN 2858 DO i = nxl g, nxrg2859 DO j = nys g, nyng2862 DO i = nxl, nxr 2863 DO j = nys, nyn 2860 2864 DO k = nzb, nzt+1 2861 2865 local_pf(i,j,k) = rad_lw_out(k,j,i) … … 2864 2868 ENDDO 2865 2869 ELSE 2866 DO i = nxl g, nxrg2867 DO j = nys g, nyng2870 DO i = nxl, nxr 2871 DO j = nys, nyn 2868 2872 DO k = nzb, nzt+1 2869 2873 local_pf(i,j,k) = rad_lw_out_av(k,j,i) … … 2876 2880 CASE ( 'rad_lw_cs_hr_xy', 'rad_lw_cs_hr_xz', 'rad_lw_cs_hr_yz' ) 2877 2881 IF ( av == 0 ) THEN 2878 DO i = nxl g, nxrg2879 DO j = nys g, nyng2882 DO i = nxl, nxr 2883 DO j = nys, nyn 2880 2884 DO k = nzb, nzt+1 2881 2885 local_pf(i,j,k) = rad_lw_cs_hr(k,j,i) … … 2884 2888 ENDDO 2885 2889 ELSE 2886 DO i = nxl g, nxrg2887 DO j = nys g, nyng2890 DO i = nxl, nxr 2891 DO j = nys, nyn 2888 2892 DO k = nzb, nzt+1 2889 2893 local_pf(i,j,k) = rad_lw_cs_hr_av(k,j,i) … … 2896 2900 CASE ( 'rad_lw_hr_xy', 'rad_lw_hr_xz', 'rad_lw_hr_yz' ) 2897 2901 IF ( av == 0 ) THEN 2898 DO i = nxl g, nxrg2899 DO j = nys g, nyng2902 DO i = nxl, nxr 2903 DO j = nys, nyn 2900 2904 DO k = nzb, nzt+1 2901 2905 local_pf(i,j,k) = rad_lw_hr(k,j,i) … … 2904 2908 ENDDO 2905 2909 ELSE 2906 DO i = nxl g, nxrg2907 DO j = nys g, nyng2910 DO i = nxl, nxr 2911 DO j = nys, nyn 2908 2912 DO k = nzb, nzt+1 2909 2913 local_pf(i,j,k) = rad_lw_hr_av(k,j,i) … … 2916 2920 CASE ( 'rad_sw_in_xy', 'rad_sw_in_xz', 'rad_sw_in_yz' ) 2917 2921 IF ( av == 0 ) THEN 2918 DO i = nxl g, nxrg2919 DO j = nys g, nyng2922 DO i = nxl, nxr 2923 DO j = nys, nyn 2920 2924 DO k = nzb, nzt+1 2921 2925 local_pf(i,j,k) = rad_sw_in(k,j,i) … … 2924 2928 ENDDO 2925 2929 ELSE 2926 DO i = nxl g, nxrg2927 DO j = nys g, nyng2930 DO i = nxl, nxr 2931 DO j = nys, nyn 2928 2932 DO k = nzb, nzt+1 2929 2933 local_pf(i,j,k) = rad_sw_in_av(k,j,i) … … 2936 2940 CASE ( 'rad_sw_out_xy', 'rad_sw_out_xz', 'rad_sw_out_yz' ) 2937 2941 IF ( av == 0 ) THEN 2938 DO i = nxl g, nxrg2939 DO j = nys g, nyng2942 DO i = nxl, nxr 2943 DO j = nys, nyn 2940 2944 DO k = nzb, nzt+1 2941 2945 local_pf(i,j,k) = rad_sw_out(k,j,i) … … 2944 2948 ENDDO 2945 2949 ELSE 2946 DO i = nxl g, nxrg2947 DO j = nys g, nyng2950 DO i = nxl, nxr 2951 DO j = nys, nyn 2948 2952 DO k = nzb, nzt+1 2949 2953 local_pf(i,j,k) = rad_sw_out_av(k,j,i) … … 2956 2960 CASE ( 'rad_sw_cs_hr_xy', 'rad_sw_cs_hr_xz', 'rad_sw_cs_hr_yz' ) 2957 2961 IF ( av == 0 ) THEN 2958 DO i = nxl g, nxrg2959 DO j = nys g, nyng2962 DO i = nxl, nxr 2963 DO j = nys, nyn 2960 2964 DO k = nzb, nzt+1 2961 2965 local_pf(i,j,k) = rad_sw_cs_hr(k,j,i) … … 2964 2968 ENDDO 2965 2969 ELSE 2966 DO i = nxl g, nxrg2967 DO j = nys g, nyng2970 DO i = nxl, nxr 2971 DO j = nys, nyn 2968 2972 DO k = nzb, nzt+1 2969 2973 local_pf(i,j,k) = rad_sw_cs_hr_av(k,j,i) … … 2976 2980 CASE ( 'rad_sw_hr_xy', 'rad_sw_hr_xz', 'rad_sw_hr_yz' ) 2977 2981 IF ( av == 0 ) THEN 2978 DO i = nxl g, nxrg2979 DO j = nys g, nyng2982 DO i = nxl, nxr 2983 DO j = nys, nyn 2980 2984 DO k = nzb, nzt+1 2981 2985 local_pf(i,j,k) = rad_sw_hr(k,j,i) … … 2984 2988 ENDDO 2985 2989 ELSE 2986 DO i = nxl g, nxrg2987 DO j = nys g, nyng2990 DO i = nxl, nxr 2991 DO j = nys, nyn 2988 2992 DO k = nzb, nzt+1 2989 2993 local_pf(i,j,k) = rad_sw_hr_av(k,j,i) … … 3028 3032 LOGICAL :: found !< 3029 3033 3030 REAL(sp), DIMENSION(nxl g:nxrg,nysg:nyng,nzb:nzt+1) :: local_pf !<3034 REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) :: local_pf !< 3031 3035 3032 3036 … … 3038 3042 CASE ( 'rad_sw_in' ) 3039 3043 IF ( av == 0 ) THEN 3040 DO i = nxl g, nxrg3041 DO j = nys g, nyng3044 DO i = nxl, nxr 3045 DO j = nys, nyn 3042 3046 DO k = nzb, nzt+1 3043 3047 local_pf(i,j,k) = rad_sw_in(k,j,i) … … 3046 3050 ENDDO 3047 3051 ELSE 3048 DO i = nxl g, nxrg3049 DO j = nys g, nyng3052 DO i = nxl, nxr 3053 DO j = nys, nyn 3050 3054 DO k = nzb, nzt+1 3051 3055 local_pf(i,j,k) = rad_sw_in_av(k,j,i) … … 3057 3061 CASE ( 'rad_sw_out' ) 3058 3062 IF ( av == 0 ) THEN 3059 DO i = nxl g, nxrg3060 DO j = nys g, nyng3063 DO i = nxl, nxr 3064 DO j = nys, nyn 3061 3065 DO k = nzb, nzt+1 3062 3066 local_pf(i,j,k) = rad_sw_out(k,j,i) … … 3065 3069 ENDDO 3066 3070 ELSE 3067 DO i = nxl g, nxrg3068 DO j = nys g, nyng3071 DO i = nxl, nxr 3072 DO j = nys, nyn 3069 3073 DO k = nzb, nzt+1 3070 3074 local_pf(i,j,k) = rad_sw_out_av(k,j,i) … … 3076 3080 CASE ( 'rad_sw_cs_hr' ) 3077 3081 IF ( av == 0 ) THEN 3078 DO i = nxl g, nxrg3079 DO j = nys g, nyng3082 DO i = nxl, nxr 3083 DO j = nys, nyn 3080 3084 DO k = nzb, nzt+1 3081 3085 local_pf(i,j,k) = rad_sw_cs_hr(k,j,i) … … 3084 3088 ENDDO 3085 3089 ELSE 3086 DO i = nxl g, nxrg3087 DO j = nys g, nyng3090 DO i = nxl, nxr 3091 DO j = nys, nyn 3088 3092 DO k = nzb, nzt+1 3089 3093 local_pf(i,j,k) = rad_sw_cs_hr_av(k,j,i) … … 3095 3099 CASE ( 'rad_sw_hr' ) 3096 3100 IF ( av == 0 ) THEN 3097 DO i = nxl g, nxrg3098 DO j = nys g, nyng3101 DO i = nxl, nxr 3102 DO j = nys, nyn 3099 3103 DO k = nzb, nzt+1 3100 3104 local_pf(i,j,k) = rad_sw_hr(k,j,i) … … 3103 3107 ENDDO 3104 3108 ELSE 3105 DO i = nxl g, nxrg3106 DO j = nys g, nyng3109 DO i = nxl, nxr 3110 DO j = nys, nyn 3107 3111 DO k = nzb, nzt+1 3108 3112 local_pf(i,j,k) = rad_sw_hr_av(k,j,i) … … 3114 3118 CASE ( 'rad_lw_in' ) 3115 3119 IF ( av == 0 ) THEN 3116 DO i = nxl g, nxrg3117 DO j = nys g, nyng3120 DO i = nxl, nxr 3121 DO j = nys, nyn 3118 3122 DO k = nzb, nzt+1 3119 3123 local_pf(i,j,k) = rad_lw_in(k,j,i) … … 3122 3126 ENDDO 3123 3127 ELSE 3124 DO i = nxl g, nxrg3125 DO j = nys g, nyng3128 DO i = nxl, nxr 3129 DO j = nys, nyn 3126 3130 DO k = nzb, nzt+1 3127 3131 local_pf(i,j,k) = rad_lw_in_av(k,j,i) … … 3133 3137 CASE ( 'rad_lw_out' ) 3134 3138 IF ( av == 0 ) THEN 3135 DO i = nxl g, nxrg3136 DO j = nys g, nyng3139 DO i = nxl, nxr 3140 DO j = nys, nyn 3137 3141 DO k = nzb, nzt+1 3138 3142 local_pf(i,j,k) = rad_lw_out(k,j,i) … … 3141 3145 ENDDO 3142 3146 ELSE 3143 DO i = nxl g, nxrg3144 DO j = nys g, nyng3147 DO i = nxl, nxr 3148 DO j = nys, nyn 3145 3149 DO k = nzb, nzt+1 3146 3150 local_pf(i,j,k) = rad_lw_out_av(k,j,i) … … 3152 3156 CASE ( 'rad_lw_cs_hr' ) 3153 3157 IF ( av == 0 ) THEN 3154 DO i = nxl g, nxrg3155 DO j = nys g, nyng3158 DO i = nxl, nxr 3159 DO j = nys, nyn 3156 3160 DO k = nzb, nzt+1 3157 3161 local_pf(i,j,k) = rad_lw_cs_hr(k,j,i) … … 3160 3164 ENDDO 3161 3165 ELSE 3162 DO i = nxl g, nxrg3163 DO j = nys g, nyng3166 DO i = nxl, nxr 3167 DO j = nys, nyn 3164 3168 DO k = nzb, nzt+1 3165 3169 local_pf(i,j,k) = rad_lw_cs_hr_av(k,j,i) … … 3171 3175 CASE ( 'rad_lw_hr' ) 3172 3176 IF ( av == 0 ) THEN 3173 DO i = nxl g, nxrg3174 DO j = nys g, nyng3177 DO i = nxl, nxr 3178 DO j = nys, nyn 3175 3179 DO k = nzb, nzt+1 3176 3180 local_pf(i,j,k) = rad_lw_hr(k,j,i) … … 3179 3183 ENDDO 3180 3184 ELSE 3181 DO i = nxl g, nxrg3182 DO j = nys g, nyng3185 DO i = nxl, nxr 3186 DO j = nys, nyn 3183 3187 DO k = nzb, nzt+1 3184 3188 local_pf(i,j,k) = rad_lw_hr_av(k,j,i)
Note: See TracChangeset
for help on using the changeset viewer.