Changeset 26751


Ignore:
Timestamp:
12/28/21 07:13:59 (3 years ago)
Author:
vverjans
Message:

CHG: modified number of spin-up steps for SMBautoregression and frontalforcingsrignotautoregression

Location:
issm/trunk-jpl/src
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/modules/FrontalForcingsx/FrontalForcingsx.cpp

    r26615 r26751  
    4646   femmodel->parameters->FindParam(&phi,&M,&Nphi,FrontalForcingsPhiEnum);  _assert_(M==numbasins); _assert_(Nphi==arorder);
    4747
    48    /*AR model spin-up with 0 noise to initialize ThermalforcingValuesAutoregressionEnum*/
    49         int nspin{2*arorder+5};
     48   /*AR model spin-up with 0 noise to initialize ThermalforcingValuesAutoregressionEnum (688 = log(0.001)/log(0.99): decaying time of inluence of phi[0]=0.99 to 0.001 of beta_0*/
     49        int nspin = 688;
    5050   for(Object* &object:femmodel->elements->objects){
    5151      Element* element      = xDynamicCast<Element*>(object); //generate element object

  • issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp

    r26615 r26751  
    164164        femmodel->parameters->FindParam(&phi,&M,&Nphi,SmbPhiEnum);  _assert_(M==numbasins); _assert_(Nphi==arorder);
    165165       
    166         /*AR model spin-up with 0 noise to initialize SmbValuesAutoregressionEnum*/
    167         int nspin{2*arorder+5};
     166        /*AR model spin-up with 0 noise to initialize SmbValuesAutoregressionEnum (688 = log(0.001)/log(0.99): decaying time of inluence of phi[0]=0.99 to 0.001 of beta_0*/
     167        int nspin = 688;
    168168        for(Object* &object:femmodel->elements->objects){
    169169      Element* element      = xDynamicCast<Element*>(object); //generate element object

  • issm/trunk-jpl/src/m/classes/SMBautoregression.m

    r26526 r26751  
    7979                        fielddisplay(self,'num_basins','number of different basins [unitless]');
    8080                        fielddisplay(self,'basin_id','basin number assigned to each element [unitless]');
    81                         fielddisplay(self,'beta0','basin-specific intercept values [m ice eq./yr]');
     81                        fielddisplay(self,'beta0','basin-specific intercept values [m ice eq./yr] (if beta_1==0 mean=beta_0/(1-sum(phi)))');
    8282                        fielddisplay(self,'beta1','basin-specific trend values [m ice eq. yr^(-2)]');
    8383                        fielddisplay(self,'ar_order','order of the autoregressive model [unitless]');

  • issm/trunk-jpl/src/m/classes/SMBautoregression.py

    r26553 r26751  
    3838        s += '{}\n'.format(fielddisplay(self, 'num_basins', 'number of different basins [unitless]'))
    3939        s += '{}\n'.format(fielddisplay(self, 'basin_id', 'basin number assigned to each element [unitless]'))
    40         s += '{}\n'.format(fielddisplay(self, 'beta0', 'basin-specific intercept values [m ice eq./yr]'))
     40        s += '{}\n'.format(fielddisplay(self, 'beta0', 'basin-specific intercept values [m ice eq./yr] (if beta_1==0 mean=beta_0/(1-sum(phi)))'))
    4141        s += '{}\n'.format(fielddisplay(self, 'beta1', 'basin-specific trend values [m ice eq. yr^(-2)]'))
    4242        s += '{}\n'.format(fielddisplay(self, 'ar_order', 'order of the autoregressive model [unitless]'))

  • issm/trunk-jpl/src/m/classes/frontalforcingsrignotautoregression.m

    r26553 r26751  
    6666         fielddisplay(self,'basin_id','basin number assigned to each element [unitless]');
    6767         fielddisplay(self,'subglacial_discharge','sum of subglacial discharge for each basin [m/d]');
    68          fielddisplay(self,'beta0','basin-specific intercept values [∘C]');
     68         fielddisplay(self,'beta0','basin-specific intercept values [∘C] (if beta_1==0 mean=beta_0/(1-sum(phi)))');
    6969         fielddisplay(self,'beta1','basin-specific trend values [∘C yr^(-1)]');
    7070         fielddisplay(self,'ar_order','order of the autoregressive model [unitless]');

  • issm/trunk-jpl/src/m/classes/frontalforcingsrignotautoregression.py

    r26553 r26751  
    3636        s += '{}\n'.format(fielddisplay(self, 'basin_id', 'basin number assigned to each element [unitless]'))
    3737        s += '{}\n'.format(fielddisplay(self, 'subglacial_discharge', 'sum of subglacial discharge for each basin [m/d]'))
    38         s += '{}\n'.format(fielddisplay(self, 'beta0', 'basin-specific intercept values [°C]'))
     38        s += '{}\n'.format(fielddisplay(self, 'beta0', 'basin-specific intercept values [°C] (if beta_1==0 mean=beta_0/(1-sum(phi)))'))
    3939        s += '{}\n'.format(fielddisplay(self, 'beta1', 'basin-specific trend values [°C yr^(-1)]'))
    4040        s += '{}\n'.format(fielddisplay(self, 'ar_order', 'order of the autoregressive model [unitless]'))
Note: See TracChangeset for help on using the changeset viewer.