Changeset 195 for palm/trunk/SCRIPTS/NCL/cross_sections.ncl
 Timestamp:
 Sep 2, 2008 8:54:58 AM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SCRIPTS/NCL/cross_sections.ncl
r194 r195 172 172 end do 173 173 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 crosssections") 177 if (checkxz)then 178 print("Select another input file or xz crosssections") 174 if (xyc .EQ. 1)then 175 if (.not. checkxy) then 176 print(" ") 177 print("Your input file doesn't have values for xy crosssections") 178 if (checkxz)then 179 print("Select another input file or xz crosssections") 180 else 181 print("Select another input file or yz crosssections") 182 end if 183 print(" ") 184 exit 179 185 else 180 print("Select another input file or yz crosssections") 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 crosssections") 192 if (checkxy)then 193 print("Select another input file or xy crosssections") 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 crosssections") 195 if (checkxy)then 196 print("Select another input file or xy crosssections") 197 else 198 print("Select another input file or yz crosssections") 199 end if 200 print(" ") 201 exit 194 202 else 195 print("Select another input file or yz crosssections") 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 crosssections") 207 if (checkxy)then 208 print("Select another input file or xy crosssections") 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 crosssections") 212 if (checkxy)then 213 print("Select another input file or xy crosssections") 214 else 215 print("Select another input file or xz crosssections") 216 end if 217 print(" ") 218 exit 209 219 else 210 print("Select another input file or xz crosssections") 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 218 224 end if 219 225 else … … 523 529 end do 524 530 end if 525 531 526 532 ; **************************************************** 527 533 ; set up ranges of x, y and zcoordinates … … 530 536 if (xs .EQ. 1.d) then 531 537 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 532 544 else 533 545 if (delta_x .EQ. 1) then … … 563 575 if (ys .EQ. 1.d) then 564 576 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 565 583 else 566 584 if (delta_y .EQ. 1) then … … 596 614 if (zs .EQ. 1) then 597 615 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 598 621 else 599 622 if (delta_z .EQ. 1) then … … 628 651 if (xe .EQ. 1) then 629 652 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 630 659 else 631 660 if (delta_x .EQ. 1) then … … 679 708 if (ye .EQ. 1) then 680 709 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 681 716 else 682 717 if (delta_y .EQ. 1) then … … 730 765 if (ze .EQ. 1) then 731 766 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 732 773 else 733 774 if (delta_z .EQ. 1) then … … 809 850 810 851 if( shape .eq. 1 ) then 811 cs_res@vpWidthF = (xexs)/(yeys) 812 cs_res@vpHeightF = 1 852 if (xyc .EQ. 1)then 853 cs_res@vpWidthF = (xexs)/(yeys) 854 cs_res@vpHeightF = 1 855 end if 856 if (xzc .EQ. 1)then 857 cs_res@vpWidthF = (xexs)/(delta_x*(zezs)) 858 cs_res@vpHeightF = 1 859 end if 860 if (yzc .EQ. 1)then 861 cs_res@vpWidthF = (yeys)/(delta_y*(zezs)) 862 cs_res@vpHeightF = 1 863 end if 813 864 end if 814 865
Note: See TracChangeset
for help on using the changeset viewer.