Changeset 1218 for trunk/NEMO/OPA_SRC/SBC/sbcice_lim_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/OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (9 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -6 +6 @@
    !! Sea-Ice model  :  LIM 2.0 Sea ice model time-stepping
    !!======================================================================
-   !! History :  9.0   !  06-06  (G. Madec)  from icestp_2.F90
+   !! History :  1.0   !  06-2006  (G. Madec)  from icestp_2.F90
+   !!            3.0   !  08-2008  (S. Masson, E. .... ) coupled interface
    !!----------------------------------------------------------------------
 #if defined key_lim2
@@ -23 +24 @@
    USE ice_oce         ! ice variables
    USE dom_ice_2
-   USE cpl_oce
 
    USE sbc_oce         ! Surface boundary condition: ocean fields
@@ -29 +29 @@
    USE sbcblk_core     ! Surface boundary condition: CORE bulk
    USE sbcblk_clio     ! Surface boundary condition: CLIO bulk
+   USE sbccpl          ! Surface boundary condition: coupled interface
    USE albedo
 
@@ -66 +67 @@
 CONTAINS
 
-   SUBROUTINE sbc_ice_lim_2( kt, kblk )
+   SUBROUTINE sbc_ice_lim_2( kt, ksbc )
       !!---------------------------------------------------------------------
       !!                  ***  ROUTINE sbc_ice_lim_2  ***
@@ -87 +88 @@
       !!---------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt      ! ocean time step
-      INTEGER, INTENT(in) ::   kblk    ! type of bulk (=3 CLIO, =4 CORE)
+      INTEGER, INTENT(in) ::   ksbc    ! type of sbc ( =3 CLIO bulk ; =4 CORE bulk ; =5 coupled )
       !!
       INTEGER  ::   ji, jj   ! dummy loop indices
@@ -131 +132 @@
          zhsnif(:,:,1) = hsnif(:,:)
 
-         ! ... ice albedo
+         ! ... ice albedo (clear sky and overcast sky)
          CALL albedo_ice( zsist, zhicif, zhsnif, alb_ice_cs, alb_ice_os )
 
@@ -147 +148 @@
          !     - fr2_i0     ! 2nd fraction of qsr penetration in ice     [%]
          !
-         SELECT CASE( kblk )
+         SELECT CASE( ksbc )
          CASE( 3 )           ! CLIO bulk formulation
-            CALL blk_ice_clio( zsist , alb_ice_cs , alb_ice_os ,                                 &
-               &                               utaui_ice , vtaui_ice  , zqns_ice   , zqsr_ice,   &
-               &                               zqla_ice  , zdqns_ice  , zdqla_ice  ,             &
-               &                               tprecip   , sprecip    ,                          &
-               &                               fr1_i0    , fr2_i0     , cl_grid  )
+            CALL blk_ice_clio( zsist, alb_ice_cs, alb_ice_os ,                         &
+               &                      utaui_ice , vtaui_ice  , zqns_ice  , zqsr_ice,   &
+               &                      zqla_ice  , zdqns_ice  , zdqla_ice ,             &
+               &                      tprecip   , sprecip    ,                         &
+               &                      fr1_i0    , fr2_i0     , cl_grid  )
 
          CASE( 4 )           ! CORE bulk formulation
-            CALL blk_ice_core( zsist , ui_ice , vi_ice   , alb_ice_cs ,                         &
-               &                               utaui_ice , vtaui_ice  , zqns_ice  , zqsr_ice,   &
-               &                               zqla_ice  , zdqns_ice  , zdqla_ice ,             &
-               &                               tprecip   , sprecip    ,                         &
-               &                               fr1_i0    , fr2_i0     , cl_grid)
+            CALL blk_ice_core( zsist, ui_ice    , vi_ice     , alb_ice_cs,             &
+               &                      utaui_ice , vtaui_ice  , zqns_ice  , zqsr_ice,   &
+               &                      zqla_ice  , zdqns_ice  , zdqla_ice ,             &
+               &                      tprecip   , sprecip    ,                         &
+               &                      fr1_i0    , fr2_i0     , cl_grid  )
+         CASE( 5 )           ! Coupled formulation : atmosphere-ice stress only (fluxes provided after ice dynamics)
+            CALL sbc_cpl_ice_tau( utaui_ice , vtaui_ice )
          END SELECT
 
@@ -186 +189 @@
             IF( ln_limdmp )             CALL lim_dmp_2      ( kt )           ! Ice damping 
          ENDIF
+#if defined key_coupled
+         IF( ksbc == 5    )             CALL sbc_cpl_ice_flx( frld, alb_ice_cs , sst_m, sist,   &
+      &                                                             qns_tot, qns_ice,   &
+      &                                                             qsr_tot, qsr_ice,   &
+      &                                                             emp_tot, emp_ice, sprecip )
+#endif
                                         CALL lim_thd_2      ( kt )      ! Ice thermodynamics 
                                         CALL lim_sbc_2      ( kt )      ! Ice/Ocean Mass & Heat fluxes 
@@ -202 +211 @@
    !!----------------------------------------------------------------------
 CONTAINS
-   SUBROUTINE sbc_ice_lim_2 ( kt, kblk )     ! Dummy routine
-      WRITE(*,*) 'sbc_ice_lim_2: You should not have seen this print! error?', kt, kblk
+   SUBROUTINE sbc_ice_lim_2 ( kt, ksbc )     ! Dummy routine
+      WRITE(*,*) 'sbc_ice_lim_2: You should not have seen this print! error?', kt, ksbc
    END SUBROUTINE sbc_ice_lim_2
 #endif