Changeset 991 for branches/dev_003_CPL/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

Timestamp:

2008-05-23T17:55:55+02:00 (16 years ago)

Author:

smasson

Message:

dev_003_CPL: preliminary draft (not working), see ticket #155

File:

• branches/dev_003_CPL/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (8 diffs)

branches/dev_003_CPL/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -12 +12 @@
    !!   'key_lim2' :                                  LIM 2.0 sea-ice model
    !!----------------------------------------------------------------------
+   !!----------------------------------------------------------------------
    !!   sbc_ice_lim_2  : sea-ice model time-stepping and
    !!                    update ocean sbc over ice-covered area
    !!----------------------------------------------------------------------
    USE oce             ! ocean dynamics and tracers
-   USE c1d             ! 1d configuration
    USE dom_oce         ! ocean space and time domain
    USE ice_2
@@ -46 +46 @@
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
-   USE ocfzpt          ! ocean freezing point
 
    IMPLICIT NONE
@@ -59 +58 @@
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/SBC  3.0 , LOCEAN-IPSL (2008) 
-   !! $Id: $
+   !!   OPA 9.0 , LOCEAN-IPSL (2006) 
+   !! $ Id: $
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -90 +89 @@
       !!
       INTEGER  ::   ji, jj   ! dummy loop indices
-      REAL(wp) ::   zinda    
       REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_os   ! albedo of the ice under overcast sky
       REAL(wp), DIMENSION(jpi,jpj,1) ::   alb_ice_cs   ! albedo of ice under clear sky
@@ -126 +124 @@
 
          ! ... masked sea surface freezing temperature [Kelvin] (set to rt0 over land)
-         tfu(:,:) = tfreez( sss_m ) +  rt0 
-
-         zsist (:,:,1) = sist (:,:)
-         zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:)
+         tfu(:,:) = tfreez( sss_m ) + rt0 
+
 
          ! ... ice albedo
-         CALL albedo_ice( zsist, zhicif, zhsnif, alb_ice_cs, alb_ice_os )
-
+         
+         IF ( nsbc /= 5 ) THEN
++++ INTERFACE 2D/3D suprimer les tableaux intermediaires
++++ les mettre sous cle cpp
++++ il faudrait utiliser les variables de transfert pour tn_ice, ice/snow thickness, albedo...
+
+            zsist (:,:,1) = sist (:,:)
+            zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:)
+            CALL albedo_ice( zsist, zhicif, zhsnif, alb_ice_cs, alb_ice_os )
          ! ... Sea-ice surface boundary conditions output from bulk formulae :
          !     - utaui_ice  ! surface ice stress i-component (I-point)   [N/m2]
@@ -154 +157 @@
                &                               tprecip   , sprecip    ,                          &
                &                               fr1_i0    , fr2_i0     , cl_grid  )
-
-            ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed
-            DO jj = 1, jpj
-               DO ji = 1, jpi
-                  zinda    = MAX(  0.e0, SIGN(  1.e0, -( -1.5 - freeze(ji,jj) )  )  )
-                  qsr(ji,jj) = zinda * qsr(ji,jj)
-               END DO
-            END DO
-
          CASE( 4 )           ! CORE bulk formulation
             CALL blk_ice_core( zsist , ui_ice , vi_ice   , alb_ice_cs ,                         &
@@ -175 +169 @@
          qla_ice(:,:) = zqla_ice(:,:,1)   ;   dqla_ice(:,:) = zdqla_ice(:,:,1)
 
+         ENDIF
          IF(ln_ctl) THEN         ! print mean trends (used for debugging)
             CALL prt_ctl_info( 'Ice Forcings ' )
@@ -188 +183 @@
          !  Ice model step  !
          ! ---------------- !
-                                        CALL lim_rst_opn_2  ( kt )      ! Open Ice restart file
-         IF( .NOT. lk_c1d ) THEN                                        ! Ice dynamics & transport (not in 1D case)
-                                        CALL lim_dyn_2      ( kt )           ! Ice dynamics    ( rheology/dynamics )
-                                        CALL lim_trp_2      ( kt )           ! Ice transport   ( Advection/diffusion )
-            IF( ln_limdmp )             CALL lim_dmp_2      ( kt )           ! Ice damping 
-         ENDIF
-                                        CALL lim_thd_2      ( kt )      ! Ice thermodynamics 
-                                        CALL lim_sbc_2      ( kt )      ! Ice/Ocean Mass & Heat fluxes 
+         ;                              CALL lim_rst_opn_2  ( kt )      ! Open Ice restart file
+         ;                              CALL lim_dyn_2      ( kt )      ! Ice dynamics    ( rheology/dynamics )
+         ;                              CALL lim_trp_2      ( kt )      ! Ice transport   ( Advection/diffusion )
+         IF( ln_limdmp )                CALL lim_dmp_2      ( kt )      ! Ice damping 
+         ;                              CALL lim_thd_2      ( kt )      ! Ice thermodynamics 
+         ;                              CALL lim_sbc_2      ( kt )      ! Ice/Ocean Mass & Heat fluxes 
          IF( MOD( kt+nn_fsbc-1, ninfo ) == 0 .OR.   &
             &  ntmoy == 1 )             CALL lim_dia_2      ( kt )      ! Ice Diagnostics 
-                                        CALL lim_wri_2      ( kt )      ! Ice outputs 
+         ;                              CALL lim_wri_2      ( kt )      ! Ice outputs 
          IF( lrst_ice )                 CALL lim_rst_write_2( kt )      ! Ice restart file 
          !