Ignore:
Timestamp:
2017-06-01T13:24:36+02:00 (3 years ago)
Author:
davestorkey
Message:

UKMO/branch/dev_r5518_GO6_package_FVPS: FVPS implementation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FVPS/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r6487 r8112  
    11331133         e3vw_0(:,:,jk) = e3w_1d(jk) 
    11341134      END DO 
    1135       DO jk = 1,jpk                         ! Computed as the minimum of neighbooring scale factors 
     1135      DO jk = 1,jpk                         ! Computed as the average of neighbouring scale factors 
    11361136         DO jj = 1, jpjm1 
    11371137            DO ji = 1, fs_jpim1   ! vector opt. 
    1138                e3u_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji+1,jj,jk) ) 
    1139                e3v_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji,jj+1,jk) ) 
    1140                e3uw_0(ji,jj,jk) = MIN( e3w_0(ji,jj,jk), e3w_0(ji+1,jj,jk) ) 
    1141                e3vw_0(ji,jj,jk) = MIN( e3w_0(ji,jj,jk), e3w_0(ji,jj+1,jk) ) 
     1138               e3u_0 (ji,jj,jk) = 0.5 * ( e3t_0(ji,jj,jk) + e3t_0(ji+1,jj,jk) ) 
     1139               e3v_0 (ji,jj,jk) = 0.5 * ( e3t_0(ji,jj,jk) + e3t_0(ji,jj+1,jk) ) 
     1140               e3uw_0(ji,jj,jk) = 0.5 * ( e3w_0(ji,jj,jk) + e3w_0(ji+1,jj,jk) ) 
     1141               e3vw_0(ji,jj,jk) = 0.5 * ( e3w_0(ji,jj,jk) + e3w_0(ji,jj+1,jk) ) 
    11421142            END DO 
    11431143         END DO 
     
    11731173         e3f_0(:,:,jk) = e3t_1d(jk) 
    11741174      END DO 
    1175       DO jk = 1, jpk                        ! Computed as the minimum of neighbooring V-scale factors 
     1175      DO jk = 1, jpk                        ! Computed as the average of neighbouring V-scale factors 
    11761176         DO jj = 1, jpjm1 
    11771177            DO ji = 1, fs_jpim1   ! vector opt. 
    1178                e3f_0(ji,jj,jk) = MIN( e3v_0(ji,jj,jk), e3v_0(ji+1,jj,jk) ) 
     1178               e3f_0(ji,jj,jk) = 0.5 * ( e3v_0(ji,jj,jk) + e3v_0(ji+1,jj,jk) ) 
    11791179            END DO 
    11801180         END DO 
Note: See TracChangeset for help on using the changeset viewer.