Ignore:
Timestamp:
2013-12-06T11:25:13+01:00 (7 years ago)
Author:
davestorkey
Message:

Remove OBC module at NEMO 3.6. See ticket #1189.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/SOL/solmat.F90

    r4153 r4328  
    2626   USE sol_oce         ! ocean solver 
    2727   USE phycst          ! physical constants 
    28    USE obc_oce         ! ocean open boundary conditions 
    2928   USE bdy_oce         ! unstructured open boundary conditions 
    3029   USE lbclnk          ! lateral boudary conditions 
     
    8685 
    8786#if defined key_dynspg_flt && ! defined key_bdy 
    88 #   if ! defined key_obc 
    8987 
    9088      DO jj = 2, jpjm1                      ! matrix of free surface elliptic system 
     
    103101         END DO 
    104102      END DO 
    105 #   else 
    106     IF ( Agrif_Root() ) THEN 
    107       DO jj = 2, jpjm1                      ! matrix of free surface elliptic system with open boundaries 
    108          DO ji = 2, jpim1 
    109             zcoef = z2dt * z2dt * grav * bmask(ji,jj) 
    110             !                                    ! south coefficient 
    111             IF( lp_obc_south .AND. ( jj == njs0p1 ) ) THEN 
    112                zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1)*(1.-vsmsk(ji,1)) 
    113             ELSE 
    114                zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1) 
    115             END IF 
    116             gcp(ji,jj,1) = zcoefs 
    117             ! 
    118             !                                    ! west coefficient 
    119             IF( lp_obc_west  .AND. ( ji == niw0p1 ) ) THEN 
    120                zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj)*(1.-uwmsk(jj,1)) 
    121             ELSE 
    122                zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj) 
    123             END IF 
    124             gcp(ji,jj,2) = zcoefw 
    125             ! 
    126             !                                    ! east coefficient 
    127             IF( lp_obc_east  .AND. ( ji == nie0 ) ) THEN 
    128                zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj)*(1.-uemsk(jj,1)) 
    129             ELSE 
    130                zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj) 
    131             END IF 
    132             gcp(ji,jj,3) = zcoefe 
    133             ! 
    134             !                                    ! north coefficient 
    135             IF( lp_obc_north .AND. ( jj == njn0 ) ) THEN 
    136                zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj)*(1.-vnmsk(ji,1)) 
    137             ELSE 
    138                zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj) 
    139             END IF 
    140             gcp(ji,jj,4) = zcoefn 
    141             ! 
    142             !                                    ! diagonal coefficient 
    143             gcdmat(ji,jj) = e1t(ji,jj)*e2t(ji,jj)*bmask(ji,jj)   & 
    144                &            - zcoefs -zcoefw -zcoefe -zcoefn 
    145          END DO 
    146       END DO 
    147     ELSE 
    148       DO jj = 2, jpjm1                      ! matrix of free surface elliptic system 
    149          DO ji = 2, jpim1 
    150             zcoef = z2dt * z2dt * grav * bmask(ji,jj) 
    151             zcoefs = -zcoef * hv(ji  ,jj-1) * e1v(ji  ,jj-1) / e2v(ji  ,jj-1)    ! south coefficient 
    152             zcoefw = -zcoef * hu(ji-1,jj  ) * e2u(ji-1,jj  ) / e1u(ji-1,jj  )    ! west coefficient 
    153             zcoefe = -zcoef * hu(ji  ,jj  ) * e2u(ji  ,jj  ) / e1u(ji  ,jj  )    ! east coefficient 
    154             zcoefn = -zcoef * hv(ji  ,jj  ) * e1v(ji  ,jj  ) / e2v(ji  ,jj  )    ! north coefficient 
    155             gcp(ji,jj,1) = zcoefs 
    156             gcp(ji,jj,2) = zcoefw 
    157             gcp(ji,jj,3) = zcoefe 
    158             gcp(ji,jj,4) = zcoefn 
    159             gcdmat(ji,jj) = e1t(ji,jj) * e2t(ji,jj) * bmask(ji,jj)    &          ! diagonal coefficient 
    160                &          - zcoefs -zcoefw -zcoefe -zcoefn 
    161          END DO 
    162       END DO 
    163     ENDIF 
    164 #   endif 
    165103 
    166104#  elif defined key_dynspg_flt && defined key_bdy  
Note: See TracChangeset for help on using the changeset viewer.