0001 function h_g =cieloprognostic_core_linear(m,params,inputs,analysis)
0002
0003
0004
0005 if size(m.y_s,1)==0,
0006 flag_y_s= 0;
0007 else
0008 flag_y_s= 1;
0009 end
0010
0011
0012 params.kflag=1; params.pflag=1;
0013 params.ktflag=0;
0014
0015
0016 [K_gg , p_g , dK_gg]=Emg(m.bgpdt,m.bgpdtb, m.est,m.lst,m.ept,m.mpt,m.geom3,params,inputs,analysis);
0017
0018
0019 [K_ff, K_fs] = Reducematrixfromgtof( K_gg, m.G_mn, flag_y_s );
0020 clear K_gg dK_gg
0021
0022
0023 [rp_f] = Reducerightside( rp_g, m.G_mn, rK_fs, m.y_s, flag_y_s );
0024 [clear p_g K_fs
0025
0026
0027 [h_f]=Solver(K_ff,p_f,params.solverstring);
0028 clear K_ff p_f
0029
0030
0031 h_g= Mergesolvec( h_f, m.G_mn, m.y_s, m.uset);
0032
0033 if count>2,
0034 clear h_f
0035 end