Changeset 2119

Timestamp:

2010-09-24T15:19:32+02:00 (14 years ago)

Author:

cbricaud

Message:

bugfix

Location:

branches/dev_1784_EVP/NEMO/LIM_SRC_2

Files:

• dom_ice_2.F90 (modified) (1 diff)
• limsbc_2.F90 (modified) (3 diffs)

branches/dev_1784_EVP/NEMO/LIM_SRC_2/dom_ice_2.F90

@@ -22 +22 @@
 
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   fs2cor , fcor     !: coriolis factor and coeficient
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   covrai ,          !: sine of geographic latitude
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   area   ,          !: surface of grid cell 
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   covrai            !: sine of geographic latitude
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   area              !: surface of grid cell 
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj)         ::   tms    , tmu      !: temperature and velocity points masks
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,2,2)     ::   wght              !: weight of the 4 neighbours to compute averages

branches/dev_1784_EVP/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -35 +35 @@
    PUBLIC   lim_sbc_2   ! called by sbc_ice_lim_2
 
-   REAL(wp)  ::   r1_rdtice = 1.e0 / rdt_ice   ! constant values
+   REAL(wp)  ::   r1_rdtice                    ! constant values
    REAL(wp)  ::   epsi16 = 1.e-16              !     -      -
    REAL(wp)  ::   rzero  = 0.e0                !     -      -
@@ -204 +204 @@
 #endif            
             ! salt exchanges at the ice/ocean interface
-            zfons =  ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) * zrdtir ) 
+            zfons =  ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) * r1_rdtice ) 
             !
             ! convert the salt flux from ice into a freshwater flux from ocean
@@ -215 +215 @@
       !
       IF( lk_diaar5 ) THEN
-         CALL iom_put( 'isnwmlt_cea'  ,                 rdmsnif(:,:) * zrdtir )
-         CALL iom_put( 'fsal_virt_cea',   soce_r(:,:) * rdmicif(:,:) * zrdtir )
-         CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdmicif(:,:) * zrdtir )
+         CALL iom_put( 'isnwmlt_cea'  ,                 rdmsnif(:,:) * r1_rdtice )
+         CALL iom_put( 'fsal_virt_cea',   soce_r(:,:) * rdmicif(:,:) * r1_rdtice )
+         CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdmicif(:,:) * r1_rdtice )
       ENDIF