Changeset 917


Ignore:
Timestamp:
2008-05-07T16:50:27+02:00 (13 years ago)
Author:
rblod
Message:

Correct a bug in revision 913 (ice/ocean stress computation)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limsbc.F90

    r913 r917  
    8282      REAL(wp) ::   zfrldu, zfrldv   ! lead fraction at U- & V-points 
    8383      REAL(wp) ::   zu_io , zv_io    ! 2 components of the ice-ocean velocity 
     84      REAL(wp) ::   ztglx , ztgly    ! temporary scalar 
    8485       
    8586#if defined key_coupled     
     
    279280               zmod  = SQRT( zu_io * zu_io + zv_io * zv_io )  
    280281               ! quadratic drag formulation 
    281                ztio_u(ji,jj) = rhoco * zmod * ( cangvg * zu_io - zsang * zv_io )  
    282                ztio_v(ji,jj) = rhoco * zmod * ( cangvg * zv_io + zsang * zu_io )  
     282               ztglx = rhoco * zmod * ( cangvg * zu_io - zsang * zv_io ) 
     283               ztgly = rhoco * zmod * ( cangvg * zv_io + zsang * zu_io ) 
     284               ! IMPORTANT 
     285               ! these lines are bound to prevent numerical oscillations 
     286               ! in the ice-ocean stress 
     287               ! They are physically ill-based. There is a cleaner solution 
     288               ! to try (remember discussion in Paris Gurvan) 
     289               ztio_u(ji,jj) = ztglx * exp( - zmod / 0.5 )  
     290               ztio_v(ji,jj) = ztgly * exp( - zmod / 0.5 )  
    283291               ! 
    284292            END DO 
Note: See TracChangeset for help on using the changeset viewer.