Changeset 5021


Ignore:
Timestamp:
2015-01-09T12:44:02+01:00 (6 years ago)
Author:
timgraham
Message:

2nd fix for ticket #1423

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftmx.F90

    r4990 r5021  
    427427      en_tmx(:,:) = - rn_tfe * rn_me * ( zem2(:,:) * 1.25 + zek1(:,:) ) * ssmask(:,:) 
    428428 
     429!============ 
     430!TG: Bug for VVL? Should this section be moved out of _init and be updated at every timestep? 
    429431      ! Vertical structure (az_tmx) 
    430432      DO jj = 1, jpj                ! part independent of the level 
    431433         DO ji = 1, jpi 
    432             zhdep(ji,jj) = fsdepw(ji,jj,mbkt(ji,jj)+1)       ! depth of the ocean 
     434            zhdep(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1)       ! depth of the ocean 
    433435            zfact(ji,jj) = rau0 * rn_htmx * ( 1. - EXP( -zhdep(ji,jj) / rn_htmx ) ) 
    434436            IF( zfact(ji,jj) /= 0 )   zfact(ji,jj) = en_tmx(ji,jj) / zfact(ji,jj) 
     
    438440         DO jj = 1, jpj 
    439441            DO ji = 1, jpi 
    440                az_tmx(ji,jj,jk) = zfact(ji,jj) * EXP( -( zhdep(ji,jj)-fsdepw(ji,jj,jk) ) / rn_htmx ) * tmask(ji,jj,jk) 
    441             END DO 
    442          END DO 
    443       END DO 
     442               az_tmx(ji,jj,jk) = zfact(ji,jj) * EXP( -( zhdep(ji,jj)-gdepw_0(ji,jj,jk) ) / rn_htmx ) * tmask(ji,jj,jk) 
     443            END DO 
     444         END DO 
     445      END DO 
     446!=========== 
    444447 
    445448      IF( nprint == 1 .AND. lwp ) THEN 
Note: See TracChangeset for help on using the changeset viewer.