Changeset 8454
- Timestamp:
- 05/31/11 15:08:44 (14 years ago)
- Location:
- issm/trunk/src/m/model
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/parameterization/parametercontrolB.m
r8298 r8454 39 39 end 40 40 41 41 42 %cm_min 42 cm_min=getfieldvalue(options,'cm_min',1); 43 if (length(cm_min)~=1) 44 md.cm_min=1; 43 cm_min=getfieldvalue(options,'cm_min',paterson(273.15+5)*ones(md.numberofnodes,1)); 44 if (length(cm_min)==1) 45 md.cm_min=cm_min*ones(md.numberofnodes,1); 46 elseif (length(cm_min)==md.numberofnodes) 47 md.cm_min=cm_min; 45 48 else 46 49 md.cm_min=cm_min; … … 48 51 49 52 %cm_max 50 cm_max=getfieldvalue(options,'cm_max',10^10); 51 if (length(cm_max)~=1) 52 md.cm_max=10^10; 53 cm_max=getfieldvalue(options,'cm_max',paterson(273.15-70)*ones(md.numberofnodes,1)); 54 if (length(cm_max)==1) 55 md.cm_max=cm_max*ones(md.numberofnodes,1); 56 elseif (length(cm_max)==md.numberofnodes) 57 md.cm_max=cm_max; 53 58 else 54 59 md.cm_max=cm_max; -
issm/trunk/src/m/model/plot/plot_basaldrag.m
r6340 r8454 30 30 31 31 %compute basal drag 32 drag=(m d.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag_coefficient).^2.*ub.^s/1000;32 drag=(max(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed),0)).^r.*(md.drag_coefficient).^2.*ub.^s/1000; 33 33 34 34 %Figure out if this is a Section plot
Note:
See TracChangeset
for help on using the changeset viewer.