Changeset 24480


Ignore:
Timestamp:
12/18/19 13:59:09 (5 years ago)
Author:
tsantos
Message:

BUG: wrong computation of basal drag in x and in y directions

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/plot/plot_basaldrag.m

    r17686 r24480  
    2020r=averaging(md,md.friction.q./md.friction.p,0);
    2121
    22 %compute horizontal velocity
     22ub_mag=sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts;
     23drag_mag=(max(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base),0)).^r.*(md.friction.coefficient).^2.*ub_mag.^s/1000;
     24sig=1;
     25
     26%compute horizontal velocity
    2327if strcmpi(type,'basal_drag')
    2428        ub=sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts;
     29        title_str='Basal drag [kPa]';
    2530elseif strcmpi(type,'basal_dragx')
    2631        ub=md.initialization.vx/md.constants.yts;
     32        sig=-1;
     33        title_str='Basal drag - x direction [kPa]';
    2734elseif strcmpi(type,'basal_dragy')
    2835        ub=md.initialization.vy/md.constants.yts;
     36        sig=-1;
     37        title_str='Basal drag - y direction [kPa]';
    2938end
    3039
    3140%compute basal drag
    32 drag=(max(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base),0)).^r.*(md.friction.coefficient).^2.*ub.^s/1000;
     41drag=sig*drag_mag.*ub./ub_mag;
    3342
    3443%Figure out if this is a Section plot
     
    4756
    4857        %apply options
    49         options=addfielddefault(options,'title','Basal drag [kPa]');
     58        options=addfielddefault(options,'title',title_str);     
    5059        options=addfielddefault(options,'view',2);
    5160        applyoptions(md,basal_drag,options);
Note: See TracChangeset for help on using the changeset viewer.