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 12622 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdfosm.F90 – NEMO

Ignore:
Timestamp:
2020-03-27T20:55:44+01:00 (4 years ago)
Author:
techene
Message:

all: add e3 substitute (sometimes it requires to add ze3t/u/v/w) and limit precompiled files lines to about 130 character

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdfosm.F90

    r12616 r12622  
    596596                     zwb_ent(ji,jj) = 0._wp 
    597597                  ENDIF 
    598                   inhml =   & 
    599                      & MAX( INT( zari * zhbl(ji,jj) / e3t(ji,jj,ibld(ji,jj),Kmm) ),  & 
    600                      &      1 ) 
     598                  inhml = MAX( INT( zari * zhbl(ji,jj)   & 
     599                     &              / e3t(ji,jj,ibld(ji,jj),Kmm) ), 1 ) 
    601600                  imld(ji,jj) = MAX( ibld(ji,jj) - inhml, 1) 
    602601                  zhml(ji,jj) = gdepw(ji,jj,imld(ji,jj),Kmm) 
     
    612611                  IF ( zdb_bl(ji,jj) > 0._wp ) THEN 
    613612               ! pycnocline thickness set by stratification - use same relationship as for neutral conditions. 
    614                      zari = MIN( 4.5 * ( zvstr(ji,jj)**2 ) & 
    615                        & / ( zdb_bl(ji,jj) * zhbl(ji,jj) ) + 0.01  , 0.2 ) 
    616                      inhml =   & 
    617                         & MAX( INT( zari * zhbl(ji,jj) / e3t(ji,jj,ibld(ji,jj),Kmm) ),   & 
    618                         &      1 ) 
     613                     zari = MIN( 4.5 * ( zvstr(ji,jj)**2 )   & 
     614                        & / ( zdb_bl(ji,jj) * zhbl(ji,jj) ) + 0.01  , 0.2 ) 
     615                     inhml = MAX( INT( zari * zhbl(ji,jj)   & 
     616                        &             / e3t(ji,jj,ibld(ji,jj),Kmm) ), 1 ) 
    619617                     imld(ji,jj) = MAX( ibld(ji,jj) - inhml, 1) 
    620618                     zhml(ji,jj) = gdepw(ji,jj,imld(ji,jj),Kmm) 
Note: See TracChangeset for help on using the changeset viewer.