Changeset 624 for palm


Ignore:
Timestamp:
Dec 10, 2010 11:46:30 AM (14 years ago)
Author:
heinze
Message:

Calculation of q*2 added

Location:
palm/trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r594 r624  
    697697                                                                        .OR. vNam(varn) .EQ. "w"+dq+"e:dz" .OR. isStrSubset(vNam(varn), "w"+dq+"e:dz_"))then
    698698                                                                        z_u = f_att->$vNam(varn+1)$
    699                                                                         break                                                                       
     699                                                                        break   
     700                                                                     else
     701                                                                        if(vNam(varn) .EQ. "qs2" .OR. isStrSubset(vNam(varn), "qs2_")\
     702                                                                           .OR. vNam(varn) .EQ. "q*2" .OR. isStrSubset(vNam(varn), "q*2_"))then
     703                                                                           z_u = f_att->$vNam(varn+1)$
     704                                                                           break 
     705                                                                       end if               
    700706                                                                     end if
    701707                                                                  end if
  • palm/trunk/SOURCE/flow_statistics.f90

    r623 r624  
    44! Current revisions:
    55! -----------------
    6 !
     6! Calculation of q*2 added
    77!
    88! Former revisions:
     
    373373                sums_l(k,33,tn) = sums_l(k,33,tn) + &
    374374                                  ( pt(k,j,i)-hom(k,1,4,sr) )**2 * rmask(j,i,sr)
     375
     376                IF ( humidity )  THEN
     377                   sums_l(k,70,tn) = sums_l(k,70,tn) + &
     378                                  ( q(k,j,i)-hom(k,1,41,sr) )**2 * rmask(j,i,sr)
     379                ENDIF
    375380!
    376381!--             Higher moments
Note: See TracChangeset for help on using the changeset viewer.