Changeset 161 for palm/trunk/SCRIPTS/NCL
- Timestamp:
- Apr 11, 2008 10:54:53 AM (17 years ago)
- Location:
- palm/trunk/SCRIPTS/NCL
- Files:
-
- 5 edited
-
.ncl_preferences (modified) (3 diffs)
-
UserGuide_for_NCLscripts.pdf (modified) (previous)
-
cross_sections.ncl (modified) (23 diffs)
-
profiles.ncl (modified) (20 diffs)
-
timeseries.ncl (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/.ncl_preferences
r157 r161 19 19 ___no_lines__________[1;2;...]__________________________________default = 2_______________[ALL] 20 20 2 21 ___'var'_____________[ E,dt,umax,;...]___________________________default = variables_______[ALL]21 ___'var'_____________[,E,dt,umax,;...]__________________________default = variables_______[ALL] 22 22 variables 23 23 ___combine___________[0] or [1]_________________________________default = 0_______________[PR] 24 24 0 25 ___number_comb_______[ 0] or [1]_________________________________default = 0_______________[PR]25 ___number_comb_______[2] or [3]_________________________________default = 0_______________[PR] 26 26 0 27 ___'c_var'___________[ umax,vmax,wmax,;...]______________________default = c_variables_____[PR]27 ___'c_var'___________[,us2,vs2,ws2,;...]________________________default = c_variables_____[PR] 28 28 c_variables 29 29 ___dash______________[0] or [1]_________________________________default = 0_______________[PR] … … 33 33 ___min_z_____________[0;10;...]_________________________________default = 0_______________[PR] 34 34 0 35 ___max_z_____________[100;1000;...]_____________________________default = max(z_ pr)_______[PR]36 max(z_ pr)35 ___max_z_____________[100;1000;...]_____________________________default = max(z_u)________[PR] 36 max(z_u) 37 37 ___'sort'____________[time] or [layer]__________________________default = time____________[CS] 38 38 time … … 63 63 ___vector____________[0] or [1]_________________________________default = 0_______________[CS] 64 64 0 65 ___'vec1'____________[ u,;v,;...]________________________________default = vec1____________[CS]65 ___'vec1'____________[,u,;,v,;...]______________________________default = vec1____________[CS] 66 66 vec1 67 ___'vec2'____________[ u,;v,;...]________________________________default = vec2____________[CS]67 ___'vec2'____________[,u,;,v,;...]______________________________default = vec2____________[CS] 68 68 vec2 69 __'plotvec'__________[ u,;v,;...]________________________________default = plotvec_________[CS]69 __'plotvec'__________[,u,;,v,;...]______________________________default = plotvec_________[CS] 70 70 plotvec 71 71 ___ref_mag___________[0.1;0.02;...]_____________________________default = 0.05____________[CS] 72 72 0.05 73 ___over______________[0] or [1]_________________________________default = 0_______________[TS+PR] 74 0 73 75 ___end of file ".ncl_preferences"_____________________________________________________________ -
palm/trunk/SCRIPTS/NCL/cross_sections.ncl
r157 r161 12 12 13 13 if (isfilepresent("~/.ncl_preferences")) then 14 parameter = asciiread("~/.ncl_preferences",7 3,"string")14 parameter = asciiread("~/.ncl_preferences",75,"string") 15 15 delete(parameter@_FillValue) 16 16 else … … 208 208 print(" ") 209 209 dim = dimsizes(vNam) 210 210 if (dim .EQ. 0) then 211 print(" ") 212 print("There are no data on file") 213 print(" ") 214 end if 215 211 216 ; *************************************************** 212 217 ; set up recourses … … 259 264 end if 260 265 261 ; ***************************************************262 ; open workstation(s)263 ; ***************************************************264 265 wks_ps = gsn_open_wks(format_out,file_out)266 gsn_define_colormap(wks_ps,"rainbow+white")267 268 plot=new((/no_columns*no_lines/),graphic)269 270 page = 0271 272 nc1 = no_columns273 nl1 = no_lines274 275 print("plots sorted by '"+sort+"'")276 print(" ")277 278 266 ; ********************************************* 279 267 ; set up of start_time_step and end_time_step 280 268 ; ********************************************* 281 269 282 t = f->time270 t = f->time 283 271 nt = dimsizes(t) 284 272 … … 288 276 289 277 if ( .not. isvar("start_time_step") ) then 290 start_time_step = 1278 start_time_step = 0 291 279 if (parameter(13) .NE. "1") then 292 if (parameter(13) .L E. "0")280 if (parameter(13) .LT. "1") 293 281 print(" ") 294 282 print("Begin with time step 1") … … 302 290 exit 303 291 end if 304 start_time_step = stringtointeger(parameter(13)) 292 start_time_step = stringtointeger(parameter(13))-1 305 293 end if 306 294 else … … 317 305 exit 318 306 end if 307 start_time_step = start_time_step - 1 319 308 end if 320 309 … … 324 313 325 314 if ( .not. isvar("end_time_step") ) then 326 end_time_step = nt 315 end_time_step = nt-1 327 316 if (parameter(15) .NE. "nt") then 328 317 if (parameter(15) .LE. "0") … … 344 333 exit 345 334 end if 346 end_time_step = stringtointeger(parameter(15)) 335 end_time_step = stringtointeger(parameter(15))-1 347 336 end if 348 337 else … … 365 354 exit 366 355 end if 367 end if 356 end_time_step = end_time_step-1 357 end if 358 359 no_time=(end_time_step-start_time_step)+1 368 360 369 361 ; **************************************************** … … 374 366 if ( .not. isvar("plotvec") ) then 375 367 plotvec = parameter(69) 376 if (parameter(69) .EQ. "plotvec") then377 print(" ")378 print("Please indicate the variable where the vector plot shall overlay ('plotvec')")379 print(" ")380 exit381 end if382 368 end if 383 369 if ( .not. isvar("vec1") ) then … … 409 395 ; **************************************************** 410 396 411 do varn= 0,dim-1412 413 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" ) then397 do varn=dim-1,0,1 398 399 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or. vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or. vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or. vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x" .or. vNam(varn) .eq. "y" .or. vNam(varn) .eq. "zu_xy" .or. vNam(varn) .eq. "zw_xy" .or. vNam(varn) .eq. "zu1_xy" .or. vNam(varn) .eq. "ind_z_xy" .or. vNam(varn) .eq. "y_xz" .or. vNam(varn) .eq. "yv_xz" .or. vNam(varn) .eq. "ind_y_xz" .or. vNam(varn) .eq. "x_yz" .or. vNam(varn) .eq. "xu_yz" .or. vNam(varn) .eq. "ind_x_yz") then 414 400 check = False 415 end if 416 if ( vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi") then 417 check = False 418 end if 401 else 402 check = True 403 end if 419 404 if ( isvar("var") ) then 420 check = isStrSubset( var, vNam(varn)+",")405 check = isStrSubset( var,","+vNam(varn)+",") 421 406 end if 422 407 if (parameter(21) .NE. "variables") then 423 408 var = parameter(21) 424 check = isStrSubset( var, vNam(varn)+"," )409 check = isStrSubset( var,","+vNam(varn)+"," ) 425 410 end if 426 411 427 412 if(check) then 428 413 print(vNam(varn)) 429 414 data_all = f->$vNam(varn)$ 430 415 … … 536 521 end do 537 522 523 if (xyc .EQ. 1) then 524 no_layer = (ze-zs)+1 525 end if 526 if (xzc .EQ. 1) then 527 no_layer = (ye-ys)+1 528 end if 529 if (yzc .EQ. 1) then 530 no_layer = (xe-xs)+1 531 end if 532 538 533 data = new((/dim,nt,(ze-zs)+1,(ye-ys)+1,(xe-xs)+1/),float) 539 534 … … 570 565 v1=0 571 566 v2=0 572 573 do varn=0,dim-1 567 no_var=0 568 n=0 569 570 do varn=dim-1,0,1 574 571 575 572 data_all = f->$vNam(varn)$ 576 573 577 574 if (vector .EQ. 1) then 578 check_vec1 = isStrSubset( vec1, vNam(varn)+",")579 check_vec2 = isStrSubset( vec2, vNam(varn)+",")575 check_vec1 = isStrSubset( vec1,","+vNam(varn)+",") 576 check_vec2 = isStrSubset( vec2,","+vNam(varn)+",") 580 577 end if 581 578 582 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" ) then579 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or. vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or. vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or. vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x".or. vNam(varn) .eq. "y" .or. vNam(varn) .eq. "zu_xy" .or. vNam(varn) .eq. "zw_xy" .or. vNam(varn) .eq. "zu1_xy" .or. vNam(varn) .eq. "ind_z_xy" .or. vNam(varn) .eq. "y_xz" .or. vNam(varn) .eq. "yv_xz" .or. vNam(varn) .eq. "ind_y_xz" .or. vNam(varn) .eq. "x_yz" .or. vNam(varn) .eq. "xu_yz" .or. vNam(varn) .eq. "ind_x_yz") then 583 580 check = False 584 end if 585 if ( vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi") then 586 check = False 587 end if 581 else 582 check = True 583 end if 588 584 if ( isvar("var") ) then 589 check = isStrSubset( var,vNam(varn)+",") 590 end if 591 585 check = isStrSubset( var,","+vNam(varn)+",") 586 end if 592 587 if (parameter(21) .NE. "variables") then 593 588 var = parameter(21) 594 check = isStrSubset( var, vNam(varn)+"," )589 check = isStrSubset( var,","+vNam(varn)+"," ) 595 590 end if 596 591 597 592 if(check) then 598 593 594 no_var=no_var+1 595 599 596 if (vector .EQ. 1) then 600 if ( vNam(varn)+"," .EQ. vec1) then597 if (","+vNam(varn)+"," .EQ. vec1) then 601 598 v1=v1+1 602 599 end if 603 if ( vNam(varn)+"," .EQ. vec2) then600 if (","+vNam(varn)+"," .EQ. vec2) then 604 601 v2=v2+1 605 602 end if … … 608 605 data(varn,:,:,:,:)=data_all(0:nt-1,zs:ze,ys:ye,xs:xe) 609 606 if (check_vec1) then 610 print("in vec1")611 607 vect1=data_all 612 608 end if 613 609 if (check_vec2) then 614 print("in vec2")615 610 vect2=data_all 616 611 end if … … 626 621 delete(data_all) 627 622 628 end do 623 end do 624 625 plot_panel=new((/no_time*no_layer*no_var/),graphic) 629 626 630 627 if (vector .EQ. 1) then … … 645 642 646 643 ; *************************************************** 644 ; open workstation(s) 645 ; *************************************************** 646 647 wks_ps = gsn_open_wks(format_out,file_out) 648 gsn_define_colormap(wks_ps,"rainbow+white") 649 650 plot=new((/no_time*no_layer*no_var/),graphic) 651 652 page = 0 653 n=0 654 655 print("plots sorted by '"+sort+"'") 656 print(" ") 657 658 ; *************************************************** 647 659 ; create plots 648 660 ; *************************************************** … … 650 662 check = True 651 663 652 do varn= 0,dim-1664 do varn=dim-1,0,1 653 665 654 666 if (vector .EQ. 1) then 655 check_vecp = isStrSubset( plotvec, vNam(varn)+",")667 check_vecp = isStrSubset( plotvec,","+vNam(varn)+",") 656 668 end if 657 669 658 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" ) then 659 check = False 660 end if 661 if ( vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi") then 662 check = False 663 end if 670 if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or. vNam(varn) .eq. "xu" .or. vNam(varn) .eq. "zu" .or. vNam(varn) .eq. "zwwi" .or. vNam(varn) .eq. "zusi" .or. vNam(varn) .eq. "time" .or. vNam(varn) .eq. "zu_3d" .or. vNam(varn) .eq. "zw_3d" .or. vNam(varn) .eq. "x".or. vNam(varn) .eq. "y" .or. vNam(varn) .eq. "zu_xy" .or. vNam(varn) .eq. "zw_xy" .or. vNam(varn) .eq. "zu1_xy" .or. vNam(varn) .eq. "ind_z_xy" .or. vNam(varn) .eq. "y_xz" .or. vNam(varn) .eq. "yv_xz" .or. vNam(varn) .eq. "ind_y_xz" .or. vNam(varn) .eq. "x_yz" .or. vNam(varn) .eq. "xu_yz" .or. vNam(varn) .eq. "ind_x_yz") then 671 check = False 672 else 673 check = True 674 end if 664 675 if ( isvar("var") ) then 665 check = isStrSubset( var, vNam(varn)+",")676 check = isStrSubset( var,","+vNam(varn)+",") 666 677 end if 667 678 if (parameter(21) .NE. "variables") then 668 679 var = parameter(21) 669 check = isStrSubset( var, vNam(varn)+"," )680 check = isStrSubset( var,","+vNam(varn)+"," ) 670 681 end if 671 682 672 683 if(check) then 673 684 674 685 ; **************************************************** 675 686 ; loops over time and layer 676 687 ; **************************************************** 677 688 678 do lo = los, loe ; lo: loop outer 679 do li = lis, lie, no_lines*no_columns ; li: loop inner 680 681 if ( sort .eq. "time" ) then 682 ta = li 683 tb = min( (/li + no_lines*no_columns - 1, end_time_step/) ) 684 if ( xyc .eq. 1 ) then 685 xa = xs 686 xb = xe 687 ya = ys 688 yb = ye 689 za = lo 690 zb = lo 691 end if 692 if ( xzc .eq. 1 ) then 693 xa = xs 694 xb = xe 695 ya = lo 696 yb = lo 697 za = zs 698 zb = ze 699 end if 700 if ( yzc .eq. 1 ) then 701 xa = lo 702 xb = lo 703 ya = ys 704 yb = ye 705 za = zs 706 zb = ze 707 end if 708 end if 709 710 if ( sort .eq. "layer" ) then 711 ta = lo 712 tb = lo 713 if ( xyc .eq. 1 ) then 714 xa = xs 715 xb = xe 716 ya = ys 717 yb = ye 718 za = li 719 zb = min( (/li + no_lines*no_columns - 1, ze/) ) 720 end if 721 if ( xzc .eq. 1 ) then 722 xa = xs 723 xb = xe 724 ya = li 725 yb = min( (/li + no_lines*no_columns - 1, ye/) ) 726 za = zs 727 zb = ze 728 end if 729 if ( yzc .eq. 1 ) then 730 xa = li 731 xb = min( (/li + no_lines*no_columns - 1, xe/) ) 732 ya = ys 733 yb = ye 734 za = zs 735 zb = ze 736 end if 737 end if 738 689 do lo = los, loe 690 do li = lis, lie 691 739 692 ; **************************************************** 740 693 ; xy cross section … … 743 696 if ( xyc .eq. 1 ) then 744 697 745 cs_res@tiXAxisString 746 cs_res@tiYAxisString 698 cs_res@tiXAxisString = "x [m]" 699 cs_res@tiYAxisString = "y [m]" 747 700 cs_res@gsnLeftString = "Plot of "+vNam(varn) 748 701 749 702 if ( sort .eq. "time" ) then 750 703 if ( data&z(lo) .eq. -1 ) then … … 753 706 level = "=" + data&z(lo) + "m" 754 707 end if 755 do n = 0,tb-ta756 cs_res@gsnCenterString = "time=" + data&t(li+n-1) +"s z"+level757 plot1 = gsn_csm_contour(wks_ps,data(varn,li+n-1,lo,:,:),cs_res)708 cs_res@gsnCenterString = "time=" + data&t(li) +"s z"+level 709 plot(n) = gsn_csm_contour(wks_ps,data(varn,li,lo,:,:),cs_res) 710 if (vector .EQ. 1) then 758 711 if (check_vecp) 759 712 vecres = True ; vector only resources 760 713 vecres@gsnDraw = False ; don't draw 761 762 763 764 714 vecres@gsnFrame = False ; don't advance frame 715 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 716 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 717 vecres@vcRefLengthF = 0.05 ; define length of vec ref 765 718 vecres@gsnRightString = " " ; turn off right string 766 719 vecres@gsnLeftString = " " ; turn off left string 767 720 vecres@tiXAxisString = " " 768 plot_vec=gsn_csm_vector(wks_ps,vect1(li+n-1,lo,:,:),vect2(li+n-1,lo,:,:),vecres) 769 overlay(plot1, plot_vec) 770 end if 771 plot(n) = plot1 772 end do 721 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 722 overlay(plot(n), plot_vec) 723 else 724 vecres = True ; vector only resources 725 vecres@gsnDraw = False ; don't draw 726 vecres@gsnFrame = False ; don't advance frame 727 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 728 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 729 vecres@vcRefLengthF = 0.05 ; define length of vec ref 730 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 731 vecres@gsnLeftString = " " ; turn off left string 732 vecres@tiXAxisString = " " 733 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 734 plot(n) = plot_vec 735 end if 736 end if 773 737 end if 774 738 775 739 if ( sort .eq. "layer" ) then 776 do n = 0, zb-za 777 if ( data&z(li+n) .eq. -1 ) then 778 level = "-average" 779 else 780 level = "=" + data&z(li+n) + "m" 740 if ( data&z(li) .eq. -1 ) then 741 level = "-average" 742 else 743 level = "=" + data&z(li) + "m" 744 end if 745 cs_res@gsnCenterString = "t=" + data&t(lo) + "s z"+ level 746 plot(n) = gsn_csm_contour(wks_ps,data(varn,lo,li,:,:),cs_res) 747 if (vector .EQ. 1) then 748 if (check_vecp) 749 vecres = True ; vector only resources 750 vecres@gsnDraw = False ; don't draw 751 vecres@gsnFrame = False ; don't advance frame 752 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 753 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 754 vecres@vcRefLengthF = 0.05 ; define length of vec ref 755 vecres@gsnRightString = " " ; turn off right string 756 vecres@gsnLeftString = " " ; turn off left string 757 vecres@tiXAxisString = " " 758 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 759 overlay(plot(n), plot_vec) 760 else 761 vecres = True ; vector only resources 762 vecres@gsnDraw = False ; don't draw 763 vecres@gsnFrame = False ; don't advance frame 764 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 765 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 766 vecres@vcRefLengthF = 0.05 ; define length of vec ref 767 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 768 vecres@gsnLeftString = " " ; turn off left string 769 vecres@tiXAxisString = " " 770 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 771 plot(n) = plot_vec 781 772 end if 782 cs_res@gsnCenterString = "t=" + data&t(lo-1) + "s z"+ level 783 plot1 = gsn_csm_contour(wks_ps,data(varn,lo-1,li+n,:,:),cs_res) 784 if (check_vecp) 785 vecres = True ; vector only resources 786 vecres@gsnDraw = False ; don't draw 787 vecres@gsnFrame = False ; don't advance frame 788 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 789 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 790 vecres@vcRefLengthF = 0.05 ; define length of vec ref 791 vecres@gsnRightString = " " ; turn off right string 792 vecres@gsnLeftString = " " ; turn off left string 793 vecres@tiXAxisString = " " 794 plot_vec=gsn_csm_vector(wks_ps,vect1(lo-1,li+n,:,:),vect2(lo-1,li+n,:,:),vecres) 795 overlay(plot1, plot_vec) 796 end if 797 plot(n) = plot1 798 end do 773 end if 799 774 end if 800 801 775 end if 802 776 … … 817 791 level = "=" + data&y(lo) + "m" 818 792 end if 819 do n = 0, tb-ta820 cs_res@gsnCenterString = "t=" + data&t(li+n-1) + "s y"+ level821 plot1 = gsn_csm_contour(wks_ps,data(varn,li+n-1,:,lo,:),cs_res)793 cs_res@gsnCenterString = "t=" + data&t(li) + "s y"+ level 794 plot(n) = gsn_csm_contour(wks_ps,data(varn,li,:,lo,:),cs_res) 795 if (vector .EQ. 1) then 822 796 if (check_vecp) 823 797 vecres = True ; vector only resources … … 830 804 vecres@gsnLeftString = " " ; turn off left string 831 805 vecres@tiXAxisString = " " 832 plot_vec=gsn_csm_vector(wks_ps,vect1(li+n-1,lo,:,:),vect2(li+n-1,lo,:,:),vecres) 833 overlay(plot1, plot_vec) 806 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 807 overlay(plot(n), plot_vec) 808 else 809 vecres = True ; vector only resources 810 vecres@gsnDraw = False ; don't draw 811 vecres@gsnFrame = False ; don't advance frame 812 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 813 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 814 vecres@vcRefLengthF = 0.05 ; define length of vec ref 815 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 816 vecres@gsnLeftString = " " ; turn off left string 817 vecres@tiXAxisString = " " 818 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 819 plot(n) = plot_vec 834 820 end if 835 plot(n) = plot1 836 end do 821 end if 837 822 end if 838 823 839 824 if ( sort .eq. "layer" ) then 840 do n = 0, yb-ya 841 if ( data&z(li+n) .eq. -1 ) then 842 level = "-average" 843 else 844 level = "=" + data&y(li+n) + "m" 825 if ( data&z(li) .eq. -1 ) then 826 level = "-average" 827 else 828 level = "=" + data&y(li) + "m" 829 end if 830 cs_res@gsnCenterString = "t=" + data&t(lo) + "s y"+ level 831 plot(n) = gsn_csm_contour(wks_ps,data(varn,lo,:,li,:),cs_res) 832 if (vector .EQ. 1) then 833 if (check_vecp) 834 vecres = True ; vector only resources 835 vecres@gsnDraw = False ; don't draw 836 vecres@gsnFrame = False ; don't advance frame 837 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 838 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 839 vecres@vcRefLengthF = 0.05 ; define length of vec ref 840 vecres@gsnRightString = " " ; turn off right string 841 vecres@gsnLeftString = " " ; turn off left string 842 vecres@tiXAxisString = " " 843 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 844 overlay(plot(n), plot_vec) 845 else 846 vecres = True ; vector only resources 847 vecres@gsnDraw = False ; don't draw 848 vecres@gsnFrame = False ; don't advance frame 849 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 850 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 851 vecres@vcRefLengthF = 0.05 ; define length of vec ref 852 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 853 vecres@gsnLeftString = " " ; turn off left string 854 vecres@tiXAxisString = " " 855 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 856 plot(n) = plot_vec 845 857 end if 846 cs_res@gsnCenterString = "t=" + data&t(lo-1) + "s y"+ level 847 plot1 = gsn_csm_contour(wks_ps,data(varn,lo-1,:,li+n,:),cs_res) 848 if (check_vecp) 849 vecres = True ; vector only resources 850 vecres@gsnDraw = False ; don't draw 851 vecres@gsnFrame = False ; don't advance frame 852 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 853 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 854 vecres@vcRefLengthF = 0.05 ; define length of vec ref 855 vecres@gsnRightString = " " ; turn off right string 856 vecres@gsnLeftString = " " ; turn off left string 857 vecres@tiXAxisString = " " 858 plot_vec=gsn_csm_vector(wks_ps,vect1(li+n-1,lo,:,:),vect2(li+n-1,lo,:,:),vecres) 859 overlay(plot1, plot_vec) 860 end if 861 plot(n) = plot1 862 end do 858 end if 863 859 end if 864 865 860 end if 866 861 … … 881 876 level = "=" + data&z(lo) + "m" 882 877 end if 883 do n = 0, tb-ta884 cs_res@gsnCenterString = "t=" + data&t(li+n-1) + "s x"+ level885 plot1 = gsn_csm_contour(wks_ps,data(varn,li+n-1,:,:,lo),cs_res)878 cs_res@gsnCenterString = "t=" + data&t(li) + "s x"+ level 879 plot(n) = gsn_csm_contour(wks_ps,data(varn,li,:,:,lo),cs_res) 880 if (vector .EQ. 1) then 886 881 if (check_vecp) 887 882 vecres = True ; vector only resources … … 894 889 vecres@gsnLeftString = " " ; turn off left string 895 890 vecres@tiXAxisString = " " 896 plot_vec=gsn_csm_vector(wks_ps,vect1(li+n-1,lo,:,:),vect2(li+n-1,lo,:,:),vecres) 897 overlay(plot1, plot_vec) 891 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 892 overlay(plot(n), plot_vec) 893 else 894 vecres = True ; vector only resources 895 vecres@gsnDraw = False ; don't draw 896 vecres@gsnFrame = False ; don't advance frame 897 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 898 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 899 vecres@vcRefLengthF = 0.05 ; define length of vec ref 900 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 901 vecres@gsnLeftString = " " ; turn off left string 902 vecres@tiXAxisString = " " 903 plot_vec=gsn_csm_vector(wks_ps,vect1(li,lo,:,:),vect2(li,lo,:,:),vecres) 904 plot(n) = plot_vec 898 905 end if 899 plot(n) = plot1 900 end do 906 end if 901 907 end if 902 908 903 909 if ( sort .eq. "layer" ) then 904 do n = 0, xb-xa 905 if ( data&z(li+n) .eq. -1 ) then 906 level = "-average" 907 else 908 level = "=" + data&z(li+n) + "m" 910 if ( data&z(li) .eq. -1 ) then 911 level = "-average" 912 else 913 level = "=" + data&z(li) + "m" 914 end if 915 cs_res@gsnCenterString = "t=" + data&t(lo) + "s x"+ level 916 plot(n) = gsn_csm_contour(wks_ps,data(varn,lo,:,:,li),cs_res) 917 if (vector .EQ. 1) then 918 if (check_vecp) 919 vecres = True ; vector only resources 920 vecres@gsnDraw = False ; don't draw 921 vecres@gsnFrame = False ; don't advance frame 922 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 923 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 924 vecres@vcRefLengthF = 0.05 ; define length of vec ref 925 vecres@gsnRightString = " " ; turn off right string 926 vecres@gsnLeftString = " " ; turn off left string 927 vecres@tiXAxisString = " " 928 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 929 overlay(plot(n), plot_vec) 930 else 931 vecres = True ; vector only resources 932 vecres@gsnDraw = False ; don't draw 933 vecres@gsnFrame = False ; don't advance frame 934 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 935 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 936 vecres@vcRefLengthF = 0.05 ; define length of vec ref 937 vecres@gsnRightString = "Vector plot of "+vec1+" and "+vec2 938 vecres@gsnLeftString = " " ; turn off left string 939 vecres@tiXAxisString = " " 940 plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres) 941 plot(n) = plot_vec 909 942 end if 910 cs_res@gsnCenterString = "t=" + data&t(lo-1) + "s x"+ level 911 plot1 = gsn_csm_contour(wks_ps,data(varn,lo-1,:,:,li+n),cs_res) 912 if (check_vecp) 913 vecres = True ; vector only resources 914 vecres@gsnDraw = False ; don't draw 915 vecres@gsnFrame = False ; don't advance frame 916 vecres@vcGlyphStyle = "CurlyVector" ; curly vectors 917 vecres@vcRefMagnitudeF = ref_mag ; define vector ref mag 918 vecres@vcRefLengthF = 0.05 ; define length of vec ref 919 vecres@gsnRightString = " " ; turn off right string 920 vecres@gsnLeftString = " " ; turn off left string 921 vecres@tiXAxisString = " " 922 plot_vec=gsn_csm_vector(wks_ps,vect1(li+n-1,lo,:,:),vect2(li+n-1,lo,:,:),vecres) 923 overlay(plot1, plot_vec) 924 end if 925 plot(n) = plot1 926 end do 943 end if 927 944 end if 928 929 end if 930 931 ; *************************************************** 932 ; merge plots onto one page 933 ; *************************************************** 934 935 gsn_panel(wks_ps, plot(0:n-1),(/no_lines,no_columns/),cs_resP) 936 937 ; To hold no_lines and no_columns constant, it must be defined again: 938 939 no_lines = nl1 940 no_columns = nc1 941 945 end if 946 n=n+1 942 947 end do 943 948 end do 944 949 end if 945 end do 950 end do 951 952 ; *************************************************** 953 ; merge plots onto one page 954 ; *************************************************** 955 956 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 957 gsn_panel(wks_ps,plot(0:(no_time*no_layer*no_var)-1),(/no_var,no_layer*no_time/),cs_resP) 958 print(" ") 959 print("Outputs to .eps or .epsi have only one frame") 960 print(" ") 961 else 962 do np = 0,no_layer*no_time*no_var-1,no_lines*no_columns 963 if ( np + no_lines*no_columns .gt. (no_layer*no_time*no_var)-1) then 964 gsn_panel(wks_ps, plot(np:(no_layer*no_time*no_var)-1),(/no_lines,no_columns/),cs_resP) 965 else 966 gsn_panel(wks_ps, plot(np:np+no_lines*no_columns-1),(/no_lines,no_columns/),cs_resP) 967 end if 968 end do 969 end if 946 970 947 971 print(" ") -
palm/trunk/SCRIPTS/NCL/profiles.ncl
r157 r161 10 10 11 11 if (isfilepresent("~/.ncl_preferences")) then 12 parameter = asciiread("~/.ncl_preferences",7 3,"string")12 parameter = asciiread("~/.ncl_preferences",75,"string") 13 13 delete(parameter@_FillValue) 14 14 else … … 71 71 if (combine .EQ. 1) then 72 72 if( .not. isvar("c_var") ) then 73 c_var=" " 74 if (parameter(27) .NE. "c_variables") then 73 if (parameter(27) .EQ. "c_variables") then 74 print(" ") 75 print("Please select variables for overlaying ('c_var')") 76 print(" ") 77 exit 78 else 75 79 c_var=parameter(27) 76 80 end if … … 101 105 end if 102 106 end if 103 107 if ( .not. isvar("over") ) then ; switches overlaying plots on 108 over = 0 109 if (stringtointeger(parameter(73)) .NE. 0) then 110 over = stringtointeger(parameter(73)) 111 if (stringtointeger(parameter(73)) .NE. 1) then 112 print(" ") 113 print("Please set 'over' to 0 or 1") 114 print(" ") 115 exit 116 end if 117 end if 118 end if 119 104 120 ; *************************************************** 105 121 ; open input file … … 108 124 f=addfile( file_in,"r") 109 125 110 vNam 126 vNam = getfilevarnames(f) 111 127 print(" ") 112 128 print("Variable on netCDF file: " + vNam) 113 129 print(" ") 114 dim = dimsizes(vNam) 115 z_pr = f->zpt 116 dimz = dimsizes(z_pr) 130 dim = dimsizes(vNam) 131 if (dim .EQ. 0) then 132 print(" ") 133 print("There are no data on file") 134 print(" ") 135 end if 136 137 do varn = dim-1,0,1 138 if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then 139 varn=varn-1 140 continue 141 end if 142 if (vNam(varn) .EQ. "u") 143 z_u = f->zu 144 break 145 else 146 if (vNam(varn) .EQ. "v")then 147 z_u = f->zv 148 break 149 else 150 if(vNam(varn) .EQ. "pt")then 151 z_u = f->zpt 152 break 153 else 154 if(vNam(varn) .EQ. "vpt")then 155 z_u = f->zvpt 156 break 157 else 158 if(vNam(varn) .EQ. "lpt")then 159 z_u = f->zlpt 160 break 161 else 162 if(vNam(varn) .EQ. "q")then 163 z_u = f->zq 164 break 165 else 166 if(vNam(varn) .EQ. "qv")then 167 z_u = f->zqv 168 break 169 else 170 if(vNam(varn) .EQ. "ql")then 171 z_u = f->zql 172 break 173 else 174 if(vNam(varn) .EQ. "rho")then 175 z_u = f->zrho 176 break 177 else 178 if(vNam(varn) .EQ. "s")then 179 z_u = f->zs 180 break 181 else 182 if(vNam(varn) .EQ. "sa")then 183 z_u = f->zsa 184 break 185 else 186 if(vNam(varn) .EQ. "e")then 187 z_u = f->ze 188 break 189 else 190 if(vNam(varn) .EQ. "es")then 191 z_u = f->zes 192 break 193 else 194 if(vNam(varn) .EQ. "km")then 195 z_u = f->zkm 196 break 197 else 198 if(vNam(varn) .EQ. "kh")then 199 z_u = f->zkh 200 break 201 else 202 if(vNam(varn) .EQ. "l")then 203 z_u = f->zl 204 break 205 else 206 if(vNam(varn) .EQ. "us2")then 207 z_u = f->zus2 208 break 209 else 210 if(vNam(varn) .EQ. "vs2")then 211 z_u = f->zvs2 212 break 213 else 214 if(vNam(varn) .EQ. "pts2")then 215 z_u = f->zpts2 216 break 217 else 218 if(vNam(varn) .EQ. "wsususodz")then 219 z_u = f->zwsususodz 220 break 221 else 222 if(vNam(varn) .EQ. "wspsodz")then 223 z_u = f->zwspsodz 224 break 225 else 226 if(vNam(varn) .EQ. "wpeodz")then 227 z_u = f->zwpeodz 228 break 229 end if 230 end if 231 end if 232 end if 233 end if 234 end if 235 end if 236 end if 237 end if 238 end if 239 end if 240 end if 241 end if 242 end if 243 end if 244 end if 245 end if 246 end if 247 end if 248 end if 249 end if 250 end if 251 varn=varn-1 252 end do 253 254 do varn=dim-1,0,1 255 if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then 256 varn=varn-1 257 continue 258 end if 259 if (vNam(varn) .EQ. "w") 260 z_w = f->zw 261 break 262 else 263 if (vNam(varn) .EQ. "wpup")then 264 z_w = f->zwpup 265 break 266 else 267 if(vNam(varn) .EQ. "wsus")then 268 z_w = f->zwsus 269 break 270 else 271 if(vNam(varn) .EQ. "wu")then 272 z_w = f->zwu 273 break 274 else 275 if(vNam(varn) .EQ. "wpvp")then 276 z_w = f->zwpvp 277 break 278 else 279 if(vNam(varn) .EQ. "wsvs")then 280 z_w = f->zwsvs 281 break 282 else 283 if(vNam(varn) .EQ. "wv")then 284 z_w = f->zwv 285 break 286 else 287 if(vNam(varn) .EQ. "wptpp")then 288 z_w = f->zwptpp 289 break 290 else 291 if(vNam(varn) .EQ. "wspts")then 292 z_w = f->zwspts 293 break 294 else 295 if(vNam(varn) .EQ. "wpt")then 296 z_w = f->zwpt 297 break 298 else 299 if(vNam(varn) .EQ. "wsptsBC")then 300 z_w = f->zwsptsBC 301 break 302 else 303 if(vNam(varn) .EQ. "wptBC")then 304 z_w = f->zwptBC 305 break 306 else 307 if(vNam(varn) .EQ. "wpvptp")then 308 z_w = f->zwpvptp 309 break 310 else 311 if(vNam(varn) .EQ. "wsvpts")then 312 z_w = f->zwsvpts 313 break 314 else 315 if(vNam(varn) .EQ. "wvpt")then 316 z_w = f->zwvpt 317 break 318 else 319 if(vNam(varn) .EQ. "wpqp")then 320 z_w = f->zwpqp 321 break 322 else 323 if(vNam(varn) .EQ. "wsqs")then 324 z_w = f->zwsqs 325 break 326 else 327 if(vNam(varn) .EQ. "wq")then 328 z_w = f->zwq 329 break 330 else 331 if(vNam(varn) .EQ. "wpqvp")then 332 z_w = f->zwpqvp 333 break 334 else 335 if(vNam(varn) .EQ. "wsqvs")then 336 z_w = f->zwsqvs 337 break 338 else 339 if(vNam(varn) .EQ. "wqv")then 340 z_w = f->zwqv 341 break 342 else 343 if(vNam(varn) .EQ. "wpsp")then 344 z_w = f->zwpsp 345 break 346 else 347 if(vNam(varn) .EQ. "wsss")then 348 z_w = f->zwsss 349 break 350 else 351 if(vNam(varn) .EQ. "ws")then 352 z_w = f->zws 353 break 354 else 355 if(vNam(varn) .EQ. "wpsap")then 356 z_w = f->zwpsap 357 break 358 else 359 if(vNam(varn) .EQ. "wssas")then 360 z_w = f->zwssas 361 break 362 else 363 if(vNam(varn) .EQ. "wsa")then 364 z_w = f->zwsa 365 break 366 else 367 if(vNam(varn) .EQ. "wses")then 368 z_w = f->zwses 369 break 370 else 371 if(vNam(varn) .EQ. "ws2")then 372 z_w = f->zws2 373 break 374 else 375 if(vNam(varn) .EQ. "ws3")then 376 z_w = f->zws3 377 break 378 else 379 if(vNam(varn) .EQ. "Sw")then 380 z_w = f->zSw 381 break 382 else 383 if(vNam(varn) .EQ. "ws2pts")then 384 z_w = f->zws2pts 385 break 386 else 387 if(vNam(varn) .EQ. "wspts2")then 388 z_w = f->zwspts2 389 break 390 end if 391 end if 392 end if 393 end if 394 end if 395 end if 396 end if 397 end if 398 end if 399 end if 400 end if 401 end if 402 end if 403 end if 404 end if 405 end if 406 end if 407 end if 408 end if 409 end if 410 end if 411 end if 412 end if 413 end if 414 end if 415 end if 416 end if 417 end if 418 end if 419 end if 420 end if 421 end if 422 end if 423 varn=varn-1 424 end do 425 if ( .not. isvar("z_u") ) then 426 zu = 0 427 else 428 zu = 1 429 end if 430 if ( .not. isvar("z_w") ) then 431 zw = 0 432 else 433 zw = 1 434 end if 435 436 if (zu .EQ. 0 .AND. zw .EQ. 0) then 437 do varn=0,dim-1 438 if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then 439 check = False 440 else 441 if (.not. isvar("var")) then 442 check = True 443 if (parameter(21) .NE. "variables") then 444 var=parameter(21) 445 check = isStrSubset( var,","+vNam(varn)+"," ) 446 end if 447 else 448 check = isStrSubset( var,","+vNam(varn)+"," ) 449 end if 450 end if 451 if (check) 452 z = f->$vNam(varn+1)$ 453 dimz = dimsizes(z) 454 break 455 end if 456 end do 457 end if 458 459 dimz = dimsizes(z_w) 460 dimz = dimsizes(z_u) 117 461 t_all = f->time 118 462 nt = dimsizes(t_all) 119 463 120 464 ; **************************************************** 121 465 ; start of time step and different types of mistakes that could be done … … 125 469 start_time_step = 1 126 470 if (parameter(13) .NE. "1") then 127 if (parameter(13) .LE. " 0")128 print(" ") 129 print("Begin with time step 1")471 if (parameter(13) .LE. "1") 472 print(" ") 473 print("Begin at least with time step 2") 130 474 print(" ") 131 475 exit … … 137 481 exit 138 482 end if 139 start_time_step = stringtointeger(parameter(13)) 483 start_time_step = stringtointeger(parameter(13))-1 140 484 end if 141 485 else 142 486 if (start_time_step .LE. 0) 143 487 print(" ") 144 print("Begin with time step 1")488 print("Begin at least with time step 1") 145 489 print(" ") 146 490 exit … … 152 496 exit 153 497 end if 498 start_time_step = start_time_step - 1 154 499 end if 155 500 … … 200 545 exit 201 546 end if 547 end_time_step = end_time_step - 1 202 548 end if 203 549 … … 206 552 ; **************************************************** 207 553 208 legend_label=new(nt -1,double)554 legend_label=new(nt,double) 209 555 do p=start_time_step,end_time_step 210 556 legend_label(p-start_time_step)=t_all(p) … … 224 570 res@gsnFrame = False 225 571 res@gsnPaperOrientation = "portrait" 226 ; res@gsnmaximize = True227 572 res@gsnPaperWidth = 8.27 228 573 res@gsnPaperHeight = 11.69 … … 239 584 res@pmLegendSide = "Top" 240 585 res@xyExplicitLegendLabels = legend_label 241 res@pmLegendParallelPosF = 1. 4242 res@pmLegendOrthogonalPosF = -1.0 85586 res@pmLegendParallelPosF = 1.15 587 res@pmLegendOrthogonalPosF = -1.0 243 588 res@pmLegendWidthF = 0.12 244 589 res@pmLegendHeightF = 0.3 245 590 res@lgLabelFontHeightF = .02 246 ; res@XBLabelConstantSpacingF = 1.0 247 591 res@txFontHeightF = 0.02 592 res@tiXAxisFontHeightF = 0.02 593 res@tiYAxisFontHeightF = 0.02 594 res@tiXAxisString = " " 595 res@tiYAxisString = "Height [z]" 596 248 597 if ( dash .eq. 0 ) then 249 598 res@xyMonoDashPattern = True … … 257 606 258 607 ; *************************************************** 259 ; set up colors and recourses for combined plot260 ; ***************************************************261 262 if (combine .EQ. 1) then263 264 if (.not. isvar("number_comb")) then265 m=0266 m=stringtointeger(parameter(25))267 if(m .EQ. 0) then268 print(" ")269 print("Please indicate the number of variables you would like to combine ('number_comb')")270 print(" ")271 exit272 end if273 end if274 275 ores=True276 277 colors=new(m*(nt-1),integer)278 do j=0,m-1279 colors(j+j*(nt-2):(j+1)*(nt-1)-1) = ispan(2,237,235/np)280 end do281 ores@xyLineColors = colors282 283 dash_oplot=new(m*(nt-1),integer)284 do j=0,m-1285 dash_oplot(j+j*(nt-2):(j+1)*(nt-1)-1)=j286 end do287 ores@xyDashPatterns = dash_oplot288 289 ores = True290 ores@gsnDraw = False291 ores@gsnFrame = False292 ores@gsnPaperOrientation = "portrait"293 ; ores@gsnmaximize = True294 ores@gsnPaperWidth = 8.27295 ores@gsnPaperHeight = 11.69296 ores@gsnPaperMargin = 0.79297 ores@txFont = "helvetica"298 ores@tiMainFont = "helvetica"299 ores@tiXAxisFont = "helvetica"300 ores@tiYAxisFont = "helvetica"301 ores@tmXBLabelFont = "helvetica"302 ores@tmYLLabelFont = "helvetica"303 ores@lgLabelFont = "helvetica"304 ores@tmLabelAutoStride = True305 ores@pmLegendDisplayMode = "Always"306 ores@pmLegendSide = "Top"307 ores@pmLegendParallelPosF = 1.4308 ores@pmLegendOrthogonalPosF = -1.085309 ores@pmLegendWidthF = 0.15310 ores@pmLegendHeightF = 0.60311 ores@lgLabelFontHeightF = .02312 ; ores@XBLabelConstantSpacingF = 1.0313 314 end if315 316 ; ***************************************************317 608 ; set up graphics for plot 318 609 ; *************************************************** 319 610 320 plot = new(dim,graphic) 611 plot = new(dim,graphic) 612 plot_ = new(dim,graphic) 613 614 if (combine .EQ. 1) then 615 if ( .not. isvar("number_comb") ) then 616 if (parameter(25) .EQ. "0") then 617 print(" ") 618 print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot") 619 print(" ") 620 exit 621 else 622 number_comb=stringtointeger(parameter(25)) 623 plot_o = new(number_comb,graphic) 624 end if 625 else 626 if(number_comb .EQ. 2 .OR. number_comb .EQ. 3) then 627 plot_o = new(number_comb,graphic) 628 else 629 print(" ") 630 print("Please set 'number_comb' to 2 or 3 if you would like to overlay 2 or 3 variables in one plot") 631 print(" ") 632 exit 633 end if 634 end if 635 label=new(number_comb,string) 636 color_o=new(number_comb,integer) 637 mini=new(number_comb,float) 638 maxi=new(number_comb,float) 639 end if 321 640 322 641 wks=gsn_open_wks(format_out,file_out) … … 324 643 325 644 ; *************************************************** 326 ; indicate plot number 645 ; indicate plot number 327 646 ; *************************************************** 328 647 329 648 if (combine .EQ. 1) then 330 n = 1 649 n = 1 331 650 else 332 651 n = 0 333 end if 334 652 end if 653 335 654 ; *************************************************** 336 655 ; set up minimum and maximum height … … 340 659 min_z=0 341 660 if (stringtointeger(parameter(33)) .NE. 0) then 342 if (stringtointeger(parameter(33)) .GE. max(z_ pr) ) then343 print(" ") 344 print("Minimum of height ('min_z'="+stringtointeger(parameter(33))+") is greater than available heights (="+max(z_ pr)+")")661 if (stringtointeger(parameter(33)) .GE. max(z_u) ) then 662 print(" ") 663 print("Minimum of height ('min_z'="+stringtointeger(parameter(33))+") is greater than available heights (="+max(z_u)+")") 345 664 print(" ") 346 665 exit … … 355 674 end if 356 675 else 357 if (min_z .GE. max(z_ pr) ) then676 if (min_z .GE. max(z_u) ) then 358 677 print(" ") 359 print("Minimum of height ('min_z'="+min_z+") is greater than available heights (="+max(z_ pr)+")")678 print("Minimum of height ('min_z'="+min_z+") is greater than available heights (="+max(z_u)+")") 360 679 print(" ") 361 680 exit … … 370 689 371 690 if (.not. isvar("max_z")) 372 max_z=max(z_ pr)373 if ((parameter(35)) .NE. "max(z_ pr)") then374 if (stringtofloat(parameter(35)) .GE. max(z_ pr) ) then375 print(" ") 376 print("Maximum of height ('max_z'="+parameter(35)+") is greater than available heights (="+max(z_ pr)+")")691 max_z=max(z_u) 692 if ((parameter(35)) .NE. "max(z_u)") then 693 if (stringtofloat(parameter(35)) .GE. max(z_u) ) then 694 print(" ") 695 print("Maximum of height ('max_z'="+parameter(35)+") is greater than available heights (="+max(z_u)+")") 377 696 print(" ") 378 697 exit … … 389 708 if (max_z .GE. max(z_pr) ) then 390 709 print(" ") 391 print("Maximum of height ('max_z'="+max_z+") is greater than available heights (="+max(z_ pr)+")")710 print("Maximum of height ('max_z'="+max_z+") is greater than available heights (="+max(z_u)+")") 392 711 print(" ") 393 712 exit … … 411 730 end if 412 731 end do 413 414 if (combine .EQ. 1) then 415 data_o=new((/m,end_time_step,dimz/),float) 416 data_o!0 = "e" 417 legend_label_oplot=new(m*(nt-1),string) 418 mini=new(m,float) 419 maxi=new(m,float) 420 e=-1 421 end if 422 423 do varn = 0, dim-1 424 425 if ( isStrSubset( vNam(varn), "NORM") .or. isStrSubset( vNam(varn), "time") ) 732 733 data = new((/dim,(end_time_step-start_time_step)+1,dimz/),float) 734 data_0 = new((/(end_time_step-start_time_step)+1,dimz/),float) 735 data_0 = 0.0 736 t = new((/(end_time_step-start_time_step)+1,dimz/),float) 737 t = 0.0 738 unit = new(dim,string) 739 740 if (over .EQ. 1) then 741 plot_u = gsn_csm_xy(wks,t,data_0(:,:),res) 742 miniu = 1.E27 743 maxiu =-1.E27 744 plot_v = gsn_csm_xy(wks,t,data_0(:,:),res) 745 miniv = 1.E27 746 maxiv =-1.E27 747 plot_w = gsn_csm_xy(wks,t,data_0(:,:),res) 748 miniw = 1.E27 749 maxiw =-1.E27 750 plot_pt = gsn_csm_xy(wks,t,data_0(:,:),res) 751 minipt = 1.E27 752 maxipt =-1.E27 753 plot_vpt = gsn_csm_xy(wks,t,data_0(:,:),res) 754 minivpt = 1.E27 755 maxivpt =-1.E27 756 plot_lpt = gsn_csm_xy(wks,t,data_0(:,:),res) 757 minilpt = 1.E27 758 maxilpt =-1.E27 759 plot_q = gsn_csm_xy(wks,t,data_0(:,:),res) 760 miniq = 1.E27 761 maxiq =-1.E27 762 plot_qv = gsn_csm_xy(wks,t,data_0(:,:),res) 763 miniqv = 1.E27 764 maxiqv =-1.E27 765 plot_ql = gsn_csm_xy(wks,t,data_0(:,:),res) 766 miniql = 1.E27 767 maxiql =-1.E27 768 plot_rho = gsn_csm_xy(wks,t,data_0(:,:),res) 769 plot_s = gsn_csm_xy(wks,t,data_0(:,:),res) 770 plot_sa = gsn_csm_xy(wks,t,data_0(:,:),res) 771 plot_e = gsn_csm_xy(wks,t,data_0(:,:),res) 772 minie = 1.E27 773 maxie =-1.E27 774 plot_es = gsn_csm_xy(wks,t,data_0(:,:),res) 775 minie = 1.E27 776 maxie =-1.E27 777 plot_km = gsn_csm_xy(wks,t,data_0(:,:),res) 778 minie = 1.E27 779 maxie =-1.E27 780 plot_kh = gsn_csm_xy(wks,t,data_0(:,:),res) 781 minie = 1.E27 782 maxie =-1.E27 783 plot_l = gsn_csm_xy(wks,t,data_0(:,:),res) 784 plot_wpup = gsn_csm_xy(wks,t,data_0(:,:),res) 785 miniwpup = 1.E27 786 maxiwpup =-1.E27 787 plot_wsus = gsn_csm_xy(wks,t,data_0(:,:),res) 788 miniwsus = 1.E27 789 maxiwsus =-1.E27 790 plot_wu = gsn_csm_xy(wks,t,data_0(:,:),res) 791 miniwu = 1.E27 792 maxiwu =-1.E27 793 plot_wpvp = gsn_csm_xy(wks,t,data_0(:,:),res) 794 miniwpvp = 1.E27 795 maxiwpvp =-1.E27 796 plot_wsvs = gsn_csm_xy(wks,t,data_0(:,:),res) 797 miniwsvs = 1.E27 798 maxiwsvs =-1.E27 799 plot_wv = gsn_csm_xy(wks,t,data_0(:,:),res) 800 miniwv = 1.E27 801 maxiwv =-1.E27 802 plot_wpptp = gsn_csm_xy(wks,t,data_0(:,:),res) 803 miniwpptp = 1.E27 804 maxiwpptp =-1.E27 805 plot_wspts = gsn_csm_xy(wks,t,data_0(:,:),res) 806 miniwspts = 1.E27 807 maxiwspts =-1.E27 808 plot_wpt = gsn_csm_xy(wks,t,data_0(:,:),res) 809 miniwpt = 1.E27 810 maxiwpt =-1.E27 811 plot_wsptsBC = gsn_csm_xy(wks,t,data_0(:,:),res) 812 miniwsptsBC = 1.E27 813 maxiwsptsBC =-1.E27 814 plot_wptBC = gsn_csm_xy(wks,t,data_0(:,:),res) 815 miniwptBC = 1.E27 816 maxiwptBC =-1.E27 817 plot_wpvptp = gsn_csm_xy(wks,t,data_0(:,:),res) 818 miniwpvptp = 1.E27 819 maxiwpvptp =-1.E27 820 plot_wsvpts = gsn_csm_xy(wks,t,data_0(:,:),res) 821 miniwsvpts = 1.E27 822 maxiewsvpts=-1.E27 823 plot_wvpt = gsn_csm_xy(wks,t,data_0(:,:),res) 824 miniwvpt = 1.E27 825 maxiwvpt =-1.E27 826 plot_wpqp = gsn_csm_xy(wks,t,data_0(:,:),res) 827 miniwpqp = 1.E27 828 maxiwpqp =-1.E27 829 plot_wsqs = gsn_csm_xy(wks,t,data_0(:,:),res) 830 miniwsqs = 1.E27 831 maxiwsqs =-1.E27 832 plot_wq = gsn_csm_xy(wks,t,data_0(:,:),res) 833 miniwq = 1.E27 834 maxiwq =-1.E27 835 plot_wpqvp = gsn_csm_xy(wks,t,data_0(:,:),res) 836 miniwpqvp = 1.E27 837 maxiwpqvp =-1.E27 838 plot_wsqvs = gsn_csm_xy(wks,t,data_0(:,:),res) 839 miniwsqvs = 1.E27 840 maxiwsqvs =-1.E27 841 plot_wqv = gsn_csm_xy(wks,t,data_0(:,:),res) 842 miniwqv = 1.E27 843 maxiwqv =-1.E27 844 plot_wpsp = gsn_csm_xy(wks,t,data_0(:,:),res) 845 miniwpsp = 1.E27 846 maxiwpsp =-1.E27 847 plot_wsss = gsn_csm_xy(wks,t,data_0(:,:),res) 848 miniwsss = 1.E27 849 maxiwsss =-1.E27 850 plot_ws = gsn_csm_xy(wks,t,data_0(:,:),res) 851 miniws = 1.E27 852 maxiws =-1.E27 853 plot_wpsap = gsn_csm_xy(wks,t,data_0(:,:),res) 854 miniwpsap = 1.E27 855 maxiwpsap =-1.E27 856 plot_wssas = gsn_csm_xy(wks,t,data_0(:,:),res) 857 miniwssas = 1.E27 858 maxiwssas =-1.E27 859 plot_wsa = gsn_csm_xy(wks,t,data_0(:,:),res) 860 miniwsa = 1.E27 861 maxiwsa =-1.E27 862 plot_wses = gsn_csm_xy(wks,t,data_0(:,:),res) 863 plot_us2 = gsn_csm_xy(wks,t,data_0(:,:),res) 864 minius2 = 1.E27 865 maxius2 =-1.E27 866 plot_vs2 = gsn_csm_xy(wks,t,data_0(:,:),res) 867 minivs2 = 1.E27 868 maxivs2 =-1.E27 869 plot_ws2 = gsn_csm_xy(wks,t,data_0(:,:),res) 870 miniws2 = 1.E27 871 maxiws2 =-1.E27 872 plot_pts2 = gsn_csm_xy(wks,t,data_0(:,:),res) 873 plot_ws3 = gsn_csm_xy(wks,t,data_0(:,:),res) 874 plot_Sw = gsn_csm_xy(wks,t,data_0(:,:),res) 875 plot_ws2pts = gsn_csm_xy(wks,t,data_0(:,:),res) 876 plot_wspts2 = gsn_csm_xy(wks,t,data_0(:,:),res) 877 plot_wsususodz = gsn_csm_xy(wks,t,data_0(:,:),res) 878 miniwsususodz = 1.E27 879 maxiwsususodz =-1.E27 880 plot_wspsodz = gsn_csm_xy(wks,t,data_0(:,:),res) 881 miniwspsodz = 1.E27 882 maxiwspsodz =-1.E27 883 plot_wpeodz = gsn_csm_xy(wks,t,data_0(:,:),res) 884 miniwpeodz = 1.E27 885 maxiwpeodz =-1.E27 886 end if 887 888 n_o=0 889 do varn = 0,dim-1 890 891 temp = f->$vNam(varn)$ 892 893 if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then 426 894 check = False 427 895 else … … 430 898 if (parameter(21) .NE. "variables") then 431 899 var=parameter(21) 432 check = isStrSubset( var, vNam(varn)+"," )900 check = isStrSubset( var,","+vNam(varn)+"," ) 433 901 end if 434 902 else 435 check = isStrSubset( var, vNam(varn)+"," )903 check = isStrSubset( var,","+vNam(varn)+"," ) 436 904 end if 437 905 end if 906 907 if (combine .EQ. 1) then 908 com=isStrSubset(c_var,","+vNam(varn)+"," ) 909 if (com) 910 check = False 911 data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1) 912 unit(varn) = temp@units 913 if (n_o .GT. number_comb-1) then 914 print(" ") 915 print("Please set 'number_comb' to the right number of overlaying variables ('c_var')") 916 print(" ") 917 exit 918 end if 919 mini(n_o)=min(data(varn,:,:)) 920 maxi(n_o)=max(data(varn,:,:)) 921 varn=varn+1 922 n_o=n_o+1 923 end if 924 end if 925 926 if(check) then 927 928 z = f->$vNam(varn+1)$ 929 930 unit(varn) = temp@units 931 data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1) 932 933 if (over .EQ. 0) then 934 res@gsnLeftString = vNam(varn) 935 res@gsnRightString = unit(varn) 936 res@trYMinF = min_z 937 res@trYMaxF = max_z 938 res@trXMinF = min(data(varn,:,:)) 939 res@trXMaxF = max(data(varn,:,:)) 940 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 941 end if 942 943 if (vNam(varn) .EQ. "u") then 944 miniu=min(data(varn,:,:)) 945 maxiu=max(data(varn,:,:)) 946 if (over .EQ. 1) then 947 res@xyDashPattern = 0 948 plot_u = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 949 else 950 res@gsnLeftString = vNam(varn) 951 res@gsnRightString = unit(varn) 952 res@trXMinF = miniu 953 res@trXMaxF = maxiu 954 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 955 end if 956 end if 957 if (vNam(varn) .EQ. "v") then 958 miniv=min(data(varn,:,:)) 959 maxiv=max(data(varn,:,:)) 960 if (over .EQ. 1) then 961 res@xyDashPattern = 1 962 plot_v = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 963 else 964 res@gsnLeftString = vNam(varn) 965 res@gsnRightString = unit(varn) 966 res@trXMinF = miniv 967 res@trXMaxF = maxiv 968 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 969 end if 970 end if 971 if (vNam(varn) .EQ. "w") then 972 miniw=min(data(varn,:,:)) 973 maxiw=max(data(varn,:,:)) 974 if (over .EQ. 1) then 975 res@xyDashPattern = 2 976 plot_w = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 977 else 978 res@gsnLeftString = vNam(varn) 979 res@gsnRightString = unit(varn) 980 res@trXMinF = miniw 981 res@trXMaxF = maxiw 982 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 983 end if 984 end if 985 986 if (vNam(varn) .EQ. "pt") then 987 minipt=min(data(varn,:,:)) 988 maxipt=max(data(varn,:,:)) 989 if (over .EQ. 1) then 990 res@xyDashPattern = 0 991 plot_pt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 992 else 993 res@gsnLeftString = vNam(varn) 994 res@gsnRightString = unit(varn) 995 res@trXMinF = minipt 996 res@trXMaxF = maxipt 997 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 998 end if 999 end if 1000 if (vNam(varn) .EQ. "vpt") then 1001 minivpt=min(data(varn,:,:)) 1002 maxivpt=max(data(varn,:,:)) 1003 if (over .EQ. 1) then 1004 res@xyDashPattern = 1 1005 plot_vpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1006 else 1007 res@gsnLeftString = vNam(varn) 1008 res@gsnRightString = unit(varn) 1009 res@trXMinF = minivpt 1010 res@trXMaxF = maxivpt 1011 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1012 end if 1013 end if 1014 if (vNam(varn) .EQ. "lpt") then 1015 minilpt=min(data(varn,:,:)) 1016 maxilpt=max(data(varn,:,:)) 1017 if (over .EQ. 1) then 1018 res@xyDashPattern = 2 1019 plot_lpt = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1020 else 1021 res@gsnLeftString = vNam(varn) 1022 res@gsnRightString = unit(varn) 1023 res@trXMinF = minilpt 1024 res@trXMaxF = maxilpt 1025 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1026 end if 1027 end if 1028 1029 if (vNam(varn) .EQ. "q") then 1030 miniq=min(data(varn,:,:)) 1031 maxiq=max(data(varn,:,:)) 1032 if (over .EQ. 1) then 1033 res@xyDashPattern = 0 1034 plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1035 else 1036 res@gsnLeftString = vNam(varn) 1037 res@gsnRightString = unit(varn) 1038 res@trXMinF = miniq 1039 res@trXMaxF = maxiq 1040 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1041 end if 1042 end if 1043 if (vNam(varn) .EQ. "qv") then 1044 miniqv=min(data(varn,:,:)) 1045 maxiqv=max(data(varn,:,:)) 1046 if (over .EQ. 1) then 1047 res@xyDashPattern = 1 1048 plot_qv = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1049 else 1050 res@gsnLeftString = vNam(varn) 1051 res@gsnRightString = unit(varn) 1052 res@trXMinF = miniqv 1053 res@trXMaxF = maxiqv 1054 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1055 end if 1056 end if 1057 if (vNam(varn) .EQ. "ql") then 1058 miniql=min(data(varn,:,:)) 1059 maxiql=max(data(varn,:,:)) 1060 if (over .EQ. 1) then 1061 res@xyDashPattern = 2 1062 plot_ql = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1063 else 1064 res@gsnLeftString = vNam(varn) 1065 res@gsnRightString = unit(varn) 1066 res@trXMinF = miniql 1067 res@trXMaxF = maxiql 1068 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1069 end if 1070 end if 1071 1072 if (vNam(varn) .EQ. "rho") then 1073 if (over .EQ. 1) then 1074 res@gsnLeftString = vNam(varn) 1075 res@gsnRightString = unit(varn) 1076 res@trXMinF = min(data(varn,:,:)) 1077 res@trXMaxF = max(data(varn,:,:)) 1078 plot_q = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1079 else 1080 res@gsnLeftString = vNam(varn) 1081 res@gsnRightString = unit(varn) 1082 res@trXMinF = min(data(varn,:,:)) 1083 res@trXMaxF = max(data(varn,:,:)) 1084 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1085 end if 1086 end if 1087 1088 if (vNam(varn) .EQ. "s") then 1089 if (over .EQ. 1) then 1090 res@gsnLeftString = vNam(varn) 1091 res@gsnRightString = unit(varn) 1092 res@trXMinF = min(data(varn,:,:)) 1093 res@trXMaxF = max(data(varn,:,:)) 1094 plot_s = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1095 else 1096 res@gsnLeftString = vNam(varn) 1097 res@gsnRightString = unit(varn) 1098 res@trXMinF = min(data(varn,:,:)) 1099 res@trXMaxF = max(data(varn,:,:)) 1100 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1101 end if 1102 end if 1103 1104 if (vNam(varn) .EQ. "sa") then 1105 if (over .EQ. 1) then 1106 res@gsnLeftString = vNam(varn) 1107 res@gsnRightString = unit(varn) 1108 res@trXMinF = min(data(varn,:,:)) 1109 res@trXMaxF = max(data(varn,:,:)) 1110 plot_sa = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1111 else 1112 res@gsnLeftString = vNam(varn) 1113 res@gsnRightString = unit(varn) 1114 res@trXMinF = min(data(varn,:,:)) 1115 res@trXMaxF = max(data(varn,:,:)) 1116 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1117 end if 1118 end if 1119 1120 if (vNam(varn) .EQ. "e") then 1121 minie=min(data(varn,:,:)) 1122 maxie=max(data(varn,:,:)) 1123 if (over .EQ. 1) then 1124 res@xyDashPattern = 0 1125 plot_e = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1126 else 1127 res@gsnLeftString = vNam(varn) 1128 res@gsnRightString = unit(varn) 1129 res@trXMinF = minie 1130 res@trXMaxF = maxie 1131 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1132 end if 1133 end if 1134 if (vNam(varn) .EQ. "es") then 1135 minies=min(data(varn,:,:)) 1136 maxies=max(data(varn,:,:)) 1137 if (over .EQ. 1) then 1138 res@xyDashPattern = 1 1139 plot_es = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1140 else 1141 res@gsnLeftString = vNam(varn) 1142 res@gsnRightString = unit(varn) 1143 res@trXMinF = minies 1144 res@trXMaxF = maxies 1145 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1146 end if 1147 end if 1148 1149 if (vNam(varn) .EQ. "km") then 1150 minikm=min(data(varn,:,:)) 1151 maxikm=max(data(varn,:,:)) 1152 if (over .EQ. 1) then 1153 res@xyDashPattern = 0 1154 plot_km = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1155 else 1156 res@gsnLeftString = vNam(varn) 1157 res@gsnRightString = unit(varn) 1158 res@trXMinF = minikm 1159 res@trXMaxF = maxikm 1160 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1161 end if 1162 end if 1163 if (vNam(varn) .EQ. "kh") then 1164 minikh=min(data(varn,:,:)) 1165 maxikh=max(data(varn,:,:)) 1166 if (over .EQ. 1) then 1167 res@xyDashPattern = 1 1168 plot_kh = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1169 else 1170 res@gsnLeftString = vNam(varn) 1171 res@gsnRightString = unit(varn) 1172 res@trXMinF = minikh 1173 res@trXMaxF = maxikh 1174 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1175 end if 1176 end if 1177 1178 if (vNam(varn) .EQ. "l") then 1179 if (over .EQ. 1) then 1180 res@gsnLeftString = vNam(varn) 1181 res@gsnRightString = unit(varn) 1182 res@trXMinF = min(data(varn,:,:)) 1183 res@trXMaxF = max(data(varn,:,:)) 1184 plot_l = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1185 else 1186 res@gsnLeftString = vNam(varn) 1187 res@gsnRightString = unit(varn) 1188 res@trXMinF = min(data(varn,:,:)) 1189 res@trXMaxF = max(data(varn,:,:)) 1190 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1191 end if 1192 end if 1193 1194 if (vNam(varn) .EQ. "wpup") then 1195 miniwpup=min(data(varn,:,:)) 1196 maxiwpup=max(data(varn,:,:)) 1197 if (over .EQ. 1) then 1198 res@xyDashPattern = 0 1199 plot_wpup = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1200 else 1201 res@gsnLeftString = vNam(varn) 1202 res@gsnRightString = unit(varn) 1203 res@trXMinF = miniwpup 1204 res@trXMaxF = maxiwpup 1205 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1206 end if 1207 end if 1208 if (vNam(varn) .EQ. "wsus") then 1209 miniwsus=min(data(varn,:,:)) 1210 maxiwsus=max(data(varn,:,:)) 1211 if (over .EQ. 1) then 1212 res@xyDashPattern = 1 1213 plot_wsus = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1214 else 1215 res@gsnLeftString = vNam(varn) 1216 res@gsnRightString = unit(varn) 1217 res@trXMinF = miniwsus 1218 res@trXMaxF = maxiwsus 1219 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1220 end if 1221 end if 1222 if (vNam(varn) .EQ. "wu") then 1223 miniwu=min(data(varn,:,:)) 1224 maxiwu=max(data(varn,:,:)) 1225 if (over .EQ. 1) then 1226 res@xyDashPattern = 2 1227 plot_wu = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1228 else 1229 res@gsnLeftString = vNam(varn) 1230 res@gsnRightString = unit(varn) 1231 res@trXMinF = miniwu 1232 res@trXMaxF = maxiwu 1233 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1234 end if 1235 end if 1236 1237 if (vNam(varn) .EQ. "wpvp") then 1238 miniwpvp=min(data(varn,:,:)) 1239 maxiwpvp=max(data(varn,:,:)) 1240 if (over .EQ. 1) then 1241 res@xyDashPattern = 0 1242 plot_wpvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1243 else 1244 res@gsnLeftString = vNam(varn) 1245 res@gsnRightString = unit(varn) 1246 res@trXMinF = miniwpvp 1247 res@trXMaxF = maxiwpvp 1248 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1249 end if 1250 end if 1251 if (vNam(varn) .EQ. "wsvs") then 1252 miniwsvs=min(data(varn,:,:)) 1253 maxiwsvs=max(data(varn,:,:)) 1254 if (over .EQ. 1) then 1255 res@xyDashPattern = 1 1256 plot_wsvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1257 else 1258 res@gsnLeftString = vNam(varn) 1259 res@gsnRightString = unit(varn) 1260 res@trXMinF = miniwsvs 1261 res@trXMaxF = maxiwsvs 1262 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1263 end if 1264 end if 1265 if (vNam(varn) .EQ. "wv") then 1266 miniwv=min(data(varn,:,:)) 1267 maxiwv=max(data(varn,:,:)) 1268 if (over .EQ. 1) then 1269 res@xyDashPattern = 2 1270 plot_wv = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1271 else 1272 res@gsnLeftString = vNam(varn) 1273 res@gsnRightString = unit(varn) 1274 res@trXMinF = miniwv 1275 res@trXMaxF = maxiwv 1276 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1277 end if 1278 end if 1279 1280 if (vNam(varn) .EQ. "wpptp") then 1281 miniwpptp=min(data(varn,:,:)) 1282 maxiwpptp=max(data(varn,:,:)) 1283 if (over .EQ. 1) then 1284 res@xyDashPattern = 0 1285 plot_wpptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1286 else 1287 res@gsnLeftString = vNam(varn) 1288 res@gsnRightString = unit(varn) 1289 res@trXMinF = miniwpptp 1290 res@trXMaxF = maxiwpptp 1291 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1292 end if 1293 end if 1294 if (vNam(varn) .EQ. "wspts") then 1295 miniwspts=min(data(varn,:,:)) 1296 maxiwspts=max(data(varn,:,:)) 1297 if (over .EQ. 1) then 1298 res@xyDashPattern = 1 1299 plot_wspts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1300 else 1301 res@gsnLeftString = vNam(varn) 1302 res@gsnRightString = unit(varn) 1303 res@trXMinF = miniwspts 1304 res@trXMaxF = maxiwspts 1305 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1306 end if 1307 end if 1308 if (vNam(varn) .EQ. "wpt") then 1309 miniwpt=min(data(varn,:,:)) 1310 maxiwpt=max(data(varn,:,:)) 1311 if (over .EQ. 1) then 1312 res@xyDashPattern = 2 1313 plot_wpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1314 else 1315 res@gsnLeftString = vNam(varn) 1316 res@gsnRightString = unit(varn) 1317 res@trXMinF = miniwpt 1318 res@trXMaxF = maxiwpt 1319 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1320 end if 1321 end if 1322 1323 if (vNam(varn) .EQ. "wsptsBC") then 1324 miniwsptsBC=min(data(varn,:,:)) 1325 maxiwsptsBC=max(data(varn,:,:)) 1326 if (over .EQ. 1) then 1327 res@xyDashPattern = 0 1328 plot_wsptsBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1329 else 1330 res@gsnLeftString = vNam(varn) 1331 res@gsnRightString = unit(varn) 1332 res@trXMinF = miniwsptsBC 1333 res@trXMaxF = maxiwsptsBC 1334 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1335 end if 1336 end if 1337 if (vNam(varn) .EQ. "wptBC") then 1338 miniwptBC=min(data(varn,:,:)) 1339 maxiwptBC=max(data(varn,:,:)) 1340 if (over .EQ. 1) then 1341 res@xyDashPattern = 1 1342 plot_wptBC = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1343 else 1344 res@gsnLeftString = vNam(varn) 1345 res@gsnRightString = unit(varn) 1346 res@trXMinF = miniwptBC 1347 res@trXMaxF = maxiwptBC 1348 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1349 end if 1350 end if 1351 1352 if (vNam(varn) .EQ. "wpvptp") then 1353 miniwpvptp=min(data(varn,:,:)) 1354 maxiwpvptp=max(data(varn,:,:)) 1355 if (over .EQ. 1) then 1356 res@xyDashPattern = 0 1357 plot_wpvptp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1358 else 1359 res@gsnLeftString = vNam(varn) 1360 res@gsnRightString = unit(varn) 1361 res@trXMinF = miniwpvptp 1362 res@trXMaxF = maxiwpvptp 1363 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1364 end if 1365 end if 1366 if (vNam(varn) .EQ. "wsvpts") then 1367 miniwsvpts=min(data(varn,:,:)) 1368 maxiwsvpts=max(data(varn,:,:)) 1369 if (over .EQ. 1) then 1370 res@xyDashPattern = 1 1371 plot_wsvpts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1372 else 1373 res@gsnLeftString = vNam(varn) 1374 res@gsnRightString = unit(varn) 1375 res@trXMinF = miniwsvpts 1376 res@trXMaxF = maxiwsvpts 1377 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1378 end if 1379 end if 1380 if (vNam(varn) .EQ. "wvpt") then 1381 miniwvpt=min(data(varn,:,:)) 1382 maxiwvpt=max(data(varn,:,:)) 1383 if (over .EQ. 1) then 1384 res@xyDashPattern = 2 1385 plot_wvpt = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1386 else 1387 res@gsnLeftString = vNam(varn) 1388 res@gsnRightString = unit(varn) 1389 res@trXMinF = miniwvpt 1390 res@trXMaxF = maxiwvpt 1391 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1392 end if 1393 end if 1394 1395 if (vNam(varn) .EQ. "wpqp") then 1396 miniwpqp=min(data(varn,:,:)) 1397 maxiwpqp=max(data(varn,:,:)) 1398 if (over .EQ. 1) then 1399 res@xyDashPattern = 0 1400 plot_wpqp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1401 else 1402 res@gsnLeftString = vNam(varn) 1403 res@gsnRightString = unit(varn) 1404 res@trXMinF = miniwpqp 1405 res@trXMaxF = maxiwpqp 1406 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1407 end if 1408 end if 1409 if (vNam(varn) .EQ. "wsqs") then 1410 miniwsqs=min(data(varn,:,:)) 1411 maxiwsqs=max(data(varn,:,:)) 1412 if (over .EQ. 1) then 1413 res@xyDashPattern = 1 1414 plot_wsqs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1415 else 1416 res@gsnLeftString = vNam(varn) 1417 res@gsnRightString = unit(varn) 1418 res@trXMinF = miniwsqs 1419 res@trXMaxF = maxiwsqs 1420 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1421 end if 1422 end if 1423 if (vNam(varn) .EQ. "wq") then 1424 miniwq=min(data(varn,:,:)) 1425 maxiwq=max(data(varn,:,:)) 1426 if (over .EQ. 1) then 1427 res@xyDashPattern = 2 1428 plot_wq = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1429 else 1430 res@gsnLeftString = vNam(varn) 1431 res@gsnRightString = unit(varn) 1432 res@trXMinF = miniwq 1433 res@trXMaxF = maxiwq 1434 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1435 end if 1436 end if 1437 1438 if (vNam(varn) .EQ. "wpqvp") then 1439 miniwpqvp=min(data(varn,:,:)) 1440 maxiwpqvp=max(data(varn,:,:)) 1441 if (over .EQ. 1) then 1442 res@xyDashPattern = 0 1443 plot_wpqvp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1444 else 1445 res@gsnLeftString = vNam(varn) 1446 res@gsnRightString = unit(varn) 1447 res@trXMinF = miniwpqvp 1448 res@trXMaxF = maxiwpqvp 1449 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1450 end if 1451 end if 1452 if (vNam(varn) .EQ. "wsqvs") then 1453 miniwsqvs=min(data(varn,:,:)) 1454 maxiwsqvs=max(data(varn,:,:)) 1455 if (over .EQ. 1) then 1456 res@xyDashPattern = 1 1457 plot_wsqvs = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1458 else 1459 res@gsnLeftString = vNam(varn) 1460 res@gsnRightString = unit(varn) 1461 res@trXMinF = miniwsqvs 1462 res@trXMaxF = maxiwsqvs 1463 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1464 end if 1465 end if 1466 if (vNam(varn) .EQ. "wqv") then 1467 miniwqv=min(data(varn,:,:)) 1468 maxiwqv=max(data(varn,:,:)) 1469 if (over .EQ. 1) then 1470 res@xyDashPattern = 2 1471 plot_wqv = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1472 else 1473 res@gsnLeftString = vNam(varn) 1474 res@gsnRightString = unit(varn) 1475 res@trXMinF = miniwqv 1476 res@trXMaxF = maxiwqv 1477 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1478 end if 1479 end if 1480 1481 if (vNam(varn) .EQ. "wpsp") then 1482 miniwpsp=min(data(varn,:,:)) 1483 maxiwpsp=max(data(varn,:,:)) 1484 if (over .EQ. 1) then 1485 res@xyDashPattern = 0 1486 plot_wpsp = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1487 else 1488 res@gsnLeftString = vNam(varn) 1489 res@gsnRightString = unit(varn) 1490 res@trXMinF = miniwpsp 1491 res@trXMaxF = maxiwpsp 1492 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1493 end if 1494 end if 1495 if (vNam(varn) .EQ. "wsss") then 1496 miniwsss=min(data(varn,:,:)) 1497 maxiwsss=max(data(varn,:,:)) 1498 if (over .EQ. 1) then 1499 res@xyDashPattern = 1 1500 plot_wsss = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1501 else 1502 res@gsnLeftString = vNam(varn) 1503 res@gsnRightString = unit(varn) 1504 res@trXMinF = miniwsss 1505 res@trXMaxF = maxiwsss 1506 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1507 end if 1508 end if 1509 if (vNam(varn) .EQ. "ws") then 1510 miniws=min(data(varn,:,:)) 1511 maxiws=max(data(varn,:,:)) 1512 if (over .EQ. 1) then 1513 res@xyDashPattern = 2 1514 plot_ws = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1515 else 1516 res@gsnLeftString = vNam(varn) 1517 res@gsnRightString = unit(varn) 1518 res@trXMinF = miniws 1519 res@trXMaxF = maxiws 1520 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1521 end if 1522 end if 1523 1524 if (vNam(varn) .EQ. "wpsap") then 1525 miniwpsap=min(data(varn,:,:)) 1526 maxiwpsap=max(data(varn,:,:)) 1527 if (over .EQ. 1) then 1528 res@xyDashPattern = 0 1529 plot_wpsap = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1530 else 1531 res@gsnLeftString = vNam(varn) 1532 res@gsnRightString = unit(varn) 1533 res@trXMinF = miniwpsap 1534 res@trXMaxF = maxiwpsap 1535 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1536 end if 1537 end if 1538 if (vNam(varn) .EQ. "wssas") then 1539 miniwssas=min(data(varn,:,:)) 1540 maxiwssas=max(data(varn,:,:)) 1541 if (over .EQ. 1) then 1542 res@xyDashPattern = 1 1543 plot_wssas = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1544 else 1545 res@gsnLeftString = vNam(varn) 1546 res@gsnRightString = unit(varn) 1547 res@trXMinF = miniwssas 1548 res@trXMaxF = maxiwssas 1549 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1550 end if 1551 end if 1552 if (vNam(varn) .EQ. "wsa") then 1553 miniwsa=min(data(varn,:,:)) 1554 maxiwsa=max(data(varn,:,:)) 1555 if (over .EQ. 1) then 1556 res@xyDashPattern = 2 1557 plot_wsa = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1558 else 1559 res@gsnLeftString = vNam(varn) 1560 res@gsnRightString = unit(varn) 1561 res@trXMinF = miniwsa 1562 res@trXMaxF = maxiwsa 1563 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1564 end if 1565 end if 1566 1567 if (vNam(varn) .EQ. "wses") then 1568 if (over .EQ. 1) then 1569 res@gsnLeftString = vNam(varn) 1570 res@gsnRightString = unit(varn) 1571 res@trXMinF = min(data(varn,:,:)) 1572 res@trXMaxF = max(data(varn,:,:)) 1573 plot_wses = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1574 else 1575 res@gsnLeftString = vNam(varn) 1576 res@gsnRightString = unit(varn) 1577 res@trXMinF = min(data(varn,:,:)) 1578 res@trXMaxF = max(data(varn,:,:)) 1579 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1580 end if 1581 end if 1582 1583 if (vNam(varn) .EQ. "us2") then 1584 minius2=min(data(varn,:,:)) 1585 maxius2=max(data(varn,:,:)) 1586 if (over .EQ. 1) then 1587 res@xyDashPattern = 0 1588 plot_us2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1589 else 1590 res@gsnLeftString = vNam(varn) 1591 res@gsnRightString = unit(varn) 1592 res@trXMinF = minius2 1593 res@trXMaxF = maxius2 1594 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1595 end if 1596 end if 1597 if (vNam(varn) .EQ. "vs2") then 1598 minivs2=min(data(varn,:,:)) 1599 maxivs2=max(data(varn,:,:)) 1600 if (over .EQ. 1) then 1601 res@xyDashPattern = 1 1602 plot_vs2 = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1603 else 1604 res@gsnLeftString = vNam(varn) 1605 res@gsnRightString = unit(varn) 1606 res@trXMinF = minivs2 1607 res@trXMaxF = maxivs2 1608 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1609 end if 1610 end if 1611 if (vNam(varn) .EQ. "ws2") then 1612 miniws2=min(data(varn,:,:)) 1613 maxiws2=max(data(varn,:,:)) 1614 if (over .EQ. 1) then 1615 res@xyDashPattern = 2 1616 plot_ws2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1617 else 1618 res@gsnLeftString = vNam(varn) 1619 res@gsnRightString = unit(varn) 1620 res@trXMinF = miniws2 1621 res@trXMaxF = maxiws2 1622 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1623 end if 1624 end if 1625 1626 if (vNam(varn) .EQ. "pts2") then 1627 minipts2=min(data(varn,:,:)) 1628 maxipts2=max(data(varn,:,:)) 1629 if (over .EQ. 1) then 1630 res@gsnLeftString = vNam(varn) 1631 res@gsnRightString = unit(varn) 1632 res@trXMinF = minipts2 1633 res@trXMaxF = maxipts2 1634 plot_pts2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1635 else 1636 res@gsnLeftString = vNam(varn) 1637 res@gsnRightString = unit(varn) 1638 res@trXMinF = minipts2 1639 res@trXMaxF = maxipts2 1640 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1641 end if 1642 end if 1643 1644 if (vNam(varn) .EQ. "ws3") then 1645 miniws3=min(data(varn,:,:)) 1646 maxiws3=max(data(varn,:,:)) 1647 if (over .EQ. 1) then 1648 res@gsnLeftString = vNam(varn) 1649 res@gsnRightString = unit(varn) 1650 res@trXMinF = miniws3 1651 res@trXMaxF = maxiws3 1652 plot_ws3 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1653 else 1654 res@gsnLeftString = vNam(varn) 1655 res@gsnRightString = unit(varn) 1656 res@trXMinF = miniws3 1657 res@trXMaxF = maxiws3 1658 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1659 end if 1660 end if 1661 1662 if (vNam(varn) .EQ. "Sw") then 1663 miniSw=min(data(varn,:,:)) 1664 maxiSw=max(data(varn,:,:)) 1665 if (over .EQ. 1) then 1666 res@gsnLeftString = vNam(varn) 1667 res@gsnRightString = unit(varn) 1668 res@trXMinF = miniSw 1669 res@trXMaxF = maxiSw 1670 plot_Sw = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1671 else 1672 res@gsnLeftString = vNam(varn) 1673 res@gsnRightString = unit(varn) 1674 res@trXMinF = miniSw 1675 res@trXMaxF = maxiSw 1676 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1677 end if 1678 end if 1679 1680 if (vNam(varn) .EQ. "ws2pts") then 1681 miniws2pts=min(data(varn,:,:)) 1682 maxiws2pts=max(data(varn,:,:)) 1683 if (over .EQ. 1) then 1684 res@gsnLeftString = vNam(varn) 1685 res@gsnRightString = unit(varn) 1686 res@trXMinF = miniws2pts 1687 res@trXMaxF = maxiws2pts 1688 plot_ws2pts = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1689 else 1690 res@gsnLeftString = vNam(varn) 1691 res@gsnRightString = unit(varn) 1692 res@trXMinF = miniws2pts 1693 res@trXMaxF = maxiws2pts 1694 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1695 end if 1696 end if 1697 1698 if (vNam(varn) .EQ. "wspts2") then 1699 miniwspts2=min(data(varn,:,:)) 1700 maxiwspts2=max(data(varn,:,:)) 1701 if (over .EQ. 1) then 1702 res@gsnLeftString = vNam(varn) 1703 res@gsnRightString = unit(varn) 1704 res@trXMinF = miniwspts2 1705 res@trXMaxF = maxiwspts2 1706 plot_wspts2 = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1707 else 1708 res@gsnLeftString = vNam(varn) 1709 res@gsnRightString = unit(varn) 1710 res@trXMinF = miniwspts2 1711 res@trXMaxF = maxiwspts2 1712 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1713 end if 1714 end if 1715 1716 if (vNam(varn) .EQ. "wsususodz") then 1717 miniwsususodz=min(data(varn,:,:)) 1718 maxiwsususodz=max(data(varn,:,:)) 1719 if (over .EQ. 1) then 1720 res@xyDashPattern = 0 1721 plot_wsususodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1722 else 1723 res@gsnLeftString = vNam(varn) 1724 res@gsnRightString = unit(varn) 1725 res@trXMinF = miniwsususodz 1726 res@trXMaxF = maxiwsususodz 1727 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1728 end if 1729 end if 1730 if (vNam(varn) .EQ. "wspsodz") then 1731 miniwspsodz=min(data(varn,:,:)) 1732 maxiwspsodz=max(data(varn,:,:)) 1733 if (over .EQ. 1) then 1734 res@xyDashPattern = 1 1735 plot_wspsodz = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1736 else 1737 res@gsnLeftString = vNam(varn) 1738 res@gsnRightString = unit(varn) 1739 res@trXMinF = miniwspsodz 1740 res@trXMaxF = maxiwspsodz 1741 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1742 end if 1743 end if 1744 if (vNam(varn) .EQ. "wpeodz") then 1745 miniwpeodz=min(data(varn,:,:)) 1746 maxiwpeodz=max(data(varn,:,:)) 1747 if (over .EQ. 1) then 1748 res@xyDashPattern = 2 1749 plot_wpeodz = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1750 else 1751 res@gsnLeftString = vNam(varn) 1752 res@gsnRightString = unit(varn) 1753 res@trXMinF = miniwpeodz 1754 res@trXMaxF = maxiwpeodz 1755 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1756 end if 1757 end if 1758 if (over .EQ. 0) then 1759 n=n+1 1760 end if 1761 varn=varn+1 1762 end if 1763 delete(temp) 1764 end do 1765 1766 if (over .EQ. 1 ) then 1767 1768 overlay(plot_u,plot_v) 1769 overlay(plot_u,plot_w) 1770 u=0 1771 overlay(plot_pt,plot_vpt) 1772 overlay(plot_pt,plot_lpt) 1773 pt=0 1774 overlay(plot_q,plot_qv) 1775 overlay(plot_q,plot_ql) 1776 q=0 1777 overlay(plot_e,plot_es) 1778 e=0 1779 overlay(plot_km,plot_kh) 1780 km=0 1781 overlay(plot_wpup,plot_wsus) 1782 overlay(plot_wpup,plot_wu) 1783 wpup=0 1784 overlay(plot_wpvp,plot_wsvs) 1785 overlay(plot_wpvp,plot_wv) 1786 wpvp=0 1787 overlay(plot_wpptp,plot_wspts) 1788 overlay(plot_wpptp,plot_wpt) 1789 wpptp=0 1790 overlay(plot_wsptsBC,plot_wptBC) 1791 wsptsBC=0 1792 overlay(plot_wpvptp,plot_wsvpts) 1793 overlay(plot_wpvptp,plot_wvpt) 1794 wpvptp=0 1795 overlay(plot_wpqp,plot_wsqs) 1796 overlay(plot_wpqp,plot_wq) 1797 wpqp=0 1798 overlay(plot_wpqvp,plot_wsqvs) 1799 overlay(plot_wpqvp,plot_wqv) 1800 wpqvp=0 1801 overlay(plot_wpsp,plot_wsss) 1802 overlay(plot_wpsp,plot_ws) 1803 wpsp=0 1804 overlay(plot_wpsap,plot_wssas) 1805 overlay(plot_wpsap,plot_wsa) 1806 wpsap=0 1807 overlay(plot_us2,plot_vs2) 1808 overlay(plot_us2,plot_ws2) 1809 us2=0 1810 overlay(plot_wsususodz,plot_wspsodz) 1811 overlay(plot_wsususodz,plot_wpeodz) 1812 wsususodz=0 1813 1814 end if 1815 1816 if (over .EQ. 1) then 1817 1818 do varn = 0,dim-1 1819 1820 if ( isStrSubset( vNam(varn), "time") .OR. isStrSubset( vNam(varn), "NORM")) then 1821 check = False 1822 else 1823 if (.not. isvar("var")) then 1824 check = True 1825 if (parameter(21) .NE. "variables") then 1826 var=parameter(21) 1827 check = isStrSubset( var,","+vNam(varn)+"," ) 1828 end if 1829 else 1830 check = isStrSubset( var,","+vNam(varn)+"," ) 1831 end if 1832 end if 1833 1834 if (check) 1835 1836 z = f->$vNam(varn+1)$ 1837 1838 ;res@gsnLeftString = vNam(varn) 1839 res@gsnRightString = unit(varn) 1840 res@trYMinF = min_z 1841 res@trYMaxF = max_z 1842 res@trXMinF = min(data(varn,:,:)) 1843 res@trXMaxF = max(data(varn,:,:)) 1844 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 1845 1846 if (vNam(varn) .EQ. "u" .OR. vNam(varn) .EQ. "v" .OR. vNam(varn) .EQ. "w") then 1847 if (u .EQ. 0) then 1848 res@gsnLeftString = "u, v and w" 1849 res@gsnRightString = unit(varn) 1850 res@trXMinF = min((/miniu,miniv,miniw/)) 1851 res@trXMaxF = max((/maxiu,maxiv,maxiw/)) 1852 if (vNam(varn) .EQ. "w") then 1853 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1854 else 1855 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1856 end if 1857 1858 ; *************************************************** 1859 ; legend for combined plot 1860 ; *************************************************** 1861 1862 lgres = True 1863 lgMonoDashIndex = False 1864 lgres@lgLabelFont = "helvetica" 1865 lgres@lgLabelFontHeightF = .1 1866 lgres@vpWidthF = 0.12 1867 lgres@vpHeightF = 0.1 1868 lgres@lgDashIndexes = (/0,1,2/) 1869 lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres) 1870 1871 amres = True 1872 amres@amParallelPosF = 0.65 1873 amres@amOrthogonalPosF = -0.2 1874 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 1875 overlay(plot(n),plot_u) 1876 u=1 1877 else 1878 varn=varn+1 1879 continue 1880 end if 1881 end if 438 1882 439 if (combine .EQ. 1) then 440 if (c_var .EQ. " ") then 441 print(" ") 442 print("Please indicate the variables you would like to combine ('c_var')") 443 print(" ") 444 exit 445 end if 446 com=isStrSubset(c_var,vNam(varn)+",") 447 if (com) then 448 e=e+1 1883 if (vNam(varn) .EQ. "pt" .OR. vNam(varn) .EQ. "vpt" .OR. vNam(varn) .EQ. "lpt") then 1884 if (pt .EQ. 0) then 1885 res@gsnLeftString = "pt, vpt and lpt" 1886 res@gsnRightString = unit(varn) 1887 res@trXMinF = min((/minipt,minivpt,minilpt/)) 1888 res@trXMaxF = max((/maxipt,maxivpt,maxilpt/)) 1889 1890 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1891 1892 ; *************************************************** 1893 ; legend for combined plot 1894 ; *************************************************** 1895 1896 lgres = True 1897 lgMonoDashIndex = False 1898 lgres@lgLabelFont = "helvetica" 1899 lgres@lgLabelFontHeightF = .1 1900 lgres@vpWidthF = 0.12 1901 lgres@vpHeightF = 0.1 1902 lgres@lgDashIndexes = (/0,1,2/) 1903 lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres) 1904 1905 amres = True 1906 amres@amParallelPosF = 0.65 1907 amres@amOrthogonalPosF = -0.2 1908 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 1909 overlay(plot(n),plot_pt) 1910 pt=1 1911 else 1912 varn=varn+1 1913 continue 1914 end if 1915 end if 1916 if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then 1917 if (q .EQ. 0) then 1918 res@gsnLeftString = "q, qv and ql" 1919 res@gsnRightString = unit(varn) 1920 res@trXMinF = min((/miniq,miniqv,miniql/)) 1921 res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) 1922 1923 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1924 1925 ; *************************************************** 1926 ; legend for combined plot 1927 ; *************************************************** 1928 1929 lgres = True 1930 lgMonoDashIndex = False 1931 lgres@lgLabelFont = "helvetica" 1932 lgres@lgLabelFontHeightF = .1 1933 lgres@vpWidthF = 0.12 1934 lgres@vpHeightF = 0.1 1935 lgres@lgDashIndexes = (/0,1,2/) 1936 lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres) 1937 1938 amres = True 1939 amres@amParallelPosF = 0.65 1940 amres@amOrthogonalPosF = -0.2 1941 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 1942 overlay(plot(n),plot_q) 1943 q=1 1944 else 1945 varn=varn+1 1946 continue 1947 end if 1948 end if 1949 if (vNam(varn) .EQ. "rho") then 1950 plot(n) = plot_rho 1951 end if 1952 if (vNam(varn) .EQ. "s") then 1953 plot(n) = plot_s 1954 end if 1955 if (vNam(varn) .EQ. "sa") then 1956 plot(n) = plot_sa 1957 end if 1958 if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then 1959 if (e .EQ. 0) then 1960 res@gsnLeftString = "e and es" 1961 res@gsnRightString = unit(varn) 1962 res@trXMinF = min((/minie,minies/)) 1963 res@trXMaxF = max((/maxie,maxies/)) 1964 1965 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1966 1967 ; *************************************************** 1968 ; legend for combined plot 1969 ; *************************************************** 1970 1971 lgres = True 1972 lgMonoDashIndex = False 1973 lgres@lgLabelFont = "helvetica" 1974 lgres@lgLabelFontHeightF = .1 1975 lgres@vpWidthF = 0.12 1976 lgres@vpHeightF = 0.1 1977 lgres@lgDashIndexes = (/0,1,2/) 1978 lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres) 1979 1980 amres = True 1981 amres@amParallelPosF = 0.65 1982 amres@amOrthogonalPosF = -0.2 1983 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 1984 overlay(plot(n),plot_e) 1985 e=1 1986 else 1987 varn=varn+1 1988 continue 1989 end if 1990 end if 1991 if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then 1992 if (km .EQ. 0) then 1993 res@gsnLeftString = "km and kh" 1994 res@gsnRightString = unit(varn) 1995 res@trXMinF = min((/minikm,minikh/)) 1996 res@trXMaxF = max((/maxikm,maxikh/)) 1997 1998 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1999 2000 ; *************************************************** 2001 ; legend for combined plot 2002 ; *************************************************** 2003 2004 lgres = True 2005 lgMonoDashIndex = False 2006 lgres@lgLabelFont = "helvetica" 2007 lgres@lgLabelFontHeightF = .1 2008 lgres@vpWidthF = 0.12 2009 lgres@vpHeightF = 0.1 2010 lgres@lgDashIndexes = (/0,1,2/) 2011 lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres) 2012 2013 amres = True 2014 amres@amParallelPosF = 0.65 2015 amres@amOrthogonalPosF = -0.2 2016 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2017 overlay(plot(n),plot_km) 2018 km=1 2019 else 2020 varn=varn+1 2021 continue 2022 end if 2023 end if 2024 if (vNam(varn) .EQ. "l") then 2025 plot(n) = plot_l 2026 end if 2027 if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then 2028 if (wpup .EQ. 0) then 2029 res@gsnLeftString = "wpup, wsus and wu" 2030 res@gsnRightString = unit(varn) 2031 res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) 2032 res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) 2033 2034 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2035 2036 ; *************************************************** 2037 ; legend for combined plot 2038 ; *************************************************** 2039 2040 lgres = True 2041 lgMonoDashIndex = False 2042 lgres@lgLabelFont = "helvetica" 2043 lgres@lgLabelFontHeightF = .1 2044 lgres@vpWidthF = 0.12 2045 lgres@vpHeightF = 0.1 2046 lgres@lgDashIndexes = (/0,1,2/) 2047 lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres) 2048 2049 amres = True 2050 amres@amParallelPosF = 0.65 2051 amres@amOrthogonalPosF = -0.2 2052 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2053 overlay(plot(n),plot_wpup) 2054 wpup=1 2055 else 2056 varn=varn+1 2057 continue 2058 end if 2059 end if 2060 if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then 2061 if (wpvp .EQ. 0) then 2062 res@gsnLeftString = "wpvp, wsus and wv" 2063 res@gsnRightString = unit(varn) 2064 res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) 2065 res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) 2066 2067 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2068 2069 ; *************************************************** 2070 ; legend for combined plot 2071 ; *************************************************** 2072 2073 lgres = True 2074 lgMonoDashIndex = False 2075 lgres@lgLabelFont = "helvetica" 2076 lgres@lgLabelFontHeightF = .1 2077 lgres@vpWidthF = 0.12 2078 lgres@vpHeightF = 0.1 2079 lgres@lgDashIndexes = (/0,1,2/) 2080 lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres) 2081 2082 amres = True 2083 amres@amParallelPosF = 0.65 2084 amres@amOrthogonalPosF = -0.2 2085 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2086 overlay(plot(n),plot_wpvp) 2087 wpup=0 2088 else 2089 varn=varn+1 2090 continue 2091 end if 2092 end if 2093 if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then 2094 if (wpptp .EQ. 0) then 2095 res@gsnLeftString = "wpptp, wspts and wv" 2096 res@gsnRightString = unit(varn) 2097 res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) 2098 res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) 2099 2100 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2101 2102 ; *************************************************** 2103 ; legend for combined plot 2104 ; *************************************************** 2105 2106 lgres = True 2107 lgMonoDashIndex = False 2108 lgres@lgLabelFont = "helvetica" 2109 lgres@lgLabelFontHeightF = .1 2110 lgres@vpWidthF = 0.12 2111 lgres@vpHeightF = 0.1 2112 lgres@lgDashIndexes = (/0,1,2/) 2113 lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres) 2114 2115 amres = True 2116 amres@amParallelPosF = 0.65 2117 amres@amOrthogonalPosF = -0.2 2118 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2119 overlay(plot(n),plot_wpptp) 2120 wpptp=1 2121 else 2122 varn=varn+1 2123 continue 2124 end if 2125 end if 2126 if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then 2127 if (wsptsBC .EQ. 0) then 2128 res@gsnLeftString = "wsptsBC and wptBC" 2129 res@gsnRightString = unit(varn) 2130 res@trXMinF = min((/miniwsptsBC,miniwptBC/)) 2131 res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) 2132 2133 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2134 2135 ; *************************************************** 2136 ; legend for combined plot 2137 ; *************************************************** 2138 2139 lgres = True 2140 lgMonoDashIndex = False 2141 lgres@lgLabelFont = "helvetica" 2142 lgres@lgLabelFontHeightF = .1 2143 lgres@vpWidthF = 0.12 2144 lgres@vpHeightF = 0.1 2145 lgres@lgDashIndexes = (/0,1,2/) 2146 lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres) 2147 2148 amres = True 2149 amres@amParallelPosF = 0.65 2150 amres@amOrthogonalPosF = -0.2 2151 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2152 overlay(plot(n),plot_wsptsBC) 2153 wsptsBC=1 2154 else 2155 varn=varn+1 2156 continue 2157 end if 2158 end if 2159 if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then 2160 if (wpvptp .EQ. 0) then 2161 res@gsnLeftString = "wpvptp, wsvpts and wv" 2162 res@gsnRightString = unit(varn) 2163 res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) 2164 res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) 2165 2166 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2167 2168 ; *************************************************** 2169 ; legend for combined plot 2170 ; *************************************************** 2171 2172 lgres = True 2173 lgMonoDashIndex = False 2174 lgres@lgLabelFont = "helvetica" 2175 lgres@lgLabelFontHeightF = .1 2176 lgres@vpWidthF = 0.12 2177 lgres@vpHeightF = 0.1 2178 lgres@lgDashIndexes = (/0,1,2/) 2179 lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres) 2180 2181 amres = True 2182 amres@amParallelPosF = 0.65 2183 amres@amOrthogonalPosF = -0.2 2184 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2185 overlay(plot(n),plot_wpvptp) 2186 wpvptp=1 2187 else 2188 varn=varn+1 2189 continue 2190 end if 2191 end if 2192 if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then 2193 if (wpqp .EQ. 0) then 2194 res@gsnLeftString = "wpqp, wsqs and wq" 2195 res@gsnRightString = unit(varn) 2196 res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) 2197 res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) 2198 2199 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2200 2201 ; *************************************************** 2202 ; legend for combined plot 2203 ; *************************************************** 2204 2205 lgres = True 2206 lgMonoDashIndex = False 2207 lgres@lgLabelFont = "helvetica" 2208 lgres@lgLabelFontHeightF = .1 2209 lgres@vpWidthF = 0.12 2210 lgres@vpHeightF = 0.1 2211 lgres@lgDashIndexes = (/0,1,2/) 2212 lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres) 2213 2214 amres = True 2215 amres@amParallelPosF = 0.65 2216 amres@amOrthogonalPosF = -0.2 2217 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2218 overlay(plot(n),plot_wpqp) 2219 wpqp=1 2220 else 2221 varn=varn+1 2222 continue 2223 end if 2224 end if 2225 if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then 2226 if (wpqvp .EQ. 0) then 2227 res@gsnLeftString = "wpqvp, wsqvs and wqv" 2228 res@gsnRightString = unit(varn) 2229 res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) 2230 res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) 2231 2232 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2233 2234 ; *************************************************** 2235 ; legend for combined plot 2236 ; *************************************************** 2237 2238 lgres = True 2239 lgMonoDashIndex = False 2240 lgres@lgLabelFont = "helvetica" 2241 lgres@lgLabelFontHeightF = .1 2242 lgres@vpWidthF = 0.12 2243 lgres@vpHeightF = 0.1 2244 lgres@lgDashIndexes = (/0,1,2/) 2245 lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres) 2246 2247 amres = True 2248 amres@amParallelPosF = 0.65 2249 amres@amOrthogonalPosF = -0.2 2250 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2251 overlay(plot(n),plot_wpqvp) 2252 wpqvp=1 2253 else 2254 varn=varn+1 2255 continue 2256 end if 2257 end if 2258 if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then 2259 if (wpsp .EQ. 0) then 2260 res@gsnLeftString = "wpsp, wsss and ws" 2261 res@gsnRightString = unit(varn) 2262 res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) 2263 res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) 2264 2265 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2266 2267 ; *************************************************** 2268 ; legend for combined plot 2269 ; *************************************************** 2270 2271 lgres = True 2272 lgMonoDashIndex = False 2273 lgres@lgLabelFont = "helvetica" 2274 lgres@lgLabelFontHeightF = .1 2275 lgres@vpWidthF = 0.12 2276 lgres@vpHeightF = 0.1 2277 lgres@lgDashIndexes = (/0,1,2/) 2278 lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres) 2279 2280 amres = True 2281 amres@amParallelPosF = 0.65 2282 amres@amOrthogonalPosF = -0.2 2283 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2284 overlay(plot(n),plot_wpsp) 2285 wpsp=1 2286 else 2287 varn=varn+1 2288 continue 2289 end if 2290 end if 2291 if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then 2292 if (wpsap .EQ. 0) then 2293 res@gsnLeftString = "wpsap, wssas and wsa" 2294 res@gsnRightString = unit(varn) 2295 res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) 2296 res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) 2297 2298 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2299 2300 ; *************************************************** 2301 ; legend for combined plot 2302 ; *************************************************** 2303 2304 lgres = True 2305 lgMonoDashIndex = False 2306 lgres@lgLabelFont = "helvetica" 2307 lgres@lgLabelFontHeightF = .1 2308 lgres@vpWidthF = 0.12 2309 lgres@vpHeightF = 0.1 2310 lgres@lgDashIndexes = (/0,1,2/) 2311 lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres) 2312 2313 amres = True 2314 amres@amParallelPosF = 0.65 2315 amres@amOrthogonalPosF = -0.2 2316 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2317 overlay(plot(n),plot_wpsap) 2318 wpsap=1 2319 else 2320 varn=varn+1 2321 continue 2322 end if 2323 end if 2324 if (vNam(varn) .EQ. "wses") then 2325 plot(n) = plot_wses 2326 end if 2327 if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then 2328 if (us2 .EQ. 0) then 2329 res@gsnLeftString = "us2, vs2 and ws2" 2330 res@gsnRightString = unit(varn) 2331 res@trXMinF = min((/minius2,minivs2,miniws2/)) 2332 res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) 2333 2334 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2335 2336 ; *************************************************** 2337 ; legend for combined plot 2338 ; *************************************************** 2339 2340 lgres = True 2341 lgMonoDashIndex = False 2342 lgres@lgLabelFont = "helvetica" 2343 lgres@lgLabelFontHeightF = .1 2344 lgres@vpWidthF = 0.12 2345 lgres@vpHeightF = 0.1 2346 lgres@lgDashIndexes = (/0,1,2/) 2347 lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres) 2348 2349 amres = True 2350 amres@amParallelPosF = 0.65 2351 amres@amOrthogonalPosF = -0.2 2352 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2353 overlay(plot(n),plot_us2) 2354 us2=1 2355 else 2356 varn=varn+1 2357 continue 2358 end if 2359 end if 2360 if (vNam(varn) .EQ. "pts2") then 2361 plot(n) = plot_pts2 2362 end if 2363 if (vNam(varn) .EQ. "ws3") then 2364 plot(n) = plot_ws3 2365 end if 2366 if (vNam(varn) .EQ. "Sw") then 2367 plot(n) = plot_Sw 2368 end if 2369 if (vNam(varn) .EQ. "ws2pts") then 2370 plot(n) = plot_ws2pts 2371 end if 2372 if (vNam(varn) .EQ. "wspts2") then 2373 plot(n) = plot_wspts2 2374 end if 2375 if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then 2376 if (wsususodz .EQ. 0) then 2377 res@gsnLeftString = "wsususodz, wspsodz and ws2" 2378 res@gsnRightString = unit(varn) 2379 res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) 2380 res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) 2381 2382 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 2383 2384 ; *************************************************** 2385 ; legend for combined plot 2386 ; *************************************************** 2387 2388 lgres = True 2389 lgMonoDashIndex = False 2390 lgres@lgLabelFont = "helvetica" 2391 lgres@lgLabelFontHeightF = .1 2392 lgres@vpWidthF = 0.12 2393 lgres@vpHeightF = 0.1 2394 lgres@lgDashIndexes = (/0,1,2/) 2395 lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres) 2396 2397 amres = True 2398 amres@amParallelPosF = 0.65 2399 amres@amOrthogonalPosF = -0.2 2400 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2401 overlay(plot(n),plot_wsususodz) 2402 wsususodz=1 2403 else 2404 varn=varn+1 2405 continue 2406 end if 2407 end if 2408 n=n+1 2409 varn=varn+1 2410 end if 2411 end do 2412 end if 2413 2414 if (combine .EQ. 1) then 2415 2416 check=True 2417 n_o=0 2418 do varn = 0,dim-1 2419 2420 if ( isStrSubset( vNam(varn), "NORM") .or. isStrSubset( vNam(varn), "time") ) 449 2421 check = False 450 temp = f->$vNam(varn)$ 451 data_o!1 = temp!0 452 data_o!2 = temp!1 453 data_o(e,:,:) = temp(start_time_step:end_time_step,:) 454 mini(e) = min(data_o(e,:,:)) 455 maxi(e) = max(data_o(e,:,:)) 456 457 do j=e,e 458 do p=start_time_step,end_time_step 459 legend_label_oplot(e*end_time_step+p-start_time_step)=" "+flt2string(doubletofloat(t_all(p)))+" "+vNam(varn) 460 end do 461 end do 462 463 delete(temp) 464 varn = varn + 1 2422 else 2423 if (.not. isvar("var")) then 2424 check = True 2425 if (parameter(21) .NE. "variables") then 2426 var=parameter(21) 2427 check = isStrSubset( var,","+vNam(varn)+"," ) 2428 end if 2429 else 2430 check = isStrSubset( var,","+vNam(varn)+"," ) 2431 end if 2432 end if 2433 2434 if(check) then 2435 2436 z = f->$vNam(varn+1)$ 2437 com=isStrSubset( c_var,","+vNam(varn)+"," ) 2438 2439 if (com) 2440 2441 if (n_o .EQ. 1) then 2442 res@xyDashPattern = 1 2443 2444 else 2445 if (n_o .EQ. 2) then 2446 res@xyDashPattern = 2 2447 else 2448 res@xyDashPattern = 0 2449 res@gsnLeftString = c_var 2450 res@gsnRightString = unit(varn) 2451 res@trXMinF = min(mini) 2452 res@trXMaxF = max(maxi) 2453 end if 2454 end if 2455 label(n_o)=vNam(varn) 2456 color_o(n_o)=237 2457 plot_o(n_o)=gsn_csm_xy(wks,data(varn,:,:),z,res) 2458 n_o=n_o+1 2459 end if 2460 varn=varn+1 2461 end if 2462 end do 2463 2464 if(number_comb .EQ. 2)then 2465 overlay(plot_o(0),plot_o(1)) 2466 end if 2467 if(number_comb .EQ. 3)then 2468 overlay(plot_o(0),plot_o(1)) 2469 overlay(plot_o(0),plot_o(2)) 2470 end if 2471 2472 ; *************************************************** 2473 ; legend for combined plot 2474 ; *************************************************** 2475 2476 lgres = True 2477 lgMonoDashIndex = False 2478 lgres@lgDashIndexes = (/0,1,2/) 2479 lgres@lgLabelFont = "helvetica" 2480 lgres@lgLabelFontHeightF = .1 2481 lgres@vpWidthF = 0.12 2482 lgres@vpHeightF = 0.1 2483 2484 lbid = gsn_create_legend(wks,number_comb,label,lgres) 2485 2486 amres = True 2487 amres@amParallelPosF = 0.65 2488 amres@amOrthogonalPosF = -0.2 2489 annoid1 = gsn_add_annotation(plot_o(0),lbid,amres) 465 2490 466 end if 467 end if 468 469 if(check) then 470 471 temp = f->$vNam(varn)$ 472 data = temp(start_time_step:end_time_step,:) 473 testx= data*100 474 print(" plot of " + vNam(varn)) 475 data!1 = "z" 476 477 res@gsnLeftString = "PROFILE plot of" 478 res@gsnRightString = vNam(varn) 479 res@txFontHeightF = 0.02 480 res@tiXAxisFontHeightF = 0.02 481 res@tiYAxisFontHeightF = 0.02 482 res@tiXAxisString = " " 483 res@tiYAxisString = "Height [z]" 484 res@trXMinF = min(data(:,:)) 485 res@trXMaxF = max(data(:,:)) 486 res@trYMinF = min_z 487 res@trYMaxF = max_z 488 489 plot(n) = gsn_csm_xy(wks,data,z_pr,res) 490 491 delete(temp) 492 delete(data) 493 n = n + 1 494 varn = varn + 1 495 496 end if 497 498 end do 499 500 if (combine .EQ. 1) then 501 data_all=new((/m*end_time_step,dimz/),float) 502 do j=0,e 503 data_all(start_time_step-1+j*end_time_step:(j+1)*end_time_step-1,:) = data_o(j,:,:) 504 end do 505 506 ores@gsnLeftString = "combined PROFILE plot of" 507 ores@gsnRightString = c_var 508 ores@txFontHeightF = 0.02 509 ores@tiXAxisFontHeightF = 0.02 510 ores@tiYAxisFontHeightF = 0.02 511 ores@tiXAxisString = " " 512 ores@tiYAxisString = "Height [z]" 513 ores@trXMinF = min(mini(:)) 514 ores@trXMaxF = max(maxi(:)) 515 ores@trYMinF = min_z 516 ores@trYMaxF = max_z 517 ores@xyExplicitLegendLabels = legend_label_oplot 518 519 print(" ") 520 print("combined plot of " + c_var) 521 522 plot(0) = gsn_csm_xy(wks,data_all,z_pr,ores) 523 2491 plot(0) = plot_o(0) 2492 524 2493 end if 525 2494 … … 527 2496 ; merge plots onto one page 528 2497 ; *************************************************** 529 530 do i = 0, n-1, no_lines*no_columns 531 if( (i+no_lines*no_columns) .gt. (n-1)) then 532 gsn_panel(wks,plot(i:n-1),(/no_lines,no_columns/),resP) 533 else 534 gsn_panel(wks,plot(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP) 535 end if 2498 print(n) 2499 do m=0,n-1 2500 plot_(m)=plot(n-1-m) 536 2501 end do 2502 2503 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 2504 gsn_panel(wks,plot_(0:n-1),(/n,1/),resP) 2505 else 2506 do i = 0,n-1, no_lines*no_columns 2507 if( (i+no_lines*no_columns) .gt. (n-1)) then 2508 gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP) 2509 else 2510 gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP) 2511 end if 2512 end do 2513 end if 537 2514 538 2515 print(" ") -
palm/trunk/SCRIPTS/NCL/timeseries.ncl
r157 r161 11 11 12 12 if (isfilepresent("~/.ncl_preferences")) then 13 parameter = asciiread("~/.ncl_preferences",7 3,"string")13 parameter = asciiread("~/.ncl_preferences",75,"string") 14 14 delete(parameter@_FillValue) 15 15 else … … 23 23 ; set up default parameter values and strings if not assigned in prompt or parameter list 24 24 ; *************************************************** 25 25 26 26 if ( .not. isvar("file_in") ) then ; path+name of input file 27 27 if (parameter(7) .EQ. "input file") then … … 60 60 if ( .not. isvar("var") ) then ; variable name 61 61 check = True 62 end if 63 if ( .not. isvar("over") ) then ; switches overlaying plots on 64 over = 0 65 if (stringtointeger(parameter(73)) .NE. 0) then 66 over = stringtointeger(parameter(73)) 67 if (stringtointeger(parameter(73)) .NE. 1) then 68 print(" ") 69 print("Please set 'over' to 0 or 1") 70 print(" ") 71 exit 72 end if 73 end if 62 74 end if 63 75 … … 68 80 f = addfile(file_in , "r" ) 69 81 70 vNam = getfilevarnames(f) 82 vNam = getfilevarnames(f) 71 83 print(" ") 72 84 print("Variable on netCDF file: " + vNam) 73 85 print(" ") 74 dim = dimsizes(vNam) 86 dim = dimsizes(vNam) 87 if (dim .EQ. 0) then 88 print(" ") 89 print("There are no data on file") 90 print(" ") 91 end if 92 75 93 t_all = f->time 76 nt 94 nt = dimsizes(t_all) 77 95 78 96 ; **************************************************** … … 81 99 82 100 if ( .not. isvar("start_time_step") ) then 83 start_time_step = 184 if (parameter(13) . LE. "1") then85 if (parameter(13) . EQ. "0")101 start_time_step = 0 102 if (parameter(13) .NE. "1") then 103 if (parameter(13) .LE. "0") 86 104 print(" ") 87 105 print("Begin with time step 1") … … 95 113 exit 96 114 end if 97 start_time_step = stringtointeger(parameter(13)) 115 start_time_step = stringtointeger(parameter(13))-1 98 116 end if 99 117 else … … 110 128 exit 111 129 end if 130 start_time_step = start_time_step - 1 112 131 end if 113 132 … … 137 156 exit 138 157 end if 139 end_time_step = stringtointeger(parameter(15)) 158 end_time_step = stringtointeger(parameter(15))-1 140 159 end if 141 160 else … … 153 172 end if 154 173 if (end_time_step .LT. start_time_step) 155 print(" ") 156 print("'end_time_step' = "+end_time_step +" is lower than 'start_time_step' = "+start_time_step) 157 print(" ") 158 exit 159 end if 174 print(" ") 175 print("'end_time_step' = "+end_time_step +" is lower than 'start_time_step' = "+start_time_step) 176 print(" ") 177 exit 178 end if 179 end_time_step = end_time_step - 1 160 180 end if 161 181 … … 181 201 res@tiYAxisFont = "helvetica" 182 202 res@tmXBLabelFont = "helvetica" 183 res@tmYLLabelFont = "helvetica" 184 res@xyLineColors = (/237/) 203 res@tmYLLabelFont = "helvetica" 204 res@xyLineColors = (/237/) 205 206 res@lgLabelFontHeightF = .02 207 185 208 resP = True 186 209 resP@txFont = "helvetica" … … 191 214 res@vpWidthF=4 192 215 216 txres = True 193 217 194 218 ; *************************************************** … … 198 222 wks_ps = gsn_open_wks(format_out,file_out) 199 223 gsn_define_colormap(wks_ps,"rainbow+white") 200 plot_ps=new(dim,graphic) 224 plot_ps=new(dim,graphic) 225 201 226 n=0 202 203 do varn = 0, dim-1 227 minE=1.E27 228 maxE=-1.E27 229 minus=1.E27 230 maxus=-1.E27 231 minu=1.E27 232 maxu=-1.E27 233 minz=1.E27 234 maxz=-1.E27 235 minw=1.E27 236 maxw=-1.E27 237 minp=1.E27 238 maxp=-1.E27 239 mins=1.E27 240 maxs=-1.E27 241 242 data = new((/dim,(end_time_step-start_time_step)+1/),float) 243 unit = new(dim,string) 244 data_0 = new((end_time_step-start_time_step)+1,float) 245 data_0 = 0.0 246 mini = new(dim,float) 247 maxi = new(dim,float) 248 249 if (over .EQ. 1) then 250 plot_E = gsn_csm_xy(wks_ps,t,data_0(:),res) 251 plot_Es = gsn_csm_xy(wks_ps,t,data_0(:),res) 252 plot_us = gsn_csm_xy(wks_ps,t,data_0(:),res) 253 plot_ws = gsn_csm_xy(wks_ps,t,data_0(:),res) 254 plot_umax = gsn_csm_xy(wks_ps,t,data_0(:),res) 255 plot_vmax = gsn_csm_xy(wks_ps,t,data_0(:),res) 256 plot_wmax = gsn_csm_xy(wks_ps,t,data_0(:),res) 257 plot_z_i_wpt = gsn_csm_xy(wks_ps,t,data_0(:),res) 258 plot_z_i_pt = gsn_csm_xy(wks_ps,t,data_0(:),res) 259 plot_wpptp0 = gsn_csm_xy(wks_ps,t,data_0(:),res) 260 plot_wpptp = gsn_csm_xy(wks_ps,t,data_0(:),res) 261 plot_wpt = gsn_csm_xy(wks_ps,t,data_0(:),res) 262 plot_pt_0_ = gsn_csm_xy(wks_ps,t,data_0(:),res) 263 plot_pt_zp_ = gsn_csm_xy(wks_ps,t,data_0(:),res) 264 plot_splptx = gsn_csm_xy(wks_ps,t,data_0(:),res) 265 plot_splpty = gsn_csm_xy(wks_ps,t,data_0(:),res) 266 plot_splptz = gsn_csm_xy(wks_ps,t,data_0(:),res) 267 end if 268 269 do varn = dim-1,0,1 270 271 if( isStrSubset (vNam(varn), "time") ) 272 check = False 273 else 274 check = True 275 end if 204 276 if( isvar("var") ) then 205 check = isStrSubset( var, vNam(varn)+"," )277 check = isStrSubset( var,","+vNam(varn)+"," ) 206 278 end if 207 279 if (parameter(21) .NE. "variables") then 208 280 var = parameter(21) 209 check = isStrSubset( var, vNam(varn)+"," )281 check = isStrSubset( var,","+vNam(varn)+"," ) 210 282 end if 211 283 284 285 if(check) then 286 287 data_all = f ->$vNam(varn)$ 288 unit(varn) = data_all@units 289 290 data(varn,:)=data_all(start_time_step:end_time_step) 291 292 if (over .EQ. 1) then 293 294 mini(varn) = min(data(varn,:)) 295 maxi(varn) = max(data(varn,:)) 296 297 if (vNam(varn) .EQ. "E" .OR. vNam(varn) .EQ. "Es") then 298 if (mini(varn) .EQ. maxi(varn)) then 299 if (min(data(varn,:)) .EQ. 0)then 300 mini(varn)= mini(varn)-1. 301 maxi(varn)= maxi(varn)+1. 302 end if 303 if (min(data(varn,:)) .LT. 0)then 304 mini(varn)= mini(varn)-1.+(mini(varn))/2 305 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 306 end if 307 if (min(data(varn,:)) .GT. 0)then 308 mini(varn)= mini(varn)-1.-(mini(varn))/2 309 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 310 end if 311 end if 312 minE=min((/minE,mini(varn)/)) 313 maxE=max((/maxE,maxi(varn)/)) 314 end if 315 316 if (vNam(varn) .EQ. "us" .OR. vNam(varn) .EQ. "ws") then 317 if (mini(varn) .EQ. maxi(varn)) then 318 if (min(data(varn,:)) .EQ. 0)then 319 mini(varn)= mini(varn)-1. 320 maxi(varn)= maxi(varn)+1. 321 end if 322 if (min(data(varn,:)) .LT. 0)then 323 mini(varn)= mini(varn)-1.+(mini(varn))/2 324 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 325 end if 326 if (min(data(varn,:)) .GT. 0)then 327 mini(varn)= mini(varn)-1.-(mini(varn))/2 328 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 329 end if 330 end if 331 minus=min((/minus,mini(varn)/)) 332 maxus=max((/maxus,maxi(varn)/)) 333 end if 334 335 if (vNam(varn) .EQ. "umax" .OR. vNam(varn) .EQ. "vmax" .OR. vNam(varn) .EQ. "wmax") then 336 if (mini(varn) .EQ. maxi(varn)) then 337 if (mini(varn) .EQ. 0)then 338 mini(varn)= mini(varn)-1. 339 maxi(varn)= maxi(varn)+1. 340 end if 341 if (mini(varn) .LT. 0)then 342 mini(varn)= mini(varn)-1.+(mini(varn))/2 343 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 344 end if 345 if (mini(varn) .GT. 0)then 346 mini(varn)= mini(varn)-1.-(mini(varn))/2 347 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 348 end if 349 end if 350 minu=min((/minu,mini(varn)/)) 351 maxu=max((/maxu,maxi(varn)/)) 352 end if 353 354 if (vNam(varn) .EQ. "z_i_wpt" .OR. vNam(varn) .EQ. "z_i_pt") then 355 if (mini(varn) .EQ. maxi(varn)) then 356 if (min(data(varn,:)) .EQ. 0)then 357 mini(varn)= mini(varn)-1. 358 maxi(varn)= maxi(varn)+1. 359 end if 360 if (min(data(varn,:)) .LT. 0)then 361 mini(varn)= mini(varn)-1.+(mini(varn))/2 362 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 363 end if 364 if (min(data(varn,:)) .GT. 0)then 365 mini(varn)= mini(varn)-1.-(mini(varn))/2 366 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 367 end if 368 end if 369 minz=min((/minz,mini(varn)/)) 370 maxz=max((/maxz,maxi(varn)/)) 371 end if 372 373 if (vNam(varn) .EQ. "wpptp0" .OR. vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wpt") then 374 if (mini(varn) .EQ. maxi(varn)) then 375 if (min(data(varn,:)) .EQ. 0)then 376 mini(varn)= mini(varn)-1. 377 maxi(varn)= maxi(varn)+1. 378 end if 379 if (min(data(varn,:)) .LT. 0)then 380 mini(varn)= mini(varn)-1.+(mini(varn))/2 381 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 382 end if 383 if (min(data(varn,:)) .GT. 0)then 384 mini(varn)= mini(varn)-1.-(mini(varn))/2 385 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 386 end if 387 end if 388 minw=min((/minw,mini(varn)/)) 389 maxw=max((/maxw,maxi(varn)/)) 390 end if 391 392 if (vNam(varn) .EQ. "pt_0_" .OR. vNam(varn) .EQ. "pt_zp_") then 393 if (mini(varn) .EQ. maxi(varn)) then 394 if (min(data(varn,:)) .EQ. 0)then 395 mini(varn)= mini(varn)-1. 396 maxi(varn)= maxi(varn)+1. 397 end if 398 if (min(data(varn,:)) .LT. 0)then 399 mini(varn)= mini(varn)-1.+(mini(varn))/2 400 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 401 end if 402 if (min(data(varn,:)) .GT. 0)then 403 mini(varn)= mini(varn)-1.-(mini(varn))/2 404 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 405 end if 406 end if 407 minp=min((/minp,mini(varn)/)) 408 maxp=max((/maxp,maxi(varn)/)) 409 end if 410 411 if (vNam(varn) .EQ. "splptx" .OR. vNam(varn) .EQ. "splpty" .OR. vNam(varn) .EQ. "splptz") then 412 if (mini(varn) .EQ. maxi(varn)) then 413 if (min(data(varn,:)) .EQ. 0)then 414 mini(varn)= mini(varn)-1. 415 maxi(varn)= maxi(varn)+1. 416 end if 417 if (min(data(varn,:)) .LT. 0)then 418 mini(varn)= mini(varn)-1.+(mini(varn))/2 419 maxi(varn)= maxi(varn)+1.-(maxi(varn))/2 420 end if 421 if (min(data(varn,:)) .GT. 0)then 422 mini(varn)= mini(varn)-1.-(mini(varn))/2 423 maxi(varn)= maxi(varn)+1.+(maxi(varn))/2 424 end if 425 end if 426 mins=min((/mins,mini(varn)/)) 427 maxs=max((/maxs,maxi(varn)/)) 428 end if 429 430 end if 431 end if 432 end do 433 434 do varn = dim-1,0,1 435 212 436 if( isStrSubset (vNam(varn), "time") ) 213 437 check = False 214 end if 438 else 439 check = True 440 end if 441 if( isvar("var") ) then 442 check = isStrSubset( var,","+vNam(varn)+"," ) 443 end if 444 if (parameter(21) .NE. "variables") then 445 var = parameter(21) 446 check = isStrSubset( var,","+vNam(varn)+"," ) 447 end if 215 448 216 449 if(check) then 450 451 if (over .EQ. 1) then 217 452 218 n=n+1 219 data = f ->$vNam(varn)$ 220 print("plot of " + vNam(varn)) 221 222 res@gsnRightString = vNam(varn) 223 res@tiXAxisString = " time [s] " 224 res@tiXAxisFontHeightF = 0.07 225 res@txFontHeightF = 0.07 226 res@tiYAxisFontHeightF = 0.07 227 228 plot_ps(n) = gsn_csm_xy(wks_ps,t,data(start_time_step:end_time_step),res) 453 res@gsnLeftString = "overlayed plot" 454 res@gsnRightString = unit(varn) 455 res@tiXAxisString = " time [s] " 456 res@tiYAxisString = " " 457 res@tiXAxisFontHeightF = 0.07 458 res@txFontHeightF = 0.07 459 res@tiYAxisFontHeightF = 0.07 460 res@xyLabelMode = "Custom" 461 462 463 if (vNam(varn) .EQ. "E") 464 E=0 465 res@xyLineColors = (/237/) 466 res@xyExplicitLabels = "E" 467 res@xyLineLabelFontHeightF = 0.05 468 res@xyLineLabelFontColor = 237 469 res@trYMaxF = minE 470 res@trYMinF = maxE 471 plot_E = gsn_csm_xy(wks_ps,t,data(varn,:),res) 472 end if 473 if (vNam(varn) .EQ. "Es") 474 Es=0 475 res@xyLineColors = (/144/) 476 res@xyExplicitLabels = "Es" 477 res@xyLineLabelFontHeightF = 0.05 478 res@xyLineLabelFontColor = 144 479 plot_Es = gsn_csm_xy(wks_ps,t,data(varn,:),res) 480 end if 481 482 if (vNam(varn) .EQ. "us") 483 us=0 484 res@xyLineColors = (/237/) 485 res@xyExplicitLabels = "us" 486 res@xyLineLabelFontHeightF = 0.05 487 res@xyLineLabelFontColor = 237 488 res@trYMaxF = minus 489 res@trYMinF = maxus 490 plot_us = gsn_csm_xy(wks_ps,t,data(varn,:),res) 491 end if 492 if (vNam(varn) .EQ. "ws") 493 ws=0 494 res@xyLineColors = (/144/) 495 res@xyExplicitLabels = "ws" 496 res@xyLineLabelFontHeightF = 0.05 497 res@xyLineLabelFontColor = 144 498 plot_ws = gsn_csm_xy(wks_ps,t,data(varn,:),res) 499 end if 500 501 if (vNam(varn) .EQ. "umax") 502 u=0 503 res@xyLineColors = (/237/) 504 res@xyExplicitLabels = "umax" 505 res@xyLineLabelFontHeightF = 0.05 506 res@xyLineLabelFontColor = 237 507 res@trYMaxF = minu 508 res@trYMinF = maxu 509 plot_umax = gsn_csm_xy(wks_ps,t,data(varn,:),res) 510 end if 511 if (vNam(varn) .EQ. "vmax") 512 v=0 513 res@xyLineColors = (/144/) 514 res@xyExplicitLabels = "vmax" 515 res@xyLineLabelFontHeightF = 0.05 516 res@xyLineLabelFontColor = 144 517 plot_vmax = gsn_csm_xy(wks_ps,t,data(varn,:),res) 518 end if 519 if (vNam(varn) .EQ. "wmax") 520 w=0 521 res@xyLineColors = (/80/) 522 res@xyExplicitLabels = "wmax" 523 res@xyLineLabelFontHeightF = 0.05 524 res@xyLineLabelFontColor = 80 525 plot_wmax = gsn_csm_xy(wks_ps,t,data(varn,:),res) 526 end if 229 527 230 end if 528 if (vNam(varn) .EQ. "z_i_wpt") 529 zw=0 530 res@xyLineColors = (/237/) 531 res@xyExplicitLabels = "z_i_wpt" 532 res@xyLineLabelFontHeightF = 0.05 533 res@xyLineLabelFontColor = 237 534 res@trYMaxF = minz 535 res@trYMinF = maxz 536 plot_z_i_wpt = gsn_csm_xy(wks_ps,t,data(varn,:),res) 537 end if 538 if (vNam(varn) .EQ. "z_i_pt") 539 z=0 540 res@xyLineColors = (/144/) 541 res@xyExplicitLabels = "z_i_pt" 542 res@xyLineLabelFontHeightF = 0.05 543 res@xyLineLabelFontColor = 144 544 plot_z_i_pt = gsn_csm_xy(wks_ps,t,data(varn,:),res) 545 end if 546 547 if (vNam(varn) .EQ. "wpptp0") 548 w0=0 549 res@xyLineColors = (/237/) 550 res@xyExplicitLabels = "wpptp0" 551 res@xyLineLabelFontHeightF = 0.05 552 res@xyLineLabelFontColor = 237 553 res@trYMaxF = minw 554 res@trYMinF = maxw 555 plot_wpptp0 = gsn_csm_xy(wks_ps,t,data(varn,:),res) 556 end if 557 if (vNam(varn) .EQ. "wpptp") 558 wp=0 559 res@xyLineColors = (/144/) 560 res@xyExplicitLabels = "wpptp" 561 res@xyLineLabelFontHeightF = 0.05 562 res@xyLineLabelFontColor = 144 563 plot_wpptp = gsn_csm_xy(wks_ps,t,data(varn,:),res) 564 end if 565 if (vNam(varn) .EQ. "wpt") 566 wt=0 567 res@xyLineColors = (/80/) 568 res@xyExplicitLabels = "wpt" 569 res@xyLineLabelFontHeightF = 0.05 570 res@xyLineLabelFontColor = 80 571 plot_wpt = gsn_csm_xy(wks_ps,t,data(varn,:),res) 572 end if 573 574 if (vNam(varn) .EQ. "pt_0_") 575 p=0 576 res@xyLineColors = (/237/) 577 res@xyExplicitLabels = "pt_0_" 578 res@xyLineLabelFontHeightF = 0.05 579 res@xyLineLabelFontColor = 237 580 res@trYMaxF = minp 581 res@trYMinF = maxp 582 plot_pt_0_ = gsn_csm_xy(wks_ps,t,data(varn,:),res) 583 end if 584 if (vNam(varn) .EQ. "pt_zp_") 585 pz=0 586 res@xyLineColors = (/144/) 587 res@xyExplicitLabels = "pt_zp_" 588 res@xyLineLabelFontHeightF = 0.05 589 res@xyLineLabelFontColor = 144 590 plot_pt_zp_ = gsn_csm_xy(wks_ps,t,data(varn,:),res) 591 end if 592 593 if (vNam(varn) .EQ. "splptx") 594 x=0 595 res@xyLineColors = (/237/) 596 res@xyExplicitLabels = "splptx" 597 res@xyLineLabelFontHeightF = 0.05 598 res@xyLineLabelFontColor = 237 599 res@trYMaxF = mins 600 res@trYMinF = maxs 601 plot_splptx = gsn_csm_xy(wks_ps,t,data(varn,:),res) 602 end if 603 if (vNam(varn) .EQ. "splpty") 604 y=0 605 res@xyLineColors = (/144/) 606 res@xyExplicitLabels = "splpty" 607 res@xyLineLabelFontHeightF = 0.05 608 res@xyLineLabelFontColor = 144 609 plot_splpty = gsn_csm_xy(wks_ps,t,data(varn,:),res) 610 end if 611 if (vNam(varn) .EQ. "splptz") 612 z=0 613 res@xyLineColors = (/80/) 614 res@xyExplicitLabels = "splptz" 615 res@xyLineLabelFontHeightF = 0.05 616 res@xyLineLabelFontColor = 80 617 plot_splptz = gsn_csm_xy(wks_ps,t,data(varn,:),res) 618 end if 619 620 end if 621 end if 622 end do 623 624 do varn = dim-1,0,1 231 625 626 if( isStrSubset (vNam(varn), "time") ) 627 check = False 628 else 629 check = True 630 end if 631 if( isvar("var") ) then 632 check = isStrSubset( var,","+vNam(varn)+"," ) 633 end if 634 if (parameter(21) .NE. "variables") then 635 var = parameter(21) 636 check = isStrSubset( var,","+vNam(varn)+"," ) 637 end if 638 639 if(check) then 640 641 if (over .EQ. 1) then 642 643 if (vNam(varn) .EQ. "E" .AND. Es .NE. 1) then 644 E=1 645 overlay(plot_E,plot_Es) 646 n=n+1 647 plot_ps(n) = plot_E 648 end if 649 if (vNam(varn) .EQ. "Es" .AND. E .NE. 1) then 650 Es=1 651 overlay(plot_E,plot_Es) 652 n=n+1 653 plot_ps(n) = plot_E 654 end if 655 656 if (vNam(varn) .EQ. "us" .AND. ws .NE. 1) then 657 us=1 658 overlay(plot_us,plot_ws) 659 n=n+1 660 plot_ps(n) = plot_us 661 end if 662 if (vNam(varn) .EQ. "ws" .AND. us .NE. 1) then 663 ws=1 664 overlay(plot_us,plot_ws) 665 n=n+1 666 plot_ps(n) = plot_us 667 end if 668 669 if (vNam(varn) .EQ. "umax" .AND. v .NE. 1) 670 if (w .NE. 1) then 671 u=1 672 overlay(plot_umax,plot_vmax) 673 overlay(plot_umax,plot_wmax) 674 n=n+1 675 plot_ps(n) = plot_umax 676 end if 677 end if 678 if (vNam(varn) .EQ. "vmax" .AND. u .NE. 1) 679 if (w .NE. 1) then 680 v=1 681 overlay(plot_umax,plot_vmax) 682 overlay(plot_umax,plot_wmax) 683 n=n+1 684 plot_ps(n) = plot_umax 685 end if 686 end if 687 if (vNam(varn) .EQ. "wmax" .AND. v .NE. 1) 688 if(u .NE. 1) then 689 w=1 690 overlay(plot_umax,plot_vmax) 691 overlay(plot_umax,plot_wmax) 692 n=n+1 693 plot_ps(n) = plot_umax 694 end if 695 end if 696 697 if (vNam(varn) .EQ. "z_i_wpt" .AND. z .NE. 1) then 698 zw=1 699 overlay(plot_z_i_wpt,plot_z_i_pt) 700 n=n+1 701 plot_ps(n) = plot_z_i_wpt 702 end if 703 if (vNam(varn) .EQ. "z_i_pt" .AND. zw .NE. 1) then 704 z=1 705 overlay(plot_z_i_wpt,plot_z_i_pt) 706 n=n+1 707 plot_ps(n) = plot_z_i_wpt 708 end if 709 710 if (vNam(varn) .EQ. "wpptp0" .AND. wp .NE. 1) 711 if (wt .NE. 1) then 712 w0=1 713 overlay(plot_wpptp0,plot_wpptp) 714 overlay(plot_wpptp0,plot_wpt) 715 n=n+1 716 plot_ps(n) = plot_wpptp0 717 end if 718 end if 719 if (vNam(varn) .EQ. "wpptp" .AND. w0 .NE. 1) 720 if (wt .NE. 1) then 721 wp=1 722 overlay(plot_wpptp0,plot_wpptp) 723 overlay(plot_wpptp0,plot_wpt) 724 n=n+1 725 plot_ps(n) = plot_wpptp0 726 end if 727 end if 728 if (vNam(varn) .EQ. "wpt" .AND. wp .NE. 1) 729 if (w0 .NE. 1) then 730 wt=1 731 overlay(plot_wpptp0,plot_wpptp) 732 overlay(plot_wpptp0,plot_wpt) 733 n=n+1 734 plot_ps(n) = plot_wpptp0 735 end if 736 end if 737 738 if (vNam(varn) .EQ. "pt_0_" .AND. pz .NE. 1) then 739 p=1 740 overlay(plot_pt_0_,plot_pt_zp_) 741 n=n+1 742 plot_ps(n) = plot_pt_0_ 743 end if 744 if (vNam(varn) .EQ. "pt_zp_" .AND. p .NE. 1) then 745 pz=1 746 overlay(plot_pt_0_,plot_pt_zp_) 747 n=n+1 748 plot_ps(n) = plot_pt_0_ 749 end if 750 751 if (vNam(varn) .EQ. "splptx" .AND. y .NE. 1) 752 if (z .NE.1 ) then 753 x=1 754 overlay(plot_splptx,plot_splpty) 755 overlay(plot_splptx,plot_splptz) 756 n=n+1 757 plot_ps(n) = plot_splptx 758 end if 759 end if 760 if (vNam(varn) .EQ. "splpty" .AND. x .NE. 1) 761 if(z .NE.1 ) then 762 y=1 763 overlay(plot_splptx,plot_splpty) 764 overlay(plot_splptx,plot_splptz) 765 n=n+1 766 plot_ps(n) = plot_splptx 767 end if 768 end if 769 if (vNam(varn) .EQ. "splptz" .AND. y .NE. 1) 770 if(x .NE.1 ) then 771 z=1 772 overlay(plot_splptx,plot_splpty) 773 overlay(plot_splptx,plot_splptz) 774 n=n+1 775 plot_ps(n) = plot_splptx 776 end if 777 end if 778 779 if(vNam(varn) .NE. "splptz" .AND. vNam(varn) .NE. "splpty" .AND. vNam(varn) .NE. "splptx" .AND. vNam(varn) .NE. "pt_zp_" .AND. vNam(varn) .NE. "pt_0_" .AND. vNam(varn) .NE. "wpt" .AND. vNam(varn) .NE. "wpptp" .AND. vNam(varn) .NE. "wpptp0" .AND. vNam(varn) .NE. "z_i_pt" .AND. vNam(varn) .NE. "z_i_wpt" .AND. vNam(varn) .NE. "wmax" .AND. vNam(varn) .NE. "vmax" .AND. vNam(varn) .NE. "umax" .AND. vNam(varn) .NE. "ws" .AND. vNam(varn) .NE. "us" .AND. vNam(varn) .NE. "Es" .AND. vNam(varn) .NE. "E") then 780 n=n+1 781 res@xyLineColors = (/237/) 782 res@xyLabelMode = False 783 res@gsnLeftString = vNam(varn) 784 res@gsnRightString = unit(varn) 785 res@trYMaxF = min(data(varn,:)) 786 res@trYMinF = max(data(varn,:)) 787 if (min(data(varn,:)) .EQ. max(data(varn,:))) then 788 if (min(data(varn,:)) .EQ. 0)then 789 res@trYMaxF = min(data(varn,:))-1. 790 res@trYMinF = max(data(varn,:))+1. 791 end if 792 if (min(data(varn,:)) .LT. 0)then 793 res@trYMaxF = min(data(varn,:))+(min(data(varn,:)))/2 794 res@trYMinF = max(data(varn,:))-(max(data(varn,:)))/2 795 end if 796 if (min(data(varn,:)) .GT. 0)then 797 res@trYMaxF = min(data(varn,:))-(min(data(varn,:)))/2 798 res@trYMinF = max(data(varn,:))+(max(data(varn,:)))/2 799 end if 800 end if 801 plot_ps(n) = gsn_csm_xy(wks_ps,t,data(varn,:),res) 802 end if 803 804 else 805 806 print("plot of " + vNam(varn)) 807 808 n=n+1 809 res@xyLineColors = (/237/) 810 res@gsnLeftString = vNam(varn) 811 res@gsnRightString = unit(varn) 812 res@tiXAxisString = " time [s] " 813 res@tiYAxisString = " " 814 res@tiXAxisFontHeightF = 0.07 815 res@txFontHeightF = 0.07 816 res@tiYAxisFontHeightF = 0.07 817 res@trYMaxF = min(data(varn,:)) 818 res@trYMinF = max(data(varn,:)) 819 if (min(data(varn,:)) .EQ. max(data(varn,:))) then 820 if (min(data(varn,:)) .EQ. 0)then 821 res@trYMaxF = min(data(varn,:))-1. 822 res@trYMinF = max(data(varn,:))+1. 823 end if 824 if (min(data(varn,:)) .LT. 0)then 825 res@trYMaxF = min(data(varn,:))+(min(data(varn,:)))/2 826 res@trYMinF = max(data(varn,:))-(max(data(varn,:)))/2 827 end if 828 if (min(data(varn,:)) .GT. 0)then 829 res@trYMaxF = min(data(varn,:))-(min(data(varn,:)))/2 830 res@trYMinF = max(data(varn,:))+(max(data(varn,:)))/2 831 end if 832 end if 833 plot_ps(n) = gsn_csm_xy(wks_ps,t,data(varn,:),res) 834 835 end if 836 end if 232 837 end do 233 838 234 839 ; *************************************************** 235 840 ; merge plots onto one page 236 841 ; *************************************************** 237 238 do np = 1,n,no_lines*no_columns 239 240 if ( np + no_lines*no_columns .gt. n) then 241 gsn_panel(wks_ps, plot_ps(np:n),(/no_lines,no_columns/),resP) 242 else 243 gsn_panel(wks_ps, plot_ps(np:np+no_lines*no_columns-1),(/no_lines,no_columns/),resP) 244 end if 245 246 end do 842 843 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 844 gsn_panel(wks_ps,plot_ps(1:n),(/n,1/),resP) 845 else 846 do np = 1,n,no_lines*no_columns 847 if ( np + no_lines*no_columns .gt. n) then 848 gsn_panel(wks_ps, plot_ps(np:n),(/no_lines,no_columns/),resP) 849 else 850 gsn_panel(wks_ps, plot_ps(np:np+no_lines*no_columns-1),(/no_lines,no_columns/),resP) 851 end if 852 end do 853 end if 247 854 248 855 print(" ")
Note: See TracChangeset
for help on using the changeset viewer.
