Ignore:
Timestamp:
Sep 15, 2011 12:03:51 PM (11 years ago)
Author:
heinze
Message:

Bugfixes for NCL scripts timeseries.ncl and cross_sections.ncl concerning usage of norm_t and vector plot

File:
1 edited

Legend:

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

    r585 r758  
    14151415   end if
    14161416
    1417    check = True
    1418    v1=0
    1419    v2=0
    1420    no_var=0
    1421    n=0
    1422    no_zu1=0
     1417   check  = True
     1418   v1     = 0
     1419   v2     = 0
     1420   no_var = 0
     1421   n      = 0
     1422   no_zu1 = 0
    14231423
    14241424   ;****************************************************
     
    15471547   
    15481548   do varn=dim-1,0,1   
    1549    
    1550       ;if (vector .EQ. 1) then   
    1551       ;   check_vec1 = isStrSubset( vec1,","+vNam(varn)+",")
    1552       ;   check_vec2 = isStrSubset( vec2,","+vNam(varn)+",")
    1553       ;end if
    15541549           
    15551550      if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or.          \
     
    16301625
    16311626   if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
    1632       plot=new((/no_time*no_layer*no_var - no_time*(no_layer-1)*no_zu1 \
    1633                                          + no_time*no_layer/),graphic)
     1627      plot=new((/no_time*no_layer/),graphic)
    16341628   else
    16351629      plot=new((/no_time*no_layer*no_var - no_time*(no_layer-1)*no_zu1/)\
     
    18391833                     cs_res@gsnCenterString = "t=" + \
    18401834                       decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level
    1841                      plot(n) = gsn_csm_contour(wks_ps,\
     1835
     1836                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     1837                        ;nothing to be done
     1838                     else
     1839                         plot(n) = gsn_csm_contour(wks_ps,\
    18421840                                           data(varn,li,lo-los,:,:),cs_res)
     1841                     end if
     1842
     1843
    18431844                     if (vector .EQ. 1 .AND. check_vecp) then
    18441845                        vecres                 = True           ; vector only
     
    18971898                     cs_res@gsnCenterString = "t=" + \
    18981899                           decimalPlaces(t_all(lo)/3600,2,True) + "h  z"+ level
    1899                      
    1900                      plot(n) = gsn_csm_contour(wks_ps,\
     1900
     1901                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     1902                        ;nothing to be done
     1903                     else
     1904                        plot(n) = gsn_csm_contour(wks_ps,\
    19011905                                            data(varn,lo,li-lis,:,:),cs_res)
     1906                     end if
     1907
    19021908                     if (vector .EQ. 1 .AND. check_vecp) then
    19031909                        vecres                 = True           ; vector only
     
    19411947                     cs_res@gsnCenterString = "t=" + \
    19421948                           decimalPlaces(t_all(li)/3600,2,True) + "h  y"+ level
    1943                      plot(n) = gsn_csm_contour(wks_ps,\
     1949
     1950                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     1951                        ;nothing to be done
     1952                     else
     1953                         plot(n) = gsn_csm_contour(wks_ps,\
    19441954                                              data(varn,li,:,lo-los,:),cs_res)
     1955                     end if
     1956
    19451957                     if (vector .EQ. 1 .AND. check_vecp) then
    19461958                        vecres                 = True           ; vector only
     
    19751987                     cs_res@gsnCenterString = "t=" + \
    19761988                           decimalPlaces(t_all(lo)/3600,2,True) + "h  y"+ level
    1977                      plot(n) = gsn_csm_contour(wks_ps,\
    1978                                               data(varn,lo,:,li-lis,:),cs_res)
     1989
     1990                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     1991                        ;nothing to be done
     1992                     else
     1993                         plot(n) = gsn_csm_contour(wks_ps,\
     1994                                              data(varn,lo,:,li-lis,:),cs_res) 
     1995                     end if
     1996
    19791997                     if (vector .EQ. 1 .AND. check_vecp) then
    19801998                        vecres                 = True           ; vector only
     
    20182036                     cs_res@gsnCenterString = "t=" + \
    20192037                          decimalPlaces(t_all(li)/3600,2,True) + "h  x"+ level
    2020                      plot(n) = gsn_csm_contour(wks_ps,\
     2038
     2039                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     2040                        ;nothing to be done
     2041                     else
     2042                        plot(n) = gsn_csm_contour(wks_ps,\
    20212043                                          data(varn,li,:,:,lo-los),cs_res)
     2044                     end if
     2045
    20222046                     if (vector .EQ. 1 .AND. check_vecp) then
    20232047                        vecres                 = True           ; vector only
     
    20522076                     cs_res@gsnCenterString = "t=" + \
    20532077                           decimalPlaces(t_all(lo)/3600,2,True) + "h  x"+ level
    2054                      plot(n) = gsn_csm_contour(wks_ps,\
     2078
     2079                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
     2080                        ;nothing to be done
     2081                     else
     2082                        plot(n) = gsn_csm_contour(wks_ps,\
    20552083                                          data(varn,lo,:,:,li-lis),cs_res)
     2084                     end if
     2085
    20562086                     if (vector .EQ. 1 .AND. check_vecp)then
    20572087                        vecres                 = True           ; vector only
Note: See TracChangeset for help on using the changeset viewer.