function B = cvdtherm2eqndyn(filterflag)% cvdtherm2eqndyn.m A nonlinear equation model object constructor, derived%                   from the neqmodel class. Called as%%          A = cvdtherm2eqndynmodel_param_script % set dimensionless parameter values for this modelif nargin == 0 filterflag = 0; endtfinal = 60;nt = 100;t = quadgrid('slab',nt,'t',[0 tfinal]);Dt = linearoperator(t,'d','t');Tw0 = 1;Tq0 = 1;Tw = scalarfield(t,Tw0);Tq = scalarfield(t,Tq0);u = mask(t,'lineseg',[0 tfinal/2]);if filterflag % filtered approach    u = bpsf(u,1);endvar   = assocarray({'Tw' Tw 'Tq' Tq});param = assocarray({ 'e_wq' e_wq 'e_ws' e_ws 'k_w' k_w ...                     'q_w' q_w 'e_qw' e_qw 'e_qs' e_qs ...                     'nu_q' nu_q 'q_q' q_q 'u' u ...                     'F' F 'alf' alf 'Ta' Ta 'Dt' Dt ...                     'Tw0' Tw0 'Tq0' Tq0 ...                   });A = naemodel('cvdtherm2eqndyn',var,param);B = struct([]);B = class(B,'cvdtherm2eqndyn',A);