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 6020 for branches/UKMO/icebergs_restart_single_file/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 – NEMO

Ignore:
Timestamp:
2015-12-08T12:39:53+01:00 (8 years ago)
Author:
timgraham
Message:

Merged with head of trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/icebergs_restart_single_file/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r6019 r6020  
    219219         &  ppsur == pp_to_be_computed           ) THEN 
    220220         ! 
     221#if defined key_agrif 
     222         za1  = (  ppdzmin - pphmax / FLOAT(jpkdta-1)  )                                                   & 
     223            & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpkdta-1) * (  LOG( COSH( (jpkdta - ppkth) / ppacr) )& 
     224            &                                                      - LOG( COSH( ( 1  - ppkth) / ppacr) )  )  ) 
     225#else 
    221226         za1  = (  ppdzmin - pphmax / FLOAT(jpkm1)  )                                                      & 
    222227            & / ( TANH((1-ppkth)/ppacr) - ppacr/FLOAT(jpk-1) * (  LOG( COSH( (jpk - ppkth) / ppacr) )      & 
    223228            &                                                   - LOG( COSH( ( 1  - ppkth) / ppacr) )  )  ) 
     229#endif 
    224230         za0  = ppdzmin - za1 *              TANH( (1-ppkth) / ppacr ) 
    225231         zsur =   - za0 - za1 * ppacr * LOG( COSH( (1-ppkth) / ppacr )  ) 
     
    236242              WRITE(numout,*) '            Uniform grid with ',jpk-1,' layers' 
    237243              WRITE(numout,*) '            Total depth    :', zhmax 
     244#if defined key_agrif 
     245              WRITE(numout,*) '            Layer thickness:', zhmax/(jpkdta-1) 
     246#else 
    238247              WRITE(numout,*) '            Layer thickness:', zhmax/(jpk-1) 
     248#endif 
    239249         ELSE 
    240250            IF( ppa1 == 0._wp .AND. ppa0 == 0._wp .AND. ppsur == 0._wp ) THEN 
     
    260270      ! Reference z-coordinate (depth - scale factor at T- and W-points) 
    261271      ! ====================== 
    262       IF( ppkth == 0._wp ) THEN            !  uniform vertical grid        
     272      IF( ppkth == 0._wp ) THEN            !  uniform vertical grid  
     273#if defined key_agrif 
     274         za1 = zhmax / FLOAT(jpkdta-1)  
     275#else 
    263276         za1 = zhmax / FLOAT(jpk-1)  
     277#endif 
    264278         DO jk = 1, jpk 
    265279            zw = FLOAT( jk ) 
     
    487501            IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    488502               !                                             ! ===================== 
    489                IF( nn_cla == 0 ) THEN 
    490                   ii0 = 140   ;   ii1 = 140                  ! Gibraltar Strait open  
    491                   ij0 = 102   ;   ij1 = 102                  ! (Thomson, Ocean Modelling, 1995) 
    492                   DO ji = mi0(ii0), mi1(ii1) 
    493                      DO jj = mj0(ij0), mj1(ij1) 
    494                         mbathy(ji,jj) = 15 
    495                      END DO 
     503               ! 
     504               ii0 = 140   ;   ii1 = 140                  ! Gibraltar Strait open  
     505               ij0 = 102   ;   ij1 = 102                  ! (Thomson, Ocean Modelling, 1995) 
     506               DO ji = mi0(ii0), mi1(ii1) 
     507                  DO jj = mj0(ij0), mj1(ij1) 
     508                     mbathy(ji,jj) = 15 
    496509                  END DO 
    497                   IF(lwp) WRITE(numout,*) 
    498                   IF(lwp) WRITE(numout,*) '      orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 
    499                   ! 
    500                   ii0 = 160   ;   ii1 = 160                  ! Bab el mandeb Strait open 
    501                   ij0 = 88    ;   ij1 = 88                   ! (Thomson, Ocean Modelling, 1995) 
    502                   DO ji = mi0(ii0), mi1(ii1) 
    503                      DO jj = mj0(ij0), mj1(ij1) 
    504                         mbathy(ji,jj) = 12 
    505                      END DO 
     510               END DO 
     511               IF(lwp) WRITE(numout,*) 
     512               IF(lwp) WRITE(numout,*) '      orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 
     513               ! 
     514               ii0 = 160   ;   ii1 = 160                  ! Bab el mandeb Strait open 
     515               ij0 = 88    ;   ij1 = 88                   ! (Thomson, Ocean Modelling, 1995) 
     516               DO ji = mi0(ii0), mi1(ii1) 
     517                  DO jj = mj0(ij0), mj1(ij1) 
     518                     mbathy(ji,jj) = 12 
    506519                  END DO 
    507                   IF(lwp) WRITE(numout,*) 
    508                   IF(lwp) WRITE(numout,*) '      orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    509                ENDIF 
     520               END DO 
     521               IF(lwp) WRITE(numout,*) 
     522               IF(lwp) WRITE(numout,*) '      orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    510523               ! 
    511524            ENDIF 
     
    531544            !        
    532545            IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    533                ! 
    534               IF( nn_cla == 0 ) THEN 
    535                  ii0 = 140   ;   ii1 = 140                   ! Gibraltar Strait open  
    536                  ij0 = 102   ;   ij1 = 102                   ! (Thomson, Ocean Modelling, 1995) 
    537                  DO ji = mi0(ii0), mi1(ii1) 
    538                     DO jj = mj0(ij0), mj1(ij1) 
    539                        bathy(ji,jj) = 284._wp 
    540                     END DO 
     546            ! 
     547              ii0 = 140   ;   ii1 = 140                   ! Gibraltar Strait open  
     548              ij0 = 102   ;   ij1 = 102                   ! (Thomson, Ocean Modelling, 1995) 
     549              DO ji = mi0(ii0), mi1(ii1) 
     550                 DO jj = mj0(ij0), mj1(ij1) 
     551                    bathy(ji,jj) = 284._wp 
    541552                 END DO 
    542                  IF(lwp) WRITE(numout,*)      
    543                  IF(lwp) WRITE(numout,*) '      orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 
    544                  ! 
    545                  ii0 = 160   ;   ii1 = 160                   ! Bab el mandeb Strait open 
    546                  ij0 = 88    ;   ij1 = 88                    ! (Thomson, Ocean Modelling, 1995) 
    547                  DO ji = mi0(ii0), mi1(ii1) 
    548                     DO jj = mj0(ij0), mj1(ij1) 
    549                        bathy(ji,jj) = 137._wp 
    550                     END DO 
     553               END DO 
     554              IF(lwp) WRITE(numout,*)      
     555              IF(lwp) WRITE(numout,*) '      orca_r2: Gibraltar strait open at i=',ii0,' j=',ij0 
     556              ! 
     557              ii0 = 160   ;   ii1 = 160                   ! Bab el mandeb Strait open 
     558              ij0 = 88    ;   ij1 = 88                    ! (Thomson, Ocean Modelling, 1995) 
     559               DO ji = mi0(ii0), mi1(ii1) 
     560                 DO jj = mj0(ij0), mj1(ij1) 
     561                    bathy(ji,jj) = 137._wp 
    551562                 END DO 
    552                  IF(lwp) WRITE(numout,*) 
    553                  IF(lwp) WRITE(numout,*) '             orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    554               ENDIF 
     563              END DO 
     564              IF(lwp) WRITE(numout,*) 
     565               IF(lwp) WRITE(numout,*) '             orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    555566              ! 
    556567           ENDIF 
Note: See TracChangeset for help on using the changeset viewer.