Changeset 162 for palm/trunk/SCRIPTS/NCL/profiles.ncl
- Timestamp:
- Apr 30, 2008 1:41:13 PM (16 years ago)
- File:
-
- 1 edited
-
palm/trunk/SCRIPTS/NCL/profiles.ncl (modified) (75 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/profiles.ncl
r161 r162 107 107 if ( .not. isvar("over") ) then ; switches overlaying plots on 108 108 over = 0 109 if (stringtointeger(parameter( 73)) .NE. 0) then110 over = stringtointeger(parameter( 73))111 if (stringtointeger(parameter( 73)) .NE. 1) then109 if (stringtointeger(parameter(37)) .NE. 0) then 110 over = stringtointeger(parameter(37)) 111 if (stringtointeger(parameter(37)) .NE. 1) then 112 112 print(" ") 113 113 print("Please set 'over' to 0 or 1") … … 456 456 end do 457 457 end if 458 459 dimz = dimsizes(z_w) 460 dimz = dimsizes(z_u) 458 if (isvar("z_u") ) then 459 dimz = dimsizes(z_u) 460 else 461 if (isvar("z_w"))then 462 dimz = dimsizes(z_w) 463 end if 464 end if 461 465 t_all = f->time 462 466 nt = dimsizes(t_all) 467 delta_t=t_all(nt-1)/nt 463 468 464 469 ; **************************************************** … … 467 472 468 473 if ( .not. isvar("start_time_step") ) then 469 start_time_step = 1470 if (parameter(13) .NE. " 1") then471 if ( parameter(13) .LE. "1")474 start_time_step=t_all(1)/3600 475 if (parameter(13) .NE. "t(0)") then 476 if (stringtodouble(parameter(13)) .GT. t_all(nt-1)/3600) 472 477 print(" ") 473 print("Begin at least with time step 2") 478 print("'start_time_step' = "+ parameter(13) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h") 479 print(" ") 480 print("Please select another 'start_time_step'") 474 481 print(" ") 475 482 exit 476 483 end if 477 if (stringto integer(parameter(13)) .GE. nt)484 if (stringtofloat(parameter(13)) .LE. t_all(0)/3600) 478 485 print(" ") 479 print("'start_time_step' = "+ parameter(13) +" is greater than available time steps = " + (nt-1)) 486 print("'start_time_step' = "+ parameter(13) +"h is equal or lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h") 487 print(" ") 488 print("Please start profiles at least with second time step = "+t_all(1)+"s = "+t_all(1)/3600+"h") 480 489 print(" ") 481 490 exit 482 491 end if 483 start_time_step = stringtointeger(parameter(13))-1492 start_time_step=stringtodouble(parameter(13)) 484 493 end if 485 494 else 486 if (start_time_step . LE.0)495 if (start_time_step .GT. t_all(nt-1)/3600) 487 496 print(" ") 488 print("Begin at least with time step 1") 497 print("'start_time_step' = "+ start_time_step +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h") 498 print(" ") 499 print("Please select another 'start_time_step'") 489 500 print(" ") 490 501 exit 491 502 end if 492 if (start_time_step . GE. nt)503 if (start_time_step .LE. t_all(0)/3600) 493 504 print(" ") 494 print("'start_time_step' = "+ start_time_step +" is greater than available time steps = " + (nt-1)) 505 print("'start_time_step' = "+ start_time_step +"h is equal or lower than first time step = " + t_all(0)+"s = "+t_all(0)/3600+"h") 506 print(" ") 507 print("Please start profiles at least with second time step = "+t_all(1)+"s = "+t_all(1)/3600+"h") 495 508 print(" ") 496 509 exit 497 510 end if 498 start_time_step = start_time_step - 1499 511 end if 500 512 start_time_step = start_time_step*3600 513 if (start_time_step .GT. t_all(0) .AND. start_time_step .LT. t_all(1)+delta_t/2)then 514 st=1 515 end if 516 do i=2,nt-1 517 if (start_time_step .GE. t_all(i)-delta_t/2 .AND. start_time_step .LT. t_all(i)+delta_t/2)then 518 st=i 519 break 520 end if 521 end do 522 501 523 ; **************************************************** 502 524 ; end of time step and different types of mistakes that could be done … … 504 526 505 527 if ( .not. isvar("end_time_step") ) then 506 end_time_step = nt-1 507 if (parameter(15) .NE. "nt") then 508 if (parameter(15) .LE. "0") 509 print(" ") 510 print("'end_time_step' = "+parameter(15)+ " is too small; 'end_time_step' should be at least 1 ") 511 print(" ") 512 exit 513 end if 514 if (stringtointeger(parameter(15)) .GE. nt) 528 end_time_step = t_all(nt-1)/3600 529 if (parameter(15) .NE. "t(end)") then 530 if (stringtodouble(parameter(15)) .GT. t_all(nt-1)/3600) 515 531 print(" ") 516 print("'end_time_step' = "+ parameter(15) +" is greater than available time steps = " + (nt-1)) 532 print("'end_time_step' = "+ parameter(15) +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h") 533 print(" ") 534 print("Please select another 'end_time_step'") 517 535 print(" ") 518 536 exit 519 537 end if 520 if (stringto integer(parameter(15)) .LT. stringtointeger(parameter(13)))538 if (stringtodouble(parameter(15)) .LT. start_time_step/3600) 521 539 print(" ") 522 print("'end_time_step' = "+ parameter(15) +" is lower than 'start_time_step' = "+parameter(13)) 540 print("'end_time_step' = "+ parameter(15) +"h is lower than 'start_time_step' = "+parameter(13)+"h") 541 print(" ") 542 print("Please select another 'start_time_step' or 'end_time_step'") 523 543 print(" ") 524 544 exit 525 545 end if 526 end_time_step = stringto integer(parameter(15))546 end_time_step = stringtodouble(parameter(15)) 527 547 end if 528 548 else 529 if (end_time_step . LE.0)549 if (end_time_step .GT. t_all(nt-1)/3600) 530 550 print(" ") 531 print("'end_time_step' = "+end_time_step+ " is too small; 'end_time_step' should be at least 1 ") 551 print("'end_time_step' = "+ end_time_step +"h is greater than last time step = " + t_all(nt-1)+"s = "+t_all(nt-1)/3600+"h") 552 print(" ") 553 print("Please select another 'end_time_step'") 532 554 print(" ") 533 555 exit 534 556 end if 535 if (end_time_step . GE. nt)557 if (end_time_step .LT. start_time_step/3600) 536 558 print(" ") 537 print("'end_time_step' = "+ end_time_step +" is greater than available time steps = "+(nt-1)) 559 print("'end_time_step' = "+ end_time_step +"h is lower than 'start_time_step' = "+start_time_step+"h") 560 print(" ") 561 print("Please select another 'start_time_step' or 'end_time_step'") 538 562 print(" ") 539 563 exit 540 564 end if 541 if (end_time_step .LT. start_time_step) 542 print(" ") 543 print("'end_time_step' = "+end_time_step +" is lower than 'start_time_step' = "+start_time_step) 544 print(" ") 545 exit 565 end if 566 end_time_step = end_time_step*3600 567 if (end_time_step .GT. t_all(0) .AND. end_time_step .LT. t_all(1)+delta_t/2)then 568 et=1 569 end if 570 do i=2,nt-1 571 if (end_time_step .GE. t_all(i)-delta_t/2 .AND. end_time_step .LT. t_all(i)+delta_t/2)then 572 et=i 573 break 546 574 end if 547 end_time_step = end_time_step - 1 548 end if 575 end do 576 577 delete(start_time_step) 578 start_time_step=round(st,3) 579 delete(end_time_step) 580 end_time_step=round(et,3) 549 581 550 582 ; **************************************************** … … 557 589 end do 558 590 559 np = end_time_step-start_time_step 591 np = end_time_step-start_time_step+1 560 592 if ( black .eq. 0 ) then 561 593 res=True … … 887 919 888 920 n_o=0 921 count_var=0 889 922 do varn = 0,dim-1 890 923 … … 908 941 com=isStrSubset(c_var,","+vNam(varn)+"," ) 909 942 if (com) 910 check = False911 943 data(varn,:,:) = temp(start_time_step:end_time_step,0:dimz-1) 912 944 unit(varn) = temp@units … … 919 951 mini(n_o)=min(data(varn,:,:)) 920 952 maxi(n_o)=max(data(varn,:,:)) 921 varn=varn+1922 953 n_o=n_o+1 923 954 end if … … 925 956 926 957 if(check) then 927 958 count_var=count_var+1 959 928 960 z = f->$vNam(varn+1)$ 929 961 … … 935 967 res@gsnRightString = unit(varn) 936 968 res@trYMinF = min_z 937 res@trYMaxF = max_z 938 res@trXMinF = min(data(varn,:,:)) 939 res@trXMaxF = max(data(varn,:,:)) 969 res@trYMaxF = max_z 970 if (.not. isvar("xs")) then 971 if (parameter(63) .NE. "x0") then 972 res@trXMinF = stringtofloat(parameter(63)) 973 else 974 res@trXMinF = min(data(varn,:,:)) 975 end if 976 else 977 res@trXMinF = xs 978 end if 979 if (.not. isvar("xe")) then 980 if (parameter(65) .NE. "xdim") then 981 res@trXMaxF = stringtofloat(parameter(65)) 982 else 983 res@trXMaxF = max(data(varn,:,:)) 984 end if 985 else 986 res@trXMaxF = xe 987 end if 940 988 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 941 989 end if … … 950 998 res@gsnLeftString = vNam(varn) 951 999 res@gsnRightString = unit(varn) 952 res@trXMinF = miniu 953 res@trXMaxF = maxiu 1000 if (.not. isvar("xs")) then 1001 if (parameter(63) .NE. "x0") then 1002 res@trXMinF = stringtofloat(parameter(63)) 1003 else 1004 res@trXMinF = miniu 1005 end if 1006 else 1007 res@trXMinF = xs 1008 end if 1009 if (.not. isvar("xe")) then 1010 if (parameter(65) .NE. "xdim") then 1011 res@trXMaxF = stringtofloat(parameter(65)) 1012 else 1013 res@trXMaxF = maxiu 1014 end if 1015 else 1016 res@trXMaxF = xe 1017 end if 954 1018 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 955 1019 end if … … 964 1028 res@gsnLeftString = vNam(varn) 965 1029 res@gsnRightString = unit(varn) 966 res@trXMinF = miniv 967 res@trXMaxF = maxiv 1030 if (.not. isvar("xs")) then 1031 if (parameter(63) .NE. "x0") then 1032 res@trXMinF = stringtofloat(parameter(63)) 1033 else 1034 res@trXMinF = miniv 1035 end if 1036 else 1037 res@trXMinF = xs 1038 end if 1039 if (.not. isvar("xe")) then 1040 if (parameter(65) .NE. "xdim") then 1041 res@trXMaxF = stringtofloat(parameter(65)) 1042 else 1043 res@trXMaxF = maxiv 1044 end if 1045 else 1046 res@trXMaxF = xe 1047 end if 968 1048 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 969 1049 end if … … 978 1058 res@gsnLeftString = vNam(varn) 979 1059 res@gsnRightString = unit(varn) 980 res@trXMinF = miniw 981 res@trXMaxF = maxiw 1060 if (.not. isvar("xs")) then 1061 if (parameter(63) .NE. "x0") then 1062 res@trXMinF = stringtofloat(parameter(63)) 1063 else 1064 res@trXMinF = miniw 1065 end if 1066 else 1067 res@trXMinF = xs 1068 end if 1069 if (.not. isvar("xe")) then 1070 if (parameter(65) .NE. "xdim") then 1071 res@trXMaxF = stringtofloat(parameter(65)) 1072 else 1073 res@trXMaxF = maxiw 1074 end if 1075 else 1076 res@trXMaxF = xe 1077 end if 982 1078 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 983 1079 end if … … 993 1089 res@gsnLeftString = vNam(varn) 994 1090 res@gsnRightString = unit(varn) 995 res@trXMinF = minipt 996 res@trXMaxF = maxipt 1091 if (.not. isvar("xs")) then 1092 if (parameter(63) .NE. "x0") then 1093 res@trXMinF = stringtofloat(parameter(63)) 1094 else 1095 res@trXMinF = minipt 1096 end if 1097 else 1098 res@trXMinF = xs 1099 end if 1100 if (.not. isvar("xe")) then 1101 if (parameter(65) .NE. "xdim") then 1102 res@trXMaxF = stringtofloat(parameter(65)) 1103 else 1104 res@trXMaxF = maxipt 1105 end if 1106 else 1107 res@trXMaxF = xe 1108 end if 997 1109 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 998 1110 end if … … 1007 1119 res@gsnLeftString = vNam(varn) 1008 1120 res@gsnRightString = unit(varn) 1009 res@trXMinF = minivpt 1010 res@trXMaxF = maxivpt 1121 if (.not. isvar("xs")) then 1122 if (parameter(63) .NE. "x0") then 1123 res@trXMinF = stringtofloat(parameter(63)) 1124 else 1125 res@trXMinF = minivpt 1126 end if 1127 else 1128 res@trXMinF = xs 1129 end if 1130 if (.not. isvar("xe")) then 1131 if (parameter(65) .NE. "xdim") then 1132 res@trXMaxF = stringtofloat(parameter(65)) 1133 else 1134 res@trXMaxF = maxivpt 1135 end if 1136 else 1137 res@trXMaxF = xe 1138 end if 1011 1139 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1012 1140 end if … … 1021 1149 res@gsnLeftString = vNam(varn) 1022 1150 res@gsnRightString = unit(varn) 1023 res@trXMinF = minilpt 1024 res@trXMaxF = maxilpt 1151 if (.not. isvar("xs")) then 1152 if (parameter(63) .NE. "x0") then 1153 res@trXMinF = stringtofloat(parameter(63)) 1154 else 1155 res@trXMinF = minilpt 1156 end if 1157 else 1158 res@trXMinF = xs 1159 end if 1160 if (.not. isvar("xe")) then 1161 if (parameter(65) .NE. "xdim") then 1162 res@trXMaxF = stringtofloat(parameter(65)) 1163 else 1164 res@trXMaxF = maxilpt 1165 end if 1166 else 1167 res@trXMaxF = xe 1168 end if 1025 1169 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1026 1170 end if … … 1036 1180 res@gsnLeftString = vNam(varn) 1037 1181 res@gsnRightString = unit(varn) 1038 res@trXMinF = miniq 1039 res@trXMaxF = maxiq 1182 if (.not. isvar("xs")) then 1183 if (parameter(63) .NE. "x0") then 1184 res@trXMinF = stringtofloat(parameter(63)) 1185 else 1186 res@trXMinF = minilq 1187 end if 1188 else 1189 res@trXMinF = xs 1190 end if 1191 if (.not. isvar("xe")) then 1192 if (parameter(65) .NE. "xdim") then 1193 res@trXMaxF = stringtofloat(parameter(65)) 1194 else 1195 res@trXMaxF = maxilq 1196 end if 1197 else 1198 res@trXMaxF = xe 1199 end if 1040 1200 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1041 1201 end if … … 1050 1210 res@gsnLeftString = vNam(varn) 1051 1211 res@gsnRightString = unit(varn) 1052 res@trXMinF = miniqv 1053 res@trXMaxF = maxiqv 1212 if (.not. isvar("xs")) then 1213 if (parameter(63) .NE. "x0") then 1214 res@trXMinF = stringtofloat(parameter(63)) 1215 else 1216 res@trXMinF = minilqv 1217 end if 1218 else 1219 res@trXMinF = xs 1220 end if 1221 if (.not. isvar("xe")) then 1222 if (parameter(65) .NE. "xdim") then 1223 res@trXMaxF = stringtofloat(parameter(65)) 1224 else 1225 res@trXMaxF = maxilqv 1226 end if 1227 else 1228 res@trXMaxF = xe 1229 end if 1054 1230 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1055 1231 end if … … 1064 1240 res@gsnLeftString = vNam(varn) 1065 1241 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,:,:)) 1242 if (.not. isvar("xs")) then 1243 if (parameter(63) .NE. "x0") then 1244 res@trXMinF = stringtofloat(parameter(63)) 1245 else 1246 res@trXMinF = miniql 1247 end if 1248 else 1249 res@trXMinF = xs 1250 end if 1251 if (.not. isvar("xe")) then 1252 if (parameter(65) .NE. "xdim") then 1253 res@trXMaxF = stringtofloat(parameter(65)) 1254 else 1255 res@trXMaxF = maxiql 1256 end if 1257 else 1258 res@trXMaxF = xe 1259 end if 1116 1260 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1117 1261 end if … … 1127 1271 res@gsnLeftString = vNam(varn) 1128 1272 res@gsnRightString = unit(varn) 1129 res@trXMinF = minie 1130 res@trXMaxF = maxie 1273 if (.not. isvar("xs")) then 1274 if (parameter(63) .NE. "x0") then 1275 res@trXMinF = stringtofloat(parameter(63)) 1276 else 1277 res@trXMinF = minie 1278 end if 1279 else 1280 res@trXMinF = xs 1281 end if 1282 if (.not. isvar("xe")) then 1283 if (parameter(65) .NE. "xdim") then 1284 res@trXMaxF = stringtofloat(parameter(65)) 1285 else 1286 res@trXMaxF = maxie 1287 end if 1288 else 1289 res@trXMaxF = xe 1290 end if 1131 1291 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1132 1292 end if … … 1141 1301 res@gsnLeftString = vNam(varn) 1142 1302 res@gsnRightString = unit(varn) 1143 res@trXMinF = minies 1144 res@trXMaxF = maxies 1303 if (.not. isvar("xs")) then 1304 if (parameter(63) .NE. "x0") then 1305 res@trXMinF = stringtofloat(parameter(63)) 1306 else 1307 res@trXMinF = minies 1308 end if 1309 else 1310 res@trXMinF = xs 1311 end if 1312 if (.not. isvar("xe")) then 1313 if (parameter(65) .NE. "xdim") then 1314 res@trXMaxF = stringtofloat(parameter(65)) 1315 else 1316 res@trXMaxF = maxies 1317 end if 1318 else 1319 res@trXMaxF = xe 1320 end if 1145 1321 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1146 1322 end if … … 1156 1332 res@gsnLeftString = vNam(varn) 1157 1333 res@gsnRightString = unit(varn) 1158 res@trXMinF = minikm 1159 res@trXMaxF = maxikm 1334 if (.not. isvar("xs")) then 1335 if (parameter(63) .NE. "x0") then 1336 res@trXMinF = stringtofloat(parameter(63)) 1337 else 1338 res@trXMinF = minikm 1339 end if 1340 else 1341 res@trXMinF = xs 1342 end if 1343 if (.not. isvar("xe")) then 1344 if (parameter(65) .NE. "xdim") then 1345 res@trXMaxF = stringtofloat(parameter(65)) 1346 else 1347 res@trXMaxF = maxikm 1348 end if 1349 else 1350 res@trXMaxF = xe 1351 end if 1160 1352 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1161 1353 end if … … 1170 1362 res@gsnLeftString = vNam(varn) 1171 1363 res@gsnRightString = unit(varn) 1172 res@trXMinF = minikh 1173 res@trXMaxF = maxikh 1364 if (.not. isvar("xs")) then 1365 if (parameter(63) .NE. "x0") then 1366 res@trXMinF = stringtofloat(parameter(63)) 1367 else 1368 res@trXMinF = minikh 1369 end if 1370 else 1371 res@trXMinF = xs 1372 end if 1373 if (.not. isvar("xe")) then 1374 if (parameter(65) .NE. "xdim") then 1375 res@trXMaxF = stringtofloat(parameter(65)) 1376 else 1377 res@trXMaxF = maxikh 1378 end if 1379 else 1380 res@trXMaxF = xe 1381 end if 1174 1382 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1175 end if1176 end if1177 1178 if (vNam(varn) .EQ. "l") then1179 if (over .EQ. 1) then1180 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 else1186 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 1383 end if 1192 1384 end if … … 1201 1393 res@gsnLeftString = vNam(varn) 1202 1394 res@gsnRightString = unit(varn) 1203 res@trXMinF = miniwpup 1204 res@trXMaxF = maxiwpup 1395 if (.not. isvar("xs")) then 1396 if (parameter(63) .NE. "x0") then 1397 res@trXMinF = stringtofloat(parameter(63)) 1398 else 1399 res@trXMinF = miniwpup 1400 end if 1401 else 1402 res@trXMinF = xs 1403 end if 1404 if (.not. isvar("xe")) then 1405 if (parameter(65) .NE. "xdim") then 1406 res@trXMaxF = stringtofloat(parameter(65)) 1407 else 1408 res@trXMaxF = maxiwpup 1409 end if 1410 else 1411 res@trXMaxF = xe 1412 end if 1205 1413 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1206 1414 end if … … 1215 1423 res@gsnLeftString = vNam(varn) 1216 1424 res@gsnRightString = unit(varn) 1217 res@trXMinF = miniwsus 1218 res@trXMaxF = maxiwsus 1425 if (.not. isvar("xs")) then 1426 if (parameter(63) .NE. "x0") then 1427 res@trXMinF = stringtofloat(parameter(63)) 1428 else 1429 res@trXMinF = miniwsus 1430 end if 1431 else 1432 res@trXMinF = xs 1433 end if 1434 if (.not. isvar("xe")) then 1435 if (parameter(65) .NE. "xdim") then 1436 res@trXMaxF = stringtofloat(parameter(65)) 1437 else 1438 res@trXMaxF = maxiwsus 1439 end if 1440 else 1441 res@trXMaxF = xe 1442 end if 1219 1443 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1220 1444 end if … … 1229 1453 res@gsnLeftString = vNam(varn) 1230 1454 res@gsnRightString = unit(varn) 1231 res@trXMinF = miniwu 1232 res@trXMaxF = maxiwu 1455 if (.not. isvar("xs")) then 1456 if (parameter(63) .NE. "x0") then 1457 res@trXMinF = stringtofloat(parameter(63)) 1458 else 1459 res@trXMinF = miniwu 1460 end if 1461 else 1462 res@trXMinF = xs 1463 end if 1464 if (.not. isvar("xe")) then 1465 if (parameter(65) .NE. "xdim") then 1466 res@trXMaxF = stringtofloat(parameter(65)) 1467 else 1468 res@trXMaxF = maxiwu 1469 end if 1470 else 1471 res@trXMaxF = xe 1472 end if 1233 1473 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1234 1474 end if … … 1244 1484 res@gsnLeftString = vNam(varn) 1245 1485 res@gsnRightString = unit(varn) 1246 res@trXMinF = miniwpvp 1247 res@trXMaxF = maxiwpvp 1486 if (.not. isvar("xs")) then 1487 if (parameter(63) .NE. "x0") then 1488 res@trXMinF = stringtofloat(parameter(63)) 1489 else 1490 res@trXMinF = miniwpvp 1491 end if 1492 else 1493 res@trXMinF = xs 1494 end if 1495 if (.not. isvar("xe")) then 1496 if (parameter(65) .NE. "xdim") then 1497 res@trXMaxF = stringtofloat(parameter(65)) 1498 else 1499 res@trXMaxF = maxiwpvp 1500 end if 1501 else 1502 res@trXMaxF = xe 1503 end if 1248 1504 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1249 1505 end if … … 1258 1514 res@gsnLeftString = vNam(varn) 1259 1515 res@gsnRightString = unit(varn) 1260 res@trXMinF = miniwsvs 1261 res@trXMaxF = maxiwsvs 1516 if (.not. isvar("xs")) then 1517 if (parameter(63) .NE. "x0") then 1518 res@trXMinF = stringtofloat(parameter(63)) 1519 else 1520 res@trXMinF = miniwsvs 1521 end if 1522 else 1523 res@trXMinF = xs 1524 end if 1525 if (.not. isvar("xe")) then 1526 if (parameter(65) .NE. "xdim") then 1527 res@trXMaxF = stringtofloat(parameter(65)) 1528 else 1529 res@trXMaxF = maxiwsvs 1530 end if 1531 else 1532 res@trXMaxF = xe 1533 end if 1262 1534 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1263 1535 end if … … 1272 1544 res@gsnLeftString = vNam(varn) 1273 1545 res@gsnRightString = unit(varn) 1274 res@trXMinF = miniwv 1275 res@trXMaxF = maxiwv 1546 if (.not. isvar("xs")) then 1547 if (parameter(63) .NE. "x0") then 1548 res@trXMinF = stringtofloat(parameter(63)) 1549 else 1550 res@trXMinF = miniwv 1551 end if 1552 else 1553 res@trXMinF = xs 1554 end if 1555 if (.not. isvar("xe")) then 1556 if (parameter(65) .NE. "xdim") then 1557 res@trXMaxF = stringtofloat(parameter(65)) 1558 else 1559 res@trXMaxF = maxiwv 1560 end if 1561 else 1562 res@trXMaxF = xe 1563 end if 1276 1564 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1277 1565 end if … … 1287 1575 res@gsnLeftString = vNam(varn) 1288 1576 res@gsnRightString = unit(varn) 1289 res@trXMinF = miniwpptp 1290 res@trXMaxF = maxiwpptp 1577 if (.not. isvar("xs")) then 1578 if (parameter(63) .NE. "x0") then 1579 res@trXMinF = stringtofloat(parameter(63)) 1580 else 1581 res@trXMinF = miniwpptp 1582 end if 1583 else 1584 res@trXMinF = xs 1585 end if 1586 if (.not. isvar("xe")) then 1587 if (parameter(65) .NE. "xdim") then 1588 res@trXMaxF = stringtofloat(parameter(65)) 1589 else 1590 res@trXMaxF = maxiwpptp 1591 end if 1592 else 1593 res@trXMaxF = xe 1594 end if 1291 1595 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1292 1596 end if … … 1301 1605 res@gsnLeftString = vNam(varn) 1302 1606 res@gsnRightString = unit(varn) 1303 res@trXMinF = miniwspts 1304 res@trXMaxF = maxiwspts 1607 if (.not. isvar("xs")) then 1608 if (parameter(63) .NE. "x0") then 1609 res@trXMinF = stringtofloat(parameter(63)) 1610 else 1611 res@trXMinF = miniwspts 1612 end if 1613 else 1614 res@trXMinF = xs 1615 end if 1616 if (.not. isvar("xe")) then 1617 if (parameter(65) .NE. "xdim") then 1618 res@trXMaxF = stringtofloat(parameter(65)) 1619 else 1620 res@trXMaxF = maxiwspts 1621 end if 1622 else 1623 res@trXMaxF = xe 1624 end if 1305 1625 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1306 1626 end if … … 1315 1635 res@gsnLeftString = vNam(varn) 1316 1636 res@gsnRightString = unit(varn) 1317 res@trXMinF = miniwpt 1318 res@trXMaxF = maxiwpt 1637 if (.not. isvar("xs")) then 1638 if (parameter(63) .NE. "x0") then 1639 res@trXMinF = stringtofloat(parameter(63)) 1640 else 1641 res@trXMinF = miniwpt 1642 end if 1643 else 1644 res@trXMinF = xs 1645 end if 1646 if (.not. isvar("xe")) then 1647 if (parameter(65) .NE. "xdim") then 1648 res@trXMaxF = stringtofloat(parameter(65)) 1649 else 1650 res@trXMaxF = maxiwpt 1651 end if 1652 else 1653 res@trXMaxF = xe 1654 end if 1319 1655 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1320 1656 end if … … 1330 1666 res@gsnLeftString = vNam(varn) 1331 1667 res@gsnRightString = unit(varn) 1332 res@trXMinF = miniwsptsBC 1333 res@trXMaxF = maxiwsptsBC 1668 if (.not. isvar("xs")) then 1669 if (parameter(63) .NE. "x0") then 1670 res@trXMinF = stringtofloat(parameter(63)) 1671 else 1672 res@trXMinF = miniwsptsBC 1673 end if 1674 else 1675 res@trXMinF = xs 1676 end if 1677 if (.not. isvar("xe")) then 1678 if (parameter(65) .NE. "xdim") then 1679 res@trXMaxF = stringtofloat(parameter(65)) 1680 else 1681 res@trXMaxF = maxiwsptsBC 1682 end if 1683 else 1684 res@trXMaxF = xe 1685 end if 1334 1686 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1335 1687 end if … … 1344 1696 res@gsnLeftString = vNam(varn) 1345 1697 res@gsnRightString = unit(varn) 1346 res@trXMinF = miniwptBC 1347 res@trXMaxF = maxiwptBC 1698 if (.not. isvar("xs")) then 1699 if (parameter(63) .NE. "x0") then 1700 res@trXMinF = stringtofloat(parameter(63)) 1701 else 1702 res@trXMinF = miniwptBC 1703 end if 1704 else 1705 res@trXMinF = xs 1706 end if 1707 if (.not. isvar("xe")) then 1708 if (parameter(65) .NE. "xdim") then 1709 res@trXMaxF = stringtofloat(parameter(65)) 1710 else 1711 res@trXMaxF = maxiwptBC 1712 end if 1713 else 1714 res@trXMaxF = xe 1715 end if 1348 1716 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1349 1717 end if … … 1359 1727 res@gsnLeftString = vNam(varn) 1360 1728 res@gsnRightString = unit(varn) 1361 res@trXMinF = miniwpvptp 1362 res@trXMaxF = maxiwpvptp 1729 if (.not. isvar("xs")) then 1730 if (parameter(63) .NE. "x0") then 1731 res@trXMinF = stringtofloat(parameter(63)) 1732 else 1733 res@trXMinF = miniwpvptp 1734 end if 1735 else 1736 res@trXMinF = xs 1737 end if 1738 if (.not. isvar("xe")) then 1739 if (parameter(65) .NE. "xdim") then 1740 res@trXMaxF = stringtofloat(parameter(65)) 1741 else 1742 res@trXMaxF = maxiwpvptp 1743 end if 1744 else 1745 res@trXMaxF = xe 1746 end if 1363 1747 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1364 1748 end if … … 1373 1757 res@gsnLeftString = vNam(varn) 1374 1758 res@gsnRightString = unit(varn) 1375 res@trXMinF = miniwsvpts 1376 res@trXMaxF = maxiwsvpts 1759 if (.not. isvar("xs")) then 1760 if (parameter(63) .NE. "x0") then 1761 res@trXMinF = stringtofloat(parameter(63)) 1762 else 1763 res@trXMinF = miniwsvpts 1764 end if 1765 else 1766 res@trXMinF = xs 1767 end if 1768 if (.not. isvar("xe")) then 1769 if (parameter(65) .NE. "xdim") then 1770 res@trXMaxF = stringtofloat(parameter(65)) 1771 else 1772 res@trXMaxF = maxiwsvpts 1773 end if 1774 else 1775 res@trXMaxF = xe 1776 end if 1377 1777 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1378 1778 end if … … 1387 1787 res@gsnLeftString = vNam(varn) 1388 1788 res@gsnRightString = unit(varn) 1389 res@trXMinF = miniwvpt 1390 res@trXMaxF = maxiwvpt 1789 if (.not. isvar("xs")) then 1790 if (parameter(63) .NE. "x0") then 1791 res@trXMinF = stringtofloat(parameter(63)) 1792 else 1793 res@trXMinF = miniwvpt 1794 end if 1795 else 1796 res@trXMinF = xs 1797 end if 1798 if (.not. isvar("xe")) then 1799 if (parameter(65) .NE. "xdim") then 1800 res@trXMaxF = stringtofloat(parameter(65)) 1801 else 1802 res@trXMaxF = maxiwvpt 1803 end if 1804 else 1805 res@trXMaxF = xe 1806 end if 1391 1807 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1392 1808 end if … … 1402 1818 res@gsnLeftString = vNam(varn) 1403 1819 res@gsnRightString = unit(varn) 1404 res@trXMinF = miniwpqp 1405 res@trXMaxF = maxiwpqp 1820 if (.not. isvar("xs")) then 1821 if (parameter(63) .NE. "x0") then 1822 res@trXMinF = stringtofloat(parameter(63)) 1823 else 1824 res@trXMinF = miniwpqp 1825 end if 1826 else 1827 res@trXMinF = xs 1828 end if 1829 if (.not. isvar("xe")) then 1830 if (parameter(65) .NE. "xdim") then 1831 res@trXMaxF = stringtofloat(parameter(65)) 1832 else 1833 res@trXMaxF = maxiwpqp 1834 end if 1835 else 1836 res@trXMaxF = xe 1837 end if 1406 1838 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1407 1839 end if … … 1416 1848 res@gsnLeftString = vNam(varn) 1417 1849 res@gsnRightString = unit(varn) 1418 res@trXMinF = miniwsqs 1419 res@trXMaxF = maxiwsqs 1850 if (.not. isvar("xs")) then 1851 if (parameter(63) .NE. "x0") then 1852 res@trXMinF = stringtofloat(parameter(63)) 1853 else 1854 res@trXMinF = miniwsqs 1855 end if 1856 else 1857 res@trXMinF = xs 1858 end if 1859 if (.not. isvar("xe")) then 1860 if (parameter(65) .NE. "xdim") then 1861 res@trXMaxF = stringtofloat(parameter(65)) 1862 else 1863 res@trXMaxF = maxiwsqs 1864 end if 1865 else 1866 res@trXMaxF = xe 1867 end if 1420 1868 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1421 1869 end if … … 1430 1878 res@gsnLeftString = vNam(varn) 1431 1879 res@gsnRightString = unit(varn) 1432 res@trXMinF = miniwq 1433 res@trXMaxF = maxiwq 1880 if (.not. isvar("xs")) then 1881 if (parameter(63) .NE. "x0") then 1882 res@trXMinF = stringtofloat(parameter(63)) 1883 else 1884 res@trXMinF = miniwq 1885 end if 1886 else 1887 res@trXMinF = xs 1888 end if 1889 if (.not. isvar("xe")) then 1890 if (parameter(65) .NE. "xdim") then 1891 res@trXMaxF = stringtofloat(parameter(65)) 1892 else 1893 res@trXMaxF = maxiwq 1894 end if 1895 else 1896 res@trXMaxF = xe 1897 end if 1434 1898 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1435 1899 end if … … 1445 1909 res@gsnLeftString = vNam(varn) 1446 1910 res@gsnRightString = unit(varn) 1447 res@trXMinF = miniwpqvp 1448 res@trXMaxF = maxiwpqvp 1911 if (.not. isvar("xs")) then 1912 if (parameter(63) .NE. "x0") then 1913 res@trXMinF = stringtofloat(parameter(63)) 1914 else 1915 res@trXMinF = miniwpqvp 1916 end if 1917 else 1918 res@trXMinF = xs 1919 end if 1920 if (.not. isvar("xe")) then 1921 if (parameter(65) .NE. "xdim") then 1922 res@trXMaxF = stringtofloat(parameter(65)) 1923 else 1924 res@trXMaxF = maxiwpqvp 1925 end if 1926 else 1927 res@trXMaxF = xe 1928 end if 1449 1929 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1450 1930 end if … … 1459 1939 res@gsnLeftString = vNam(varn) 1460 1940 res@gsnRightString = unit(varn) 1461 res@trXMinF = miniwsqvs 1462 res@trXMaxF = maxiwsqvs 1941 if (.not. isvar("xs")) then 1942 if (parameter(63) .NE. "x0") then 1943 res@trXMinF = stringtofloat(parameter(63)) 1944 else 1945 res@trXMinF = miniwsqvs 1946 end if 1947 else 1948 res@trXMinF = xs 1949 end if 1950 if (.not. isvar("xe")) then 1951 if (parameter(65) .NE. "xdim") then 1952 res@trXMaxF = stringtofloat(parameter(65)) 1953 else 1954 res@trXMaxF = maxiwsqvs 1955 end if 1956 else 1957 res@trXMaxF = xe 1958 end if 1463 1959 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1464 1960 end if … … 1473 1969 res@gsnLeftString = vNam(varn) 1474 1970 res@gsnRightString = unit(varn) 1475 res@trXMinF = miniwqv 1476 res@trXMaxF = maxiwqv 1971 if (.not. isvar("xs")) then 1972 if (parameter(63) .NE. "x0") then 1973 res@trXMinF = stringtofloat(parameter(63)) 1974 else 1975 res@trXMinF = miniwqv 1976 end if 1977 else 1978 res@trXMinF = xs 1979 end if 1980 if (.not. isvar("xe")) then 1981 if (parameter(65) .NE. "xdim") then 1982 res@trXMaxF = stringtofloat(parameter(65)) 1983 else 1984 res@trXMaxF = maxiwqv 1985 end if 1986 else 1987 res@trXMaxF = xe 1988 end if 1477 1989 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1478 1990 end if … … 1488 2000 res@gsnLeftString = vNam(varn) 1489 2001 res@gsnRightString = unit(varn) 1490 res@trXMinF = miniwpsp 1491 res@trXMaxF = maxiwpsp 2002 if (.not. isvar("xs")) then 2003 if (parameter(63) .NE. "x0") then 2004 res@trXMinF = stringtofloat(parameter(63)) 2005 else 2006 res@trXMinF = miniwpsp 2007 end if 2008 else 2009 res@trXMinF = xs 2010 end if 2011 if (.not. isvar("xe")) then 2012 if (parameter(65) .NE. "xdim") then 2013 res@trXMaxF = stringtofloat(parameter(65)) 2014 else 2015 res@trXMaxF = maxiwpsp 2016 end if 2017 else 2018 res@trXMaxF = xe 2019 end if 1492 2020 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1493 2021 end if … … 1502 2030 res@gsnLeftString = vNam(varn) 1503 2031 res@gsnRightString = unit(varn) 1504 res@trXMinF = miniwsss 1505 res@trXMaxF = maxiwsss 2032 if (.not. isvar("xs")) then 2033 if (parameter(63) .NE. "x0") then 2034 res@trXMinF = stringtofloat(parameter(63)) 2035 else 2036 res@trXMinF = miniwsss 2037 end if 2038 else 2039 res@trXMinF = xs 2040 end if 2041 if (.not. isvar("xe")) then 2042 if (parameter(65) .NE. "xdim") then 2043 res@trXMaxF = stringtofloat(parameter(65)) 2044 else 2045 res@trXMaxF = maxiwsss 2046 end if 2047 else 2048 res@trXMaxF = xe 2049 end if 1506 2050 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1507 2051 end if … … 1516 2060 res@gsnLeftString = vNam(varn) 1517 2061 res@gsnRightString = unit(varn) 1518 res@trXMinF = miniws 1519 res@trXMaxF = maxiws 2062 if (.not. isvar("xs")) then 2063 if (parameter(63) .NE. "x0") then 2064 res@trXMinF = stringtofloat(parameter(63)) 2065 else 2066 res@trXMinF = miniws 2067 end if 2068 else 2069 res@trXMinF = xs 2070 end if 2071 if (.not. isvar("xe")) then 2072 if (parameter(65) .NE. "xdim") then 2073 res@trXMaxF = stringtofloat(parameter(65)) 2074 else 2075 res@trXMaxF = maxiws 2076 end if 2077 else 2078 res@trXMaxF = xe 2079 end if 1520 2080 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1521 2081 end if … … 1531 2091 res@gsnLeftString = vNam(varn) 1532 2092 res@gsnRightString = unit(varn) 1533 res@trXMinF = miniwpsap 1534 res@trXMaxF = maxiwpsap 2093 if (.not. isvar("xs")) then 2094 if (parameter(63) .NE. "x0") then 2095 res@trXMinF = stringtofloat(parameter(63)) 2096 else 2097 res@trXMinF = miniwpsap 2098 end if 2099 else 2100 res@trXMinF = xs 2101 end if 2102 if (.not. isvar("xe")) then 2103 if (parameter(65) .NE. "xdim") then 2104 res@trXMaxF = stringtofloat(parameter(65)) 2105 else 2106 res@trXMaxF = maxiwpsap 2107 end if 2108 else 2109 res@trXMaxF = xe 2110 end if 1535 2111 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1536 2112 end if … … 1545 2121 res@gsnLeftString = vNam(varn) 1546 2122 res@gsnRightString = unit(varn) 1547 res@trXMinF = miniwssas 1548 res@trXMaxF = maxiwssas 2123 if (.not. isvar("xs")) then 2124 if (parameter(63) .NE. "x0") then 2125 res@trXMinF = stringtofloat(parameter(63)) 2126 else 2127 res@trXMinF = miniwssas 2128 end if 2129 else 2130 res@trXMinF = xs 2131 end if 2132 if (.not. isvar("xe")) then 2133 if (parameter(65) .NE. "xdim") then 2134 res@trXMaxF = stringtofloat(parameter(65)) 2135 else 2136 res@trXMaxF = maxiwssas 2137 end if 2138 else 2139 res@trXMaxF = xe 2140 end if 1549 2141 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1550 2142 end if … … 1559 2151 res@gsnLeftString = vNam(varn) 1560 2152 res@gsnRightString = unit(varn) 1561 res@trXMinF = miniwsa1562 res@trXMaxF = maxiwsa1563 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res)1564 end if1565 end if1566 1567 if (vNam(varn) .EQ. "wses") then1568 if (over .EQ. 1) then1569 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 else1575 res@gsnLeftString = vNam(varn)1576 res@gsnRightString = unit(varn)1577 res@trXMinF = min(data(varn,:,:))1578 res@trXMaxF = max(data(varn,:,:))2153 if (.not. isvar("xs")) then 2154 if (parameter(63) .NE. "x0") then 2155 res@trXMinF = stringtofloat(parameter(63)) 2156 else 2157 res@trXMinF = miniwsa 2158 end if 2159 else 2160 res@trXMinF = xs 2161 end if 2162 if (.not. isvar("xe")) then 2163 if (parameter(65) .NE. "xdim") then 2164 res@trXMaxF = stringtofloat(parameter(65)) 2165 else 2166 res@trXMaxF = maxiwsa 2167 end if 2168 else 2169 res@trXMaxF = xe 2170 end if 1579 2171 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1580 2172 end if … … 1590 2182 res@gsnLeftString = vNam(varn) 1591 2183 res@gsnRightString = unit(varn) 1592 res@trXMinF = minius2 1593 res@trXMaxF = maxius2 2184 if (.not. isvar("xs")) then 2185 if (parameter(63) .NE. "x0") then 2186 res@trXMinF = stringtofloat(parameter(63)) 2187 else 2188 res@trXMinF = minius2 2189 end if 2190 else 2191 res@trXMinF = xs 2192 end if 2193 if (.not. isvar("xe")) then 2194 if (parameter(65) .NE. "xdim") then 2195 res@trXMaxF = stringtofloat(parameter(65)) 2196 else 2197 res@trXMaxF = maxius2 2198 end if 2199 else 2200 res@trXMaxF = xe 2201 end if 1594 2202 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1595 2203 end if … … 1604 2212 res@gsnLeftString = vNam(varn) 1605 2213 res@gsnRightString = unit(varn) 1606 res@trXMinF = minivs2 1607 res@trXMaxF = maxivs2 2214 if (.not. isvar("xs")) then 2215 if (parameter(63) .NE. "x0") then 2216 res@trXMinF = stringtofloat(parameter(63)) 2217 else 2218 res@trXMinF = minivs2 2219 end if 2220 else 2221 res@trXMinF = xs 2222 end if 2223 if (.not. isvar("xe")) then 2224 if (parameter(65) .NE. "xdim") then 2225 res@trXMaxF = stringtofloat(parameter(65)) 2226 else 2227 res@trXMaxF = maxivs2 2228 end if 2229 else 2230 res@trXMaxF = xe 2231 end if 1608 2232 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1609 2233 end if … … 1618 2242 res@gsnLeftString = vNam(varn) 1619 2243 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 2244 if (.not. isvar("xs")) then 2245 if (parameter(63) .NE. "x0") then 2246 res@trXMinF = stringtofloat(parameter(63)) 2247 else 2248 res@trXMinF = miniws2 2249 end if 2250 else 2251 res@trXMinF = xs 2252 end if 2253 if (.not. isvar("xe")) then 2254 if (parameter(65) .NE. "xdim") then 2255 res@trXMaxF = stringtofloat(parameter(65)) 2256 else 2257 res@trXMaxF = maxiws2 2258 end if 2259 else 2260 res@trXMaxF = xe 2261 end if 1712 2262 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1713 2263 end if … … 1723 2273 res@gsnLeftString = vNam(varn) 1724 2274 res@gsnRightString = unit(varn) 1725 res@trXMinF = miniwsususodz 1726 res@trXMaxF = maxiwsususodz 2275 if (.not. isvar("xs")) then 2276 if (parameter(63) .NE. "x0") then 2277 res@trXMinF = stringtofloat(parameter(63)) 2278 else 2279 res@trXMinF = miniwsususodz 2280 end if 2281 else 2282 res@trXMinF = xs 2283 end if 2284 if (.not. isvar("xe")) then 2285 if (parameter(65) .NE. "xdim") then 2286 res@trXMaxF = stringtofloat(parameter(65)) 2287 else 2288 res@trXMaxF = maxiwsususodz 2289 end if 2290 else 2291 res@trXMaxF = xe 2292 end if 1727 2293 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1728 2294 end if … … 1737 2303 res@gsnLeftString = vNam(varn) 1738 2304 res@gsnRightString = unit(varn) 1739 res@trXMinF = miniwspsodz 1740 res@trXMaxF = maxiwspsodz 2305 if (.not. isvar("xs")) then 2306 if (parameter(63) .NE. "x0") then 2307 res@trXMinF = stringtofloat(parameter(63)) 2308 else 2309 res@trXMinF = miniwspsodz 2310 end if 2311 else 2312 res@trXMinF = xs 2313 end if 2314 if (.not. isvar("xe")) then 2315 if (parameter(65) .NE. "xdim") then 2316 res@trXMaxF = stringtofloat(parameter(65)) 2317 else 2318 res@trXMaxF = maxiwspsodz 2319 end if 2320 else 2321 res@trXMaxF = xe 2322 end if 1741 2323 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) 1742 2324 end if … … 1751 2333 res@gsnLeftString = vNam(varn) 1752 2334 res@gsnRightString = unit(varn) 1753 res@trXMinF = miniwpeodz 1754 res@trXMaxF = maxiwpeodz 2335 if (.not. isvar("xs")) then 2336 if (parameter(63) .NE. "x0") then 2337 res@trXMinF = stringtofloat(parameter(63)) 2338 else 2339 res@trXMinF = miniwpeodz 2340 end if 2341 else 2342 res@trXMinF = xs 2343 end if 2344 if (.not. isvar("xe")) then 2345 if (parameter(65) .NE. "xdim") then 2346 res@trXMaxF = stringtofloat(parameter(65)) 2347 else 2348 res@trXMaxF = maxiwpeodz 2349 end if 2350 else 2351 res@trXMaxF = xe 2352 end if 1755 2353 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) 1756 2354 end if … … 1763 2361 delete(temp) 1764 2362 end do 1765 2363 2364 if (count_var .EQ. 0) then 2365 print(" ") 2366 print("Please select a variable 'var=' or use the default value") 2367 print(" ") 2368 print("Your selection '"+var+"' does not exist on the input file") 2369 print(" ") 2370 exit 2371 end if 2372 1766 2373 if (over .EQ. 1 ) then 1767 2374 … … 1836 2443 z = f->$vNam(varn+1)$ 1837 2444 1838 ;res@gsnLeftString = vNam(varn)2445 res@gsnLeftString = vNam(varn) 1839 2446 res@gsnRightString = unit(varn) 1840 2447 res@trYMinF = min_z 1841 2448 res@trYMaxF = max_z 1842 res@trXMinF = min(data(varn,:,:)) 1843 res@trXMaxF = max(data(varn,:,:)) 2449 if (.not. isvar("xs")) then 2450 if (parameter(63) .NE. "x0") then 2451 res@trXMinF = stringtofloat(parameter(63)) 2452 else 2453 res@trXMinF = min(data(varn,:,:)) 2454 end if 2455 else 2456 res@trXMinF = xs 2457 end if 2458 if (.not. isvar("xe")) then 2459 if (parameter(65) .NE. "xdim") then 2460 res@trXMaxF = stringtofloat(parameter(65)) 2461 else 2462 res@trXMaxF = max(data(varn,:,:)) 2463 end if 2464 else 2465 res@trXMaxF = xe 2466 end if 1844 2467 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 1845 2468 … … 1848 2471 res@gsnLeftString = "u, v and w" 1849 2472 res@gsnRightString = unit(varn) 1850 res@trXMinF = min((/miniu,miniv,miniw/)) 1851 res@trXMaxF = max((/maxiu,maxiv,maxiw/)) 2473 if (.not. isvar("xs")) then 2474 if (parameter(63) .NE. "x0") then 2475 res@trXMinF = stringtofloat(parameter(63)) 2476 else 2477 res@trXMinF = min((/miniu,miniv,miniw/)) 2478 end if 2479 else 2480 res@trXMinF = xs 2481 end if 2482 if (.not. isvar("xe")) then 2483 if (parameter(65) .NE. "xdim") then 2484 res@trXMaxF = stringtofloat(parameter(65)) 2485 else 2486 res@trXMaxF = max((/maxiu,maxiv,maxiw/)) 2487 end if 2488 else 2489 res@trXMaxF = xe 2490 end if 1852 2491 if (vNam(varn) .EQ. "w") then 1853 2492 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 1885 2524 res@gsnLeftString = "pt, vpt and lpt" 1886 2525 res@gsnRightString = unit(varn) 1887 res@trXMinF = min((/minipt,minivpt,minilpt/)) 1888 res@trXMaxF = max((/maxipt,maxivpt,maxilpt/)) 2526 if (.not. isvar("xs")) then 2527 if (parameter(63) .NE. "x0") then 2528 res@trXMinF = stringtofloat(parameter(63)) 2529 else 2530 res@trXMinF = min((/minipt,minivpt,minilpt/)) 2531 end if 2532 else 2533 res@trXMinF = xs 2534 end if 2535 if (.not. isvar("xe")) then 2536 if (parameter(65) .NE. "xdim") then 2537 res@trXMaxF = stringtofloat(parameter(65)) 2538 else 2539 res@trXMaxF = max((/maxipt,maxivpt,maxilpt/)) 2540 end if 2541 else 2542 res@trXMaxF = xe 2543 end if 1889 2544 1890 2545 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) … … 1918 2573 res@gsnLeftString = "q, qv and ql" 1919 2574 res@gsnRightString = unit(varn) 1920 res@trXMinF = min((/miniq,miniqv,miniql/)) 1921 res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) 2575 if (.not. isvar("xs")) then 2576 if (parameter(63) .NE. "x0") then 2577 res@trXMinF = stringtofloat(parameter(63)) 2578 else 2579 res@trXMinF = min((/miniq,miniqv,miniql/)) 2580 end if 2581 else 2582 res@trXMinF = xs 2583 end if 2584 if (.not. isvar("xe")) then 2585 if (parameter(65) .NE. "xdim") then 2586 res@trXMaxF = stringtofloat(parameter(65)) 2587 else 2588 res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) 2589 end if 2590 else 2591 res@trXMaxF = xe 2592 end if 1922 2593 1923 2594 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) … … 1946 2617 continue 1947 2618 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 2619 end if 2620 1958 2621 if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then 1959 2622 if (e .EQ. 0) then 1960 2623 res@gsnLeftString = "e and es" 1961 2624 res@gsnRightString = unit(varn) 1962 res@trXMinF = min((/minie,minies/)) 1963 res@trXMaxF = max((/maxie,maxies/)) 2625 if (.not. isvar("xs")) then 2626 if (parameter(63) .NE. "x0") then 2627 res@trXMinF = stringtofloat(parameter(63)) 2628 else 2629 res@trXMinF = min((/minie,minies/)) 2630 end if 2631 else 2632 res@trXMinF = xs 2633 end if 2634 if (.not. isvar("xe")) then 2635 if (parameter(65) .NE. "xdim") then 2636 res@trXMaxF = stringtofloat(parameter(65)) 2637 else 2638 res@trXMaxF = max((/maxie,maxies/)) 2639 end if 2640 else 2641 res@trXMaxF = xe 2642 end if 1964 2643 1965 2644 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) … … 1993 2672 res@gsnLeftString = "km and kh" 1994 2673 res@gsnRightString = unit(varn) 1995 res@trXMinF = min((/minikm,minikh/)) 1996 res@trXMaxF = max((/maxikm,maxikh/)) 2674 if (.not. isvar("xs")) then 2675 if (parameter(63) .NE. "x0") then 2676 res@trXMinF = stringtofloat(parameter(63)) 2677 else 2678 res@trXMinF = min((/minikm,minikh/)) 2679 end if 2680 else 2681 res@trXMinF = xs 2682 end if 2683 if (.not. isvar("xe")) then 2684 if (parameter(65) .NE. "xdim") then 2685 res@trXMaxF = stringtofloat(parameter(65)) 2686 else 2687 res@trXMaxF = max((/maxikm,maxikh/)) 2688 end if 2689 else 2690 res@trXMaxF = xe 2691 end if 1997 2692 1998 2693 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_u,res) … … 2022 2717 end if 2023 2718 end if 2024 if (vNam(varn) .EQ. "l") then 2025 plot(n) = plot_l 2026 end if 2719 2027 2720 if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then 2028 2721 if (wpup .EQ. 0) then 2029 2722 res@gsnLeftString = "wpup, wsus and wu" 2030 2723 res@gsnRightString = unit(varn) 2031 res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) 2032 res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) 2724 if (.not. isvar("xs")) then 2725 if (parameter(63) .NE. "x0") then 2726 res@trXMinF = stringtofloat(parameter(63)) 2727 else 2728 res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) 2729 end if 2730 else 2731 res@trXMinF = xs 2732 end if 2733 if (.not. isvar("xe")) then 2734 if (parameter(65) .NE. "xdim") then 2735 res@trXMaxF = stringtofloat(parameter(65)) 2736 else 2737 res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) 2738 end if 2739 else 2740 res@trXMaxF = xe 2741 end if 2033 2742 2034 2743 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2062 2771 res@gsnLeftString = "wpvp, wsus and wv" 2063 2772 res@gsnRightString = unit(varn) 2064 res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) 2065 res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) 2773 if (.not. isvar("xs")) then 2774 if (parameter(63) .NE. "x0") then 2775 res@trXMinF = stringtofloat(parameter(63)) 2776 else 2777 res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) 2778 end if 2779 else 2780 res@trXMinF = xs 2781 end if 2782 if (.not. isvar("xe")) then 2783 if (parameter(65) .NE. "xdim") then 2784 res@trXMaxF = stringtofloat(parameter(65)) 2785 else 2786 res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) 2787 end if 2788 else 2789 res@trXMaxF = xe 2790 end if 2066 2791 2067 2792 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2095 2820 res@gsnLeftString = "wpptp, wspts and wv" 2096 2821 res@gsnRightString = unit(varn) 2097 res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) 2098 res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) 2822 if (.not. isvar("xs")) then 2823 if (parameter(63) .NE. "x0") then 2824 res@trXMinF = stringtofloat(parameter(63)) 2825 else 2826 res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) 2827 end if 2828 else 2829 res@trXMinF = xs 2830 end if 2831 if (.not. isvar("xe")) then 2832 if (parameter(65) .NE. "xdim") then 2833 res@trXMaxF = stringtofloat(parameter(65)) 2834 else 2835 res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) 2836 end if 2837 else 2838 res@trXMaxF = xe 2839 end if 2099 2840 2100 2841 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2128 2869 res@gsnLeftString = "wsptsBC and wptBC" 2129 2870 res@gsnRightString = unit(varn) 2130 res@trXMinF = min((/miniwsptsBC,miniwptBC/)) 2131 res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) 2871 if (.not. isvar("xs")) then 2872 if (parameter(63) .NE. "x0") then 2873 res@trXMinF = stringtofloat(parameter(63)) 2874 else 2875 res@trXMinF = min((/miniwsptsBC,miniwptBC/)) 2876 end if 2877 else 2878 res@trXMinF = xs 2879 end if 2880 if (.not. isvar("xe")) then 2881 if (parameter(65) .NE. "xdim") then 2882 res@trXMaxF = stringtofloat(parameter(65)) 2883 else 2884 res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) 2885 end if 2886 else 2887 res@trXMaxF = xe 2888 end if 2132 2889 2133 2890 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2161 2918 res@gsnLeftString = "wpvptp, wsvpts and wv" 2162 2919 res@gsnRightString = unit(varn) 2163 res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) 2164 res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) 2920 if (.not. isvar("xs")) then 2921 if (parameter(63) .NE. "x0") then 2922 res@trXMinF = stringtofloat(parameter(63)) 2923 else 2924 res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) 2925 end if 2926 else 2927 res@trXMinF = xs 2928 end if 2929 if (.not. isvar("xe")) then 2930 if (parameter(65) .NE. "xdim") then 2931 res@trXMaxF = stringtofloat(parameter(65)) 2932 else 2933 res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) 2934 end if 2935 else 2936 res@trXMaxF = xe 2937 end if 2165 2938 2166 2939 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2194 2967 res@gsnLeftString = "wpqp, wsqs and wq" 2195 2968 res@gsnRightString = unit(varn) 2196 res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) 2197 res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) 2969 if (.not. isvar("xs")) then 2970 if (parameter(63) .NE. "x0") then 2971 res@trXMinF = stringtofloat(parameter(63)) 2972 else 2973 res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) 2974 end if 2975 else 2976 res@trXMinF = xs 2977 end if 2978 if (.not. isvar("xe")) then 2979 if (parameter(65) .NE. "xdim") then 2980 res@trXMaxF = stringtofloat(parameter(65)) 2981 else 2982 res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) 2983 end if 2984 else 2985 res@trXMaxF = xe 2986 end if 2198 2987 2199 2988 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2227 3016 res@gsnLeftString = "wpqvp, wsqvs and wqv" 2228 3017 res@gsnRightString = unit(varn) 2229 res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) 2230 res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) 3018 if (.not. isvar("xs")) then 3019 if (parameter(63) .NE. "x0") then 3020 res@trXMinF = stringtofloat(parameter(63)) 3021 else 3022 res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) 3023 end if 3024 else 3025 res@trXMinF = xs 3026 end if 3027 if (.not. isvar("xe")) then 3028 if (parameter(65) .NE. "xdim") then 3029 res@trXMaxF = stringtofloat(parameter(65)) 3030 else 3031 res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) 3032 end if 3033 else 3034 res@trXMaxF = xe 3035 end if 2231 3036 2232 3037 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2260 3065 res@gsnLeftString = "wpsp, wsss and ws" 2261 3066 res@gsnRightString = unit(varn) 2262 res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) 2263 res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) 3067 if (.not. isvar("xs")) then 3068 if (parameter(63) .NE. "x0") then 3069 res@trXMinF = stringtofloat(parameter(63)) 3070 else 3071 res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) 3072 end if 3073 else 3074 res@trXMinF = xs 3075 end if 3076 if (.not. isvar("xe")) then 3077 if (parameter(65) .NE. "xdim") then 3078 res@trXMaxF = stringtofloat(parameter(65)) 3079 else 3080 res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) 3081 end if 3082 else 3083 res@trXMaxF = xe 3084 end if 2264 3085 2265 3086 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2293 3114 res@gsnLeftString = "wpsap, wssas and wsa" 2294 3115 res@gsnRightString = unit(varn) 2295 res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) 2296 res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) 3116 if (.not. isvar("xs")) then 3117 if (parameter(63) .NE. "x0") then 3118 res@trXMinF = stringtofloat(parameter(63)) 3119 else 3120 res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) 3121 end if 3122 else 3123 res@trXMinF = xs 3124 end if 3125 if (.not. isvar("xe")) then 3126 if (parameter(65) .NE. "xdim") then 3127 res@trXMaxF = stringtofloat(parameter(65)) 3128 else 3129 res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) 3130 end if 3131 else 3132 res@trXMaxF = xe 3133 end if 2297 3134 2298 3135 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2322 3159 end if 2323 3160 end if 2324 if (vNam(varn) .EQ. "wses") then 2325 plot(n) = plot_wses 2326 end if 3161 2327 3162 if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then 2328 3163 if (us2 .EQ. 0) then 2329 3164 res@gsnLeftString = "us2, vs2 and ws2" 2330 3165 res@gsnRightString = unit(varn) 2331 res@trXMinF = min((/minius2,minivs2,miniws2/)) 2332 res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) 3166 if (.not. isvar("xs")) then 3167 if (parameter(63) .NE. "x0") then 3168 res@trXMinF = stringtofloat(parameter(63)) 3169 else 3170 res@trXMinF = min((/minius2,minivs2,miniws2/)) 3171 end if 3172 else 3173 res@trXMinF = xs 3174 end if 3175 if (.not. isvar("xe")) then 3176 if (parameter(65) .NE. "xdim") then 3177 res@trXMaxF = stringtofloat(parameter(65)) 3178 else 3179 res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) 3180 end if 3181 else 3182 res@trXMaxF = xe 3183 end if 2333 3184 2334 3185 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2358 3209 end if 2359 3210 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 3211 2375 3212 if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then 2376 3213 if (wsususodz .EQ. 0) then 2377 3214 res@gsnLeftString = "wsususodz, wspsodz and ws2" 2378 3215 res@gsnRightString = unit(varn) 2379 res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) 2380 res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) 3216 if (.not. isvar("xs")) then 3217 if (parameter(63) .NE. "x0") then 3218 res@trXMinF = stringtofloat(parameter(63)) 3219 else 3220 res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) 3221 end if 3222 else 3223 res@trXMinF = xs 3224 end if 3225 if (.not. isvar("xe")) then 3226 if (parameter(65) .NE. "xdim") then 3227 res@trXMaxF = stringtofloat(parameter(65)) 3228 else 3229 res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) 3230 end if 3231 else 3232 res@trXMaxF = xe 3233 end if 2381 3234 2382 3235 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z_w,res) … … 2449 3302 res@gsnLeftString = c_var 2450 3303 res@gsnRightString = unit(varn) 2451 res@trXMinF = min(mini) 2452 res@trXMaxF = max(maxi) 3304 if (.not. isvar("xs")) then 3305 if (parameter(63) .NE. "x0") then 3306 res@trXMinF = stringtofloat(parameter(63)) 3307 else 3308 res@trXMinF = min(mini) 3309 end if 3310 else 3311 res@trXMinF = xs 3312 end if 3313 if (.not. isvar("xe")) then 3314 if (parameter(65) .NE. "xdim") then 3315 res@trXMaxF = stringtofloat(parameter(65)) 3316 else 3317 res@trXMaxF = max(maxi) 3318 end if 3319 else 3320 res@trXMaxF = xe 3321 end if 2453 3322 end if 2454 3323 end if … … 2496 3365 ; merge plots onto one page 2497 3366 ; *************************************************** 2498 print(n) 3367 2499 3368 do m=0,n-1 2500 3369 plot_(m)=plot(n-1-m)
Note: See TracChangeset
for help on using the changeset viewer.
