Parent Directory | Revision Log | Revision Graph
Links to HEAD: | (view) (download) (annotate) |
Links to checkpoint64q: | (view) (download) (annotate) |
Sticky Tag: |
add SEAICE_SIZE.h anywhere SEAICE.h is included
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 |