Opened 11 years ago

Closed 11 years ago

Last modified 3 years ago

#587 closed Enhancement (fixed)

tilt in ssh in lim2 rheology

Reported by: gm Owned by: nemo
Priority: low Milestone:
Component: LIM2 Version: release-3.2
Severity: Keywords: 2010 LIM*
Cc: Branch review:
MP ready?: Task progress:

Description

The tilt of the sea surface in the rheology is better represented by the gradient of ssh rather that its approximation using geostrophic current as it is done currently in lim2. This change has already being done in LIM3.

Here we propose the modification same modification in LIM2.

This changes the ice thickness by improving the ice velocity.

The main modification is to replace the following lines:

            ! Computation of the velocity field taking into account the ice-ice interaction.                                 
            ! Terms that are independent of the ice velocity field.
            za1ct(ji,jj) = ztagnx - zcorl(ji,jj) * v_oce(ji,jj) - zgphsx
            za2ct(ji,jj) = ztagny + zcorl(ji,jj) * u_oce(ji,jj) - zgphsy

by

            ! Gradient of the sea surface height
            zgsshx =  (   (ssh_m(ji  ,jj  ) - ssh_m(ji-1,jj  ))/e1u(ji-1,jj  )   &
               &       +  (ssh_m(ji  ,jj-1) - ssh_m(ji-1,jj-1))/e1u(ji-1,jj-1)   ) * 0.5
            zgsshy =  (   (ssh_m(ji  ,jj  ) - ssh_m(ji  ,jj-1))/e2v(ji  ,jj-1)   &
               &       +  (ssh_m(ji-1,jj  ) - ssh_m(ji-1,jj-1))/e2v(ji-1,jj-1)   ) * 0.5

            ! Computation of the velocity field taking into account the ice-ice interaction.                                 
            ! Terms that are independent of the ice velocity field.
            za1ct(ji,jj) = ztagnx - zmass(ji,jj) * grav * zgsshx - zgphsx
            za2ct(ji,jj) = ztagny - zmass(ji,jj) * grav * zgsshy - zgphsy

and add in the header:

   USE sbc_oce        ! surface boundary condition: ocean variables

The resulting routine can be found in attachment. It has been tested in a 100y long run with ORCA2_LIM

Gurvan & Yevgeny

Commit History (1)

ChangesetAuthorTimeChangeLog
1711rblod2009-11-04T17:01:10+01:00

Add tilt off sea-surface in LIM2 rheology, see ticket #587

Attachments (1)

limrhg_2_v3.2.F90 (33.2 KB) - added by gm 11 years ago.
limrhg_2.F90 for the v3.2 with the use of ssh_m

Download all attachments as: .zip

Change History (13)

Changed 11 years ago by gm

limrhg_2.F90 for the v3.2 with the use of ssh_m

comment:1 Changed 11 years ago by rblod

  • Resolution set to fixed
  • Status changed from new to closed

comment:2 Changed 5 years ago by nicolasmartin

  • Keywords 2010 3: Developments Implementation Stream features new of added

comment:3 Changed 5 years ago by nicolasmartin

  • Keywords 3: removed

comment:4 Changed 5 years ago by nicolasmartin

  • Keywords developments added; Developments removed

comment:5 Changed 5 years ago by nicolasmartin

  • Keywords of removed

comment:6 Changed 5 years ago by nicolasmartin

  • Keywords nemo_v3_2* added

comment:7 Changed 5 years ago by nicolasmartin

  • Keywords nemo_v3_2_beta added; Implementation Stream developments features removed

comment:8 Changed 5 years ago by nicolasmartin

  • Keywords new removed

comment:9 Changed 5 years ago by nicolasmartin

  • Milestone 2010 Stream 3: Developments & Implementation of new features deleted

Milestone 2010 Stream 3: Developments & Implementation of new features deleted

comment:10 Changed 5 years ago by nicolasmartin

  • Keywords nemo_v3_2_beta removed

comment:11 Changed 3 years ago by nemo

  • Keywords LIM* added

comment:12 Changed 3 years ago by nemo

  • Keywords nemo_v3_2* removed
Note: See TracTickets for help on using tickets.