Ignore:
Timestamp:
Apr 30, 2014 12:15:41 PM (7 years ago)
Author:
boeske
Message:

minor changes in profile data output of lsf tendencies, variables renamed

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/flow_statistics.f90

    r1375 r1382  
    2121! Current revisions:
    2222! -----------------
    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
    2428!
    2529! Former revisions:
     
    103107
    104108    USE arrays_3d,                                                             &
    105         ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, pt_lsa, pt_subs, q,&
    106                qc, 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, zw
     109        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
    109113       
    110114    USE cloud_parameters,                                                      &
     
    10141018
    10151019          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)
    10211028
    10221029          IF ( large_scale_subsidence .AND. use_subsidence_tendencies )  THEN
    10231030
    10241031             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)
    10301040
    10311041          ENDIF
     
    11921202
    11931203       IF ( large_scale_forcing )  THEN
    1194           hom(:,1,81,sr) = sums_ls_l(:,0)          ! pt_lsa
    1195           hom(:,1,82,sr) = sums_ls_l(:,1)          ! q_lsa
     1204          hom(:,1,81,sr) = sums_ls_l(:,0)          ! td_lsa_lpt
     1205          hom(:,1,82,sr) = sums_ls_l(:,1)          ! td_lsa_q
    11961206          IF ( use_subsidence_tendencies )  THEN
    1197              hom(:,1,83,sr) = sums_ls_l(:,2)       ! pt_subs
    1198              hom(:,1,84,sr) = sums_ls_l(:,3)       ! q_subs
     1207             hom(:,1,83,sr) = sums_ls_l(:,2)       ! td_sub_lpt
     1208             hom(:,1,84,sr) = sums_ls_l(:,3)       ! td_sub_q
    11991209          ELSE
    1200              hom(:,1,83,sr) = sums(:,83)           ! pt_subs
    1201              hom(:,1,84,sr) = sums(:,84)           ! q_subs
     1210             hom(:,1,83,sr) = sums(:,83)           ! td_sub_lpt
     1211             hom(:,1,84,sr) = sums(:,84)           ! td_sub_q
    12021212          ENDIF
    1203           hom(:,1,85,sr) = sums(:,85)              ! pt_nudge
    1204           hom(:,1,86,sr) = sums(:,86)              ! q_nudge
    1205           hom(:,1,87,sr) = sums(:,87)              ! u_nudge
    1206           hom(:,1,88,sr) = sums(:,88)              ! v_nudge
     1213          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
    12071217       ENDIF
    12081218
     
    13861396
    13871397    USE arrays_3d,                                                             &
    1388         ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, pt_lsa, pt_subs, q,&
    1389                qc, 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, zw
     1398        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
    13921402                 
    13931403       
     
    27632773
    27642774          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)
    27702783
    27712784          IF ( large_scale_subsidence .AND. use_subsidence_tendencies )  THEN
    27722785
    27732786             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)
    27792795
    27802796          ENDIF
     
    29432959
    29442960       IF ( large_scale_forcing )  THEN
    2945           hom(:,1,81,sr) = sums_ls_l(:,0)          ! pt_lsa
    2946           hom(:,1,82,sr) = sums_ls_l(:,1)          ! q_lsa
     2961          hom(:,1,81,sr) = sums_ls_l(:,0)          ! td_lsa_lpt
     2962          hom(:,1,82,sr) = sums_ls_l(:,1)          ! td_lsa_q
    29472963          IF ( use_subsidence_tendencies )  THEN
    2948              hom(:,1,83,sr) = sums_ls_l(:,2)       ! pt_subs
    2949              hom(:,1,84,sr) = sums_ls_l(:,3)       ! q_subs
     2964             hom(:,1,83,sr) = sums_ls_l(:,2)       ! td_sub_lpt
     2965             hom(:,1,84,sr) = sums_ls_l(:,3)       ! td_sub_q
    29502966          ELSE
    2951              hom(:,1,83,sr) = sums(:,83)           ! pt_subs
    2952              hom(:,1,84,sr) = sums(:,84)           ! q_subs
     2967             hom(:,1,83,sr) = sums(:,83)           ! td_sub_lpt
     2968             hom(:,1,84,sr) = sums(:,84)           ! td_sub_q
    29532969          ENDIF
    2954           hom(:,1,85,sr) = sums(:,85)              ! pt_nudge
    2955           hom(:,1,86,sr) = sums(:,86)              ! q_nudge
    2956           hom(:,1,87,sr) = sums(:,87)              ! u_nudge
    2957           hom(:,1,88,sr) = sums(:,88)              ! v_nudge
     2970          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
    29582974       END IF
    29592975
Note: See TracChangeset for help on using the changeset viewer.