Changeset 514 for palm/trunk

Timestamp:

Mar 17, 2010 9:42:33 AM (15 years ago)

Author:

heinze

Message:

Bugfix concerning the plot of xy cross sections which have only one layer (zu1_xy)

File:

• palm/trunk/SCRIPTS/NCL/cross_sections.ncl (modified) (3 diffs)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -1318 +1318 @@
             temp = f[:]->$vNam(varn)$
             data_att = f_att->$vNam(varn)$
-            if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")
-              ;these variables depend zu1_xy and that's why they have only one z-layer
+            if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"     \
+              .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
+              .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
+              .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
+              ;these variables depend on zu1_xy and that's why they have only one z-layer
               data(varn,:,0,:,:)=temp(:,0,ys:ye,xs:xe) 
               no_zu1=no_zu1+1
@@ -1600 +1603 @@
                      end if
 
-                     if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")then
+                     if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"       \
+                         .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
+                         .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
+                         .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
                          loe = 0
                          level = "=" + zu1(0) + "m"
                      end if
                     
-                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level                                
+                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level           
                      plot(n) = gsn_csm_contour(wks_ps,data(varn,li,lo-los,:,:),cs_res) 
                      if (vector .EQ. 1 .AND. check_vecp) then
@@ -1633 +1639 @@
                      end if 
        
-                     if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")then
+                     if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"       \
+                         .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
+                         .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
+                         .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
                          lie = 0
                          level = "=" + zu1(0) + "m"