Ignore:
Timestamp:
Oct 21, 2016 3:11:58 PM (5 years ago)
Author:
knoop
Message:

Renamed variable rho to rho_ocean, rho_init to rho_ocean_init and rho_av to rho_ocean_av

File:
1 edited

Legend:

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

    r2027 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    212212    USE arrays_3d,                                                             &
    213213        ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, ol, p, prho, prr, pt, q, qc, ql,&
    214                qr, qs, qsws, qswst, rho, s, sa, ss, ssws, sswst, saswsb,       &
     214               qr, qs, qsws, qswst, rho_ocean, s, sa, ss, ssws, sswst, saswsb,       &
    215215               saswst, shf, td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q,        &
    216216               time_vert, ts, tswst, u, ug, us, usws, uswst, vsws, v, vg, vpt, &
     
    926926                                        rmask(j,i,sr)
    927927                   ENDIF
    928                    sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) *            &
     928                   sums_l(k,64,tn) = sums_l(k,64,tn) + rho_ocean(k,j,i) *            &
    929929                                                       rmask(j,i,sr)
    930930                   sums_l(k,71,tn) = sums_l(k,71,tn) + prho(k,j,i) *           &
     
    14001400       hom(:,1,55,sr) = sums(:,55)     ! w*u*u*/dz
    14011401       hom(:,1,56,sr) = sums(:,56)     ! w*p*/dz
    1402        hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho )/dz
     1402       hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho_ocean )/dz
    14031403       hom(:,1,58,sr) = sums(:,58)     ! u"pt"
    14041404       hom(:,1,59,sr) = sums(:,59)     ! u*pt*
     
    14071407       hom(:,1,62,sr) = sums(:,62)     ! v*pt*
    14081408       hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
    1409        hom(:,1,64,sr) = sums(:,64)     ! rho
     1409       hom(:,1,64,sr) = sums(:,64)     ! rho_ocean
    14101410       hom(:,1,65,sr) = sums(:,65)     ! w"sa"
    14111411       hom(:,1,66,sr) = sums(:,66)     ! w*sa*
    14121412       hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
    14131413       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    1414        hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     1414       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho_ocean
    14151415       hom(:,1,70,sr) = sums(:,70)     ! q*2
    14161416       hom(:,1,71,sr) = sums(:,71)     ! prho
     
    17061706    USE arrays_3d,                                                             &
    17071707        ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, q, qc, ql, qr, qs, &
    1708                qsws, qswst, rho, s, sa, saswsb, saswst, shf, ss, ssws, sswst,  &
     1708               qsws, qswst, rho_ocean, s, sa, saswsb, saswst, shf, ss, ssws, sswst,  &
    17091709               td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q, time_vert, ts,      &
    17101710               tswst, u, ug, us, usws, uswst, vsws, v, vg, vpt, vswst, w,      &
     
    27862786          ENDIF
    27872787
    2788           !$acc parallel loop gang present( rflags_invers, rho, prho, rmask, sums_l ) create( s1, s2 )
     2788          !$acc parallel loop gang present( rflags_invers, rho_ocean, prho, rmask, sums_l ) create( s1, s2 )
    27892789          DO  k = nzb, nzt_diff
    27902790             s1 = 0
     
    27932793             DO  i = nxl, nxr
    27942794                DO  j = nys, nyn
    2795                    s1 = s1 + rho(k,j,i)  * rmask(j,i,sr) * rflags_invers(j,i,k+1)
     2795                   s1 = s1 + rho_ocean(k,j,i)  * rmask(j,i,sr) * rflags_invers(j,i,k+1)
    27962796                   s2 = s2 + prho(k,j,i) * rmask(j,i,sr) * rflags_invers(j,i,k+1)
    27972797                ENDDO
     
    34243424       hom(:,1,55,sr) = sums(:,55)     ! w*u*u*/dz
    34253425       hom(:,1,56,sr) = sums(:,56)     ! w*p*/dz
    3426        hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho )/dz
     3426       hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho_ocean )/dz
    34273427       hom(:,1,58,sr) = sums(:,58)     ! u"pt"
    34283428       hom(:,1,59,sr) = sums(:,59)     ! u*pt*
     
    34313431       hom(:,1,62,sr) = sums(:,62)     ! v*pt*
    34323432       hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
    3433        hom(:,1,64,sr) = sums(:,64)     ! rho
     3433       hom(:,1,64,sr) = sums(:,64)     ! rho_ocean
    34343434       hom(:,1,65,sr) = sums(:,65)     ! w"sa"
    34353435       hom(:,1,66,sr) = sums(:,66)     ! w*sa*
    34363436       hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
    34373437       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    3438        hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     3438       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho_ocean
    34393439       hom(:,1,70,sr) = sums(:,70)     ! q*2
    34403440       hom(:,1,71,sr) = sums(:,71)     ! prho
Note: See TracChangeset for help on using the changeset viewer.