Changeset 624 for palm

Timestamp:

Dec 10, 2010 11:46:30 AM (14 years ago)

Author:

heinze

Message:

Calculation of q*2 added

Location:

palm/trunk

Files:

• SCRIPTS/NCL/profiles.ncl (modified) (1 diff)
• SOURCE/flow_statistics.f90 (modified) (2 diffs)

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -697 +697 @@
                                                                         .OR. vNam(varn) .EQ. "w"+dq+"e:dz" .OR. isStrSubset(vNam(varn), "w"+dq+"e:dz_"))then
                                                                         z_u = f_att->$vNam(varn+1)$
-                                                                        break                                                                        
+                                                                        break   
+                                                                     else
+                                                                        if(vNam(varn) .EQ. "qs2" .OR. isStrSubset(vNam(varn), "qs2_")\
+                                                                           .OR. vNam(varn) .EQ. "q*2" .OR. isStrSubset(vNam(varn), "q*2_"))then
+                                                                           z_u = f_att->$vNam(varn+1)$
+                                                                           break  
+                                                                       end if               
                                                                      end if
                                                                   end if

palm/trunk/SOURCE/flow_statistics.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Calculation of q*2 added
 !
 ! Former revisions:
@@ -373 +373 @@
                 sums_l(k,33,tn) = sums_l(k,33,tn) + &
                                   ( pt(k,j,i)-hom(k,1,4,sr) )**2 * rmask(j,i,sr)
+
+                IF ( humidity )  THEN
+                   sums_l(k,70,tn) = sums_l(k,70,tn) + &
+                                  ( q(k,j,i)-hom(k,1,41,sr) )**2 * rmask(j,i,sr)
+                ENDIF
 !
 !--             Higher moments