Changeset 769 for palm/trunk/SCRIPTS/NCL/cross_sections.ncl

Timestamp:

Oct 26, 2011 2:20:45 PM (11 years ago)

Author:

heinze

Message:

Bugfixes in case of plot of t=0h and plot of topography zusi/zwwi possible

File:

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

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -454 +454 @@
          st=i
          break
+      else
+         st=0
       end if
    end do
-   
-   if (.not. isvar("st"))then
-      print(" ")
-      print("'start_time_step' = "+ start_time_step +"h is invalid")
-      print(" ")
-      print("Select another 'start_time_step'")
-      print(" ")
-      exit
-   end if
+
         
    ; ****************************************************
@@ -506 +500 @@
          et=i
          break
+      else
+         et=0
       end if
    end do 
-   
-   if (.not. isvar("et"))then
-      print(" ")
-      print("'end_time_step' = "+ end_time_step +"h is invalid")
-      print(" ")
-      print("Select another 'end_time_step'")
-      print(" ")
-      exit
-   end if
   
    delete(start_time_step)
@@ -1421 +1408 @@
    n      = 0
    no_zu1 = 0
+   no_topo= 0
 
    ;****************************************************
@@ -1461 +1449 @@
       if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or.          \
            vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or.          \
-           vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or.      \
            vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or.     \
            vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x".or.         \
@@ -1485 +1472 @@
          if (xyc .EQ. 1) then
             temp = f[:]->$vNam(varn)$
-            data_att = f_att->$vNam(varn)$
+            if (vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi")
+               dummy=0
+            else
+               data_att = f_att->$vNam(varn)$
+            end if
             if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"     \
               .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
@@ -1499 +1490 @@
               no_zu1=no_zu1+1
             else
-              data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe) 
+              if (vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi")
+                  ;these variables depend on x and y
+                  data(varn,0,0,:,:)=doubletofloat(temp(ys:ye,xs:xe))
+                  no_topo=no_topo+1
+               else
+                  data(varn,:,:,:,:)=temp(:,zs:ze,ys:ye,xs:xe) 
+               end if
             end if
             delete(temp)                
@@ -1527 +1524 @@
                                                0:(ze-zs),0:(ye-ys),0:(xe-xs)))
          
-         unit(varn) = data_att@units
-         delete(data_att)
+         if (vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi")
+            unit(varn) = "meters"
+          else
+            unit(varn) = data_att@units
+            delete(data_att)
+          end if
 
       end if
@@ -1550 +1551 @@
       if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or.          \
            vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or.          \
-           vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or.      \
            vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or.     \
            vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x".or.         \
@@ -1627 +1627 @@
       plot=new((/no_time*no_layer/),graphic)
    else
-      plot=new((/no_time*no_layer*no_var - no_time*(no_layer-1)*no_zu1/)\
-                                                                ,graphic)
+      plot=new((/no_time*no_layer*(no_var-no_topo) + no_topo - \
+                       no_time*(no_layer-1)*no_zu1/),graphic)
    end if
    dim_plot=dimsizes(plot)
@@ -1756 +1756 @@
       if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or.          \
            vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or.          \
-           vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or.      \
            vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or.     \
            vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x".or.         \
@@ -1829 +1828 @@
                          loe = 0
                          level = "=" + zu1(0) + "m"
+                      else
+                         loe = laye
+                     end if
+
+                     if(vNam(varn) .eq. "zwwi"  .or. \
+                        vNam(varn) .eq. "zusi") then
+                         loe = 0
+                         los = 0
+                         lie  = 0
+                         lis = 0
+                         level = ""
+                      else
+                         lis = start_time_step
+                         lie = end_time_step
+                         los = lays
+                         loe = laye
                      end if
                     
-                     cs_res@gsnCenterString = "t=" + \
-                       decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level
+                     if(vNam(varn) .eq. "zwwi"  .or. \
+                        vNam(varn) .eq. "zusi") then
+                         cs_res@gsnCenterString = ""
+                     else
+                         cs_res@gsnCenterString = "t=" + \
+                          decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level
+                     end if
 
                      if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
@@ -1894 +1914 @@
                          lie = 0
                          level = "=" + zu1(0) + "m"
+                     else
+                         lie = laye
                      end if
-                
-                     cs_res@gsnCenterString = "t=" + \
-                           decimalPlaces(t_all(lo)/3600,2,True) + "h  z"+ level
+
+                     if(vNam(varn) .eq. "zwwi"  .or. \
+                        vNam(varn) .eq. "zusi") then
+                         lie = 0
+                         lis = 0
+                         loe = 0
+                         los = 0
+                         level = ""
+                     else
+                         lis = lays
+                         lie = laye
+                         los = start_time_step 
+                         loe = end_time_step
+                     end if
+
+                     if(vNam(varn) .eq. "zwwi"  .or. \
+                        vNam(varn) .eq. "zusi") then
+                         cs_res@gsnCenterString = ""
+                     else
+                         cs_res@gsnCenterString = "t=" + \
+                          decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level
+                     end if
 
                      if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then