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 7616 for branches/2016 – NEMO

Changeset 7616 for branches/2016


Ignore:
Timestamp:
2017-01-27T18:09:47+01:00 (7 years ago)
Author:
acc
Message:

Branch 2016/dev_merge_2016. Wetting and drying updates. Includes a slightly modified MY_SRC copy of domain.F90 to enable a north-south cyclic test domain (jperio=7). Other changes include fixes in usrdef_zgr.F90 to correct metrics on the northern boundary and a new test case (#8) which uses the north-south cyclic option.

Location:
branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/MY_SRC
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/MY_SRC/usrdef_istate.F90

    r7609 r7616  
    9090            end do 
    9191            !                                     ! ==================== 
    92          CASE ( 2 )                               ! WAD 2 configuration 
     92         CASE ( 2, 8 )                            ! WAD 2 configuration 
    9393            !                                     ! ==================== 
    9494            ! 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/MY_SRC/usrdef_nam.F90

    r7547 r7616  
    9494      !                             ! Set the lateral boundary condition of the global domain 
    9595      kperio = 0                    ! WAD_TEST_CASES configuration : closed domain 
     96      IF( nn_wad_test == 8 ) kperio = 7 ! North-South cyclic test 
    9697      ! 
    9798      WRITE(ldtxt(ii),*) '   '                                                                          ;   ii = ii + 1 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/WAD_TEST_CASES/MY_SRC/usrdef_zgr.F90

    r7609 r7616  
    112112               zht(:,mj0(jpjglo):mj1(jpjglo)) = -4._wp 
    113113               !                                     ! ==================== 
    114             CASE ( 2, 3 )                            ! WAD 2 or 3  configuration 
     114            CASE ( 2, 3, 8 )                         ! WAD 2 or 3  configuration 
    115115               !                                     ! ==================== 
    116116               ! 
     
    125125               zht(mi0(1):mi1(1),:) = -4._wp 
    126126               zht(mi0(jpiglo):mi1(jpiglo),:) = -4._wp 
    127                zht(:,mj0(1):mj1(1)) = -4._wp 
    128                zht(:,mj0(jpjglo):mj1(jpjglo)) = -4._wp 
     127               IF( nn_cfg /= 8 ) THEN 
     128                  zht(:,mj0(1):mj1(1)) = -4._wp 
     129                  zht(:,mj0(jpjglo):mj1(jpjglo)) = -4._wp 
     130               ENDIF 
    129131               !                                     ! ==================== 
    130132            CASE ( 4 )                               ! WAD 4 configuration 
     
    305307                  zwet = MAX( zhu(ji,jj), rn_wdmin1 ) 
    306308                  pe3u (ji,jj,jk) = zwet * z1_jpkm1 
     309                  pe3uw(ji,jj,jk) = zwet * z1_jpkm1 
    307310                  pe3f (ji,jj,jk) = zwet * z1_jpkm1 
    308                   pe3uw(ji,jj,jk) = zwet * z1_jpkm1 
    309311                  zwet = MAX( zhv(ji,jj), rn_wdmin1 ) 
     312                  pe3v (ji,jj,jk) = zwet * z1_jpkm1 
    310313                  pe3vw(ji,jj,jk) = zwet * z1_jpkm1 
    311                   pe3v (ji,jj,jk) = zwet * z1_jpkm1 
    312314              END DO       
    313315           END DO       
    314316         END DO       
     317         CALL lbc_lnk( pdept, 'T', 1. ) 
     318         CALL lbc_lnk( pdepw, 'T', 1. ) 
     319         CALL lbc_lnk( pe3t , 'T', 1. ) 
     320         CALL lbc_lnk( pe3w , 'T', 1. ) 
    315321         CALL lbc_lnk( pe3u , 'U', 1. ) 
    316322         CALL lbc_lnk( pe3uw, 'U', 1. ) 
     323         CALL lbc_lnk( pe3f , 'F', 1. ) 
    317324         CALL lbc_lnk( pe3v , 'V', 1. ) 
    318325         CALL lbc_lnk( pe3vw, 'V', 1. ) 
    319          CALL lbc_lnk( pe3f , 'F', 1. ) 
    320326         WHERE( pe3t (:,:,:) == 0._wp )   pe3t (:,:,:) = 1._wp 
    321327         WHERE( pe3u (:,:,:) == 0._wp )   pe3u (:,:,:) = 1._wp 
Note: See TracChangeset for help on using the changeset viewer.