Changeset 8308
- Timestamp:
- 05/16/11 15:59:04 (14 years ago)
- Location:
- issm/trunk/test
- Files:
-
- 12 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Miscellaneous/Bump/runme.m
r8011 r8308 5 5 md.cluster=astrid; 6 6 md.cluster.np=13; 7 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);7 md=solve(md,DiagnosticSolutionEnum); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk/test/Miscellaneous/GJM_test1/test1.m
r8011 r8308 4 4 md=setelementstype(md,'macayeal','all'); 5 5 md.cluster=none; 6 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);6 md=solve(md,DiagnosticSolutionEnum); -
issm/trunk/test/Miscellaneous/Regularization/runme.m
r8307 r8308 15 15 16 16 md.cluster=eval(oshostname); 17 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);17 md=solve(md,DiagnosticSolutionEnum); 18 18 md.vx_obs=PatchToVec(md.results.DiagnosticSolution.Vx); 19 19 md.vy_obs=PatchToVec(md.results.DiagnosticSolution.Vy); … … 32 32 md.weights=1*ones(md.numberofnodes,1); 33 33 md.eps_cm=NaN; 34 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);34 md=solve(md,DiagnosticSolutionEnum); 35 35 md.results.DiagnosticSolution_NO=md.results.DiagnosticSolution; 36 36 end … … 47 47 md.weights=1*ones(md.numberofnodes,1); 48 48 md.eps_cm=NaN; 49 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);49 md=solve(md,DiagnosticSolutionEnum); 50 50 md.results.DiagnosticSolution_YES1=md.results.DiagnosticSolution; 51 51 end … … 62 62 md.weights=1*ones(md.numberofnodes,1); 63 63 md.eps_cm=NaN; 64 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);64 md=solve(md,DiagnosticSolutionEnum); 65 65 md.results.DiagnosticSolution_YES2=md.results.DiagnosticSolution; 66 66 end … … 77 77 md.weights=1*ones(md.numberofnodes,1); 78 78 md.eps_cm=NaN; 79 md=solve(md, 'analysis_type',DiagnosticSolutionEnum);79 md=solve(md,DiagnosticSolutionEnum); 80 80 md.results.DiagnosticSolution_YES3=md.results.DiagnosticSolution; 81 81 end -
issm/trunk/test/Miscellaneous/connectivity/runme.m
r2553 r8308 49 49 if isparallel 50 50 md.cluster='wilkes'; 51 md=solve(md,' analysis_type','diagnostic_horiz');51 md=solve(md,'diagnostic_horiz'); 52 52 else 53 md=solve(md,' analysis_type','diagnostic');53 md=solve(md,'diagnostic'); 54 54 end 55 55 t2=clock; -
issm/trunk/test/Par/79North.par
r8305 r8308 43 43 %Boundary conditions: 44 44 md=SetMarineIceSheetBC(md); 45 pos=find(md. gridonboundary);45 pos=find(md.nodeonboundary); 46 46 md.spcthickness(pos,1)=1; 47 47 md.spcthickness(pos,2)=md.thickness(pos); -
issm/trunk/test/Par/ISMIPA.par
r8305 r8308 21 21 22 22 disp(' boundary conditions for diagnostic model'); 23 %Create gridon boundary fist (because we cannot use mesh)23 %Create node on boundary fist (because we cannot use mesh) 24 24 md=SetIceSheetBC(md); 25 25 -
issm/trunk/test/Par/ISMIPB.par
r8305 r8308 21 21 22 22 disp(' boundary conditions for diagnostic model'); 23 %Create gridon boundary fist (because we cannot use mesh)23 %Create node on boundary fist (because we cannot use mesh) 24 24 md=SetIceSheetBC(md); 25 25 md=SetParallel(md,8); -
issm/trunk/test/Par/ISMIPC.par
r8305 r8308 22 22 23 23 disp(' boundary conditions for diagnostic model: '); 24 %Create gridon boundary fist (because wi can not use mesh)24 %Create node on boundary fist (because wi can not use mesh) 25 25 md=SetIceSheetBC(md); 26 26 md=SetParallel(md,8); -
issm/trunk/test/Par/ISMIPD.par
r8305 r8308 21 21 22 22 disp(' boundary conditions for diagnostic model: '); 23 %Create gridon boundary fist (because wi can not use mesh)23 %Create node on boundary fist (because wi can not use mesh) 24 24 md=SetIceSheetBC(md); 25 25 -
issm/trunk/test/Par/ISMIPE.par
r8305 r8308 30 30 31 31 disp(' boundary conditions for diagnostic model: '); 32 %Create gridon boundary fist (because wi can not use mesh)32 %Create node on boundary fist (because wi can not use mesh) 33 33 md=SetIceSheetBC(md); 34 34 -
issm/trunk/test/Par/ISMIPF.par
r8305 r8308 21 21 22 22 disp(' boundary conditions for diagnostic model'); 23 %Create gridon boundary fist (because we cannot use mesh)23 %Create node on boundary fist (because we cannot use mesh) 24 24 md=SetIceSheetBC(md); 25 25 md.spcvelocity(:,4)=100*ones(md.numberofnodes,1); -
issm/trunk/test/Par/RoundSheetShelf.par
r8305 r8308 13 13 md.bed=-md.rho_ice/md.rho_water*md.thickness; 14 14 15 pos=find(md. gridonicesheet);15 pos=find(md.nodeonicesheet); 16 16 md.bed(pos)=md.bed(pos)-300*(radius(pos)-(rad-shelfextent))/(rad-shelfextent); 17 17 md.surface=md.bed+md.thickness; … … 48 48 md.accumulation_rate=-10*ones(md.numberofnodes,1); 49 49 md.melting_rate=zeros(md.numberofnodes,1); 50 pos=find(md. gridoniceshelf);md.melting_rate(pos)=10;50 pos=find(md.nodeoniceshelf);md.melting_rate(pos)=10; 51 51 52 52 %Friction … … 72 72 md.gl_migration=AgressiveMigrationEnum; 73 73 md.bathymetry=md.bed; 74 pos=find(md. gridoniceshelf); md.bathymetry(pos)=md.bed(pos)-900;74 pos=find(md.nodeoniceshelf); md.bathymetry(pos)=md.bed(pos)-900; 75 75 76 76 %Deal with boundary conditions:
Note:
See TracChangeset
for help on using the changeset viewer.