![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 27251 for issm/trunk-jpl/src/m/classes/stochasticforcing.m
- Timestamp:
- 08/31/22 06:32:32 (3 years ago)
- File:
-
- 1 edited
-
issm/trunk-jpl/src/m/classes/stochasticforcing.m (modified) (10 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/stochasticforcing.m
r27165 r27251 52 52 for field=self.fields 53 53 %Checking agreement of classes 54 if(contains(field,'SMBa utoregression'))54 if(contains(field,'SMBarma')) 55 55 mdname = structstoch.mdnames(find(strcmp(structstoch.fields,char(field)))); 56 56 if~(isequal(class(md.smb),char(mdname))) … … 88 88 end 89 89 end 90 if(contains(field,'BasalforcingsDeepwaterMeltingRate Autoregression'))90 if(contains(field,'BasalforcingsDeepwaterMeltingRatearma')) 91 91 mdname = structstoch.mdnames(find(strcmp(structstoch.fields,char(field)))); 92 92 if~(isequal(class(md.basalforcings),char(mdname))) … … 111 111 end 112 112 %Checking for specific dimensions 113 if ~(strcmp(field,'SMBa utoregression') || strcmp(field,'FrontalForcingsRignotAutoregression') || strcmp(field,'BasalforcingsDeepwaterMeltingRateAutoregression'))113 if ~(strcmp(field,'SMBarma') || strcmp(field,'FrontalForcingsRignotarma') || strcmp(field,'BasalforcingsDeepwaterMeltingRatearma')) 114 114 checkdefaults = true; %field with non-specific dimensionality 115 115 end … … 118 118 %Retrieve all the field dimensionalities 119 119 dimensions = self.defaultdimension*ones(1,num_fields); 120 indSMBar = -1; %about to check for index of SMBautoregression121 indTFar = -1; %about to check for index of FrontalForcingsRignotAutoregression122 indBDWar = -1; %about to check for index of BasalforcingsDeepwaterMeltingRateAutoregression123 if any(contains(self.fields,'SMBa utoregression'))124 indSMBar = find(contains(self.fields,'SMBautoregression')); %index of SMBar, now check for consistency with other artimesteps125 dimensions(indSMBar ) = md.smb.num_basins;126 if(md.smb.ar _timestep<self.stochastictimestep)127 error('SMBa utoregressioncannot have a timestep shorter than stochastictimestep');128 end 129 end 130 if any(contains(self.fields,'FrontalForcingsRignot Autoregression'))131 indTFar = find(contains(self.fields,'FrontalForcingsRignotAutoregression')); %index of TFar, now check for consistency with other artimesteps132 dimensions(indTFar ) = md.frontalforcings.num_basins;133 if(md.frontalforcings.ar _timestep<self.stochastictimestep)134 error('FrontalForcingsRignot Autoregressioncannot have a timestep shorter than stochastictimestep');135 end 136 end 137 if any(contains(self.fields,'BasalforcingsDeepwaterMeltingRate Autoregression'))138 indBDWar = find(contains(self.fields,'BasalforcingsDeepwaterMeltingRateAutoregression')); %index of BDWar, now check for consistency with other artimesteps139 dimensions(indBDWar ) = md.basalforcings.num_basins;140 if(md.basalforcings.ar _timestep<self.stochastictimestep)141 error('BasalforcingsDeepwaterMeltingRate Autoregressioncannot have a timestep shorter than stochastictimestep');120 indSMBarma = -1; %about to check for index of SMBarma 121 indTFarma = -1; %about to check for index of FrontalForcingsRignotarma 122 indBDWarma = -1; %about to check for index of BasalforcingsDeepwaterMeltingRatearma 123 if any(contains(self.fields,'SMBarma')) 124 indSMBarma = find(contains(self.fields,'SMBarma')); %index of SMBarma, now check for consistency with other arma timesteps 125 dimensions(indSMBarma) = md.smb.num_basins; 126 if(md.smb.arma_timestep<self.stochastictimestep) 127 error('SMBarma cannot have a timestep shorter than stochastictimestep'); 128 end 129 end 130 if any(contains(self.fields,'FrontalForcingsRignotarma')) 131 indTFarma = find(contains(self.fields,'FrontalForcingsRignotarma')); %index of TFarma, now check for consistency with other arma timesteps 132 dimensions(indTFarma) = md.frontalforcings.num_basins; 133 if(md.frontalforcings.arma_timestep<self.stochastictimestep) 134 error('FrontalForcingsRignotarma cannot have a timestep shorter than stochastictimestep'); 135 end 136 end 137 if any(contains(self.fields,'BasalforcingsDeepwaterMeltingRatearma')) 138 indBDWarma = find(contains(self.fields,'BasalforcingsDeepwaterMeltingRatearma')); %index of BDWarma, now check for consistency with other arma timesteps 139 dimensions(indBDWarma) = md.basalforcings.num_basins; 140 if(md.basalforcings.arma_timestep<self.stochastictimestep) 141 error('BasalforcingsDeepwaterMeltingRatearma cannot have a timestep shorter than stochastictimestep'); 142 142 end 143 143 end 144 144 size_tot = sum(dimensions); 145 145 146 if(indSMBar ~=-1 && indTFar~=-1) %both autoregressive models are used: check autoregressivetime step consistency147 if(md.smb.ar _timestep~=md.frontalforcings.ar_timestep)148 crossentries = reshape(self.covariance(1+sum(dimensions(1:indSMBar -1)):sum(dimensions(1:indSMBar)),1+sum(dimensions(1:indTFar-1)):sum(dimensions(1:indTFar))),1,[]);146 if(indSMBarma~=-1 && indTFarma~=-1) %both ARMA models are used: check ARMA time step consistency 147 if(md.smb.arma_timestep~=md.frontalforcings.arma_timestep) 148 crossentries = reshape(self.covariance(1+sum(dimensions(1:indSMBarma-1)):sum(dimensions(1:indSMBarma)),1+sum(dimensions(1:indTFarma-1)):sum(dimensions(1:indTFarma))),1,[]); 149 149 if any(crossentries~=0) 150 error('SMBa utoregression and FrontalForcingsRignotAutoregression have different ar_timestep and non-zero covariance');151 end 152 end 153 end 154 if(indSMBar ~=-1 && indBDWar~=-1) %both autoregressive models are used: check autoregressivetime step consistency155 if(md.smb.ar _timestep~=md.basalforcings.ar_timestep)156 crossentries = reshape(self.covariance(1+sum(dimensions(1:indSMBar -1)):sum(dimensions(1:indSMBar)),1+sum(dimensions(1:indBDWar-1)):sum(dimensions(1:indBDWar))),1,[]);150 error('SMBarma and FrontalForcingsRignotarma have different arma_timestep and non-zero covariance'); 151 end 152 end 153 end 154 if(indSMBarma~=-1 && indBDWarma~=-1) %both ARMA models are used: check ARMA time step consistency 155 if(md.smb.arma_timestep~=md.basalforcings.arma_timestep) 156 crossentries = reshape(self.covariance(1+sum(dimensions(1:indSMBarma-1)):sum(dimensions(1:indSMBarma)),1+sum(dimensions(1:indBDWarma-1)):sum(dimensions(1:indBDWarma))),1,[]); 157 157 if any(crossentries~=0) 158 error('SMBa utoregression and BasalforcingsDeepwaterMeltingRateAutoregression have different ar_timestep and non-zero covariance');159 end 160 end 161 end 162 if(indTFar ~=-1 && indBDWar~=-1) %both autoregressive models are used: check autoregressivetime step consistency163 if(md.frontalforcings.ar _timestep~=md.basalforcings.ar_timestep)164 crossentries = reshape(self.covariance(1+sum(dimensions(1:indTFar -1)):sum(dimensions(1:indTFar)),1+sum(dimensions(1:indBDWar-1)):sum(dimensions(1:indBDWar))),1,[]);158 error('SMBarma and BasalforcingsDeepwaterMeltingRatearma have different arma_timestep and non-zero covariance'); 159 end 160 end 161 end 162 if(indTFarma~=-1 && indBDWarma~=-1) %both ARMA models are used: check ARMA time step consistency 163 if(md.frontalforcings.arma_timestep~=md.basalforcings.arma_timestep) 164 crossentries = reshape(self.covariance(1+sum(dimensions(1:indTFarma-1)):sum(dimensions(1:indTFarma)),1+sum(dimensions(1:indBDWarma-1)):sum(dimensions(1:indBDWarma))),1,[]); 165 165 if any(crossentries~=0) 166 error('FrontalForcingsRignot Autoregression and BasalforcingsDeepwaterMeltingRateAutoregression have different ar_timestep and non-zero covariance');166 error('FrontalForcingsRignotarma and BasalforcingsDeepwaterMeltingRatearma have different arma_timestep and non-zero covariance'); 167 167 end 168 168 end … … 181 181 disp(sprintf(' stochasticforcing parameters:')); 182 182 fielddisplay(self,'isstochasticforcing','is stochasticity activated?'); 183 fielddisplay(self,'fields','fields with stochasticity applied, ex: [{''SMBa utoregression''}], or [{''SMBforcing''},{''DefaultCalving''}]');183 fielddisplay(self,'fields','fields with stochasticity applied, ex: [{''SMBarma''}], or [{''SMBforcing''},{''DefaultCalving''}]'); 184 184 fielddisplay(self,'defaultdimension','dimensionality of the noise terms (does not apply to fields with their specific dimension)'); 185 185 fielddisplay(self,'default_id','id of each element for partitioning of the noise terms (does not apply to fields with their specific partition)'); … … 212 212 for field=self.fields 213 213 %Checking for specific dimensions 214 if(strcmp(field,'SMBa utoregression'))214 if(strcmp(field,'SMBarma')) 215 215 dimensions(ind) = md.smb.num_basins; 216 216 end 217 if(strcmp(field,'FrontalForcingsRignot Autoregression'))217 if(strcmp(field,'FrontalForcingsRignotarma')) 218 218 dimensions(ind) = md.frontalforcings.num_basins; 219 219 end 220 if(strcmp(field,'BasalforcingsDeepwaterMeltingRate Autoregression'))220 if(strcmp(field,'BasalforcingsDeepwaterMeltingRatearma')) 221 221 dimensions(ind) = md.basalforcings.num_basins; 222 222 end … … 225 225 226 226 %Scaling covariance matrix (scale column-by-column and row-by-row) 227 scaledfields = {'BasalforcingsDeepwaterMeltingRate Autoregression','BasalforcingsSpatialDeepwaterMeltingRate','DefaultCalving','FloatingMeltRate','SMBautoregression','SMBforcing'}; %list of fields that need scaling *1/yts227 scaledfields = {'BasalforcingsDeepwaterMeltingRatearma','BasalforcingsSpatialDeepwaterMeltingRate','DefaultCalving','FloatingMeltRate','SMBarma','SMBforcing'}; %list of fields that need scaling *1/yts 228 228 tempcovariance = self.covariance; %copy of covariance to avoid writing back in member variable 229 229 for i=1:num_fields … … 266 266 % supported and corresponding md names 267 267 structure.fields = {... 268 'BasalforcingsDeepwaterMeltingRate Autoregression',...268 'BasalforcingsDeepwaterMeltingRatearma',... 269 269 'BasalforcingsSpatialDeepwaterMeltingRate',... 270 270 'DefaultCalving',... … … 273 273 'FrictionCoulombWaterPressure',... 274 274 'FrictionSchoofWaterPressure',... 275 'FrontalForcingsRignot Autoregression',...276 'SMBa utoregression',...275 'FrontalForcingsRignotarma',... 276 'SMBarma',... 277 277 'SMBforcing' 278 278 }; 279 279 structure.mdnames = {... 280 ' autoregressionlinearbasalforcings',...280 'linearbasalforcingsarma',... 281 281 'spatiallinearbasalforcings',... 282 282 'calving',... … … 285 285 'frictioncoulomb',... 286 286 'frictionschoof',... 287 'frontalforcingsrignota utoregression',...288 'SMBa utoregression',...287 'frontalforcingsrignotarma',... 288 'SMBarma',... 289 289 'SMBforcing' 290 290 };
Note:
See TracChangeset
for help on using the changeset viewer.
