Log of /MITgcm/pkg/seaice/seaice_ocean_stress.F

---

Revision 1.10 - (view) (download) (annotate) - [select for diffs]
Thu Mar 8 11:21:34 2007 UTC (17 years, 1 month ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58w_post, checkpoint58x_post, checkpoint58y_post
Changes since 1.9: +16 -1 lines
Diff to previous 1.9

if ALLOW_ATM_WIND is undefined, use the wind stress directly to drive
the ice model

---

Revision 1.9 - (view) (download) (annotate) - [select for diffs]
Wed Jun 14 16:08:36 2006 UTC (17 years, 10 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58i_post, checkpoint58j_post, checkpoint58k_post, checkpoint58l_post, checkpoint58m_post, checkpoint58n_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58r_post, checkpoint58s_post, checkpoint58t_post, checkpoint58u_post, checkpoint58v_post
Changes since 1.8: +9 -3 lines
Diff to previous 1.8

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

---

Revision 1.8 - (view) (download) (annotate) - [select for diffs]
Mon Mar 20 21:36:11 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58c_post, checkpoint58d_post, checkpoint58e_post, checkpoint58f_post, checkpoint58g_post, checkpoint58h_post
Changes since 1.7: +4 -1 lines
Diff to previous 1.7

  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.

---

Revision 1.7 - (view) (download) (annotate) - [select for diffs]
Fri Mar 17 15:53:38 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
Changes since 1.6: +5 -5 lines
Diff to previous 1.6

- fix factor 0.5 for pressure term in seaice_lsr
- small fix for ice-ocean stress (Hibler and Bryan)

---

Revision 1.6 - (view) (download) (annotate) - [select for diffs]
Thu Mar 16 14:25:40 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
Changes since 1.5: +27 -19 lines
Diff to previous 1.5

fix the stress one more time (and fix turning angles to work on both
hemispheres)

---

Revision 1.5 - (view) (download) (annotate) - [select for diffs]
Wed Mar 15 19:49:04 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
Changes since 1.4: +191 -5 lines
Diff to previous 1.4

  - 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)

---

Revision 1.4 - (view) (download) (annotate) - [select for diffs]
Thu Mar 9 20:22:40 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58b_post
Changes since 1.3: +6 -5 lines
Diff to previous 1.3

o another fix for the ocean stress: use weighted average with weights
  derived from AREA interpolated to U and V points

---

Revision 1.3 - (view) (download) (annotate) - [select for diffs]
Mon Mar 6 20:53:25 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
Changes since 1.2: +7 -54 lines
Diff to previous 1.2

o fix ice-ocean stresses

---

Revision 1.2 - (view) (download) (annotate) - [select for diffs]
Mon Mar 6 13:27:40 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN
Changes since 1.1: +3 -3 lines
Diff to previous 1.1

small bug fix

---

Revision 1.1 - (view) (download) (annotate) - [select for diffs]
Mon Mar 6 13:17:37 2006 UTC (18 years, 1 month ago) by mlosch
Branch: MAIN

 - 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.