Changeset 195 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Sep 2, 2008 8:54:58 AM (16 years ago)
Author:
steinfeld
Message:

Bug fix in the NCL scripts

File:
1 edited

Legend:

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

    r194 r195  
    172172   end do
    173173   if (.not. check_3d) then
    174       if (xyc .EQ. 1 .AND. .not. checkxy) then
    175          print(" ")
    176          print("Your input file doesn't have values for xy cross-sections")
    177          if (checkxz)then
    178             print("Select another input file or xz cross-sections")
     174      if (xyc .EQ. 1)then
     175         if (.not. checkxy) then
     176            print(" ")
     177            print("Your input file doesn't have values for xy cross-sections")
     178            if (checkxz)then
     179               print("Select another input file or xz cross-sections")
     180            else
     181               print("Select another input file or yz cross-sections")
     182            end if
     183            print(" ")
     184            exit
    179185         else
    180             print("Select another input file or yz cross-sections")
    181          end if
    182          print(" ")
    183          exit
    184       else
    185          print(" ")
    186          print("Your input file contains xy data")
    187          print(" ")
    188       end if
    189       if (xzc .EQ. 1 .AND. .not. checkxz) then
    190          print(" ")
    191          print("Your input file doesn't have values for xz cross-sections")
    192          if (checkxy)then
    193             print("Select another input file or xy cross-sections")
     186            print(" ")
     187            print("Your input file contains xy data")
     188            print(" ")
     189         end if
     190      end if
     191      if (xzc .EQ. 1)then
     192         if (.not. checkxz) then
     193            print(" ")
     194            print("Your input file doesn't have values for xz cross-sections")
     195            if (checkxy)then
     196               print("Select another input file or xy cross-sections")
     197            else
     198               print("Select another input file or yz cross-sections")
     199            end if
     200            print(" ")
     201            exit
    194202         else
    195             print("Select another input file or yz cross-sections")
    196          end if
    197          print(" ")
    198          exit
    199       else
    200          print(" ")
    201          print("Your input file contains xz data")
    202          print(" ")
    203       end if
    204       if (yzc .EQ. 1 .AND. .not. checkyz) then
    205          print(" ")
    206          print("Your input file doesn't have values for yz cross-sections")
    207          if (checkxy)then
    208             print("Select another input file or xy cross-sections")
     203            print(" ")
     204            print("Your input file contains xz data")
     205            print(" ")
     206         end if
     207      end if
     208      if (yzc .EQ. 1)then
     209         if (.not. checkyz) then
     210            print(" ")
     211            print("Your input file doesn't have values for yz cross-sections")
     212            if (checkxy)then
     213               print("Select another input file or xy cross-sections")
     214            else
     215               print("Select another input file or xz cross-sections")
     216            end if
     217            print(" ")
     218            exit
    209219         else
    210             print("Select another input file or xz cross-sections")
    211          end if
    212          print(" ")
    213          exit
    214       else
    215          print(" ")
    216          print("Your input file contains yz data")
    217          print(" ")
     220            print(" ")
     221            print("Your input file contains yz data")
     222            print(" ")
     223         end if
    218224      end if
    219225   else
     
    523529      end do
    524530   end if
    525   
     531 
    526532   ; ****************************************************
    527533   ; set up ranges of x-, y- and z-coordinates
     
    530536   if (xs .EQ. -1.d) then               
    531537      xs = x_d(0)
     538      if (delta_x .EQ. -1) then
     539         print(" ")
     540         print("You cannot choose a start value for x, there are preseted layers for x")
     541         print(" ")
     542         xstart=0
     543      end if
    532544   else
    533545      if (delta_x .EQ. -1) then
     
    563575   if (ys .EQ. -1.d) then   
    564576      ys = y_d(0)
     577      if (delta_y .EQ. -1) then
     578         print(" ")
     579         print("You cannot choose a start value for y, there are preseted layers for y")
     580         print(" ")
     581         ystart=0
     582      end if
    565583   else
    566584      if (delta_y .EQ. -1) then
     
    596614   if (zs .EQ. -1) then                         
    597615      zs = 0
     616      if (delta_z .EQ. -1) then
     617         print(" ")
     618         print("You cannot choose a start value for z, there are preseted layers for z")
     619         print(" ")
     620      end if
    598621   else
    599622      if (delta_z .EQ. -1) then
     
    628651   if (xe .EQ. -1) then   
    629652      xe = x_d(xdim)
     653      if (delta_x .EQ. -1) then
     654         print(" ")
     655         print("You cannot choose an end value for x, there are preseted layers for x")
     656         print(" ")
     657         xend=xdim
     658      end if
    630659   else
    631660      if (delta_x .EQ. -1) then
     
    679708   if (ye .EQ. -1) then 
    680709      ye = y_d(ydim)
     710      if (delta_y .EQ. -1) then
     711         print(" ")
     712         print("You cannot choose an end value for y, there are preseted layers for y")
     713         print(" ")
     714         yend=ydim
     715      end if
    681716   else
    682717      if (delta_y .EQ. -1) then
     
    730765   if (ze .EQ. -1) then 
    731766      ze = zdim
     767      if (delta_z .EQ. -1) then
     768         print(" ")
     769         print("You cannot choose an end value for z, there are preseted layers for z")
     770         print(" ")
     771         ze = zdim
     772      end if
    732773   else
    733774      if (delta_z .EQ. -1) then
     
    809850   
    810851   if( shape .eq. 1 ) then
    811       cs_res@vpWidthF    = (xe-xs)/(ye-ys)       
    812       cs_res@vpHeightF   = 1
     852      if (xyc .EQ. 1)then
     853         cs_res@vpWidthF    = (xe-xs)/(ye-ys)       
     854         cs_res@vpHeightF   = 1
     855      end if
     856      if (xzc .EQ. 1)then
     857         cs_res@vpWidthF    = (xe-xs)/(delta_x*(ze-zs))       
     858         cs_res@vpHeightF   = 1
     859      end if
     860      if (yzc .EQ. 1)then
     861         cs_res@vpWidthF    = (ye-ys)/(delta_y*(ze-zs))       
     862         cs_res@vpHeightF   = 1
     863      end if
    813864   end if
    814865   
Note: See TracChangeset for help on using the changeset viewer.