Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint61t: | (view) (download) (annotate) |
Sticky Tag: |
Adding diagnostics SHIForcT and SHIForcS, the ice shelf forcing for T and S.
fix bug with SHELFICEconserve flag, prepare for real fresh water flux (not yet) this simplifies the code considerable, but does not have too large an effect and the model solutions
- convert potential temperature into in-situ temperature. This changes the isomip tests => update tests
- rearrange computation of tLoc, sLoc, pLoc to make code more efficient
- sort out freshwater flux; convert freshwater fluxes to mass fluxes (kg/m^2/s), does not break tests, because in verification/isomip, rhoConstFresh=rhoConst (by accident)
clean up recent check-in with too comments only meant for myself
how embarrassing: total mixup of ice and sea water densities is resolved
clean up a little
add option for convservative form of Hellmer&Olbers (1989) thermodynamics - requires a little reorganization, which affects the testreport results of isomip.htd (only 7 digits of cg2d agree) => update experiment The results are extremely sensitive to details of the algorithm (e.g., order in which sums and differences are taken, etc) so that replacing a term with a (theoretically) exactly equal other term causes the big difference; possibly, because the solution of a quadratic equation is involed. I am not happy about this, but the present formulation appears to be a little more robust. - add a few comments
Change EmPmR units: from m/s to kg/m2/s
replace "recip_rhoConst*horiVertRatio" by mass2rUnit and "rhoConst*recip_horiVertRatio" by rUnit2mass
add a few diagnostics
add a simple boundary layer scheme to reduce noise, off by default for now
replace a '<' by a '.LT.'
- add more sophisticated thermodynamics of Hellmer (1989), this is now the default - give a few constants more intuitive names
- separate shelfice load anomaly from pload (breaks with time dependent forcing), => introduce constant field shelficeLoadAnomaly. Its default is 0., but it may be computed more cleverly from (unknown) t- and s-profiles (tRef, sRef) and actual EOS. For now this has to be done offline. A good approximation of the pressure load anomaly is necessary to avoid large initial adjustment processes underneath deep-reaching shelfice. - small fix in shelfice_thermodynamics
o add new package shelfice
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.
ViewVC Help | |
Powered by ViewVC 1.1.22 |