Changeset 918 for trunk/NEMO/OPA_SRC/SBC/sbcmod.F90

Timestamp:

2008-05-09T12:00:09+02:00 (16 years ago)

Author:

rblod

Message:

Change ice_ocean stress computation, step II, see ticket #128

File:

• trunk/NEMO/OPA_SRC/SBC/sbcmod.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -55 +55 @@
    INTEGER , PUBLIC ::   nn_ice      = 0         !: flag on ice in the surface boundary condition (=0/1/2)
    INTEGER , PUBLIC ::   nn_fwb      = 0         !: type of FreshWater Budget control (=0/1/2)
+   INTEGER          ::   nn_ico_cpl = 0    ! ice-ocean coupling indicator:         !!gm   ===>> to be put in namelist
+   !                                       !  = 0   LIM-3 old case
+   !                                       !  = 1   stresses computed using now ocean velocity
+   !                                       !  = 2   combination of 0 and 1 cases
+
    NAMELIST/namsbc/ nn_fsbc, ln_ana, ln_flx, ln_blk_clio, ln_blk_core, ln_cpl,   &
-      &             nn_ice , ln_dm2dc, ln_rnf, ln_ssr, nn_fwb
+      &             nn_ice , ln_dm2dc, ln_rnf, ln_ssr, nn_fwb, nn_ico_cpl
 
    INTEGER ::   nsbc   ! type of surface boundary condition (deduced from namsbc informations)
@@ -116 +121 @@
          WRITE(numout,*) '           Misc. options of sbc : '
          WRITE(numout,*) '              ice management in the sbc (=0/1/2)         nn_ice      = ', nn_ice 
+         WRITE(numout,*) '              ice-ocean stress computation (=0/1/2)      nn_ico_cpl  = ', nn_ico_cpl
          WRITE(numout,*) '              daily mean to diurnal cycle qsr            ln_dm2dc    = ', ln_dm2dc 
          WRITE(numout,*) '              runoff / runoff mouths                     ln_rnf      = ', ln_rnf
@@ -228 +234 @@
       
       SELECT CASE( nn_ice )                                     ! Update sbc over ice-covered areas
-      CASE(  1 )   ;         CALL sbc_ice_if ( kt )                  ! Ice-cover climatology ("Ice-if" model)
-         !                                                           ! (update heat and freshwater fluxes)
-      CASE(  2 )   ;         CALL sbc_ice_lim_2( kt, nsbc )          ! LIM 2.0 ice model
-         !                                                           ! (update heat and freshwater fluxes)
-      CASE(  3 )   ;         CALL sbc_ice_lim  ( kt, nsbc )          ! LIM 3.0 ice model
-      END SELECT                                                     ! (update all fluxes using bulk + LIM)
+      CASE(  1 )   ;         CALL sbc_ice_if ( kt )                     ! Ice-cover climatology ("Ice-if" model)
+         !                                                              ! (update heat and freshwater fluxes)
+      CASE(  2 )   ;         CALL sbc_ice_lim_2( kt, nsbc )             ! LIM 2.0 ice model
+         !                                                              ! (update heat and freshwater fluxes)
+      CASE(  3 )   ;         CALL sbc_ice_lim  ( kt, nsbc, nn_ico_cpl)          ! LIM 3.0 ice model
+      END SELECT                                                        ! (update all fluxes using bulk + LIM)
 
       IF( ln_ssr         )   CALL sbc_ssr( kt )                 ! add SST/SSS damping term