Open Hub
Open Hub

 Settings |  Report Duplicate

1
I Use This!
Very Low Activity

Commits : Listings

Analyzed 14 days ago. based on code collected 14 days ago.
Jan 29, 2024 — Jan 29, 2025
Commit Message Contributor Files Modified Lines Added Lines Removed Code Location Date
Most work done on Test K. Basically works; try 'pismv -test K -d HTt -Mx 6 -My 6 -Mz 401 -Mbz 101 -y 135000 -maxdt 100.0' and note that the exact solution first has z=0 ice hit pressure-melting at about 133465 years. Still needs proper error reporting and also documentation in User's Manual.
bueler
as Ed Bueler
More... over 17 years ago
Added a tool (pism/test/series.py) for turning PISM standard output into time series in a NetCDF file.
bueler
as Ed Bueler
More... over 17 years ago
Several improvements: 1. Improved transfer of SIA velocities from staggered to regular grid. 2. Viewers for basal velocity (both -d and -matv). 3. Clarified save of ubar,vbar in preparation for SSA calculations. 4. Cleaned up IceRYANModel.
bueler
as Ed Bueler
More... over 17 years ago
Fixed -constant_nu option so it should work.
bueler
as Ed Bueler
More... over 17 years ago
Bug fix: dbdt term removed from incompressibility; it is incorrect when there is no mass balance. For now it is removed; it must be fixed up for long term.
bueler
as Ed Bueler
More... over 17 years ago
Fixed up initialization in legacy bootstrap. Removed initSounding() as unnecessary.
bueler
as Ed Bueler
More... over 17 years ago
I CAUGHT A BIG BUG. Or at least a long-lived bug. It was in putTempAtDepth() in iMIOnetcdf.cc and it meant that when using thermal bedrock in realistic circumstances we were making big errors (in the whole thermocoupled flow). It is now clear that much of what we diagnosed as surface roughness was this error.
bueler
as Ed Bueler
More... over 17 years ago
Fixing up -bif_legacy. Also split long iMutil.cc; now iMreport.cc contains methods related to runtime outputs to standard out, esp. summary().
bueler
as Ed Bueler
More... over 17 years ago
Cleaned up and completed the Matlab-readable views; nearly totally matches runtime graphical viewers user interface.
bueler
as Ed Bueler
More... over 17 years ago
Renovation of diagnostic viewers code. Added IceModel::tn[], an array of names and titles. Add -matv option. Now the diagnostic viewers and the Matlab output more or less match in the user interface.
bueler
as Ed Bueler
More... over 17 years ago
Work on documentation. Work on variable names (when written/read from NetCDF or matlab) and short titles (e.g. on viewers). Should greatly ease conversion to CISM/GLIMMER naming conventions.
bueler
as Ed Bueler
More... over 17 years ago
Working on documentation (i.e. comments which doxygen turns into pism/doc/refman.pdf and a source browser). Need to back up now as I am working outside and around kids ...
bueler
as Ed Bueler
More... over 17 years ago
EISMINT II experiments I,J,K,L in better shape. Experiment I should be a basis for plastic till demonstration but still not getting it quite right; why? At least all signs of asymmetry are gone.
bueler
as Ed Bueler
More... over 17 years ago
Test J now converges at O(dx^2) ON THE DOT (power is 2.0002 experimentally). THERE WAS NO FLAW IN JED'S macayeal equations code. Rather, test J truely uses a periodic grid and thus rescale() was changed to optionally accomodate that. Also L changed to 300 km in test J.
bueler
as Ed Bueler
More... over 17 years ago
Preliminary completion of test J. Test J is a linearized periodic non-constant thickness ice shelf. (By 'linearized' and 'non-constant thickness' I mean that nu(x,y) H(x,y) = nu_0 H_0 even though H(x,y) is not constant.) The real significance is that the SSA differential operator (i.e. the left side of the 'MacAyeal' equations) is so well understood that one can know its exact eigenvalues, and thus apply iterative linear algebra with confidence.
bueler
as Ed Bueler
More... over 17 years ago
Finally rewrote the vertical integral that gets vertical velocity from incompressibility. Better documentation of same. Now back to test J.
bueler
as Ed Bueler
More... over 17 years ago
Work on implementing exact solution test J, which is a Fourier series solution to the SSA equations for an ice shelf in a linearized (and constant coefficient) case with periodic boundary conditions.
bueler
as Ed Bueler
More... over 17 years ago
Forgot to add src/base/iMsia.cc. Also fiddled with Matlab output from SSA solution. Working on bug #9836 by building new nontrivial exact ice shelf solution (found by Fourier series soln). Fiddles with Makefile; no significant changes though.
bueler
as Ed Bueler
More... over 17 years ago
In search of bug #9836. The pervasive 'MacAyeal equations' language has been replaced by 'shallow shelf approximation' (SSA) which is Weis-Greve-Hutter language. I have added comments to iMssa.cc (formerly iMmacayeal.cc) to indicate what continuum equations correspond to what routines.
bueler
as Ed Bueler
More... over 17 years ago
In writing subsection 3.3 of the User's Manual a asymmetry of the numerical solution of the MacAyeal equations (actually Schoof plastic till case) was found. See bug #9836. Note all of EISMINT II is now in pism, including unpublished experiments.
bueler
as Ed Bueler
More... over 17 years ago
Fixed a serious bug: EISMINT II experiments ABCDEF were all sliding on the base because of the renaming basal() --> basalVelocity(). This is fixed; in particular the main example in section 2 (Getting Started) now produces the usual spokes and not more extreme ones. Also: trivial changes to document build and on-going User's Manual improvements.
bueler
as Ed Bueler
More... over 17 years ago
Start on using doxygen to automatically generate continuum model and numerical methods documentation (PDF); will also generate a class structure browser in (HTML). Read doc/README. Added EISMINT II heterodox experiments I,J,K,L; note experiment I will be used to document setting of mask and use of -mv and -plastic options. Renamed test/5k to test/green5k; these are utilities for reading NSIDC and PARCA Greenland data into NetCDF.
bueler
as Ed Bueler
More... over 17 years ago
Added some directory structure to the test/ directory. Completed the redo of EISMINT-Ross: it is now parallel and at arbitrary resolution. See section 8.2 of the manual as it is pretty cool. Additional modifications to the manual.
bueler
as Ed Bueler
More... over 17 years ago
So I did it: Directory structure added to src/ directory. (scons is hosed again; sorry) Also pismd -ross can now compute chi^2 against RIGGS data, and a hardness tuning script is included.
bueler
as Ed Bueler
More... over 17 years ago
Further work on pismd -ross and documentation thereof. Generalized the forcing.[hh|cc] mechanism; added class Data1D of which IceSheetForcing is now a derived class. Modified constants DEFAULT_MAX_HMELT and a coefficient in excessToFromBasalMeltLayer() to represent a more conservative, slower production-of-basal-melt-water. (How should this mechanism be constrained, anyway?) Additional revisions to manual.
bueler
as Ed Bueler
More... over 17 years ago
More work on pismd -ross and related matters. Documentation thereof.
bueler
as Ed Bueler
More... over 17 years ago
More work on EISMINT-Ross and the general mechanism for diagnostic computations of velocity. Basically works under pismd; see comment at start of src/pismd.cc.
bueler
as Ed Bueler
More... over 17 years ago
1. New PISM executable "pismd"; "d" for diagnostic. Purpose is to compute and save full velocity field from given geometry and temperature fields. May be used either with -if or -bif. Saves state with full 3D velocity field; see #3 below.
bueler
as Ed Bueler
More... over 17 years ago
fixed forcing bug with "size_t" data type. added accumulation to regridding (since it was forgotten). Other minor changes
Nathan Shemonski
More... over 17 years ago
Oops, minor bug fix.
Jed Brown
More... over 17 years ago