Changeset 353 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Jul 9, 2009 1:59:55 PM (15 years ago)
Author:
heinze
Message:

small adjustments of plotting area

File:
1 edited

Legend:

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

    r351 r353  
    15241524   res@xyDashPattern = 1*nof
    15251525
    1526    
     1526   over_remind = False
     1527   if ( over .eq. 1)then
     1528     over_remind = True
     1529   end if
    15271530   
    15281531   do varn = 0,dim-1
     
    16911694            end do           
    16921695         end if
     1696
     1697         if (isStrSubset(data@long_name," SR " ) ) then
     1698           over = 0
     1699         end if
    16931700         
    16941701         if (nof .EQ. 0) then
     
    16981705            z=z/norm_z
    16991706         end if
     1707
     1708         delta_z = z(1) - z(0)
     1709
     1710         max_z_int=doubletoint(max_z/delta_z)
     1711         min_z_int=doubletoint(min_z/delta_z)
     1712
     1713         if(max_z_int .eq. min_z_int)
     1714             print(" ")
     1715             print("Please increase 'max_z' or decrease 'min_z' so that there are")
     1716             print("at least two layers for the z-axis to plot")
     1717             print(" ")
     1718             exit
     1719         end if
     1720
     1721         if(min_z_int .gt. max_z_int)
     1722            print(" ")
     1723            print("'min_z' is greater than 'max_z',")
     1724            print("please change this")
     1725            print(" ")
     1726            exit
     1727         end if
     1728
     1729         if(max_z_int .gt. dimz-1)
     1730            max_z_int = dimz-1
     1731         end if
     1732
     1733         if(min_z_int .lt. 0)
     1734            min_z_int = 0
     1735         end if         
     1736         
    17001737         
    1701          if (abs(min(data(varn,:,:))) .GT. 10)then
    1702             min_value = abs(0.01*min(data(varn,:,:)))
    1703             max_value = abs(0.01*max(data(varn,:,:)))
     1738         if (abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 10)then
     1739            min_value = abs(0.001*min(data(varn,:,min_z_int:max_z_int)))
     1740            max_value = abs(0.001*max(data(varn,:,min_z_int:max_z_int)))
    17041741         else
    1705             if (abs(min(data(varn,:,:))) .LT. 0.01 .AND. abs(max(data(varn,:,:))) .GT. 0.01)then
    1706                min_value = abs(0.1*max(data(varn,:,:)))
    1707                max_value = abs(0.1*max(data(varn,:,:)))
     1742            if (abs(min(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND. \
     1743                abs(max(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then
     1744               min_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
     1745               max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
    17081746            else
    1709                if (abs(max(data(varn,:,:))) .LT. 0.01 .AND. abs(min(data(varn,:,:))) .GT. 0.01)then
    1710                   min_value = abs(0.1*min(data(varn,:,:)))
    1711                   max_value = abs(0.1*min(data(varn,:,:)))
     1747               if (abs(max(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND.\
     1748                   abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then
     1749                  min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
     1750                  max_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
    17121751               else
    1713                   min_value = abs(0.1*min(data(varn,:,:)))
    1714                   max_value = abs(0.1*max(data(varn,:,:)))
     1752                  min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
     1753                  max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
    17151754               end if
    17161755            end if
    17171756         end if
    1718          if (min(data(varn,:,:)) .EQ. 0 .AND. max(data(varn,:,:)) .EQ. 0)then
     1757         if (min(data(varn,:,min_z_int:max_z_int)) .EQ. 0 .AND. \
     1758             max(data(varn,:,min_z_int:max_z_int)) .EQ. 0)then
    17191759            min_value = 0.1
    17201760            max_value = 0.1
     
    17281768            res@trYMaxF            = max_z
    17291769            if (xs .EQ. -1) then
    1730                res@trXMinF            = min(data(varn,:,:))-min_value
     1770               res@trXMinF            = min(data(varn,:,min_z_int:max_z_int))-min_value
    17311771            else
    17321772               res@trXMinF            = xs
    17331773            end if
    17341774            if (xe .EQ. -1) then
    1735                res@trXMaxF            = max(data(varn,:,:))+max_value
     1775               res@trXMaxF            = max(data(varn,:,min_z_int:max_z_int))+max_value
    17361776            else
    17371777               res@trXMaxF            = xe 
    17381778            end if         
    17391779            plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 
    1740          end if 
    1741 
    1742             if (vNam(varn) .EQ. "u") then
    1743                miniu=min(data(varn,:,:))-min_value
    1744                maxiu=max(data(varn,:,:))+max_value
     1780         end if
     1781 
     1782     
     1783         if (vNam(varn) .EQ. "u" ) then
     1784               miniu=min(data(varn,:,min_z_int:max_z_int))-min_value
     1785               maxiu=max(data(varn,:,min_z_int:max_z_int))+max_value
    17451786               if (over .EQ. 1) then
    17461787                  res@xyDashPattern  = 0
     
    17631804               end if
    17641805            end if
    1765             if (vNam(varn) .EQ. "v") then
    1766                miniv=min(data(varn,:,:))-min_value
    1767                maxiv=max(data(varn,:,:))+max_value
     1806            if (vNam(varn) .EQ. "v") then 
     1807               miniv=min(data(varn,:,min_z_int:max_z_int))-min_value
     1808               maxiv=max(data(varn,:,min_z_int:max_z_int))+max_value
    17681809               if (over .EQ. 1) then
    17691810                  res@xyMonoDashPattern = True
     
    17881829            end if
    17891830            if (vNam(varn) .EQ. "w") then
    1790                miniw=min(data(varn,:,:))-min_value
    1791                maxiw=max(data(varn,:,:))+max_value
     1831               miniw=min(data(varn,:,min_z_int:max_z_int))-min_value
     1832               maxiw=max(data(varn,:,min_z_int:max_z_int))+max_value
    17921833               if (over .EQ. 1) then
    17931834                  res@xyDashPattern = 2
     
    18121853
    18131854            if (vNam(varn) .EQ. "pt") then
    1814                minipt=min(data(varn,:,:))-min_value
    1815                maxipt=max(data(varn,:,:))+max_value
     1855               minipt=min(data(varn,:,min_z_int:max_z_int))-min_value
     1856               maxipt=max(data(varn,:,min_z_int:max_z_int))+max_value
    18161857               if (over .EQ. 1) then
    18171858                  res@xyDashPattern  = 0
     
    18351876            end if
    18361877            if (vNam(varn) .EQ. "vpt") then
    1837                minivpt=min(data(varn,:,:))-min_value
    1838                maxivpt=max(data(varn,:,:))+max_value
     1878               minivpt=min(data(varn,:,min_z_int:max_z_int))-min_value
     1879               maxivpt=max(data(varn,:,min_z_int:max_z_int))+max_value
    18391880               if (over .EQ. 1) then
    18401881                  res@xyDashPattern  = 1
     
    18581899            end if
    18591900            if (vNam(varn) .EQ. "lpt") then
    1860                minilpt=min(data(varn,:,:))-min_value
    1861                maxilpt=max(data(varn,:,:))+max_value
     1901               minilpt=min(data(varn,:,min_z_int:max_z_int))-min_value
     1902               maxilpt=max(data(varn,:,min_z_int:max_z_int))+max_value
    18621903               if (over .EQ. 1) then
    18631904                  res@xyDashPattern  = 2
     
    18821923
    18831924            if (vNam(varn) .EQ. "q") then
    1884                miniq=min(data(varn,:,:))-min_value
    1885                maxiq=max(data(varn,:,:))+max_value
     1925               miniq=min(data(varn,:,min_z_int:max_z_int))-min_value
     1926               maxiq=max(data(varn,:,min_z_int:max_z_int))+max_value
    18861927               if (over .EQ. 1) then
    18871928                  res@xyDashPattern  = 0
     
    19051946            end if
    19061947            if (vNam(varn) .EQ. "qv") then
    1907                miniqv=min(data(varn,:,:))-min_value
    1908                maxiqv=max(data(varn,:,:))+max_value
     1948               miniqv=min(data(varn,:,min_z_int:max_z_int))-min_value
     1949               maxiqv=max(data(varn,:,min_z_int:max_z_int))+max_value
    19091950               if (over .EQ. 1) then
    19101951                  res@xyDashPattern  = 1
     
    19281969            end if
    19291970            if (vNam(varn) .EQ. "ql") then
    1930                miniql=min(data(varn,:,:))-min_value
    1931                maxiql=max(data(varn,:,:))+max_value
     1971               miniql=min(data(varn,:,min_z_int:max_z_int))-min_value
     1972               maxiql=max(data(varn,:,min_z_int:max_z_int))+max_value
    19321973               if (over .EQ. 1) then
    19331974                  res@xyDashPattern  = 2
     
    19521993
    19531994            if (vNam(varn) .EQ. "e") then
    1954                minie=min(data(varn,:,:))-min_value
    1955                maxie=max(data(varn,:,:))+max_value
     1995               minie=min(data(varn,:,min_z_int:max_z_int))-min_value
     1996               maxie=max(data(varn,:,min_z_int:max_z_int))+max_value
    19561997               if (over .EQ. 1) then
    19571998                  res@xyDashPattern  = 0
     
    19752016            end if
    19762017            if (vNam(varn) .EQ. "es") then
    1977                minies=min(data(varn,:,:))-min_value
    1978                maxies=max(data(varn,:,:))+max_value
     2018               minies=min(data(varn,:,min_z_int:max_z_int))-min_value
     2019               maxies=max(data(varn,:,min_z_int:max_z_int))+max_value
    19792020               if (over .EQ. 1) then
    19802021                  res@xyDashPattern  = 1
     
    19992040
    20002041            if (vNam(varn) .EQ. "km") then
    2001                minikm=min(data(varn,:,:))-min_value
    2002                maxikm=max(data(varn,:,:))+max_value
     2042               minikm=min(data(varn,:,min_z_int:max_z_int))-min_value
     2043               maxikm=max(data(varn,:,min_z_int:max_z_int))+max_value
    20032044               if (over .EQ. 1) then
    20042045                  res@xyDashPattern  = 0
     
    20222063            end if
    20232064            if (vNam(varn) .EQ. "kh") then
    2024                minikh=min(data(varn,:,:))-min_value
    2025                maxikh=max(data(varn,:,:))+max_value
     2065               minikh=min(data(varn,:,min_z_int:max_z_int))-min_value
     2066               maxikh=max(data(varn,:,min_z_int:max_z_int))+max_value
    20262067               if (over .EQ. 1) then
    20272068                  res@xyDashPattern  = 1
     
    20462087
    20472088            if (vNam(varn) .EQ. "wpup") then
    2048                miniwpup=min(data(varn,:,:))-min_value
    2049                maxiwpup=max(data(varn,:,:))+max_value
     2089               miniwpup=min(data(varn,:,min_z_int:max_z_int))-min_value
     2090               maxiwpup=max(data(varn,:,min_z_int:max_z_int))+max_value
    20502091               if (over .EQ. 1) then
    20512092                  res@xyDashPattern  = 0
     
    20692110            end if
    20702111            if (vNam(varn) .EQ. "wsus") then
    2071                miniwsus=min(data(varn,:,:))-min_value
    2072                maxiwsus=max(data(varn,:,:))+max_value
     2112               miniwsus=min(data(varn,:,min_z_int:max_z_int))-min_value
     2113               maxiwsus=max(data(varn,:,min_z_int:max_z_int))+max_value
    20732114               if (over .EQ. 1) then
    20742115                  res@xyDashPattern  = 1
     
    20922133            end if
    20932134            if (vNam(varn) .EQ. "wu") then
    2094                miniwu=min(data(varn,:,:))-min_value
    2095                maxiwu=max(data(varn,:,:))+max_value
     2135               miniwu=min(data(varn,:,min_z_int:max_z_int))-min_value
     2136               maxiwu=max(data(varn,:,min_z_int:max_z_int))+max_value
    20962137               if (over .EQ. 1) then
    20972138                  res@xyDashPattern  = 2
     
    21162157
    21172158            if (vNam(varn) .EQ. "wpvp") then
    2118                miniwpvp=min(data(varn,:,:))-min_value
    2119                maxiwpvp=max(data(varn,:,:))+max_value
     2159               miniwpvp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2160               maxiwpvp=max(data(varn,:,min_z_int:max_z_int))+max_value
    21202161               if (over .EQ. 1) then
    21212162                  res@xyDashPattern  = 0
     
    21392180            end if
    21402181            if (vNam(varn) .EQ. "wsvs") then
    2141                miniwsvs=min(data(varn,:,:))-min_value
    2142                maxiwsvs=max(data(varn,:,:))+max_value
     2182               miniwsvs=min(data(varn,:,min_z_int:max_z_int))-min_value
     2183               maxiwsvs=max(data(varn,:,min_z_int:max_z_int))+max_value
    21432184               if (over .EQ. 1) then
    21442185                  res@xyDashPattern  = 1
     
    21622203            end if
    21632204            if (vNam(varn) .EQ. "wv") then
    2164                miniwv=min(data(varn,:,:))-min_value
    2165                maxiwv=max(data(varn,:,:))+max_value
     2205               miniwv=min(data(varn,:,min_z_int:max_z_int))-min_value
     2206               maxiwv=max(data(varn,:,min_z_int:max_z_int))+max_value
    21662207               if (over .EQ. 1) then
    21672208                  res@xyDashPattern  = 2
     
    21862227
    21872228            if (vNam(varn) .EQ. "wpptp") then
    2188                miniwpptp=min(data(varn,:,:))-min_value
    2189                maxiwpptp=max(data(varn,:,:))+max_value
     2229               miniwpptp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2230               maxiwpptp=max(data(varn,:,min_z_int:max_z_int))+max_value
    21902231               if (over .EQ. 1) then
    21912232                  res@xyDashPattern  = 0
     
    22092250            end if
    22102251            if (vNam(varn) .EQ. "wspts") then
    2211                miniwspts=min(data(varn,:,:))-min_value
    2212                maxiwspts=max(data(varn,:,:))+max_value
     2252               miniwspts=min(data(varn,:,min_z_int:max_z_int))-min_value
     2253               maxiwspts=max(data(varn,:,min_z_int:max_z_int))+max_value
    22132254               if (over .EQ. 1) then
    22142255                  res@xyDashPattern  = 1
     
    22322273            end if
    22332274            if (vNam(varn) .EQ. "wpt") then       
    2234                miniwpt=min(data(varn,:,:))-min_value
    2235                maxiwpt=max(data(varn,:,:))+max_value
     2275               miniwpt=min(data(varn,:,min_z_int:max_z_int))-min_value
     2276               maxiwpt=max(data(varn,:,min_z_int:max_z_int))+max_value
    22362277               if (over .EQ. 1) then
    22372278                  res@xyDashPattern  = 2
     
    22562297
    22572298            if (vNam(varn) .EQ. "wsptsBC") then
    2258                miniwsptsBC=min(data(varn,:,:))-min_value
    2259                maxiwsptsBC=max(data(varn,:,:))+max_value
     2299               miniwsptsBC=min(data(varn,:,min_z_int:max_z_int))-min_value
     2300               maxiwsptsBC=max(data(varn,:,min_z_int:max_z_int))+max_value
    22602301               if (over .EQ. 1) then
    22612302                  res@xyDashPattern  = 0
     
    22792320            end if             
    22802321            if (vNam(varn) .EQ. "wptBC") then
    2281                miniwptBC=min(data(varn,:,:))-min_value
    2282                maxiwptBC=max(data(varn,:,:))+max_value
     2322               miniwptBC=min(data(varn,:,min_z_int:max_z_int))-min_value
     2323               maxiwptBC=max(data(varn,:,min_z_int:max_z_int))+max_value
    22832324               if (over .EQ. 1) then
    22842325                  res@xyDashPattern  = 1
     
    23032344
    23042345            if (vNam(varn) .EQ. "wpvptp") then
    2305                miniwpvptp=min(data(varn,:,:))-min_value
    2306                maxiwpvptp=max(data(varn,:,:))+max_value
     2346               miniwpvptp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2347               maxiwpvptp=max(data(varn,:,min_z_int:max_z_int))+max_value
    23072348               if (over .EQ. 1) then
    23082349                  res@xyDashPattern  = 0
     
    23262367            end if
    23272368            if (vNam(varn) .EQ. "wsvpts") then
    2328                miniwsvpts=min(data(varn,:,:))-min_value
    2329                maxiwsvpts=max(data(varn,:,:))+max_value
     2369               miniwsvpts=min(data(varn,:,min_z_int:max_z_int))-min_value
     2370               maxiwsvpts=max(data(varn,:,min_z_int:max_z_int))+max_value
    23302371               if (over .EQ. 1) then
    23312372                  res@xyDashPattern  = 1
     
    23492390            end if
    23502391            if (vNam(varn) .EQ. "wvpt") then
    2351                miniwvpt=min(data(varn,:,:))-min_value
    2352                maxiwvpt=max(data(varn,:,:))+max_value
     2392               miniwvpt=min(data(varn,:,min_z_int:max_z_int))-min_value
     2393               maxiwvpt=max(data(varn,:,min_z_int:max_z_int))+max_value
    23532394               if (over .EQ. 1) then
    23542395                  res@xyDashPattern  = 2
     
    23732414
    23742415            if (vNam(varn) .EQ. "wpqp") then
    2375                miniwpqp=min(data(varn,:,:))-min_value
    2376                maxiwpqp=max(data(varn,:,:))+max_value
     2416               miniwpqp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2417               maxiwpqp=max(data(varn,:,min_z_int:max_z_int))+max_value
    23772418               if (over .EQ. 1) then
    23782419                  res@xyDashPattern  = 0
     
    23962437            end if
    23972438            if (vNam(varn) .EQ. "wsqs") then
    2398                miniwsqs=min(data(varn,:,:))-min_value
    2399                maxiwsqs=max(data(varn,:,:))+max_value
     2439               miniwsqs=min(data(varn,:,min_z_int:max_z_int))-min_value
     2440               maxiwsqs=max(data(varn,:,min_z_int:max_z_int))+max_value
    24002441               if (over .EQ. 1) then
    24012442                  res@xyDashPattern  = 1
     
    24192460            end if
    24202461            if (vNam(varn) .EQ. "wq") then
    2421                miniwq=min(data(varn,:,:))-min_value
    2422                maxiwq=max(data(varn,:,:))+max_value
     2462               miniwq=min(data(varn,:,min_z_int:max_z_int))-min_value
     2463               maxiwq=max(data(varn,:,min_z_int:max_z_int))+max_value
    24232464               if (over .EQ. 1) then
    24242465                  res@xyDashPattern  = 2
     
    24432484
    24442485            if (vNam(varn) .EQ. "wpqvp") then
    2445                miniwpqvp=min(data(varn,:,:))-min_value
    2446                maxiwpqvp=max(data(varn,:,:))+max_value
     2486               miniwpqvp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2487               maxiwpqvp=max(data(varn,:,min_z_int:max_z_int))+max_value
    24472488               if (over .EQ. 1) then
    24482489                  res@xyDashPattern  = 0
     
    24662507            end if
    24672508            if (vNam(varn) .EQ. "wsqvs") then
    2468                miniwsqvs=min(data(varn,:,:))-min_value
    2469                maxiwsqvs=max(data(varn,:,:))+max_value
     2509               miniwsqvs=min(data(varn,:,min_z_int:max_z_int))-min_value
     2510               maxiwsqvs=max(data(varn,:,min_z_int:max_z_int))+max_value
    24702511               if (over .EQ. 1) then
    24712512                  res@xyDashPattern  = 1
     
    24892530            end if
    24902531            if (vNam(varn) .EQ. "wqv") then
    2491                miniwqv=min(data(varn,:,:))-min_value
    2492                maxiwqv=max(data(varn,:,:))+max_value
     2532               miniwqv=min(data(varn,:,min_z_int:max_z_int))-min_value
     2533               maxiwqv=max(data(varn,:,min_z_int:max_z_int))+max_value
    24932534               if (over .EQ. 1) then
    24942535                  res@xyDashPattern  = 2
     
    25132554
    25142555            if (vNam(varn) .EQ. "wpsp") then
    2515                miniwpsp=min(data(varn,:,:))-min_value
    2516                maxiwpsp=max(data(varn,:,:))+max_value
     2556               miniwpsp=min(data(varn,:,min_z_int:max_z_int))-min_value
     2557               maxiwpsp=max(data(varn,:,min_z_int:max_z_int))+max_value
    25172558               if (over .EQ. 1) then
    25182559                  res@xyDashPattern  = 0
     
    25362577            end if
    25372578            if (vNam(varn) .EQ. "wsss") then
    2538                miniwsss=min(data(varn,:,:))-min_value
    2539                maxiwsss=max(data(varn,:,:))+max_value
     2579               miniwsss=min(data(varn,:,min_z_int:max_z_int))-min_value
     2580               maxiwsss=max(data(varn,:,min_z_int:max_z_int))+max_value
    25402581               if (over .EQ. 1) then
    25412582                  res@xyDashPattern  = 1
     
    25592600            end if
    25602601            if (vNam(varn) .EQ. "ws") then
    2561                miniws=min(data(varn,:,:))-min_value
    2562                maxiws=max(data(varn,:,:))+max_value
     2602               miniws=min(data(varn,:,min_z_int:max_z_int))-min_value
     2603               maxiws=max(data(varn,:,min_z_int:max_z_int))+max_value
    25632604               if (over .EQ. 1) then
    25642605                  res@xyDashPattern  = 2
     
    25832624
    25842625            if (vNam(varn) .EQ. "wpsap") then
    2585                miniwpsap=min(data(varn,:,:))-min_value
    2586                maxiwpsap=max(data(varn,:,:))+max_value
     2626               miniwpsap=min(data(varn,:,min_z_int:max_z_int))-min_value
     2627               maxiwpsap=max(data(varn,:,min_z_int:max_z_int))+max_value
    25872628               if (over .EQ. 1) then
    25882629                  res@xyDashPattern  = 0
     
    26062647            end if
    26072648            if (vNam(varn) .EQ. "wssas") then
    2608                miniwssas=min(data(varn,:,:))-min_value
    2609                maxiwssas=max(data(varn,:,:))+max_value
     2649               miniwssas=min(data(varn,:,min_z_int:max_z_int))-min_value
     2650               maxiwssas=max(data(varn,:,min_z_int:max_z_int))+max_value
    26102651               if (over .EQ. 1) then
    26112652                  res@xyDashPattern  = 1
     
    26292670            end if
    26302671            if (vNam(varn) .EQ. "wsa") then
    2631                miniwsa=min(data(varn,:,:))-min_value
    2632                maxiwsa=max(data(varn,:,:))+max_value
     2672               miniwsa=min(data(varn,:,min_z_int:max_z_int))-min_value
     2673               maxiwsa=max(data(varn,:,min_z_int:max_z_int))+max_value
    26332674               if (over .EQ. 1) then
    26342675                  res@xyDashPattern  = 2
     
    26532694
    26542695            if (vNam(varn) .EQ. "us2") then
    2655                minius2=min(data(varn,:,:))-min_value
    2656                maxius2=max(data(varn,:,:))+max_value
     2696               minius2=min(data(varn,:,min_z_int:max_z_int))-min_value
     2697               maxius2=max(data(varn,:,min_z_int:max_z_int))+max_value
    26572698               if (over .EQ. 1) then
    26582699                  res@xyDashPattern  = 0
     
    26762717            end if
    26772718            if (vNam(varn) .EQ. "vs2") then
    2678                minivs2=min(data(varn,:,:))-min_value
    2679                maxivs2=max(data(varn,:,:))+max_value
     2719               minivs2=min(data(varn,:,min_z_int:max_z_int))-min_value
     2720               maxivs2=max(data(varn,:,min_z_int:max_z_int))+max_value
    26802721               if (over .EQ. 1) then
    26812722                  res@xyDashPattern  = 1
     
    26992740            end if
    27002741            if (vNam(varn) .EQ. "ws2") then
    2701                miniws2=min(data(varn,:,:))-min_value
    2702                maxiws2=max(data(varn,:,:))+max_value
     2742               miniws2=min(data(varn,:,min_z_int:max_z_int))-min_value
     2743               maxiws2=max(data(varn,:,min_z_int:max_z_int))+max_value
    27032744               if (over .EQ. 1) then
    27042745                  res@xyDashPattern  = 2
     
    27232764
    27242765            if (vNam(varn) .EQ. "wsususodz") then
    2725                miniwsususodz=min(data(varn,:,:))-min_value
    2726                maxiwsususodz=max(data(varn,:,:))+max_value
     2766               miniwsususodz=min(data(varn,:,min_z_int:max_z_int))-min_value
     2767               maxiwsususodz=max(data(varn,:,min_z_int:max_z_int))+max_value
    27272768               if (over .EQ. 1) then
    27282769                  res@xyDashPattern  = 0
     
    27462787            end if
    27472788            if (vNam(varn) .EQ. "wspsodz") then
    2748                miniwspsodz=min(data(varn,:,:))-min_value
    2749                maxiwspsodz=max(data(varn,:,:))+max_value
     2789               miniwspsodz=min(data(varn,:,min_z_int:max_z_int))-min_value
     2790               maxiwspsodz=max(data(varn,:,min_z_int:max_z_int))+max_value
    27502791               if (over .EQ. 1) then
    27512792                  res@xyDashPattern  = 1
     
    27692810            end if
    27702811            if (vNam(varn) .EQ. "wpeodz") then
    2771                miniwpeodz=min(data(varn,:,:))-min_value
    2772                maxiwpeodz=max(data(varn,:,:))+max_value
     2812               miniwpeodz=min(data(varn,:,min_z_int:max_z_int))-min_value
     2813               maxiwpeodz=max(data(varn,:,min_z_int:max_z_int))+max_value
    27732814               if (over .EQ. 1) then
    27742815                  res@xyDashPattern  = 2
     
    28162857   end do
    28172858   ;#########ENDE DO LOOP FOR no_files GT 1#############
    2818    
     2859
     2860   if (isStrSubset(data@long_name," SR " ) .and. over_remind) then
     2861      print(" ")
     2862      print("If you have outputs of statistic regions you cannot overlay variables;")
     2863      print("'over' is set to 0" )
     2864      print(" ")
     2865      over = 0
     2866   end if
     2867
    28192868   if (count_var .EQ. 0) then
    28202869      print(" ")
     
    28242873      print(" ")
    28252874      exit       
    2826    end if
    2827 
    2828    if (isStrSubset(data@long_name," SR " ))then
    2829       print(" ")
    2830       print("If you have outputs of statistic regions you cannot overlay variables;")
    2831       print("'over' is set to 0" )
    2832       print(" ")
    2833       over = 0
    28342875   end if
    28352876   
     
    30023043            z=z/norm_z
    30033044           
    3004             if (abs(min(data(varn,:,:))) .GT. 10)then
    3005                min_value = abs(0.01*min(data(varn,:,:)))
    3006                max_value = abs(0.01*max(data(varn,:,:)))
     3045            if (abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 10)then
     3046               min_value = abs(0.01*min(data(varn,:,min_z_int:max_z_int)))
     3047               max_value = abs(0.01*max(data(varn,:,min_z_int:max_z_int)))
    30073048            else
    3008                if (abs(min(data(varn,:,:))) .LT. 0.01 .AND. abs(max(data(varn,:,:))) .GT. 0.01)then
    3009                   min_value = abs(0.1*max(data(varn,:,:)))
    3010                   max_value = abs(0.1*max(data(varn,:,:)))
     3049               if (abs(min(data(varn,:,min_z_int:max_z_int))) .LT. 0.01 .AND. \
     3050                   abs(max(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then
     3051                  min_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
     3052                  max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
    30113053               else
    3012                   if (abs(max(data(varn,:,:))) .LT. 0.01 .AND. abs(min(data(varn,:,:))) .GT. 0.01)then
    3013                      min_value = abs(0.1*min(data(varn,:,:)))
    3014                      max_value = abs(0.1*min(data(varn,:,:)))
     3054                  if (abs(max(data(varn,:,:))) .LT. 0.01 .AND. \
     3055                      abs(min(data(varn,:,min_z_int:max_z_int))) .GT. 0.01)then
     3056                     min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
     3057                     max_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
    30153058                  else
    3016                      min_value = abs(0.1*min(data(varn,:,:)))
    3017                      max_value = abs(0.1*max(data(varn,:,:)))
     3059                     min_value = abs(0.1*min(data(varn,:,min_z_int:max_z_int)))
     3060                     max_value = abs(0.1*max(data(varn,:,min_z_int:max_z_int)))
    30183061                  end if
    30193062               end if
    30203063            end if
    3021             if (min(data(varn,:,:)) .EQ. 0 .AND. max(data(varn,:,:)) .EQ. 0)then
     3064            if (min(data(varn,:,min_z_int:max_z_int)) .EQ. 0 .AND.\
     3065                max(data(varn,:,min_z_int:max_z_int)) .EQ. 0)then
    30223066               min_value = 0.1
    30233067               max_value = 0.1
     
    30323076
    30333077            if (xs .EQ. -1) then
    3034                res@trXMinF = min(data(varn,:,:))-min_value
     3078               res@trXMinF = min(data(varn,:,min_z_int:max_z_int))-min_value
    30353079            else
    30363080               res@trXMinF = xs
    30373081            end if
    30383082            if (xe .EQ. -1) then
    3039                res@trXMaxF = max(data(varn,:,:))+max_value
     3083               res@trXMaxF = max(data(varn,:,min_z_int:max_z_int))+max_value
    30403084            else
    30413085               res@trXMaxF = xe 
Note: See TracChangeset for help on using the changeset viewer.