When you increase the grid size, you have to change
the temperary size of t in trdiag.f,trdgs.f, trdgv.f
The default is 200, but we may get more than that.

1. Everything is in file saltled.f, you have to change the scheme 
manually. Several parameters can be adjusted (swit,bpp).

2. when you change the scheme, you have to change "tmpn=" acording to
the method. High order scheme uses quick as base, don't forget to change 
sharp coefficient to be 1. instead of 1.e-9 in data/datain.
shp=1e-9 originally.
