Ignore:
Timestamp:
2018-02-23T16:52:00+01:00 (2 years ago)
Author:
jchanut
Message:

Tilde coordinate update: first draft

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2018/dev_r8864_nemo_v3_6_ZTILDE/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r7522 r9353  
    7070      !!              - 1D configuration, move Coriolis, u and v at T-point 
    7171      !!---------------------------------------------------------------------- 
    72       INTEGER ::   jk          ! dummy loop argument 
     72      INTEGER ::   jj, jk      ! dummy loop argument 
    7373      INTEGER ::   iconf = 0   ! local integers 
    7474      !!---------------------------------------------------------------------- 
     
    9292      hu_0(:,:) = 0.0_wp                       ! Reference ocean depth at U-points 
    9393      hv_0(:,:) = 0.0_wp                       ! Reference ocean depth at V-points 
     94      hf_0(:,:) = 0.0_wp                       ! Reference ocean depth at F-points 
    9495      DO jk = 1, jpk 
    9596         ht_0(:,:) = ht_0(:,:) + e3t_0(:,:,jk) * tmask(:,:,jk) 
     
    9798         hv_0(:,:) = hv_0(:,:) + e3v_0(:,:,jk) * vmask(:,:,jk) 
    9899      END DO 
     100      !    
     101      DO jj=1, jpjm1 
     102         DO jk = 1, jpk 
     103            hf_0(:,jj) = hf_0(:,jj) + e3f_0(:,jj,jk) * umask(:,jj,jk)*umask(:,jj+1,jk) 
     104         END DO 
     105      END DO 
     106      CALL lbc_lnk( hf_0, 'F', 1._wp ) 
    99107      ! 
    100108      IF( lk_c1d )           CALL cor_c1d      ! 1D configuration: Coriolis set at T-point 
Note: See TracChangeset for help on using the changeset viewer.