Changeset 26047 for issm/trunk-jpl/src/m/classes/sealevelmodel.m

Timestamp:

03/08/21 17:46:39 (5 years ago)

Author:

Eric.Larour

Message:

CHG: huge commit on solid earth capability rewrite. Complete cleanup of the sea level core.
New mass transport capabilities for ocean and tws. No more giacore. GiaIvins folded into
the sea level core. Debugging of Materials.

File:

• issm/trunk-jpl/src/m/classes/sealevelmodel.m (modified) (1 diff)

issm/trunk-jpl/src/m/classes/sealevelmodel.m

@@ -79 +79 @@
                                 end
                         end
+
+                        %make sure grd is the same everywhere:
+                        for i=1:length(slm.icecaps),
+                                md= slm.icecaps{i}; 
+                                if md.solidearthsettings.isgrd~=slm.earth.solidearthsettings.isgrd),
+                                        error(sprintf('isgrd on ice cap %s is not the same as for the earth\n',md.miscellaneous.name));
+                                end
+                        end
+                        %make sure that there is no solid earth external forcing on the basins: 
+                        for i=1:length(slm.icecaps),
+                                md= slm.icecaps{i}; 
+                                if ~isempty(md.solidearth.external),
+                                        error(sprintf('cannot run external forcings on an ice sheet when runing a coupling earth/ice sheet model');
+                                end
+
+                        end
+                        %make sure that we have the rigth grd model for computing ouf sealevel patterns: 
+                        for i=1:length(slm.icecaps),
+                                md= slm.icecaps{i}; 
+                                if md.solidearth.settings.grdmodel~=0
+                                        error(sprintf('sealevelmodel checkconsistenty error message: ice sheets do not run GRD module, specify solidearth.settings.grdmodel=0 on ice cap %i',i));
+                                end
+                        end
+
 
                 end