Opened 10 years ago

Closed 10 years ago

Last modified 3 years ago

#757 closed Bug (fixed)

LIM-2 - bulk CLIO bug in surface ice stress (sbcblk_clio.F90)

Reported by: gm Owned by: nemo
Priority: low Milestone:
Component: OCE Version: release-3.2
Severity: Keywords: CLIO LIM*
Cc:

Description

In clio bulk formulae, the air-ice stress is deduced from the air-sea stress by applying a ratio cai / cao,
BUT only in case of 'C'-grid ice (LIM-3) not for LIM-2 ('B'-grid case).

change :

      SELECT CASE( cd_grid )
      CASE( 'C' )                          ! C-grid ice dynamics
         ! Change from wind speed to wind stress over OCEAN (cao is used)
         zcoef = cai / cao 
!CDIR COLLAPSE
         DO jj = 1 , jpj
            DO ji = 1, jpi
               p_taui(ji,jj) = zcoef * utau(ji,jj)
               p_tauj(ji,jj) = zcoef * vtau(ji,jj)
            END DO
         END DO
      CASE( 'B' )                          ! B-grid ice dynamics:  I-point (i.e. F-point lower-left corner)
!CDIR COLLAPSE
         DO jj = 2, jpj         ! stress from ocean U- and V-points to ice U,V point
            DO ji = 2, jpi   ! I-grid : no vector opt.
               p_taui(ji,jj) = 0.5 * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) )
               p_taui(ji,jj) = 0.5 * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) )
            END DO
         END DO
         CALL lbc_lnk( p_taui(:,:), 'I', -1. )   ;   CALL lbc_lnk( p_tauj(:,:), 'I', -1. )   ! I-point
      END SELECT

into

      SELECT CASE( cd_grid )
      CASE( 'C' )                          ! C-grid ice dynamics
      zcoef  = cai / cao                         ! Change from air-sea stress to air-ice stress
!CDIR COLLAPSE
         DO jj = 1 , jpj
            DO ji = 1, jpi
               p_taui(ji,jj) = zcoef * utau(ji,jj)
               p_tauj(ji,jj) = zcoef * vtau(ji,jj)
            END DO
         END DO
      CASE( 'B' )                          ! B-grid ice dynamics:  I-point (i.e. F-point lower-left corner)
      zcoef  = 0.5_wp * cai / cao                ! Change from air-sea stress to air-ice stress!CDIR COLLAPSE
!CDIR COLLAPSE
         DO jj = 2, jpj         ! stress from ocean U- and V-points to ice U,V point
            DO ji = 2, jpi   ! I-grid : no vector opt.
               p_taui(ji,jj) = zcoef2 * ( utau(ji-1,jj  ) + utau(ji-1,jj-1) )
               p_tauj(ji,jj) = zcoef2 * ( vtau(ji  ,jj-1) + vtau(ji-1,jj-1) )
            END DO
         END DO
         CALL lbc_lnk( p_taui(:,:), 'I', -1. )   ;   CALL lbc_lnk( p_tauj(:,:), 'I', -1. )   ! I-point
      END SELECT

Commit History (1)

ChangesetAuthorTimeChangeLog
2388smasson2010-11-15T16:24:52+01:00

ticket #757, LIM-2 - bulk CLIO bug in surface ice stress

Change History (8)

comment:1 Changed 10 years ago by gm

hups, Obviously in the modified code, zcoef2 stands for zcoef

Gurvan

comment:2 Changed 10 years ago by gm

done in v3.3 beta see Changeset 2370

=⇒>> update v3.2 remains to be done

comment:3 Changed 10 years ago by smasson

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

done for NEMO v3.2, see changeset:2388

comment:4 Changed 4 years ago by nicolasmartin

  • Keywords LIM2 added; LIM-2 removed

comment:5 Changed 4 years ago by nicolasmartin

  • Keywords nemo_v3_2* added

comment:6 Changed 4 years ago by nicolasmartin

  • Keywords and removed

comment:7 Changed 3 years ago by nemo

  • Keywords LIM* added; LIM2 removed

comment:8 Changed 3 years ago by nemo

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