Parent Directory
| 
Revision Log
| 
Revision Graph
| Links to HEAD: | (view) (download) (annotate) |
| Links to checkpoint59e: | (view) (download) (annotate) |
| Sticky Tag: |
fix a small bux: replace one Oly by Olx (should not have any effect unless Oly.NE.Olx for some reason)
Re-arrange store directives.
- make lsr boundary conditions free slip (by default, changes results of
lsr solver)
- no slip boundary conditions for lsr will follow, currently,
SEAICE_no_slip = .true., restores the old behavior (not good, but
only temporarily)
- set minimum of zeta to zero (and not 4e8) (changes results)
- remove WINDX & WINDY from C-grid dynamics version.
o Remove recomputations after latest modifs o EVP TLM and ADM gradient checks completely wrong (1.E34)
remove '#include SEAICE_FFIELDS.h'
o add new interface for atmospheric wind forcing (C-grid only) - new routine - two new field (taux, tauy for stress over ice) - removed one global field (DAIRN) for C-grid case only - small clean up - even remembered to update seaice_ad_diff.list
fix out-of-bound index problem.
if ALLOW_ATM_WIND is undefined, use the wind stress directly to drive the ice model
Change store.
seaice adjoint, part 2 (SEAICE_ALLOW_EVP).
apply the correct exchange (u,v point) to seaiceMaskU,seaiceMaskV
ice-dynamics forcing: add atmospheric and Sea-Ice loading to sea-surface slope
Avoid recomputation call to seaice_advdiff to circumvene TAF bug (1.8.42)
First stab at SEAICE_CGRID.
seaice: add an EVP solver following Hunke and Dukowicz (1997) and the documentation of CICE. Turn on by defining SEAICE_ALLOW_EVP in SEAICE_OPTIONS.h and SEAICEuseEVP=.true. in data.seaice. Works only with SEAICE_CGRID defined. Use at own risk.
fix turning angles to work on both hemispheres
add flags for clipping ice velocities and/or masking solver forcing
- added a few comments to seaice_diffusion.F
- separate C-grid and B-grid versions more cleanly (UVM and seaiceMaskU/V
are now exclusive)
- add a new ocean-ice stress coupling method (Hibler and Bryan, 1987) which
required re-ordering the code a little:
+ make DAIRN global variable (defined in SEAICE)
+ move computation of the viscosities eta and zeta into a new separate
routine (because I need to recompute them in seaice_ocean_stress.F)
o fix bug that caused problems in multiple tile/processor runs: extend loop ranges o remove velocity clipping from seaice_dynsolver
o continue fixing ice-ocean stresses
o another bug fix
- add c-grid version of the dynamic solver:
+ three new routines that are the c-grid substitute lsr, dynsolver
and ostres
+ put a few fields that were local to dynsolver into global common
blocks, so that I can move the computation of stresses etc into
seaice_lsr (saves coding but may break the adjoint; Patrick, I am
sorry!).
- replace more hardwired parameters by runtime parameters
- add ice masks that mask the rhs of the implicit solvers where there
is no ice (commented out in seaice_dynsolver, because i am not sure
if this works properly), eventually this should replace the clipping
of ice velocities in seaice_dynsolver to +/-40cm/s.
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 |