Changeset 25956 for issm/trunk-jpl/test/NightlyRun/test2021.m
- Timestamp:
- 01/27/21 16:14:34 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2021.m
r25428 r25956 1 %Test Name: SESAWsl r2 % SESAW method of solving GRD sl r1 %Test Name: SESAWslc 2 % SESAW method of solving GRD slc 3 3 % reference: Adhikari et al., 2016, GMD, https://doi.org/10.5194/gmd-9-1087-2016 4 4 … … 36 36 area_element=GetAreasSphericalTria(md.mesh.elements,md.mesh.lat,md.mesh.long,md.solidearth.planetradius); 37 37 38 % Parameters input for SESAWsl rsolver.38 % Parameters input for SESAWslc solver. 39 39 para.ocean_element = ocean_element; 40 40 para.loads_element = loads_element; … … 61 61 62 62 % solve: Rigid without rotational feedbacks. 63 disp(['Solving sesaw-sl rfor Rigid Earth WITHOUT rotational feedback...']);63 disp(['Solving sesaw-slc for Rigid Earth WITHOUT rotational feedback...']); 64 64 [eus_rigid,rsl_rigid] = SESAWslr(md.mesh.elements,md.mesh.lat,md.mesh.long,greens,para); 65 65 66 66 % solve: Rigid with rotational feedbacks. 67 67 para.rotational.flag = 1; 68 disp(['Solving sesaw-sl rfor Rigid Earth WITH rotational feedback...']);68 disp(['Solving sesaw-slc for Rigid Earth WITH rotational feedback...']); 69 69 [eus_rigid_rot,rsl_rigid_rot] = SESAWslr(md.mesh.elements,md.mesh.lat,md.mesh.long,greens,para); 70 70 71 71 % solve: Elastic with rotational feedbacks. 72 72 para.solidearth = 'elastic'; 73 disp(['Solving sesaw-sl rfor Elastic Earth WITH rotational feedback...']);73 disp(['Solving sesaw-slc for Elastic Earth WITH rotational feedback...']); 74 74 [eus_elast_rot,rsl_elast_rot] = SESAWslr(md.mesh.elements,md.mesh.lat,md.mesh.long,greens,para); 75 75 76 76 % solve: Elastic with rotational feedbacks. 77 77 para.rotational.flag = 0; 78 disp(['Solving sesaw-sl rfor Elastic Earth WITHOUT rotational feedback...']);78 disp(['Solving sesaw-slc for Elastic Earth WITHOUT rotational feedback...']); 79 79 [eus_elast,rsl_elast] = SESAWslr(md.mesh.elements,md.mesh.lat,md.mesh.long,greens,para); 80 80
Note:
See TracChangeset
for help on using the changeset viewer.