New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8595 – NEMO

Changeset 8595


Ignore:
Timestamp:
2017-10-05T11:49:56+02:00 (6 years ago)
Author:
flavoni
Message:

#1715 fix for zdfric if ln_mldw=true

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

    r7049 r8595  
    178178      !  Compute Ekman depth from wind stress forcing. 
    179179      ! ------------------------------------------------------- 
    180       zflageos = ( 0.5 + SIGN( 0.5, nn_eos - 1. ) ) * rau0 
     180      !SF  zflageos = ( 0.5 + SIGN( 0.5, nn_eos - 1. ) ) * rau0 
     181      !SF  DO jj = 1, jpj 
     182      !SF     DO ji = 1, jpi 
     183      !SF        zrhos          = rhop(ji,jj,1) + zflageos * ( 1. - tmask(ji,jj,1) ) 
     184      !SF        zustar         = SQRT( taum(ji,jj) / ( zrhos +  rsmall ) ) 
     185      !SF        ekm_dep(ji,jj) = rn_ekmfc * zustar / ( ABS( ff(ji,jj) ) + rsmall ) 
     186      !SF        ekm_dep(ji,jj) = MAX(ekm_dep(ji,jj),rn_mldmin) ! Minimun allowed 
     187      !SF        ekm_dep(ji,jj) = MIN(ekm_dep(ji,jj),rn_mldmax) ! Maximum allowed 
     188      !SF     END DO 
     189      !SF  END DO 
     190 
    181191      DO jj = 1, jpj 
    182192         DO ji = 1, jpi 
    183             zrhos          = rhop(ji,jj,1) + zflageos * ( 1. - tmask(ji,jj,1) ) 
    184             zustar         = SQRT( taum(ji,jj) / ( zrhos +  rsmall ) ) 
     193            zustar         = SQRT( taum(ji,jj) * r1_rau0 ) 
    185194            ekm_dep(ji,jj) = rn_ekmfc * zustar / ( ABS( ff(ji,jj) ) + rsmall ) 
    186195            ekm_dep(ji,jj) = MAX(ekm_dep(ji,jj),rn_mldmin) ! Minimun allowed 
Note: See TracChangeset for help on using the changeset viewer.