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 5845 for branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90 – NEMO

Ignore:
Timestamp:
2015-10-31T08:40:45+01:00 (8 years ago)
Author:
gm
Message:

#1613: vvl by default: suppression of domzgr_substitute.h90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r5132 r5845  
    5959#endif 
    6060 
    61 #  include "domzgr_substitute.h90" 
    6261   !!---------------------------------------------------------------------- 
    6362   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    289288                              DO ik = 1, jpkm1 
    290289                                 dta%u2d(ib) = dta%u2d(ib) & 
    291                        &                          + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik) 
     290                       &                          + e3u_n(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik) 
    292291                              END DO 
    293                               dta%u2d(ib) =  dta%u2d(ib) * hur(ii,ij) 
     292                              dta%u2d(ib) =  dta%u2d(ib) * r1_hu_n(ii,ij) 
    294293                           END DO 
    295294                           igrd = 3                      ! meridional velocity 
     
    300299                              DO ik = 1, jpkm1 
    301300                                 dta%v2d(ib) = dta%v2d(ib) & 
    302                        &                       + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik) 
     301                       &                       + e3v_n(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik) 
    303302                              END DO 
    304                               dta%v2d(ib) =  dta%v2d(ib) * hvr(ii,ij) 
     303                              dta%v2d(ib) =  dta%v2d(ib) * r1_hv_n(ii,ij) 
    305304                           END DO 
    306305                        ENDIF                     
     
    353352                        DO ik = 1, jpkm1 
    354353                           dta%u2d(ib) = dta%u2d(ib) & 
    355                  &                       + fse3u(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik) 
     354                 &                       + e3u_n(ii,ij,ik) * umask(ii,ij,ik) * dta%u3d(ib,ik) 
    356355                        END DO 
    357                         dta%u2d(ib) =  dta%u2d(ib) * hur(ii,ij) 
     356                        dta%u2d(ib) =  dta%u2d(ib) * r1_hu_n(ii,ij) 
    358357                        DO ik = 1, jpkm1 
    359358                           dta%u3d(ib,ik) = dta%u3d(ib,ik) - dta%u2d(ib) 
     
    367366                        DO ik = 1, jpkm1 
    368367                           dta%v2d(ib) = dta%v2d(ib) & 
    369                  &                       + fse3v(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik) 
     368                 &                       + e3v_n(ii,ij,ik) * vmask(ii,ij,ik) * dta%v3d(ib,ik) 
    370369                        END DO 
    371                         dta%v2d(ib) =  dta%v2d(ib) * hvr(ii,ij) 
     370                        dta%v2d(ib) =  dta%v2d(ib) * r1_hv_n(ii,ij) 
    372371                        DO ik = 1, jpkm1 
    373372                           dta%v3d(ib,ik) = dta%v3d(ib,ik) - dta%v2d(ib) 
     
    882881         ENDIF 
    883882#endif 
    884  
    885       ENDDO ! ib_bdy  
    886  
     883         ! 
     884      END DO ! ib_bdy  
     885      ! 
    887886      IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init') 
    888  
    889       END SUBROUTINE bdy_dta_init 
     887      ! 
     888   END SUBROUTINE bdy_dta_init 
    890889 
    891890#else 
Note: See TracChangeset for help on using the changeset viewer.