Changeset 7359


Ignore:
Timestamp:
2016-11-28T18:44:31+01:00 (4 years ago)
Author:
emanuelaclementi
Message:

#1805 updated nomenclature in 2016/dev_INGV_UKMO_2016

Location:
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r7353 r7359  
    412412         srcv(jpr_itx1)%clgrid = 'U'                  ! ocean stress used after its transformation 
    413413         srcv(jpr_ity1)%clgrid = 'V'                  ! i.e. it is always at U- & V-points for i- & j-comp. resp. 
     414      ENDIF 
    414415      ENDIF 
    415416 
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

    r7353 r7359  
    9696               zk = ABS(zsp0)/MAX(ABS(5.97_wp*ztransp),0.0000001_wp) 
    9797               ! Depth attenuation 
    98                zfac = EXP(-2.0_wp*zk*fsdept(ji,jj,jk))/(1.0_wp+8.0_wp*zk*fsdept(ji,jj,jk)) 
     98               zfac = EXP(-2.0_wp*zk*gdept_n(ji,jj,jk))/(1.0_wp+8.0_wp*zk*gdept_n(ji,jj,jk)) 
    9999               ! 
    100100               usd3dt(ji,jj,jk) = zfac * zusd2dt(ji,jj) * tmask(ji,jj,jk) 
     
    139139      ! 
    140140      DO jk = jpkm1, 1, -1                   ! integrate from the bottom the e3t * hor. divergence 
    141          wsd3d(:,:,jk) = wsd3d(:,:,jk+1) - fse3t_n(:,:,jk) * ze3hdiv(:,:,jk) 
     141         wsd3d(:,:,jk) = wsd3d(:,:,jk+1) - e3t_n(:,:,jk) * ze3hdiv(:,:,jk) 
    142142      END DO 
    143143#if defined key_bdy 
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r7351 r7359  
    108108      IF( ln_wave .AND. ln_sdw )  THEN 
    109109         DO jk = 1, jpkm1 
    110             zun(:,:,jk) = e2u(:,:) * fse3u(:,:,jk) *      & 
     110            zun(:,:,jk) = e2u(:,:) * e3u_n(:,:,jk) *      & 
    111111                        &  ( un(:,:,jk) + usd3d(:,:,jk) )                       ! eulerian transport + Stokes Drift 
    112             zvn(:,:,jk) = e1v(:,:) * fse3v(:,:,jk) *      & 
     112            zvn(:,:,jk) = e1v(:,:) * e3v_n(:,:,jk) *      & 
    113113                        &  ( vn(:,:,jk) + vsd3d(:,:,jk) ) 
    114114            zwn(:,:,jk) = e1e2t(:,:) *                    & 
     
    117117      ELSE 
    118118         DO jk = 1, jpkm1 
    119             zun(:,:,jk) = e2u  (:,:) * fse3u(:,:,jk) * un(:,:,jk)               ! eulerian transport only 
    120             zvn(:,:,jk) = e1v  (:,:) * fse3v(:,:,jk) * vn(:,:,jk) 
     119            zun(:,:,jk) = e2u  (:,:) * e3u_n(:,:,jk) * un(:,:,jk)               ! eulerian transport only 
     120            zvn(:,:,jk) = e1v  (:,:) * e3v_n(:,:,jk) * vn(:,:,jk) 
    121121            zwn(:,:,jk) = e1e2t(:,:)                 * wn(:,:,jk) 
    122122         END DO 
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfqiao.F90

    r7353 r7359  
    7070         DO jj = 1, jpjm1 
    7171            DO ji = 1, fs_jpim1 
    72                qbv(ji,jj,jk) = 1.0 * 0.353553 * swh(ji,jj) * tsd2d(ji,jj) *           & 
    73             &                  EXP(3.0 * wnum(ji,jj) *                                &                      
    74             &                  (-MIN( fsdepw(ji  ,jj  ,jk), fsdepw(ji+1,jj  ,jk),     & 
    75             &                         fsdepw(ji  ,jj+1,jk), fsdepw(ji+1,jj+1,jk))))   & 
     72               qbv(ji,jj,jk) = 1.0 * 0.353553 * swh(ji,jj) * tsd2d(ji,jj) *             & 
     73            &                  EXP(3.0 * wnum(ji,jj) *                                  &                      
     74            &                  (-MIN( gdepw_n(ji  ,jj  ,jk), gdepw_n(ji+1,jj  ,jk),     & 
     75            &                         gdepw_n(ji  ,jj+1,jk), gdepw_n(ji+1,jj+1,jk))))   & 
    7676            &                          * wmask(ji,jj,jk) 
    7777            END DO 
  • branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

    r7351 r7359  
    205205         DO jj = 2, jpjm1 
    206206            DO ji = fs_2, fs_jpim1 
    207                IF( fsdept(ji,jj,jk) < ekm_dep(ji,jj) ) THEN 
     207               IF( gdepw_n(ji,jj,jk) < ekm_dep(ji,jj) ) THEN 
    208208                  avmv(ji,jj,jk) = MAX( avmv(ji,jj,jk), rn_wvmix ) 
    209209                  avmu(ji,jj,jk) = MAX( avmu(ji,jj,jk), rn_wvmix ) 
Note: See TracChangeset for help on using the changeset viewer.