1 120 !ifile ndays. ifile=1: generate elev2D.th; =2: salt3D.th and temp3D.th; =3: uv3D.th); ndays is the # of days needed;
1 1 !total # of boundary segments that need *3D.th, followed by list of boundary segment indices
0 !0: normal; 1: more outputs for debug