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 5443 for branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90 – NEMO

Ignore:
Timestamp:
2015-06-19T17:18:00+02:00 (9 years ago)
Author:
davestorkey
Message:

Update 2015/dev_r5021_UKMO1_CICE_coupling branch to revision 5442 of the trunk.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    • Property svn:keywords set to Id
    r4996 r5443  
    8585   !!---------------------------------------------------------------------- 
    8686   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    87    !! $Header:$  
     87   !! $Id$ 
    8888   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    8989   !!---------------------------------------------------------------------- 
     
    136136            DO jj = 1, jpj 
    137137               DO ji = 1, jpi 
    138                   zcoef = ryyss * cvol(ji,jj,1)  
     138                  zcoef = ryyss * e1e2t(ji,jj) * h_rnf(ji,jj)  
    139139                  rivalk(ji,jj) =   sf_river(jr_dic)%fnow(ji,jj,1)                                    & 
    140140                     &              * 1.E3        / ( 12. * zcoef + rtrn ) 
     
    187187      INTEGER  :: ierr, ierr1, ierr2, ierr3 
    188188      INTEGER  :: ios                 ! Local integer output status for namelist read 
     189      INTEGER  :: ik50                !  last level where depth less than 50 m 
     190      INTEGER  :: isrow             ! index for ORCA1 starting row 
    189191      REAL(wp) :: zexpide, zdenitide, zmaskt 
    190192      REAL(wp) :: ztimes_dust, ztimes_riv, ztimes_ndep  
     
    216218902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampissbc in configuration namelist', lwp ) 
    217219      IF(lwm) WRITE ( numonp, nampissbc ) 
     220 
     221      IF ( ( nn_ice_tr >= 0 ) .AND. ln_ironice ) THEN 
     222         IF(lwp) THEN 
     223            WRITE(numout,*) ' ln_ironice incompatible with nn_ice_tr = ', nn_ice_tr 
     224            WRITE(numout,*) ' Specify your sea ice iron concentration in nampisice instead ' 
     225            WRITE(numout,*) ' ln_ironice is forced to .FALSE. ' 
     226            ln_ironice = .FALSE. 
     227         ENDIF 
     228      ENDIF 
    218229 
    219230      IF(lwp) THEN 
     
    247258      ENDIF 
    248259 
     260      ! set the number of level over which river runoffs are applied  
     261      ! online configuration : computed in sbcrnf 
     262      IF( lk_offline ) THEN 
     263        nk_rnf(:,:) = 1 
     264        h_rnf (:,:) = fsdept(:,:,1) 
     265      ENDIF 
     266 
    249267      ! dust input from the atmosphere 
    250268      ! ------------------------------ 
     
    358376         rivalkinput = 0._wp 
    359377      END IF  
    360  
    361378      ! nutrient input from dust 
    362379      ! ------------------------ 
     
    410427         CALL iom_close( numiron ) 
    411428         ! 
    412          DO jk = 1, 5 
     429         ik50 = 5        !  last level where depth less than 50 m 
     430         DO jk = jpkm1, 1, -1 
     431            IF( gdept_1d(jk) > 50. )  ik50 = jk - 1 
     432         END DO 
     433         IF (lwp) WRITE(numout,*) 
     434         IF (lwp) WRITE(numout,*) ' Level corresponding to 50m depth ',  ik50,' ', gdept_1d(ik50+1) 
     435         IF (lwp) WRITE(numout,*) 
     436         DO jk = 1, ik50 
    413437            DO jj = 2, jpjm1 
    414438               DO ji = fs_2, fs_jpim1 
     
    421445            END DO 
    422446         END DO 
    423          IF( cp_cfg == 'orca' .AND. jp_cfg == 2 ) THEN 
    424             ii0 = 176   ;   ii1 =  176        ! Southern Island : Kerguelen 
    425             ij0 =  37   ;   ij1 =   37  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    426             ! 
    427             ii0 = 119   ;   ii1 =  119        ! South Georgia 
    428             ij0 =  29   ;   ij1 =   29  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    429             ! 
    430             ii0 = 111   ;   ii1 =  111        ! Falklands 
    431             ij0 =  35   ;   ij1 =   35  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    432             ! 
    433             ii0 = 168   ;   ii1 =  168        ! Crozet 
    434             ij0 =  40   ;   ij1 =   40  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    435             ! 
    436             ii0 = 119   ;   ii1 =  119        ! South Orkney 
    437             ij0 =  28   ;   ij1 =   28  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    438             ! 
    439             ii0 = 140   ;   ii1 =  140        ! Bouvet Island 
    440             ij0 =  33   ;   ij1 =   33  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    441             ! 
    442             ii0 = 178   ;   ii1 =  178        ! Prince edwards 
    443             ij0 =  34   ;   ij1 =   34  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    444             ! 
    445             ii0 =  43   ;   ii1 =   43        ! Balleny islands 
    446             ij0 =  21   ;   ij1 =   21  ;   zcmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) =  0.3_wp   
    447          ENDIF 
     447         ! 
    448448         CALL lbc_lnk( zcmask , 'T', 1. )      ! lateral boundary conditions on cmask   (sign unchanged) 
     449         ! 
    449450         DO jk = 1, jpk 
    450451            DO jj = 1, jpj 
Note: See TracChangeset for help on using the changeset viewer.