Changeset 1218 for trunk/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2008-10-28T10:12:16+01:00 (16 years ago)

Author:

smasson

Message:

first implementation of the new coupling interface in the trunk, see ticket:155

File:

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

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -29 +29 @@
    USE albedo           ! albedo parameters
    USE prtctl           ! Print control
+   USE cpl_oasis3, ONLY : lk_cpl
 
    IMPLICIT NONE
@@ -85 +86 @@
       REAL(wp) ::   zutau , zvtau    ! lead fraction at U- & V-points
       REAL(wp) ::   zu_io , zv_io    ! 2 components of the ice-ocean velocity
-#if defined key_coupled    
-      REAL(wp), DIMENSION(jpi,jpj) ::   zalb     ! albedo of ice under overcast sky
-      REAL(wp), DIMENSION(jpi,jpj) ::   zalbp    ! albedo of ice under clear sky
-#endif
+! interface 2D --> 3D
+      REAL(wp), DIMENSION(jpi,jpj,1) ::   zalb     ! albedo of ice under overcast sky
+      REAL(wp), DIMENSION(jpi,jpj,1) ::   zalbp    ! albedo of ice under clear sky
+      REAL(wp), DIMENSION(jpi,jpj,1) ::   zsist    ! surface ice temperature (K)
+      REAL(wp), DIMENSION(jpi,jpj,1) ::   zhicif   ! ice thickness
+      REAL(wp), DIMENSION(jpi,jpj,1) ::   zhsnif   ! snow thickness
       REAL(wp) ::   zsang, zmod, zfm
       REAL(wp), DIMENSION(jpi,jpj) ::   ztio_u, ztio_v   ! ocean stress below sea-ice
@@ -119 +122 @@
             ifral   = ( 1  - i1mfr * ( 1 - ial ) )   
             ifrdv   = ( 1  - ifral * ( 1 - ial ) ) * iadv 
+
+!!$            zinda   = 1.0 - AINT( pfrld(ji,jj) )                   !   = 0. if pure ocean else 1. (at previous time)
+!!$
+!!$            i1mfr   = 1.0 - AINT(  frld(ji,jj) )                   !   = 0. if pure ocean else 1. (at current  time)
+!!$
+!!$            IF( phicif(ji,jj) <= 0. ) THEN   ;   ifvt = zinda      !   = 1. if (snow and no ice at previous time) else 0. ???
+!!$            ELSE                             ;   ifvt = 0.
+!!$            ENDIF
+!!$
+!!$            IF( frld(ji,jj) >= pfrld(ji,jj) ) THEN   ;   idfr = 0.  !   = 0. if lead fraction increases from previous to current
+!!$            ELSE                                     ;   idfr = 1.   
+!!$            ENDIF
+!!$
+!!$            iflt    = zinda  * (1 - i1mfr) * (1 - ifvt )    !   = 1. if ice (not only snow) at previous and pure ocean at current
+!!$
+!!$            ial     = ifvt   * i1mfr    +    ( 1 - ifvt ) * idfr
+!!$!                 snow no ice   ice         ice or nothing  lead fraction increases
+!!$!                 at previous   now           at previous
+!!$!                -> ice aera increases  ???         -> ice aera decreases ???
+!!$
+!!$            iadv    = ( 1  - i1mfr ) * zinda  
+!!$!                     pure ocean      ice at
+!!$!                     at current      previous
+!!$!                        -> = 1. if ice disapear between previous and current
+!!$
+!!$            ifral   = ( 1  - i1mfr * ( 1 - ial ) )  
+!!$!                            ice at     ???
+!!$!                            current         
+!!$!                         -> ???
+!!$ 
+!!$            ifrdv   = ( 1  - ifral * ( 1 - ial ) ) * iadv 
+!!$!                                                    ice disapear                           
+!!$
+!!$
+
             !   computation the solar flux at ocean surface
-            zqsr    = pfrld(ji,jj) * qsr(ji,jj)  + ( 1. - pfrld(ji,jj) ) * fstric(ji,jj)
+#if defined key_coupled 
+            zqsr = tqsr(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj) ) * ( 1.0 - pfrld(ji,jj) )
+#else
+            zqsr = pfrld(ji,jj) * qsr(ji,jj)  + ( 1.  - pfrld(ji,jj) ) * fstric(ji,jj)
+#endif            
             !  computation the non solar heat flux at ocean surface
             zqns    =  - ( 1. - thcm(ji,jj) ) * zqsr   &   ! part of the solar energy used in leads
@@ -145 +187 @@
          DO ji = 1, jpi
             
+#if defined key_coupled
+          zemp = emp_tot(ji,jj) - emp_ice(ji,jj) * ( 1. - pfrld(ji,jj) )    &   ! 
+             &   + rdmsnif(ji,jj) / rdt_ice                                     !  freshwaterflux due to snow melting 
+#else
+!!$            !  computing freshwater exchanges at the ice/ocean interface
+!!$            zpme = - evap(ji,jj)    *   frld(ji,jj)           &   !  evaporation over oceanic fraction
+!!$               &   + tprecip(ji,jj)                           &   !  total precipitation
+!!$               &   - sprecip(ji,jj) * ( 1. - pfrld(ji,jj) )   &   !  remov. snow precip over ice
+!!$               &   - rdmsnif(ji,jj) / rdt_ice                     !  freshwaterflux due to snow melting 
             !  computing freshwater exchanges at the ice/ocean interface
             zemp = + emp(ji,jj)     *         frld(ji,jj)      &   !  e-p budget over open ocean fraction 
@@ -151 +202 @@
                &   + rdmsnif(ji,jj) / rdt_ice                      !  freshwaterflux due to snow melting 
                !                                                   !  ice-covered fraction:
+#endif            
 
             !  computing salt exchanges at the ice/ocean interface
@@ -217 +269 @@
 
       fr_i  (:,:) = 1.0 - frld(:,:)       ! sea-ice fraction
-      tn_ice(:,:) = sist(:,:)             ! sea-ice surface temperature                      
-
-#if defined key_coupled            
-      !------------------------------------------------!
-      !    Computation of snow/ice and ocean albedo    !
-      !------------------------------------------------!
-      zalb  (:,:) = 0.e0
-      zalbp (:,:) = 0.e0
-
-      CALL albedo_ice( sist, hicif, hsnif, zalbp, zalb )
-
-      alb_ice(:,:) =  0.5 * zalbp(:,:) + 0.5 * zalb (:,:)   ! Ice albedo (mean clear and overcast skys)
-#endif
+
+      IF ( lk_cpl ) THEN           
+         ! Ice surface temperature 
+         tn_ice(:,:) = sist(:,:)          ! sea-ice surface temperature       
+         ! Computation of snow/ice and ocean albedo
+         ! INTERFACE 3D versus 2D
+         zsist (:,:,1) = sist (:,:)
+         zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:)
+         CALL albedo_ice( zsist, zhicif, zhsnif, zalbp, zalb )
+         alb_ice(:,:) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
+      ENDIF
 
       IF(ln_ctl) THEN