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 10955 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90 – NEMO

Ignore:
Timestamp:
2019-05-10T11:15:50+02:00 (5 years ago)
Author:
acc
Message:

2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps : Convert ZDF modules and all knock on effects of these conversions. SETTE tested

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90

    r10954 r10955  
    941941 ! Stable conditions 
    942942               DO jk = 2, ibld(ji,jj) 
    943                   zznd_d=gdepw_n(ji,jj,jk) / dstokes(ji,jj) 
     943                  zznd_d=gdepw(ji,jj,jk,Kmm) / dstokes(ji,jj) 
    944944                  ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + 1.5 * EXP ( -0.9 * zznd_d ) & 
    945945                       &          *                 ( 1.0 - EXP ( -4.0 * zznd_d ) ) * zsc_wth_1(ji,jj) 
     
    10661066            IF ( lconv(ji,jj) ) THEN 
    10671067               DO jk = 2, imld(ji,jj) 
    1068                   zznd_ml=gdepw_n(ji,jj,jk) / zhml(ji,jj) 
     1068                  zznd_ml=gdepw(ji,jj,jk,Kmm) / zhml(ji,jj) 
    10691069                  ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + 0.3 * zsc_wth_1(ji,jj)                & 
    10701070                       &          * ( -2.0 + 2.75 * (       ( 1.0 + 0.6 * zznd_ml**4 )      & 
     
    11881188             IF ( lconv(ji,jj) ) THEN 
    11891189               DO jk = 1, imld(ji,jj) - 1 
    1190                   znd=gdepw_n(ji,jj,jk) / zhml(ji,jj) 
     1190                  znd=gdepw(ji,jj,jk,Kmm) / zhml(ji,jj) 
    11911191                  ghamt(ji,jj,jk) = ghamt(ji,jj,jk) + zwth_ent(ji,jj) * znd 
    11921192                  ghams(ji,jj,jk) = ghams(ji,jj,jk) + zws_ent(ji,jj) * znd 
Note: See TracChangeset for help on using the changeset viewer.