Changeset 922 for issm/trunk/src/m/solutions/cielo/diagnostic.m
- Timestamp:
- 06/11/09 15:53:37 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/solutions/cielo/diagnostic.m
r738 r922 11 11 %Build all models requested for diagnostic simulation 12 12 displaystring(md.debug,'%s',['reading diagnostic horiz model data']); 13 md.analysis_type='diagnostic'; md.sub_analysis_type='horiz'; m _dh=CreateFemModel(md);13 md.analysis_type='diagnostic'; md.sub_analysis_type='horiz'; models.dh=CreateFemModel(md); 14 14 15 15 displaystring(md.debug,'\n%s',['reading diagnostic vert model data']); 16 md.analysis_type='diagnostic'; md.sub_analysis_type='vert'; m _dv=CreateFemModel(md);16 md.analysis_type='diagnostic'; md.sub_analysis_type='vert'; models.dv=CreateFemModel(md); 17 17 18 18 displaystring(md.debug,'\n%s',['reading diagnostic stokes model data']); 19 md.analysis_type='diagnostic'; md.sub_analysis_type='stokes'; m _ds=CreateFemModel(md);19 md.analysis_type='diagnostic'; md.sub_analysis_type='stokes'; models.ds=CreateFemModel(md); 20 20 21 21 displaystring(md.debug,'\n%s',['reading diagnostic hutter model data']); 22 md.analysis_type='diagnostic'; md.sub_analysis_type='hutter'; m _dhu=CreateFemModel(md);22 md.analysis_type='diagnostic'; md.sub_analysis_type='hutter'; models.dhu=CreateFemModel(md); 23 23 24 24 displaystring(md.debug,'\n%s',['reading surface and bed slope computation model data']); 25 md.analysis_type='slope_compute'; md.sub_analysis_type=''; m _sl=CreateFemModel(md);25 md.analysis_type='slope_compute'; md.sub_analysis_type=''; models.sl=CreateFemModel(md); 26 26 27 27 % figure out number of dof: just for information purposes. 28 md.dof=modelsize(m _dh,m_dv,m_ds,m_dhu,m_sl);28 md.dof=modelsize(models); 29 29 30 30 %initialize inputs 31 31 inputs=inputlist; 32 inputs=add(inputs,'velocity',m _dh.parameters.u_g,'doublevec',3,m_dh.parameters.numberofnodes);32 inputs=add(inputs,'velocity',models.dh.parameters.u_g,'doublevec',3,models.dh.parameters.numberofnodes); 33 33 34 34 %compute solution 35 [u_g,p_g]=diagnostic_core(m_dh,m_dhu,m_dv,m_ds,m_sl,inputs);35 results=diagnostic_core(models,inputs); 36 36 37 % Load results onto model38 md=loadresults(md, u_g,p_g,m_dh,m_ds,m_dhu);37 %load results onto model 38 md=loadresults(md,models,results); 39 39 40 40 %stop timing
Note:
See TracChangeset
for help on using the changeset viewer.