Changeset 1370 for trunk/NEMO/LIM_SRC_2

Timestamp:

2009-04-02T09:23:12+02:00 (15 years ago)

Author:

ctlod

Message:

bound salt exchange between Ice/Ocean? in Baltic Sea for ORCA2_LIM, see ticket: #388

File:

• trunk/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (4 diffs)

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -38 +38 @@
    REAL(wp)  ::   rzero  = 0.e0    
    REAL(wp)  ::   rone   = 1.e0
+   REAL(wp), DIMENSION(jpi,jpj)  ::   soce_r
+   REAL(wp), DIMENSION(jpi,jpj)  ::   sice_r
 
    !! * Substitutions
@@ -78 +80 @@
       INTEGER  ::   ifvt, i1mfr, idfr               ! some switches
       INTEGER  ::   iflt, ial, iadv, ifral, ifrdv
+      INTEGER  ::   ii0, ii1, ij0, ij1  ! temporary integers
       REAL(wp) ::   zqsr  , zqns     ! solar & non solar heat flux
       REAL(wp) ::   zinda            ! switch for testing the values of ice concentration
@@ -100 +103 @@
          IF(lwp) WRITE(numout,*) 'lim_sbc_2 : LIM 2.0 sea-ice - surface boundary condition'
          IF(lwp) WRITE(numout,*) '~~~~~~~~~   '
+         !
+         IF( cp_cfg == "orca"  .AND. jp_cfg == 2 ) THEN
+            !                                        ! =======================
+            !                                        !  ORCA_R2 configuration
+            !                                        ! =======================
+            soce_r(:,:) = soce
+            sice_r(:,:) = sice
+            ii0 = 145   ;   ii1 = 180        ! Baltic Sea
+            ij0 = 113   ;   ij1 = 130   ;   soce_r(mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 4.e0
+                                            sice_r(mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 2.e0
+         ENDIF
+         !
       ENDIF
 
@@ -204 +219 @@
 
             !  computing salt exchanges at the ice/ocean interface
-            zfons =  ( soce - sice ) * ( rdmicif(ji,jj) / rdt_ice ) 
+            zfons =  ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) / rdt_ice ) 
             
             !  converting the salt flux from ice to a freshwater flux from ocean