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 10787 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DOM/domvvl.F90 – NEMO

Ignore:
Timestamp:
2019-03-21T09:14:33+01:00 (5 years ago)
Author:
davestorkey
Message:

dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps: rename a few more grid variables.

File:
1 edited

Legend:

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

    r10767 r10787  
    616616         tilde_e3t_n(:,:,:) = tilde_e3t_a(:,:,:) 
    617617      ENDIF 
    618       gdept_b(:,:,:) = gdept_n(:,:,:) 
    619       gdepw_b(:,:,:) = gdepw_n(:,:,:) 
    620618 
    621619!! TO BE DELETED 
     620!!$      gdept_b(:,:,:) = gdept_n(:,:,:) 
     621!!$      gdepw_b(:,:,:) = gdepw_n(:,:,:) 
     622 
    622623!!$      e3t_n(:,:,:) = e3t_a(:,:,:) 
    623624!!$      e3u_n(:,:,:) = e3u_a(:,:,:) 
     
    643644 
    644645      ! t- and w- points depth (set the isf depth as it is in the initial step) 
    645       gdept_n(:,:,1) = 0.5_wp * e3w(:,:,1,kNnn_2lev) 
    646       gdepw_n(:,:,1) = 0.0_wp 
    647       gde3w_n(:,:,1) = gdept_n(:,:,1) - sshn(:,:) 
     646      gdept(:,:,1,kNnn_2lev) = 0.5_wp * e3w(:,:,1,kNnn_2lev) 
     647      gdepw(:,:,1,kNnn_2lev) = 0.0_wp 
     648      gde3w(:,:,1)           = gdept(:,:,1,kNnn_2lev) - sshn(:,:) 
    648649      DO jk = 2, jpk 
    649650         DO jj = 1,jpj 
     
    652653                                                                 ! 1 for jk = mikt 
    653654               zcoef = (tmask(ji,jj,jk) - wmask(ji,jj,jk)) 
    654                gdepw_n(ji,jj,jk) = gdepw_n(ji,jj,jk-1) + e3t(ji,jj,jk-1,kNnn) 
    655                gdept_n(ji,jj,jk) =    zcoef  * ( gdepw_n(ji,jj,jk  ) + 0.5 * e3w(ji,jj,jk,kNnn_2lev) )  & 
    656                    &             + (1-zcoef) * ( gdept_n(ji,jj,jk-1) +       e3w(ji,jj,jk,kNnn_2lev) )  
    657                gde3w_n(ji,jj,jk) = gdept_n(ji,jj,jk) - sshn(ji,jj) 
     655               gdepw(ji,jj,jk,kNnn_2lev) = gdepw(ji,jj,jk-1,kNnn_2lev) + e3t(ji,jj,jk-1,kNnn) 
     656               gdept(ji,jj,jk,kNnn_2lev) =    zcoef  * ( gdepw(ji,jj,jk,kNnn_2lev  ) + 0.5 * e3w(ji,jj,jk,kNnn_2lev) )  & 
     657                   &             + (1-zcoef) * ( gdept(ji,jj,jk-1,kNnn_2lev) +       e3w(ji,jj,jk,kNnn_2lev) )  
     658               gde3w(ji,jj,jk) = gdept(ji,jj,jk,kNnn_2lev) - sshn(ji,jj) 
    658659            END DO 
    659660         END DO 
     
    662663      ! Local depth and Inverse of the local depth of the water 
    663664      ! ------------------------------------------------------- 
    664       hu_n(:,:) = hu_a(:,:)   ;   r1_hu_n(:,:) = r1_hu_a(:,:) 
    665       hv_n(:,:) = hv_a(:,:)   ;   r1_hv_n(:,:) = r1_hv_a(:,:) 
    666       ! 
    667       ht_n(:,:) = e3t(:,:,1,kNnn) * tmask(:,:,1) 
     665!! TO BE DELETED 
     666!!$      hu_n(:,:) = hu_a(:,:)   ;   r1_hu_n(:,:) = r1_hu_a(:,:) 
     667!!$      hv_n(:,:) = hv_a(:,:)   ;   r1_hv_n(:,:) = r1_hv_a(:,:) 
     668!! TO BE DELETED 
     669      ! 
     670      ht(:,:) = e3t(:,:,1,kNnn) * tmask(:,:,1) 
    668671      DO jk = 2, jpkm1 
    669          ht_n(:,:) = ht_n(:,:) + e3t(:,:,jk,kNnn) * tmask(:,:,jk) 
     672         ht(:,:) = ht(:,:) + e3t(:,:,jk,kNnn) * tmask(:,:,jk) 
    670673      END DO 
    671674 
Note: See TracChangeset for help on using the changeset viewer.