Changeset 250
- Timestamp:
- Mar 2, 2009 12:18:17 PM (16 years ago)
- Location:
- palm/trunk/SCRIPTS/NCL
- Files:
-
- 3 edited
-
profiles.ncl (modified) (19 diffs)
-
spectra.ncl (modified) (1 diff)
-
timeseries.ncl (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/profiles.ncl
r236 r250 441 441 else 442 442 if(vNam(varn) .EQ. "vpt" .OR. isStrSubset(vNam(varn), "vpt_"))then 443 z_u = f ->$vNam(varn+1)$443 z_u = f_att->$vNam(varn+1)$ 444 444 break 445 445 else … … 812 812 813 813 delta_t=t_all(nt-1)/nt 814 815 816 print(nt)817 print(t_all)818 819 820 ;exit821 814 822 815 ; **************************************************** … … 946 939 res@pmLegendOrthogonalPosF = -1.0 947 940 res@pmLegendWidthF = 0.12 948 res@pmLegendHeightF = 0.0 4*(end_time_step-start_time_step+1)941 res@pmLegendHeightF = 0.05*(end_time_step-start_time_step+1) 949 942 res@lgLabelFontHeightF = font_size 950 943 res@lgTitleString = "Time [h]" … … 957 950 res@tiXAxisString = " " 958 951 if ( black .eq. 0 ) then 959 res@xyLineColors = ispan(2,237,235/np)952 res@xyLineColors = -(ispan(-237,-2,235/np)) 960 953 end if 961 954 if (norm_z .EQ. 1)then … … 972 965 res@xyMonoDashPattern = True 973 966 else 967 res@xyMonoDashPattern = False 974 968 if (no_files .GT. 1) 975 969 res@xyMonoDashPattern = True … … 979 973 end if 980 974 end if 981 975 982 976 res@tmXBMinorPerMajor = 4 983 977 res@tmYLMinorPerMajor = 4 … … 1404 1398 n_o=0 1405 1399 count_var=0 1400 1406 1401 res@xyDashPattern = 1*nof 1407 1402 1408 1403 do varn = 0,dim-1 1409 1404 … … 1579 1574 end if 1580 1575 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 1581 end if 1582 1576 end if 1577 1583 1578 if (vNam(varn) .EQ. "u") then 1584 1579 miniu=min(data(varn,:,:)) … … 1608 1603 maxiv=max(data(varn,:,:)) 1609 1604 if (over .EQ. 1) then 1605 res@xyMonoDashPattern = True 1610 1606 res@xyDashPattern = 1 1611 1607 plot_v = gsn_csm_xy(wks,data(varn,:,:),z,res) … … 1631 1627 maxiw=max(data(varn,:,:)) 1632 1628 if (over .EQ. 1) then 1633 res@xyDashPattern 1629 res@xyDashPattern = 2 1634 1630 plot_w = gsn_csm_xy(wks,data(varn,:,:),z,res) 1635 1631 else … … 2802 2798 if (prof3d .EQ. 0) then 2803 2799 if (log_z .EQ. 1) then 2804 z = f ->$vNam(varn+1)$(1:dimz-1)2805 else 2806 z = f ->$vNam(varn+1)$2800 z = f_att->$vNam(varn+1)$(1:dimz-1) 2801 else 2802 z = f_att->$vNam(varn+1)$ 2807 2803 end if 2808 2804 else … … 2835 2831 res@trYMinF = min_z 2836 2832 res@trYMaxF = max_z 2833 res@xyDashPattern = 0 2834 2837 2835 if (xs .EQ. -1) then 2838 2836 res@trXMinF = min(data(varn,:,:)) … … 2862 2860 res@trXMaxF = xe 2863 2861 end if 2862 if (vNam(varn) .EQ. "v")then 2863 res@xyDashPattern = 1 2864 end if 2865 if (vNam(varn) .EQ. "w")then 2866 res@xyDashPattern = 2 2867 end if 2864 2868 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2865 2869 … … 2872 2876 lgres@lgLabelFont = "helvetica" 2873 2877 lgres@lgLabelFontHeightF = font_size_legend 2874 lgres@vpWidthF = 0. 122875 lgres@vpHeightF = 0.1 2878 lgres@vpWidthF = 0.06 2879 lgres@vpHeightF = 0.12 2876 2880 lgres@lgDashIndexes = (/0,1,2/) 2877 2881 lbid = gsn_create_legend(wks,3,(/"u","v","w"/),lgres) … … 2906 2910 res@trXMaxF = xe 2907 2911 end if 2908 2912 if (vNam(varn) .EQ. "vpt")then 2913 res@xyDashPattern = 1 2914 end if 2915 if (vNam(varn) .EQ. "lpt")then 2916 res@xyDashPattern = 2 2917 end if 2918 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2919 2920 ; *************************************************** 2921 ; legend for combined plot 2922 ; *************************************************** 2923 2924 lgres = True 2925 lgMonoDashIndex = False 2926 lgres@lgLabelFont = "helvetica" 2927 lgres@lgLabelFontHeightF = font_size_legend 2928 lgres@vpWidthF = 0.07 2929 lgres@vpHeightF = 0.12 2930 lgres@lgDashIndexes = (/0,1,2/) 2931 lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres) 2932 2933 amres = True 2934 amres@amParallelPosF = 0.65 2935 amres@amOrthogonalPosF = -0.2 2936 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2937 overlay(plot(n),plot_pt) 2938 pt=1 2939 else 2940 if (prof3d .EQ. 0)then 2941 varn=varn+1 2942 end if 2943 continue 2944 end if 2945 end if 2946 if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then 2947 if (q .EQ. 0) then 2948 res@gsnLeftString = "q, qv and ql" 2949 res@tiXAxisString = "["+unit(varn)+"]" 2950 res@gsnRightString = " " 2951 if (xs .EQ. -1) then 2952 res@trXMinF = min((/miniq,miniqv,miniql/)) 2953 else 2954 res@trXMinF = xs 2955 end if 2956 if (xe .EQ. -1) then 2957 res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) 2958 else 2959 res@trXMaxF = xe 2960 end if 2961 if (vNam(varn) .EQ. "qv")then 2962 res@xyDashPattern = 1 2963 end if 2964 if (vNam(varn) .EQ. "ql")then 2965 res@xyDashPattern = 2 2966 end if 2967 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2968 2969 ; *************************************************** 2970 ; legend for combined plot 2971 ; *************************************************** 2972 2973 lgres = True 2974 lgMonoDashIndex = False 2975 lgres@lgLabelFont = "helvetica" 2976 lgres@lgLabelFontHeightF = font_size_legend 2977 lgres@vpWidthF = 0.06 2978 lgres@vpHeightF = 0.12 2979 lgres@lgDashIndexes = (/0,1,2/) 2980 lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres) 2981 2982 amres = True 2983 amres@amParallelPosF = 0.65 2984 amres@amOrthogonalPosF = -0.2 2985 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2986 overlay(plot(n),plot_q) 2987 q=1 2988 else 2989 if (prof3d .EQ. 0)then 2990 varn=varn+1 2991 end if 2992 continue 2993 end if 2994 end if 2995 2996 if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then 2997 if (e .EQ. 0) then 2998 res@gsnLeftString = "e and es" 2999 res@tiXAxisString = "["+unit(varn)+"]" 3000 res@gsnRightString = " " 3001 if (xs .EQ. -1) then 3002 res@trXMinF = min((/minie,minies/)) 3003 else 3004 res@trXMinF = xs 3005 end if 3006 if (xe .EQ. -1) then 3007 res@trXMaxF = max((/maxie,maxies/)) 3008 else 3009 res@trXMaxF = xe 3010 end if 3011 if (vNam(varn) .EQ. "es")then 3012 res@xyDashPattern = 1 3013 end if 3014 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3015 3016 ; *************************************************** 3017 ; legend for combined plot 3018 ; *************************************************** 3019 3020 lgres = True 3021 lgMonoDashIndex = False 3022 lgres@lgLabelFont = "helvetica" 3023 lgres@lgLabelFontHeightF = font_size_legend 3024 lgres@vpWidthF = 0.06 3025 lgres@vpHeightF = 0.08 3026 lgres@lgDashIndexes = (/0,1,2/) 3027 lbid = gsn_create_legend(wks,2,(/"e","es"/),lgres) 3028 3029 amres = True 3030 amres@amParallelPosF = 0.65 3031 amres@amOrthogonalPosF = -0.2 3032 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3033 overlay(plot(n),plot_e) 3034 e=1 3035 else 3036 if (prof3d .EQ. 0)then 3037 varn=varn+1 3038 end if 3039 continue 3040 end if 3041 end if 3042 if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then 3043 if (km .EQ. 0) then 3044 res@gsnLeftString = "km and kh" 3045 res@tiXAxisString = "["+unit(varn)+"]" 3046 res@gsnRightString = " " 3047 if (xs .EQ. -1) then 3048 res@trXMinF = min((/minikm,minikh/)) 3049 else 3050 res@trXMinF = xs 3051 end if 3052 if (xe .EQ. -1) then 3053 res@trXMaxF = max((/maxikm,maxikh/)) 3054 else 3055 res@trXMaxF = xe 3056 end if 3057 if (vNam(varn) .EQ. "kh")then 3058 res@xyDashPattern = 1 3059 end if 3060 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3061 3062 ; *************************************************** 3063 ; legend for combined plot 3064 ; *************************************************** 3065 3066 lgres = True 3067 lgMonoDashIndex = False 3068 lgres@lgLabelFont = "helvetica" 3069 lgres@lgLabelFontHeightF = font_size_legend 3070 lgres@vpWidthF = 0.06 3071 lgres@vpHeightF = 0.08 3072 lgres@lgDashIndexes = (/0,1,2/) 3073 lbid = gsn_create_legend(wks,2,(/"km","kh"/),lgres) 3074 3075 amres = True 3076 amres@amParallelPosF = 0.65 3077 amres@amOrthogonalPosF = -0.2 3078 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3079 overlay(plot(n),plot_km) 3080 km=1 3081 else 3082 if (prof3d .EQ. 0)then 3083 varn=varn+1 3084 end if 3085 continue 3086 end if 3087 end if 3088 3089 if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then 3090 if (wpup .EQ. 0) then 3091 res@gsnLeftString = "wpup, wsus and wu" 3092 res@tiXAxisString = "["+unit(varn)+"]" 3093 res@gsnRightString = " " 3094 if (xs .EQ. -1) then 3095 res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) 3096 else 3097 res@trXMinF = xs 3098 end if 3099 if (xe .EQ. -1) then 3100 res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) 3101 else 3102 res@trXMaxF = xe 3103 end if 3104 if (vNam(varn) .EQ. "wsus")then 3105 res@xyDashPattern = 1 3106 end if 3107 if (vNam(varn) .EQ. "wu")then 3108 res@xyDashPattern = 2 3109 end if 3110 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3111 3112 ; *************************************************** 3113 ; legend for combined plot 3114 ; *************************************************** 3115 3116 lgres = True 3117 lgMonoDashIndex = False 3118 lgres@lgLabelFont = "helvetica" 3119 lgres@lgLabelFontHeightF = font_size_legend 3120 lgres@vpWidthF = 0.08 3121 lgres@vpHeightF = 0.12 3122 lgres@lgDashIndexes = (/0,1,2/) 3123 lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres) 3124 3125 amres = True 3126 amres@amParallelPosF = 0.65 3127 amres@amOrthogonalPosF = -0.2 3128 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3129 overlay(plot(n),plot_wpup) 3130 wpup=1 3131 else 3132 if (prof3d .EQ. 0)then 3133 varn=varn+1 3134 end if 3135 continue 3136 end if 3137 end if 3138 if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then 3139 if (wpvp .EQ. 0) then 3140 res@gsnLeftString = "wpvp, wsvs and wv" 3141 res@tiXAxisString = "["+unit(varn)+"]" 3142 res@gsnRightString = " " 3143 if (xs .EQ. -1) then 3144 res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) 3145 else 3146 res@trXMinF = xs 3147 end if 3148 if (xe .EQ. -1) then 3149 res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) 3150 else 3151 res@trXMaxF = xe 3152 end if 3153 if (vNam(varn) .EQ. "wsvs")then 3154 res@xyDashPattern = 1 3155 end if 3156 if (vNam(varn) .EQ. "wv")then 3157 res@xyDashPattern = 2 3158 end if 3159 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3160 3161 ; *************************************************** 3162 ; legend for combined plot 3163 ; *************************************************** 3164 3165 lgres = True 3166 lgMonoDashIndex = False 3167 lgres@lgLabelFont = "helvetica" 3168 lgres@lgLabelFontHeightF = font_size_legend 3169 lgres@vpWidthF = 0.08 3170 lgres@vpHeightF = 0.12 3171 lgres@lgDashIndexes = (/0,1,2/) 3172 lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres) 3173 3174 amres = True 3175 amres@amParallelPosF = 0.65 3176 amres@amOrthogonalPosF = -0.2 3177 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3178 overlay(plot(n),plot_wpvp) 3179 wpvp=1 3180 else 3181 if (prof3d .EQ. 0)then 3182 varn=varn+1 3183 end if 3184 continue 3185 end if 3186 end if 3187 if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then 3188 if (wpptp .EQ. 0) then 3189 res@gsnLeftString = "wpptp, wspts and wpt" 3190 res@tiXAxisString = "["+unit(varn)+"]" 3191 res@gsnRightString = " " 3192 if (xs .EQ. -1) then 3193 res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) 3194 else 3195 res@trXMinF = xs 3196 end if 3197 if (xe .EQ. -1) then 3198 res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) 3199 else 3200 res@trXMaxF = xe 3201 end if 3202 if (vNam(varn) .EQ. "wspts")then 3203 res@xyDashPattern = 1 3204 end if 3205 if (vNam(varn) .EQ. "wpt")then 3206 res@xyDashPattern = 2 3207 end if 3208 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3209 3210 ; *************************************************** 3211 ; legend for combined plot 3212 ; *************************************************** 3213 3214 lgres = True 3215 lgMonoDashIndex = False 3216 lgres@lgLabelFont = "helvetica" 3217 lgres@lgLabelFontHeightF = font_size_legend 3218 lgres@vpWidthF = 0.09 3219 lgres@vpHeightF = 0.12 3220 lgres@lgDashIndexes = (/0,1,2/) 3221 lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres) 3222 3223 amres = True 3224 amres@amParallelPosF = 0.65 3225 amres@amOrthogonalPosF = -0.2 3226 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3227 overlay(plot(n),plot_wpptp) 3228 wpptp=1 3229 else 3230 if (prof3d .EQ. 0)then 3231 varn=varn+1 3232 end if 3233 continue 3234 end if 3235 end if 3236 if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then 3237 if (wsptsBC .EQ. 0) then 3238 res@gsnLeftString = "wsptsBC and wptBC" 3239 res@tiXAxisString = "["+unit(varn)+"]" 3240 res@gsnRightString = " " 3241 if (xs .EQ. -1) then 3242 res@trXMinF = min((/miniwsptsBC,miniwptBC/)) 3243 else 3244 res@trXMinF = xs 3245 end if 3246 if (xe .EQ. -1) then 3247 res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) 3248 else 3249 res@trXMaxF = xe 3250 end if 3251 if (vNam(varn) .EQ. "wptBC")then 3252 res@xyDashPattern = 1 3253 end if 3254 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3255 3256 ; *************************************************** 3257 ; legend for combined plot 3258 ; *************************************************** 3259 3260 lgres = True 3261 lgMonoDashIndex = False 3262 lgres@lgLabelFont = "helvetica" 3263 lgres@lgLabelFontHeightF = font_size_legend 3264 lgres@vpWidthF = 0.1 3265 lgres@vpHeightF = 0.12 3266 lgres@lgDashIndexes = (/0,1,2/) 3267 lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres) 3268 3269 amres = True 3270 amres@amParallelPosF = 0.65 3271 amres@amOrthogonalPosF = -0.2 3272 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3273 overlay(plot(n),plot_wsptsBC) 3274 wsptsBC=1 3275 else 3276 if (prof3d .EQ. 0)then 3277 varn=varn+1 3278 end if 3279 continue 3280 end if 3281 end if 3282 if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then 3283 if (wpvptp .EQ. 0) then 3284 res@gsnLeftString = "wpvptp, wsvpts and wvpt" 3285 res@tiXAxisString = "["+unit(varn)+"]" 3286 res@gsnRightString = " " 3287 if (xs .EQ. -1) then 3288 res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) 3289 else 3290 res@trXMinF = xs 3291 end if 3292 if (xe .EQ. -1) then 3293 res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) 3294 else 3295 res@trXMaxF = xe 3296 end if 3297 if (vNam(varn) .EQ. "wsvpts")then 3298 res@xyDashPattern = 1 3299 end if 3300 if (vNam(varn) .EQ. "wvpt")then 3301 res@xyDashPattern = 2 3302 end if 3303 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3304 3305 ; *************************************************** 3306 ; legend for combined plot 3307 ; *************************************************** 3308 3309 lgres = True 3310 lgMonoDashIndex = False 3311 lgres@lgLabelFont = "helvetica" 3312 lgres@lgLabelFontHeightF = font_size_legend 3313 lgres@vpWidthF = 0.1 3314 lgres@vpHeightF = 0.12 3315 lgres@lgDashIndexes = (/0,1,2/) 3316 lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres) 3317 3318 amres = True 3319 amres@amParallelPosF = 0.65 3320 amres@amOrthogonalPosF = -0.2 3321 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3322 overlay(plot(n),plot_wpvptp) 3323 wpvptp=1 3324 else 3325 if (prof3d .EQ. 0)then 3326 varn=varn+1 3327 end if 3328 continue 3329 end if 3330 end if 3331 if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then 3332 if (wpqp .EQ. 0) then 3333 res@gsnLeftString = "wpqp, wsqs and wq" 3334 res@tiXAxisString = "["+unit(varn)+"]" 3335 res@gsnRightString = " " 3336 if (xs .EQ. -1) then 3337 res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) 3338 else 3339 res@trXMinF = xs 3340 end if 3341 if (xe .EQ. -1) then 3342 res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) 3343 else 3344 res@trXMaxF = xe 3345 end if 3346 if (vNam(varn) .EQ. "wsqs")then 3347 res@xyDashPattern = 1 3348 end if 3349 if (vNam(varn) .EQ. "wq")then 3350 res@xyDashPattern = 2 3351 end if 3352 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3353 3354 ; *************************************************** 3355 ; legend for combined plot 3356 ; *************************************************** 3357 3358 lgres = True 3359 lgMonoDashIndex = False 3360 lgres@lgLabelFont = "helvetica" 3361 lgres@lgLabelFontHeightF = font_size_legend 3362 lgres@vpWidthF = 0.08 3363 lgres@vpHeightF = 0.12 3364 lgres@lgDashIndexes = (/0,1,2/) 3365 lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres) 3366 3367 amres = True 3368 amres@amParallelPosF = 0.65 3369 amres@amOrthogonalPosF = -0.2 3370 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3371 overlay(plot(n),plot_wpqp) 3372 wpqp=1 3373 else 3374 if (prof3d .EQ. 0)then 3375 varn=varn+1 3376 end if 3377 continue 3378 end if 3379 end if 3380 if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then 3381 if (wpqvp .EQ. 0) then 3382 res@gsnLeftString = "wpqvp, wsqvs and wqv" 3383 res@tiXAxisString = "["+unit(varn)+"]" 3384 res@gsnRightString = " " 3385 if (xs .EQ. -1) then 3386 res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) 3387 else 3388 res@trXMinF = xs 3389 end if 3390 if (xe .EQ. -1) then 3391 res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) 3392 else 3393 res@trXMaxF = xe 3394 end if 3395 if (vNam(varn) .EQ. "wsqvs")then 3396 res@xyDashPattern = 1 3397 end if 3398 if (vNam(varn) .EQ. "wqv")then 3399 res@xyDashPattern = 2 3400 end if 3401 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3402 3403 ; *************************************************** 3404 ; legend for combined plot 3405 ; *************************************************** 3406 3407 lgres = True 3408 lgMonoDashIndex = False 3409 lgres@lgLabelFont = "helvetica" 3410 lgres@lgLabelFontHeightF = font_size_legend 3411 lgres@vpWidthF = 0.09 3412 lgres@vpHeightF = 0.12 3413 lgres@lgDashIndexes = (/0,1,2/) 3414 lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres) 3415 3416 amres = True 3417 amres@amParallelPosF = 0.65 3418 amres@amOrthogonalPosF = -0.2 3419 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3420 overlay(plot(n),plot_wpqvp) 3421 wpqvp=1 3422 else 3423 if (prof3d .EQ. 0)then 3424 varn=varn+1 3425 end if 3426 continue 3427 end if 3428 end if 3429 if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then 3430 if (wpsp .EQ. 0) then 3431 res@gsnLeftString = "wpsp, wsss and ws" 3432 res@tiXAxisString = "["+unit(varn)+"]" 3433 res@gsnRightString = " " 3434 if (xs .EQ. -1) then 3435 res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) 3436 else 3437 res@trXMinF = xs 3438 end if 3439 if (xe .EQ. -1) then 3440 res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) 3441 else 3442 res@trXMaxF = xe 3443 end if 3444 if (vNam(varn) .EQ. "wsss")then 3445 res@xyDashPattern = 1 3446 end if 3447 if (vNam(varn) .EQ. "ws")then 3448 res@xyDashPattern = 2 3449 end if 3450 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3451 3452 ; *************************************************** 3453 ; legend for combined plot 3454 ; *************************************************** 3455 3456 lgres = True 3457 lgMonoDashIndex = False 3458 lgres@lgLabelFont = "helvetica" 3459 lgres@lgLabelFontHeightF = font_size_legend 3460 lgres@vpWidthF = 0.08 3461 lgres@vpHeightF = 0.12 3462 lgres@lgDashIndexes = (/0,1,2/) 3463 lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres) 3464 3465 amres = True 3466 amres@amParallelPosF = 0.65 3467 amres@amOrthogonalPosF = -0.2 3468 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3469 overlay(plot(n),plot_wpsp) 3470 wpsp=1 3471 else 3472 if (prof3d .EQ. 0)then 3473 varn=varn+1 3474 end if 3475 continue 3476 end if 3477 end if 3478 if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then 3479 if (wpsap .EQ. 0) then 3480 res@gsnLeftString = "wpsap, wssas and wsa" 3481 res@tiXAxisString = "["+unit(varn)+"]" 3482 res@gsnRightString = " " 3483 if (xs .EQ. -1) then 3484 res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) 3485 else 3486 res@trXMinF = xs 3487 end if 3488 if (xe .EQ. -1) then 3489 res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) 3490 else 3491 res@trXMaxF = xe 3492 end if 3493 if (vNam(varn) .EQ. "wssas")then 3494 res@xyDashPattern = 1 3495 end if 3496 if (vNam(varn) .EQ. "wsa")then 3497 res@xyDashPattern = 2 3498 end if 3499 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3500 3501 ; *************************************************** 3502 ; legend for combined plot 3503 ; *************************************************** 3504 3505 lgres = True 3506 lgMonoDashIndex = False 3507 lgres@lgLabelFont = "helvetica" 3508 lgres@lgLabelFontHeightF = font_size_legend 3509 lgres@vpWidthF = 0.09 3510 lgres@vpHeightF = 0.12 3511 lgres@lgDashIndexes = (/0,1,2/) 3512 lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres) 3513 3514 amres = True 3515 amres@amParallelPosF = 0.65 3516 amres@amOrthogonalPosF = -0.2 3517 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3518 overlay(plot(n),plot_wpsap) 3519 wpsap=1 3520 else 3521 if (prof3d .EQ. 0)then 3522 varn=varn+1 3523 end if 3524 continue 3525 end if 3526 end if 3527 3528 if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then 3529 if (us2 .EQ. 0) then 3530 res@gsnLeftString = "us2, vs2 and ws2" 3531 res@tiXAxisString = "["+unit(varn)+"]" 3532 res@gsnRightString = " " 3533 if (xs .EQ. -1) then 3534 res@trXMinF = min((/minius2,minivs2,miniws2/)) 3535 else 3536 res@trXMinF = xs 3537 end if 3538 if (xe .EQ. -1) then 3539 res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) 3540 else 3541 res@trXMaxF = xe 3542 end if 3543 if (vNam(varn) .EQ. "vs2")then 3544 res@xyDashPattern = 1 3545 end if 3546 if (vNam(varn) .EQ. "ws2")then 3547 res@xyDashPattern = 2 3548 end if 3549 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3550 3551 ; *************************************************** 3552 ; legend for combined plot 3553 ; *************************************************** 3554 3555 lgres = True 3556 lgMonoDashIndex = False 3557 lgres@lgLabelFont = "helvetica" 3558 lgres@lgLabelFontHeightF = font_size_legend 3559 lgres@vpWidthF = 0.07 3560 lgres@vpHeightF = 0.12 3561 lgres@lgDashIndexes = (/0,1,2/) 3562 lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres) 3563 3564 amres = True 3565 amres@amParallelPosF = 0.65 3566 amres@amOrthogonalPosF = -0.2 3567 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3568 overlay(plot(n),plot_us2) 3569 us2=1 3570 else 3571 if (prof3d .EQ. 0)then 3572 varn=varn+1 3573 end if 3574 continue 3575 end if 3576 end if 3577 3578 if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then 3579 if (wsususodz .EQ. 0) then 3580 res@gsnLeftString = "wsususodz, wspsodz and wpeodz" 3581 res@tiXAxisString = "["+unit(varn)+"]" 3582 res@gsnRightString = " " 3583 if (xs .EQ. -1) then 3584 res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) 3585 else 3586 res@trXMinF = xs 3587 end if 3588 if (xe .EQ. -1) then 3589 res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) 3590 else 3591 res@trXMaxF = xe 3592 end if 3593 if (vNam(varn) .EQ. "wspsodz")then 3594 res@xyDashPattern = 1 3595 end if 3596 if (vNam(varn) .EQ. "wpeodz")then 3597 res@xyDashPattern = 2 3598 end if 2909 3599 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2910 3600 … … 2918 3608 lgres@lgLabelFontHeightF = font_size_legend 2919 3609 lgres@vpWidthF = 0.12 2920 lgres@vpHeightF = 0.1 2921 lgres@lgDashIndexes = (/0,1,2/) 2922 lbid = gsn_create_legend(wks,3,(/"pt","vpt","lpt"/),lgres) 2923 2924 amres = True 2925 amres@amParallelPosF = 0.65 2926 amres@amOrthogonalPosF = -0.2 2927 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2928 overlay(plot(n),plot_pt) 2929 pt=1 2930 else 2931 if (prof3d .EQ. 0)then 2932 varn=varn+1 2933 end if 2934 continue 2935 end if 2936 end if 2937 if (vNam(varn) .EQ. "q" .OR. vNam(varn) .EQ. "qv" .OR. vNam(varn) .EQ. "ql") then 2938 if (q .EQ. 0) then 2939 res@gsnLeftString = "q, qv and ql" 2940 res@tiXAxisString = "["+unit(varn)+"]" 2941 res@gsnRightString = " " 2942 if (xs .EQ. -1) then 2943 res@trXMinF = min((/miniq,miniqv,miniql/)) 2944 else 2945 res@trXMinF = xs 2946 end if 2947 if (xe .EQ. -1) then 2948 res@trXMaxF = max((/maxiq,maxiqv,maxiql/)) 2949 else 2950 res@trXMaxF = xe 2951 end if 2952 2953 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2954 2955 ; *************************************************** 2956 ; legend for combined plot 2957 ; *************************************************** 2958 2959 lgres = True 2960 lgMonoDashIndex = False 2961 lgres@lgLabelFont = "helvetica" 2962 lgres@lgLabelFontHeightF = font_size_legend 2963 lgres@vpWidthF = 0.12 2964 lgres@vpHeightF = 0.1 2965 lgres@lgDashIndexes = (/0,1,2/) 2966 lbid = gsn_create_legend(wks,3,(/"q","qv","ql"/),lgres) 2967 2968 amres = True 2969 amres@amParallelPosF = 0.65 2970 amres@amOrthogonalPosF = -0.2 2971 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 2972 overlay(plot(n),plot_q) 2973 q=1 2974 else 2975 if (prof3d .EQ. 0)then 2976 varn=varn+1 2977 end if 2978 continue 2979 end if 2980 end if 2981 2982 if (vNam(varn) .EQ. "e" .OR. vNam(varn) .EQ. "es") then 2983 if (e .EQ. 0) then 2984 res@gsnLeftString = "e and es" 2985 res@tiXAxisString = "["+unit(varn)+"]" 2986 res@gsnRightString = " " 2987 if (xs .EQ. -1) then 2988 res@trXMinF = min((/minie,minies/)) 2989 else 2990 res@trXMinF = xs 2991 end if 2992 if (xe .EQ. -1) then 2993 res@trXMaxF = max((/maxie,maxies/)) 2994 else 2995 res@trXMaxF = xe 2996 end if 2997 2998 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 2999 3000 ; *************************************************** 3001 ; legend for combined plot 3002 ; *************************************************** 3003 3004 lgres = True 3005 lgMonoDashIndex = False 3006 lgres@lgLabelFont = "helvetica" 3007 lgres@lgLabelFontHeightF = font_size_legend 3008 lgres@vpWidthF = 0.12 3009 lgres@vpHeightF = 0.1 3010 lgres@lgDashIndexes = (/0,1,2/) 3011 lbid = gsn_create_legend(wks,3,(/"e","es"/),lgres) 3012 3013 amres = True 3014 amres@amParallelPosF = 0.65 3015 amres@amOrthogonalPosF = -0.2 3016 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3017 overlay(plot(n),plot_e) 3018 e=1 3019 else 3020 if (prof3d .EQ. 0)then 3021 varn=varn+1 3022 end if 3023 continue 3024 end if 3025 end if 3026 if (vNam(varn) .EQ. "km" .OR. vNam(varn) .EQ. "kh") then 3027 if (km .EQ. 0) then 3028 res@gsnLeftString = "km and kh" 3029 res@tiXAxisString = "["+unit(varn)+"]" 3030 res@gsnRightString = " " 3031 if (xs .EQ. -1) then 3032 res@trXMinF = min((/minikm,minikh/)) 3033 else 3034 res@trXMinF = xs 3035 end if 3036 if (xe .EQ. -1) then 3037 res@trXMaxF = max((/maxikm,maxikh/)) 3038 else 3039 res@trXMaxF = xe 3040 end if 3041 3042 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3043 3044 ; *************************************************** 3045 ; legend for combined plot 3046 ; *************************************************** 3047 3048 lgres = True 3049 lgMonoDashIndex = False 3050 lgres@lgLabelFont = "helvetica" 3051 lgres@lgLabelFontHeightF = font_size_legend 3052 lgres@vpWidthF = 0.12 3053 lgres@vpHeightF = 0.1 3054 lgres@lgDashIndexes = (/0,1,2/) 3055 lbid = gsn_create_legend(wks,3,(/"km","kh"/),lgres) 3056 3057 amres = True 3058 amres@amParallelPosF = 0.65 3059 amres@amOrthogonalPosF = -0.2 3060 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3061 overlay(plot(n),plot_km) 3062 km=1 3063 else 3064 if (prof3d .EQ. 0)then 3065 varn=varn+1 3066 end if 3067 continue 3068 end if 3069 end if 3070 3071 if (vNam(varn) .EQ. "wpup" .OR. vNam(varn) .EQ. "wsus" .OR. vNam(varn) .EQ. "wu") then 3072 if (wpup .EQ. 0) then 3073 res@gsnLeftString = "wpup, wsus and wu" 3074 res@tiXAxisString = "["+unit(varn)+"]" 3075 res@gsnRightString = " " 3076 if (xs .EQ. -1) then 3077 res@trXMinF = min((/miniwpup,miniwsus,miniwu/)) 3078 else 3079 res@trXMinF = xs 3080 end if 3081 if (xe .EQ. -1) then 3082 res@trXMaxF = max((/maxiwpup,maxiwsus,maxiwu/)) 3083 else 3084 res@trXMaxF = xe 3085 end if 3086 3087 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3088 3089 ; *************************************************** 3090 ; legend for combined plot 3091 ; *************************************************** 3092 3093 lgres = True 3094 lgMonoDashIndex = False 3095 lgres@lgLabelFont = "helvetica" 3096 lgres@lgLabelFontHeightF = font_size_legend 3097 lgres@vpWidthF = 0.12 3098 lgres@vpHeightF = 0.1 3099 lgres@lgDashIndexes = (/0,1,2/) 3100 lbid = gsn_create_legend(wks,3,(/"wpup","wsus","wu"/),lgres) 3101 3102 amres = True 3103 amres@amParallelPosF = 0.65 3104 amres@amOrthogonalPosF = -0.2 3105 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3106 overlay(plot(n),plot_wpup) 3107 wpup=1 3108 else 3109 if (prof3d .EQ. 0)then 3110 varn=varn+1 3111 end if 3112 continue 3113 end if 3114 end if 3115 if (vNam(varn) .EQ. "wpvp" .OR. vNam(varn) .EQ. "wsvs" .OR.vNam(varn) .EQ. "wv") then 3116 if (wpvp .EQ. 0) then 3117 res@gsnLeftString = "wpvp, wsus and wv" 3118 res@tiXAxisString = "["+unit(varn)+"]" 3119 res@gsnRightString = " " 3120 if (xs .EQ. -1) then 3121 res@trXMinF = min((/miniwpvp,miniwsvs,miniwv/)) 3122 else 3123 res@trXMinF = xs 3124 end if 3125 if (xe .EQ. -1) then 3126 res@trXMaxF = max((/maxiwpvp,maxiwsvs,maxiwv/)) 3127 else 3128 res@trXMaxF = xe 3129 end if 3130 3131 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3132 3133 ; *************************************************** 3134 ; legend for combined plot 3135 ; *************************************************** 3136 3137 lgres = True 3138 lgMonoDashIndex = False 3139 lgres@lgLabelFont = "helvetica" 3140 lgres@lgLabelFontHeightF = font_size_legend 3141 lgres@vpWidthF = 0.12 3142 lgres@vpHeightF = 0.1 3143 lgres@lgDashIndexes = (/0,1,2/) 3144 lbid = gsn_create_legend(wks,3,(/"wpvp","wsvs","wv"/),lgres) 3145 3146 amres = True 3147 amres@amParallelPosF = 0.65 3148 amres@amOrthogonalPosF = -0.2 3149 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3150 overlay(plot(n),plot_wpvp) 3151 wpup=0 3152 else 3153 if (prof3d .EQ. 0)then 3154 varn=varn+1 3155 end if 3156 continue 3157 end if 3158 end if 3159 if (vNam(varn) .EQ. "wpptp" .OR. vNam(varn) .EQ. "wspts" .OR. vNam(varn) .EQ. "wpt") then 3160 if (wpptp .EQ. 0) then 3161 res@gsnLeftString = "wpptp, wspts and wv" 3162 res@tiXAxisString = "["+unit(varn)+"]" 3163 res@gsnRightString = " " 3164 if (xs .EQ. -1) then 3165 res@trXMinF = min((/miniwpptp,miniwspts,miniwpt/)) 3166 else 3167 res@trXMinF = xs 3168 end if 3169 if (xe .EQ. -1) then 3170 res@trXMaxF = max((/maxiwpptp,maxiwspts,maxiwpt/)) 3171 else 3172 res@trXMaxF = xe 3173 end if 3174 3175 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3176 3177 ; *************************************************** 3178 ; legend for combined plot 3179 ; *************************************************** 3180 3181 lgres = True 3182 lgMonoDashIndex = False 3183 lgres@lgLabelFont = "helvetica" 3184 lgres@lgLabelFontHeightF = font_size_legend 3185 lgres@vpWidthF = 0.12 3186 lgres@vpHeightF = 0.1 3187 lgres@lgDashIndexes = (/0,1,2/) 3188 lbid = gsn_create_legend(wks,3,(/"wpptp","wspts","wpt"/),lgres) 3189 3190 amres = True 3191 amres@amParallelPosF = 0.65 3192 amres@amOrthogonalPosF = -0.2 3193 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3194 overlay(plot(n),plot_wpptp) 3195 wpptp=1 3196 else 3197 if (prof3d .EQ. 0)then 3198 varn=varn+1 3199 end if 3200 continue 3201 end if 3202 end if 3203 if (vNam(varn) .EQ. "wsptsBC" .OR. vNam(varn) .EQ. "wptBC") then 3204 if (wsptsBC .EQ. 0) then 3205 res@gsnLeftString = "wsptsBC and wptBC" 3206 res@tiXAxisString = "["+unit(varn)+"]" 3207 res@gsnRightString = " " 3208 if (xs .EQ. -1) then 3209 res@trXMinF = min((/miniwsptsBC,miniwptBC/)) 3210 else 3211 res@trXMinF = xs 3212 end if 3213 if (xe .EQ. -1) then 3214 res@trXMaxF = max((/maxiwsptsBC,maxiwptBC/)) 3215 else 3216 res@trXMaxF = xe 3217 end if 3218 3219 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3220 3221 ; *************************************************** 3222 ; legend for combined plot 3223 ; *************************************************** 3224 3225 lgres = True 3226 lgMonoDashIndex = False 3227 lgres@lgLabelFont = "helvetica" 3228 lgres@lgLabelFontHeightF = font_size_legend 3229 lgres@vpWidthF = 0.12 3230 lgres@vpHeightF = 0.1 3231 lgres@lgDashIndexes = (/0,1,2/) 3232 lbid = gsn_create_legend(wks,3,(/"wsptsBC","wptBC"/),lgres) 3233 3234 amres = True 3235 amres@amParallelPosF = 0.65 3236 amres@amOrthogonalPosF = -0.2 3237 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3238 overlay(plot(n),plot_wsptsBC) 3239 wsptsBC=1 3240 else 3241 if (prof3d .EQ. 0)then 3242 varn=varn+1 3243 end if 3244 continue 3245 end if 3246 end if 3247 if (vNam(varn) .EQ. "wpvptp" .OR. vNam(varn) .EQ. "wsvpts" .OR. vNam(varn) .EQ. "wvpt") then 3248 if (wpvptp .EQ. 0) then 3249 res@gsnLeftString = "wpvptp, wsvpts and wv" 3250 res@tiXAxisString = "["+unit(varn)+"]" 3251 res@gsnRightString = " " 3252 if (xs .EQ. -1) then 3253 res@trXMinF = min((/miniwpvptp,miniwsvpts,miniwvpt/)) 3254 else 3255 res@trXMinF = xs 3256 end if 3257 if (xe .EQ. -1) then 3258 res@trXMaxF = max((/maxiwpvptp,maxiwsvpts,maxiwvpt/)) 3259 else 3260 res@trXMaxF = xe 3261 end if 3262 3263 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3264 3265 ; *************************************************** 3266 ; legend for combined plot 3267 ; *************************************************** 3268 3269 lgres = True 3270 lgMonoDashIndex = False 3271 lgres@lgLabelFont = "helvetica" 3272 lgres@lgLabelFontHeightF = font_size_legend 3273 lgres@vpWidthF = 0.12 3274 lgres@vpHeightF = 0.1 3275 lgres@lgDashIndexes = (/0,1,2/) 3276 lbid = gsn_create_legend(wks,3,(/"wpvptp","wsvpts","wvpt"/),lgres) 3277 3278 amres = True 3279 amres@amParallelPosF = 0.65 3280 amres@amOrthogonalPosF = -0.2 3281 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3282 overlay(plot(n),plot_wpvptp) 3283 wpvptp=1 3284 else 3285 if (prof3d .EQ. 0)then 3286 varn=varn+1 3287 end if 3288 continue 3289 end if 3290 end if 3291 if (vNam(varn) .EQ. "wpqp" .OR. vNam(varn) .EQ. "wsqs" .OR. vNam(varn) .EQ. "wq") then 3292 if (wpqp .EQ. 0) then 3293 res@gsnLeftString = "wpqp, wsqs and wq" 3294 res@tiXAxisString = "["+unit(varn)+"]" 3295 res@gsnRightString = " " 3296 if (xs .EQ. -1) then 3297 res@trXMinF = min((/miniwpqp,miniwsqs,miniwq/)) 3298 else 3299 res@trXMinF = xs 3300 end if 3301 if (xe .EQ. -1) then 3302 res@trXMaxF = max((/maxiwpqp,maxiwsqs,maxiwq/)) 3303 else 3304 res@trXMaxF = xe 3305 end if 3306 3307 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3308 3309 ; *************************************************** 3310 ; legend for combined plot 3311 ; *************************************************** 3312 3313 lgres = True 3314 lgMonoDashIndex = False 3315 lgres@lgLabelFont = "helvetica" 3316 lgres@lgLabelFontHeightF = font_size_legend 3317 lgres@vpWidthF = 0.12 3318 lgres@vpHeightF = 0.1 3319 lgres@lgDashIndexes = (/0,1,2/) 3320 lbid = gsn_create_legend(wks,3,(/"wpqp","wsqs","wq"/),lgres) 3321 3322 amres = True 3323 amres@amParallelPosF = 0.65 3324 amres@amOrthogonalPosF = -0.2 3325 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3326 overlay(plot(n),plot_wpqp) 3327 wpqp=1 3328 else 3329 if (prof3d .EQ. 0)then 3330 varn=varn+1 3331 end if 3332 continue 3333 end if 3334 end if 3335 if (vNam(varn) .EQ. "wpqvp" .OR. vNam(varn) .EQ. "wsqvs" .OR. vNam(varn) .EQ. "wqv") then 3336 if (wpqvp .EQ. 0) then 3337 res@gsnLeftString = "wpqvp, wsqvs and wqv" 3338 res@tiXAxisString = "["+unit(varn)+"]" 3339 res@gsnRightString = " " 3340 if (xs .EQ. -1) then 3341 res@trXMinF = min((/miniwpqp,miniwsqvs,miniwqv/)) 3342 else 3343 res@trXMinF = xs 3344 end if 3345 if (xe .EQ. -1) then 3346 res@trXMaxF = max((/maxiwpqp,maxiwsqvs,maxiwqv/)) 3347 else 3348 res@trXMaxF = xe 3349 end if 3350 3351 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3352 3353 ; *************************************************** 3354 ; legend for combined plot 3355 ; *************************************************** 3356 3357 lgres = True 3358 lgMonoDashIndex = False 3359 lgres@lgLabelFont = "helvetica" 3360 lgres@lgLabelFontHeightF = font_size_legend 3361 lgres@vpWidthF = 0.12 3362 lgres@vpHeightF = 0.1 3363 lgres@lgDashIndexes = (/0,1,2/) 3364 lbid = gsn_create_legend(wks,3,(/"wpqvp","wsqvs","wqv"/),lgres) 3365 3366 amres = True 3367 amres@amParallelPosF = 0.65 3368 amres@amOrthogonalPosF = -0.2 3369 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3370 overlay(plot(n),plot_wpqvp) 3371 wpqvp=1 3372 else 3373 if (prof3d .EQ. 0)then 3374 varn=varn+1 3375 end if 3376 continue 3377 end if 3378 end if 3379 if (vNam(varn) .EQ. "wpsp" .OR. vNam(varn) .EQ. "wsss" .OR. vNam(varn) .EQ. "ws") then 3380 if (wpsp .EQ. 0) then 3381 res@gsnLeftString = "wpsp, wsss and ws" 3382 res@tiXAxisString = "["+unit(varn)+"]" 3383 res@gsnRightString = " " 3384 if (xs .EQ. -1) then 3385 res@trXMinF = min((/miniwpsp,miniwsss,miniws/)) 3386 else 3387 res@trXMinF = xs 3388 end if 3389 if (xe .EQ. -1) then 3390 res@trXMaxF = max((/maxiwpsp,maxiwsss,maxiws/)) 3391 else 3392 res@trXMaxF = xe 3393 end if 3394 3395 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3396 3397 ; *************************************************** 3398 ; legend for combined plot 3399 ; *************************************************** 3400 3401 lgres = True 3402 lgMonoDashIndex = False 3403 lgres@lgLabelFont = "helvetica" 3404 lgres@lgLabelFontHeightF = font_size_legend 3405 lgres@vpWidthF = 0.12 3406 lgres@vpHeightF = 0.1 3407 lgres@lgDashIndexes = (/0,1,2/) 3408 lbid = gsn_create_legend(wks,3,(/"wpsp","wsss","ws"/),lgres) 3409 3410 amres = True 3411 amres@amParallelPosF = 0.65 3412 amres@amOrthogonalPosF = -0.2 3413 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3414 overlay(plot(n),plot_wpsp) 3415 wpsp=1 3416 else 3417 if (prof3d .EQ. 0)then 3418 varn=varn+1 3419 end if 3420 continue 3421 end if 3422 end if 3423 if (vNam(varn) .EQ. "wpsap" .OR.vNam(varn) .EQ. "wssas" .OR. vNam(varn) .EQ. "wsa") then 3424 if (wpsap .EQ. 0) then 3425 res@gsnLeftString = "wpsap, wssas and wsa" 3426 res@tiXAxisString = "["+unit(varn)+"]" 3427 res@gsnRightString = " " 3428 if (xs .EQ. -1) then 3429 res@trXMinF = min((/miniwpsap,miniwssas,miniwsa/)) 3430 else 3431 res@trXMinF = xs 3432 end if 3433 if (xe .EQ. -1) then 3434 res@trXMaxF = max((/maxiwpsap,maxiwssas,maxiwsa/)) 3435 else 3436 res@trXMaxF = xe 3437 end if 3438 3439 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3440 3441 ; *************************************************** 3442 ; legend for combined plot 3443 ; *************************************************** 3444 3445 lgres = True 3446 lgMonoDashIndex = False 3447 lgres@lgLabelFont = "helvetica" 3448 lgres@lgLabelFontHeightF = font_size_legend 3449 lgres@vpWidthF = 0.12 3450 lgres@vpHeightF = 0.1 3451 lgres@lgDashIndexes = (/0,1,2/) 3452 lbid = gsn_create_legend(wks,3,(/"wpsap","wssas","wsa"/),lgres) 3453 3454 amres = True 3455 amres@amParallelPosF = 0.65 3456 amres@amOrthogonalPosF = -0.2 3457 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3458 overlay(plot(n),plot_wpsap) 3459 wpsap=1 3460 else 3461 if (prof3d .EQ. 0)then 3462 varn=varn+1 3463 end if 3464 continue 3465 end if 3466 end if 3467 3468 if (vNam(varn) .EQ. "us2" .OR. vNam(varn) .EQ. "vs2" .OR. vNam(varn) .EQ. "ws2") then 3469 if (us2 .EQ. 0) then 3470 res@gsnLeftString = "us2, vs2 and ws2" 3471 res@tiXAxisString = "["+unit(varn)+"]" 3472 res@gsnRightString = " " 3473 if (xs .EQ. -1) then 3474 res@trXMinF = min((/minius2,minivs2,miniws2/)) 3475 else 3476 res@trXMinF = xs 3477 end if 3478 if (xe .EQ. -1) then 3479 res@trXMaxF = max((/maxius2,maxivs2,maxiws2/)) 3480 else 3481 res@trXMaxF = xe 3482 end if 3483 3484 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3485 3486 ; *************************************************** 3487 ; legend for combined plot 3488 ; *************************************************** 3489 3490 lgres = True 3491 lgMonoDashIndex = False 3492 lgres@lgLabelFont = "helvetica" 3493 lgres@lgLabelFontHeightF = font_size_legend 3494 lgres@vpWidthF = 0.12 3495 lgres@vpHeightF = 0.1 3496 lgres@lgDashIndexes = (/0,1,2/) 3497 lbid = gsn_create_legend(wks,3,(/"us2","vs2","ws2"/),lgres) 3498 3499 amres = True 3500 amres@amParallelPosF = 0.65 3501 amres@amOrthogonalPosF = -0.2 3502 annoid1 = gsn_add_annotation(plot(n),lbid,amres) 3503 overlay(plot(n),plot_us2) 3504 us2=1 3505 else 3506 if (prof3d .EQ. 0)then 3507 varn=varn+1 3508 end if 3509 continue 3510 end if 3511 end if 3512 3513 if (vNam(varn) .EQ. "wsususodz" .OR. vNam(varn) .EQ. "wspsodz" .OR. vNam(varn) .EQ. "wpeodz" ) then 3514 if (wsususodz .EQ. 0) then 3515 res@gsnLeftString = "wsususodz, wspsodz and ws2" 3516 res@tiXAxisString = "["+unit(varn)+"]" 3517 res@gsnRightString = " " 3518 if (xs .EQ. -1) then 3519 res@trXMinF = min((/miniwsususodz,miniwspsodz,miniwpeodz/)) 3520 else 3521 res@trXMinF = xs 3522 end if 3523 if (xe .EQ. -1) then 3524 res@trXMaxF = max((/maxiwsususodz,maxiwspsodz,maxiwpeodz/)) 3525 else 3526 res@trXMaxF = xe 3527 end if 3528 3529 plot(n) = gsn_csm_xy(wks,data(varn,:,:),z,res) 3530 3531 ; *************************************************** 3532 ; legend for combined plot 3533 ; *************************************************** 3534 3535 lgres = True 3536 lgMonoDashIndex = False 3537 lgres@lgLabelFont = "helvetica" 3538 lgres@lgLabelFontHeightF = font_size_legend 3539 lgres@vpWidthF = 0.12 3540 lgres@vpHeightF = 0.1 3610 lgres@vpHeightF = 0.12 3541 3611 lgres@lgDashIndexes = (/0,1,2/) 3542 3612 lbid = gsn_create_legend(wks,3,(/"wsususodz","wspsodz","wpeodz"/),lgres) … … 3590 3660 if (prof3d .EQ. 0) then 3591 3661 if (log_z .EQ. 1) then 3592 z = f ->$vNam(varn+1)$(1:dimz-1)3593 else 3594 z = f ->$vNam(varn+1)$3662 z = f_att->$vNam(varn+1)$(1:dimz-1) 3663 else 3664 z = f_att->$vNam(varn+1)$ 3595 3665 end if 3596 3666 else … … 3718 3788 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 3719 3789 gsn_panel(wks,plot_,(/n,1/),resP) 3790 print(" ") 3791 print("Outputs to .eps or .epsi have only one frame") 3792 print(" ") 3720 3793 else 3721 3794 do i = 0,n-1, no_rows*no_columns … … 3723 3796 gsn_panel(wks,plot_(i:n-1),(/no_rows,no_columns/),resP) 3724 3797 else 3725 gsn_panel(wks,plot_(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP)3798 gsn_panel(wks,plot_(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP) 3726 3799 end if 3727 3800 end do -
palm/trunk/SCRIPTS/NCL/spectra.ncl
r218 r250 639 639 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 640 640 gsn_panel(wks,plot,(/n,1/),resP) 641 print(" ") 642 print("Outputs to .eps or .epsi have only one frame") 643 print(" ") 641 644 else 642 645 do i = 0,n-1, no_rows*no_columns -
palm/trunk/SCRIPTS/NCL/timeseries.ncl
r218 r250 1243 1243 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 1244 1244 gsn_panel(wks_ps,plot_ps(1:n),(/n,1/),resP) 1245 print(" ") 1246 print("Outputs to .eps or .epsi have only one frame") 1247 print(" ") 1245 1248 else 1246 1249 do np = 1,n,no_rows*no_columns
Note: See TracChangeset
for help on using the changeset viewer.
