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 6934 – NEMO

Changeset 6934


Ignore:
Timestamp:
2016-09-16T18:09:35+02:00 (8 years ago)
Author:
csanchez
Message:

Added Enda O'Dea changes for SINGV

Location:
branches/UKMO/SINGV_NEMO_5518/NEMOGCM/NEMO/OPA_SRC/DOM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/SINGV_NEMO_5518/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90

    r6933 r6934  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    34    !! $Id$ 
     34   !! $Id$  
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
     
    7979      !                                   !  workspaces 
    8080      REAL(wp), POINTER, DIMENSION(:,:  ) :: zprt, zprw  
    81       REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepu, zdepv 
     81      REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepu, zdepv, zdepf 
    8282      !!---------------------------------------------------------------------- 
    8383      ! 
     
    8585      ! 
    8686      CALL wrk_alloc( jpi, jpj, zprt, zprw ) 
    87       CALL wrk_alloc( jpi, jpj, jpk, zdepu, zdepv ) 
     87      CALL wrk_alloc( jpi, jpj, jpk, zdepu, zdepv, zdepf ) 
    8888      ! 
    8989      IF(lwp) WRITE(numout,*) 
     
    216216         CALL iom_rstput( 0, 0, inum4, 'gdepw_1d' , gdepw_1d ) 
    217217      ENDIF 
     218     
     219      !CEOD Added for gdep uvf 
     220      IF( ln_sco ) THEN                                         ! s-coordinate 
     221          
     222         DO jk = 1,jpk 
     223            DO jj = 1, jpjm1 
     224               DO ji = 1, fs_jpim1   ! vector opt. 
     225                  zdepu(ji,jj,jk) = 0.5*(gdept_0(ji,jj,jk) + gdept_0(ji+1,jj ,jk) ) 
     226                  zdepv(ji,jj,jk) = 0.5*(gdept_0(ji,jj,jk) + gdept_0(ji ,jj+1,jk) ) 
     227 
     228                  zdepf(ji,jj,jk) = 0.25*(gdept_0(ji,jj,jk) + gdept_0(ji+1 ,jj+1,jk) +gdept_0(ji+1,jj  ,jk)  + gdept_0(ji  ,jj+1,jk) )  
     229               END DO 
     230            END DO 
     231         END DO 
     232          
     233         CALL lbc_lnk( zdepu, 'U', 1. )   ;   CALL lbc_lnk( zdepv, 'V', 1. ) 
     234         CALL iom_rstput( 0, 0, inum4, 'gdepu', zdepu, ktype = jp_r4 ) 
     235         CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r4 ) 
     236         CALL iom_rstput( 0, 0, inum4, 'gdepf', zdepf, ktype = jp_r4 ) 
     237         CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r4 ) 
     238       
     239      END IF  
    218240       
    219241      IF( ln_zps ) THEN                                         ! z-coordinate - partial steps 
  • branches/UKMO/SINGV_NEMO_5518/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r6933 r6934  
    19581958         DO jj = 1, jpj 
    19591959            DO ji = 1, jpi 
    1960                ztaper = EXP( -(gphit(ji,jj)/8._wp)**2._wp ) 
     1960               ztaper = 0.0 !CEOD WARN EXP( -(gphit(ji,jj)/8._wp)**2._wp ) 
    19611961               hbatt(ji,jj) = rn_sbot_max * ztaper + hbatt(ji,jj) * ( 1._wp - ztaper ) 
    19621962            END DO 
  • branches/UKMO/SINGV_NEMO_5518/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r6933 r6934  
    211211      ! 
    212212      DO jk = 1, jpk 
    213          tsn(:,:,jk,jp_tem) = (  ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH((fsdept(:,:,jk)-80.)/30.) )   & 
    214             &                + 10. * ( 5000. - fsdept(:,:,jk) ) /5000.)  ) * tmask(:,:,jk) 
     213         tsn(:,:,jk,jp_tem) = 10.0 !CEOD (  ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH((fsdept(:,:,jk)-80.)/30.) )   & 
     214!             &                + 10. * ( 5000. - fsdept(:,:,jk) ) /5000.)  ) * tmask(:,:,jk) 
    215215         tsb(:,:,jk,jp_tem) = tsn(:,:,jk,jp_tem) 
    216216      END DO 
Note: See TracChangeset for help on using the changeset viewer.