Ignore:
Timestamp:
Sep 23, 2009 9:40:33 AM (12 years ago)
Author:
raasch
Message:

in-situ AND potential density are calculated and used in the ocean version

File:
1 edited

Legend:

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

    r343 r388  
    44! Current revisions:
    55! -----------------
     6! Vertical profiles of potential density and hydrostatic pressure are
     7! calculated.
    68! Added missing timeseries calculation of w"q"(0), moved timeseries q* to the
    79! end.
     
    531533                IF ( humidity )  THEN
    532534                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
    533                                        qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     535                                       qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv")
    534536                   IF ( cloud_physics )  THEN
    535537                      sums_l(nzt,45,tn) = sums_l(nzt,45,tn) + (           &
     
    546548                IF ( passive_scalar )  THEN
    547549                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
    548                                        qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     550                                       qswst(j,i) * rmask(j,i,sr) ! w"q" (w"qv")
    549551                ENDIF
    550552             ENDIF
     
    597599                                                       rmask(j,i,sr)
    598600                   sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) * &
     601                                                       rmask(j,i,sr)
     602                   sums_l(k,71,tn) = sums_l(k,71,tn) + prho(k,j,i) * &
    599603                                                       rmask(j,i,sr)
    600604                ENDIF
     
    644648!
    645649!--    Density at top follows Neumann condition
    646        IF ( ocean )  sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn)
     650       IF ( ocean )  THEN
     651          sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn)
     652          sums_l(nzt+1,71,tn) = sums_l(nzt,71,tn)
     653       ENDIF
    647654
    648655!
     
    897904       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
    898905       hom(:,1,70,sr) = sums(:,70)     ! q*2
     906       hom(:,1,71,sr) = sums(:,71)     ! prho
     907       hom(:,1,72,sr) = hyp * 1E-4     ! hyp in kPa
    899908
    900909       hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)
     
    916925!--    is less than 1.5 times the height where the heat flux becomes negative
    917926!--    (positive) for the first time.
    918 !--    NOTE: This criterion is still capable of improving!
    919927       z_i(1) = 0.0
    920928       first = .TRUE.
Note: See TracChangeset for help on using the changeset viewer.