- Timestamp:
- 05/20/09 09:51:48 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/test/Verification/IceShelfIceFrontM3d_2/runme.m
r489 r519 45 45 md=extrude(md,8,4); 46 46 md=setelementstype(md,'macayeal','all'); 47 if strcmpi(package,'cielo_parallel'), md.cluster='wilkes'; end48 47 if md.numberofgrids==832 49 48 load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure; … … 51 50 52 51 %compute solution 53 if (strcmpi(package,'cielo_serial') | strcmpi(package,'cielo_parallel')), 54 md=solve(md,'analysis_type',solution,'package','cielo'); 55 else 56 md=solve(md,'analysis_type',solution,'package',package); 57 end 52 [analysis_type sub_analysis_type]=testsgetanalysis(solution); 53 [md packagefinal]=testsgetpackage(md,package); 54 md=solve(md,'analysis_type',analysis_type,'sub_analysis_type',sub_analysis_type,'package',package); 58 55 59 %compute field to be checked 60 if strcmpi(solution,'diagnostic'), 61 fields={'vy','vz'}; 62 elseif strcmpi(solution,'thermalsteady'), 63 fields={'temperature','melting'}; 64 elseif strcmpi(solution,'thermaltransient'), 65 fields={'thermaltransient_results(end).temperature','thermaltransient_results(end).melting'}; 66 elseif strcmpi(solution,'prognostic'), 67 fields={'new_thickness'}; 68 elseif strcmpi(solution,'transient'), 69 fields={'transient_results(end).vy','transient_results(end).vz','transient_results(end).temperature','transient_results(end).melting','transient_results(end).thickness'}; 70 end 56 %compute fields to be checked 57 fields=testsgetfields(md.type,solution); 71 58 72 59 %load archive
Note:
See TracChangeset
for help on using the changeset viewer.