Changeset 514


Ignore:
Timestamp:
Mar 17, 2010 9:42:33 AM (14 years ago)
Author:
heinze
Message:

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

File:
1 edited

Legend:

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

    r513 r514  
    13181318            temp = f[:]->$vNam(varn)$
    13191319            data_att = f_att->$vNam(varn)$
    1320             if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")
    1321               ;these variables depend zu1_xy and that's why they have only one z-layer
     1320            if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"     \
     1321              .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
     1322              .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
     1323              .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
     1324              ;these variables depend on zu1_xy and that's why they have only one z-layer
    13221325              data(varn,:,0,:,:)=temp(:,0,ys:ye,xs:xe)
    13231326              no_zu1=no_zu1+1
     
    16001603                     end if
    16011604
    1602                      if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")then
     1605                     if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"       \
     1606                         .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
     1607                         .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
     1608                         .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
    16031609                         loe = 0
    16041610                         level = "=" + zu1(0) + "m"
    16051611                     end if
    16061612                   
    1607                      cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level                                
     1613                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level           
    16081614                     plot(n) = gsn_csm_contour(wks_ps,data(varn,li,lo-los,:,:),cs_res)
    16091615                     if (vector .EQ. 1 .AND. check_vecp) then
     
    16331639                     end if
    16341640       
    1635                      if(vNam(varn) .eq. "ts_xy" .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy")then
     1641                     if(vNam(varn) .eq. "lwps_xy" .or. vNam(varn) .eq. "pras_xy"       \
     1642                         .or. vNam(varn) .eq. "prrs_xy" .or. vNam(varn) .eq. "qsws_xy" \
     1643                         .or. vNam(varn) .eq. "shfs_xy" .or. vNam(varn) .eq. "ts_xy"   \
     1644                         .or. vNam(varn) .eq. "us_xy" .or. vNam(varn) .eq. "z0s_xy") then
    16361645                         lie = 0
    16371646                         level = "=" + zu1(0) + "m"
Note: See TracChangeset for help on using the changeset viewer.