Changeset 1928


Ignore:
Timestamp:
2010-06-11T19:07:07+02:00 (10 years ago)
Author:
cetlod
Message:

bound the salt exchange at the Ice/Ocean? interface in the Baltic Sea ; generalisation for all ORCA configurations, see ticket:388

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_2/limsbc_2.F90

    r1756 r1928  
    8080      INTEGER  ::   ifvt, i1mfr, idfr               ! some switches 
    8181      INTEGER  ::   iflt, ial, iadv, ifral, ifrdv 
    82       INTEGER  ::   ii0, ii1, ij0, ij1  ! temporary integers 
    8382      REAL(wp) ::   zrdtir           ! 1. / rdt_ice 
    8483      REAL(wp) ::   zqsr  , zqns     ! solar & non solar heat flux 
     
    109108         sice_r(:,:) = sice 
    110109         ! 
    111          IF( cp_cfg == "orca"  .AND. jp_cfg == 2 ) THEN 
    112             !                                        ! ======================= 
    113             !                                        !  ORCA_R2 configuration 
    114             !                                        ! ======================= 
    115             ii0 = 145   ;   ii1 = 180        ! Baltic Sea 
    116             ij0 = 113   ;   ij1 = 130   ;   soce_r(mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 4.e0 
    117                                             sice_r(mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 2.e0 
    118          ENDIF 
    119          ! 
    120       ENDIF 
     110         IF( cp_cfg == "orca" ) THEN 
     111           !   ocean/ice salinity in the Baltic sea  
     112           DO jj = 1, jpj 
     113              DO ji = 1, jpi 
     114                 IF( glamt(ji,jj) >= 14. .AND.  glamt(ji,jj) <= 32. .AND. gphit(ji,jj) >= 54. .AND. gphit(ji,jj) <= 66. ) THEN  
     115                   soce_r(ji,jj) = 4.e0  
     116                   sice_r(ji,jj) = 2.e0 
     117                 END IF 
     118              END DO 
     119           END DO 
     120           ! 
     121         END IF 
     122      END IF 
    121123 
    122124      !------------------------------------------! 
Note: See TracChangeset for help on using the changeset viewer.