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 13024 for utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_recompute_scalefactors.f90 – NEMO

Ignore:
Timestamp:
2020-06-03T16:26:23+02:00 (4 years ago)
Author:
rblod
Message:

First version of new nesting tools merged with domaincfg, see ticket #2129

File:
1 edited

Legend:

Unmodified
Added
Removed
  • utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/agrif_recompute_scalefactors.f90

    r12414 r13024  
    99      ! Scale factors and depth at U-, V-, UW and VW-points 
    1010      DO jk = 1, jpk                        ! initialisation to z-scale factors 
    11 !         e3u_0 (:,:,jk) = e3t_1d(jk) 
    12 !         e3v_0 (:,:,jk) = e3t_1d(jk) 
     11         e3u_0 (:,:,jk) = e3t_1d(jk) 
     12         e3v_0 (:,:,jk) = e3t_1d(jk) 
    1313         e3uw_0(:,:,jk) = e3w_1d(jk) 
    1414         e3vw_0(:,:,jk) = e3w_1d(jk) 
     
    1818         DO jj = 1, jpjm1 
    1919            DO ji = 1, jpim1   ! vector opt. 
    20 !               e3u_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji+1,jj,jk) ) 
    21 !               e3v_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji,jj+1,jk) ) 
     20               e3u_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji+1,jj,jk) ) 
     21               e3v_0 (ji,jj,jk) = MIN( e3t_0(ji,jj,jk), e3t_0(ji,jj+1,jk) ) 
    2222               e3uw_0(ji,jj,jk) = MIN( e3w_0(ji,jj,jk), e3w_0(ji+1,jj,jk) ) 
    2323               e3vw_0(ji,jj,jk) = MIN( e3w_0(ji,jj,jk), e3w_0(ji,jj+1,jk) ) 
     
    2727      IF ( ln_isfcav ) THEN 
    2828      ! (ISF) define e3uw (adapted for 2 cells in the water column) 
     29      print *,'NOT READY SINCE:' 
     30      print *,'MBATHY HAS NOT BEEN CORRECTED / UPDATED' 
     31      print *,'EVEN NOT COMPUTED IN THE CASE ln_read_cfg = .TRUE.' 
     32      STOP 
    2933         DO jj = 2, jpjm1  
    3034            DO ji = 2, jpim1   ! vector opt.  
Note: See TracChangeset for help on using the changeset viewer.