Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint64o: | (view) (download) (annotate) |
Sticky Tag: |
Fix quotes
fix a CPP flag (EVP->JFNK) (found by Jean-Michel) fix treatment of JFNK and AUTODIFF (pointed out by Jean-Michel) add a few short comments and references (on my own account)
add hooks for JFNK-solver
- fix ice strength assymptote at SEAICE_area_max<1 (no change for the standard case of SEAICE_area_max=1)
- added possibility to change ice strength/PRESS0 relationship to HEFF, according to the following run time parameters C SEAICEpresPow0 :: HEFF exponent for ice strength below SEAICEpresH0 (positive integer) C SEAICEpresPow1 :: HEFF exponent for ice strength above SEAICEpresH0 (positive integer) C SEAICEpresH0 :: HEFF threshold for ice strength (m)
- added SEAICEuseTILT run time switch : if set to false then omit the tilt term in seaice dynamics.
add SEAICE_SIZE.h anywhere SEAICE.h is included
More prepares.
Add some basic init.s to break TAF dependencies.
- fix SEAICEuseFREEDRIFT for new OBCS version (ndef OBCS_UVICE_OLD)
- add OBCS_APPLY_UVICE call (new implenetation); - disable SEAICE_maskRHS with ALLOW_AUTODIFF
- reduce memory footprint of seaice mutldim advection tapes. - move lsr iterations from seaice_dynsolver to seaice_lsr.
INTRODUCING SUBROUTINE SEAICE_FREEDRIFT o Solve ice approximate momentum equation analytically.
bug fix: the sea ice UV mask update in seaice_dynsolver.F needs to account for the land mask (as in seaice_init_varia.F)
avoiding recomputations in seaice dynamics adjoint
make maximum zeta a runtime parameter
- add CPP-flag SEAICE_ALLOW_CLIPVELS (undefined by default)
o initializations of *nm1 to break dependencies o adjust some store directives
third and step of replacing 3D versions of UICE,VICE,HEFF,AREA by 2D versions. Keep your fingers crossed.
Invoke MPSEUDOTIMESTEPS only for adjoint.
Adapt store directives to new pseudotimestep loop.
add pseudo time stepping for LSOR
Extra store.
Update tamc.h for single-prec comlev option
Add comment.
Looks like re-init. is working properly for forward, so we can restrict it to adjoint.
o Update stores after recent IF (DIFFERENT_MULTIPLE ... o re-initialise to break some dependencies (not sure if all correct) o turn some scalares into tw-dim.
removed uneeded arrays KGEO, GWATX, and GWATY from SEAICE_CGRID code
o pkg/seaice: moved DIFFERENT_MULTIPLE from seaice_model to seaice_dynsolver http://forge.csail.mit.edu/pipermail/mitgcm-devel/2008-January/003173.html
add arguments myTime, myIter & iStep (= sub-time-step) (easier for debugging) to S/R SEAICE_CALC_STRAINRATES & SEAICE_CALC_VISCOSITIES
replace hardcoded lower limit for zeta (zMin) by runtime parameter
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 |