Changeset 8354


Ignore:
Timestamp:
05/19/11 11:31:19 (14 years ago)
Author:
Mathieu Morlighem
Message:

trunk: fixed plot overlay when radarpower already present

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/plot/plot_overlay.m

    r8298 r8354  
    2626end
    2727
    28 %Get xlim and ylim (used to extract radar image)
    29 xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
    30 ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
    31 xlim2=xlim/getfieldvalue(options,'unit',1);
    32 ylim2=ylim/getfieldvalue(options,'unit',1);
    3328
    3429%radar power
    3530if ~any(isnan(md.sarxm)) & ~any(isnan(md.sarym)) & ~any(isnan(md.sarpwr)),
    3631        disp('plot_overlay info: the radar image held by the model is being used');
     32        xlim=[min(md.sarxm) max(md.sarxm)];
     33        ylim=[min(md.sarym) max(md.sarym)];
    3734else
    3835        disp('Extracting radar image...');
    39         options=addfielddefault(options,'xlim',xlim2);
    40         options=addfielddefault(options,'ylim',ylim2);
     36        %Get xlim and ylim (used to extract radar image)
     37        xlim=getfieldvalue(options,'xlim',[min(x) max(x)])/getfieldvalue(options,'unit',1);
     38        ylim=getfieldvalue(options,'ylim',[min(y) max(y)])/getfieldvalue(options,'unit',1);
     39        options=addfielddefault(options,'xlim',xlim);
     40        options=addfielddefault(options,'ylim',ylim);
    4141        md=radarpower(md,options);
    4242end
Note: See TracChangeset for help on using the changeset viewer.