Changeset 1382 for palm/trunk/SOURCE/flow_statistics.f90
- Timestamp:
- Apr 30, 2014 12:15:41 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/flow_statistics.f90
r1375 r1382 21 21 ! Current revisions: 22 22 ! ----------------- 23 ! 23 ! Renamed variables which store large scale forcing tendencies 24 ! pt_lsa -> td_lsa_lpt, pt_subs -> td_sub_lpt, 25 ! q_lsa -> td_lsa_q, q_subs -> td_sub_q, 26 ! added Neumann boundary conditions for profile data output of large scale 27 ! advection and subsidence terms at nzt+1 24 28 ! 25 29 ! Former revisions: … … 103 107 104 108 USE arrays_3d, & 105 ONLY: ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, pt_lsa, pt_subs, q,&106 q c, ql, qr, qs, qsws, qswst, q_lsa, q_subs, rho, sa, saswsb,&107 saswst, shf, time_vert, ts, tswst, u, ug, us, usws, uswst, vsws,&108 v, vg, vpt, vswst, w, w_subs, zw109 ONLY: ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, q, qc, ql, qr, qs, & 110 qsws, qswst, rho, sa, saswsb, saswst, shf, td_lsa_lpt, td_lsa_q,& 111 td_sub_lpt, td_sub_q, time_vert, ts, tswst, u, ug, us, usws, & 112 uswst, vsws, v, vg, vpt, vswst, w, w_subs, zw 109 113 110 114 USE cloud_parameters, & … … 1014 1018 1015 1019 DO k = nzb, nzt 1016 sums_ls_l(k,0) = pt_lsa(k,nt) & 1017 + fac * ( pt_lsa(k,nt+1) - pt_lsa(k,nt) ) 1018 sums_ls_l(k,1) = q_lsa(k,nt) & 1019 + fac * ( q_lsa(k,nt+1) - q_lsa(k,nt) ) 1020 ENDDO 1020 sums_ls_l(k,0) = td_lsa_lpt(k,nt) & 1021 + fac * ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) ) 1022 sums_ls_l(k,1) = td_lsa_q(k,nt) & 1023 + fac * ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) 1024 ENDDO 1025 1026 sums_ls_l(nzt+1,0) = sums_ls_l(nzt,0) 1027 sums_ls_l(nzt+1,1) = sums_ls_l(nzt,1) 1021 1028 1022 1029 IF ( large_scale_subsidence .AND. use_subsidence_tendencies ) THEN 1023 1030 1024 1031 DO k = nzb, nzt 1025 sums_ls_l(k,2) = pt_subs(k,nt) & 1026 + fac * ( pt_subs(k,nt+1) - pt_subs(k,nt) ) 1027 sums_ls_l(k,3) = q_subs(k,nt) & 1028 + fac * ( q_subs(k,nt+1) - q_subs(k,nt) ) 1029 ENDDO 1032 sums_ls_l(k,2) = td_sub_lpt(k,nt) + fac * & 1033 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) ) 1034 sums_ls_l(k,3) = td_sub_q(k,nt) + fac * & 1035 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) 1036 ENDDO 1037 1038 sums_ls_l(nzt+1,2) = sums_ls_l(nzt,2) 1039 sums_ls_l(nzt+1,3) = sums_ls_l(nzt,3) 1030 1040 1031 1041 ENDIF … … 1192 1202 1193 1203 IF ( large_scale_forcing ) THEN 1194 hom(:,1,81,sr) = sums_ls_l(:,0) ! pt_lsa1195 hom(:,1,82,sr) = sums_ls_l(:,1) ! q_lsa1204 hom(:,1,81,sr) = sums_ls_l(:,0) ! td_lsa_lpt 1205 hom(:,1,82,sr) = sums_ls_l(:,1) ! td_lsa_q 1196 1206 IF ( use_subsidence_tendencies ) THEN 1197 hom(:,1,83,sr) = sums_ls_l(:,2) ! pt_subs1198 hom(:,1,84,sr) = sums_ls_l(:,3) ! q_subs1207 hom(:,1,83,sr) = sums_ls_l(:,2) ! td_sub_lpt 1208 hom(:,1,84,sr) = sums_ls_l(:,3) ! td_sub_q 1199 1209 ELSE 1200 hom(:,1,83,sr) = sums(:,83) ! pt_subs1201 hom(:,1,84,sr) = sums(:,84) ! q_subs1210 hom(:,1,83,sr) = sums(:,83) ! td_sub_lpt 1211 hom(:,1,84,sr) = sums(:,84) ! td_sub_q 1202 1212 ENDIF 1203 hom(:,1,85,sr) = sums(:,85) ! pt_nudge1204 hom(:,1,86,sr) = sums(:,86) ! q_nudge1205 hom(:,1,87,sr) = sums(:,87) ! u_nudge1206 hom(:,1,88,sr) = sums(:,88) ! v_nudge1213 hom(:,1,85,sr) = sums(:,85) ! td_nud_lpt 1214 hom(:,1,86,sr) = sums(:,86) ! td_nud_q 1215 hom(:,1,87,sr) = sums(:,87) ! td_nud_u 1216 hom(:,1,88,sr) = sums(:,88) ! td_nud_v 1207 1217 ENDIF 1208 1218 … … 1386 1396 1387 1397 USE arrays_3d, & 1388 ONLY: ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, pt_lsa, pt_subs, q,&1389 q c, ql, qr, qs, qsws, qswst, q_lsa, q_subs, rho, sa, saswsb,&1390 saswst, shf, time_vert, ts, tswst, u, ug, us, usws, uswst, vsws,&1391 v, vg, vpt, vswst, w, w_subs, zw1398 ONLY: ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, q, qc, ql, qr, qs, & 1399 qsws, qswst, rho, sa, saswsb, saswst, shf, td_lsa_lpt, td_lsa_q,& 1400 td_sub_lpt, td_sub_q, time_vert, ts, tswst, u, ug, us, usws, & 1401 uswst, vsws, v, vg, vpt, vswst, w, w_subs, zw 1392 1402 1393 1403 … … 2763 2773 2764 2774 DO k = nzb, nzt 2765 sums_ls_l(k,0) = pt_lsa(k,nt) & 2766 + fac * ( pt_lsa(k,nt+1) - pt_lsa(k,nt) ) 2767 sums_ls_l(k,1) = q_lsa(k,nt) & 2768 + fac * ( q_lsa(k,nt+1) - q_lsa(k,nt) ) 2769 ENDDO 2775 sums_ls_l(k,0) = td_lsa_lpt(k,nt) & 2776 + fac * ( td_lsa_lpt(k,nt+1) - td_lsa_lpt(k,nt) ) 2777 sums_ls_l(k,1) = td_lsa_q(k,nt) & 2778 + fac * ( td_lsa_q(k,nt+1) - td_lsa_q(k,nt) ) 2779 ENDDO 2780 2781 sums_ls_l(nzt+1,0) = sums_ls_l(nzt,0) 2782 sums_ls_l(nzt+1,1) = sums_ls_l(nzt,1) 2770 2783 2771 2784 IF ( large_scale_subsidence .AND. use_subsidence_tendencies ) THEN 2772 2785 2773 2786 DO k = nzb, nzt 2774 sums_ls_l(k,2) = pt_subs(k,nt) & 2775 + fac * ( pt_subs(k,nt+1) - pt_subs(k,nt) ) 2776 sums_ls_l(k,3) = q_subs(k,nt) & 2777 + fac * ( q_subs(k,nt+1) - q_subs(k,nt) ) 2778 ENDDO 2787 sums_ls_l(k,2) = td_sub_lpt(k,nt) + fac * & 2788 ( td_sub_lpt(k,nt+1) - td_sub_lpt(k,nt) ) 2789 sums_ls_l(k,3) = td_sub_q(k,nt) + fac * & 2790 ( td_sub_q(k,nt+1) - td_sub_q(k,nt) ) 2791 ENDDO 2792 2793 sums_ls_l(nzt+1,2) = sums_ls_l(nzt,2) 2794 sums_ls_l(nzt+1,3) = sums_ls_l(nzt,3) 2779 2795 2780 2796 ENDIF … … 2943 2959 2944 2960 IF ( large_scale_forcing ) THEN 2945 hom(:,1,81,sr) = sums_ls_l(:,0) ! pt_lsa2946 hom(:,1,82,sr) = sums_ls_l(:,1) ! q_lsa2961 hom(:,1,81,sr) = sums_ls_l(:,0) ! td_lsa_lpt 2962 hom(:,1,82,sr) = sums_ls_l(:,1) ! td_lsa_q 2947 2963 IF ( use_subsidence_tendencies ) THEN 2948 hom(:,1,83,sr) = sums_ls_l(:,2) ! pt_subs2949 hom(:,1,84,sr) = sums_ls_l(:,3) ! q_subs2964 hom(:,1,83,sr) = sums_ls_l(:,2) ! td_sub_lpt 2965 hom(:,1,84,sr) = sums_ls_l(:,3) ! td_sub_q 2950 2966 ELSE 2951 hom(:,1,83,sr) = sums(:,83) ! pt_subs2952 hom(:,1,84,sr) = sums(:,84) ! q_subs2967 hom(:,1,83,sr) = sums(:,83) ! td_sub_lpt 2968 hom(:,1,84,sr) = sums(:,84) ! td_sub_q 2953 2969 ENDIF 2954 hom(:,1,85,sr) = sums(:,85) ! pt_nudge2955 hom(:,1,86,sr) = sums(:,86) ! q_nudge2956 hom(:,1,87,sr) = sums(:,87) ! u_nudge2957 hom(:,1,88,sr) = sums(:,88) ! v_nudge2970 hom(:,1,85,sr) = sums(:,85) ! td_nud_lpt 2971 hom(:,1,86,sr) = sums(:,86) ! td_nud_q 2972 hom(:,1,87,sr) = sums(:,87) ! td_nud_u 2973 hom(:,1,88,sr) = sums(:,88) ! td_nud_v 2958 2974 END IF 2959 2975
Note: See TracChangeset
for help on using the changeset viewer.