Parent Directory
| 
Revision Log
| 
Revision Graph
| Links to HEAD: | (view) (download) (annotate) |
| Links to checkpoint62r: | (view) (download) (annotate) |
| Sticky Tag: |
forgot to add a comment explaining this absolutely absurd code change
split a loop for better vectorization of the adjoin code
fix argument type of intrisic fct "SIGN" when REAL4_IS_SLOW is undef
cosmetic change: replace uVel(i,j,1,bi,bj) with uVel(i,j,kSrf,bi,bj) etc.
third and step of replacing 3D versions of UICE,VICE,HEFF,AREA by 2D versions. Keep your fingers crossed.
- clean up computation of Hibler+Bryan (1987) stress coupling for the case
of LSR (this change is expected to change the results slightly because now
it uses slightly difference moduli from the second last LSR solution, but
that is more consistent with the stress computations; this part of the code
is not tested in the verification experiments)
- turn strain rates eij and press into global fields within global
common block in SEAICE.h for more straighforward diagnostics
removed superfluous many metric terms of FV discretization of stress divergence for non-evp solver (this part is not tested anywhere)
removed uneeded arrays KGEO, GWATX, and GWATY from SEAICE_CGRID code
fix bug in seaice_ocean_stress (call of seaice_calc_strainrates), requires change of parameter list of seaice_calc_strainrates
add arguments myTime, myIter & iStep (= sub-time-step) (easier for debugging) to S/R SEAICE_CALC_STRAINRATES & SEAICE_CALC_VISCOSITIES
- get rid of three more exchanges seaice_calc_viscosities (required
different loop boundaries in seaice_lsr and seaice_ocean_stress)
- bug fix in seaice_ocean_stress (update lab_sea)
- reduce loop boundaries to necessary minimum also in seaice_evp
- go back to the old way of computing deltaC/Z as the adjoint likes that
better
- remove WINDX & WINDY from C-grid dynamics version.
one day we will get it right: replace ALLOW_EVP with SEAICE_ALLOW_EVP
Forward code broken: stressDivergence declared in SEAICE.h only for SEAICE_ALLOW_EVP
- modify Hibler+Bryan (1987) stress coupling so that it reuses information
of evp, complete rewrite for lsr case (latter is does not yet work
in all cases)
fix bug in evp solver, while doing that
- change (improve?) discretization of stress computations for both
cgrid-lsr and evp solver
- add a new routine to compute strain rates
- changes lab_sea results
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
if ALLOW_ATM_WIND is undefined, use the wind stress directly to drive the ice model
add parameter that allows to scale/turn off effect of seaice on wind stress (SEAICEstressFactor defaults to 1.); with SEAICEuseDynamics, SEAICEstressFactor = 0. (data.seaice) and stressReduction = 0. (data.ice), thsice gives the same results regardless of the state of useSEAICE
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 factor 0.5 for pressure term in seaice_lsr - small fix for ice-ocean stress (Hibler and Bryan)
fix the stress one more time (and fix turning angles to work on both hemispheres)
- 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 another fix for the ocean stress: use weighted average with weights derived from AREA interpolated to U and V points
o fix ice-ocean stresses
small 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 |