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

Timestamp:

2008-05-30T18:21:28+02:00 (16 years ago)

Author:

ctlod

Message:

trunk: replace freeze(:,:) variable with fr_i(:,:), use the tfreez function defined in eosbn2.F90 and remove the useless ocfzpt.F90 module, see ticket: #177

File:

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

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -20 +20 @@
    USE sbc_oce          ! surface boundary condition
    USE phycst           ! physical constants
-   USE ocfzpt           ! surface ocean freezing point
    USE ice_oce          ! sea-ice variable
    USE ice_2            ! LIM sea-ice variables
@@ -61 +60 @@
       !!              - Update 
       !!     
-      !! ** Outputs : - qsr    : sea heat flux:     solar 
-      !!              - qns    : sea heat flux: non solar
-      !!              - emp    : freshwater budget: volume flux 
-      !!              - emps   : freshwater budget: concentration/dillution 
-      !!              - utau   : sea surface i-stress (ocean referential)
-      !!              - vtau   : sea surface j-stress (ocean referential)
+      !! ** Outputs : - qsr     : sea heat flux:     solar 
+      !!              - qns     : sea heat flux: non solar
+      !!              - emp     : freshwater budget: volume flux 
+      !!              - emps    : freshwater budget: concentration/dillution 
+      !!              - utau    : sea surface i-stress (ocean referential)
+      !!              - vtau    : sea surface j-stress (ocean referential)
+      !!              - fr_i    : ice fraction
+      !!              - tn_ice  : sea-ice surface temperature
+      !!              - alb_ice : sea-ice alberdo (lk_cpl=T)
       !!
       !! References : Goosse, H. et al. 1996, Bul. Soc. Roy. Sc. Liege, 65, 87-90.
@@ -214 +216 @@
       !-----------------------------------------------!
 
-      freeze(:,:) = 1.0 - frld(:,:)       ! Sea ice cover            
-      tn_ice(:,:) = sist(:,:)             ! Ice surface temperature                      
+      fr_i  (:,:) = 1.0 - frld(:,:)       ! sea-ice fraction
+      tn_ice(:,:) = sist(:,:)             ! sea-ice surface temperature                      
 
 #if defined key_coupled            
@@ -234 +236 @@
          CALL prt_ctl(tab2d_1=utau  , clinfo1=' lim_sbc: utau   : ', mask1=umask,   &
             &         tab2d_2=vtau  , clinfo2=' vtau    : '        , mask2=vmask )
-         CALL prt_ctl(tab2d_1=freeze, clinfo1=' lim_sbc: freeze : ', tab2d_2=tn_ice, clinfo2=' tn_ice  : ')
+         CALL prt_ctl(tab2d_1=fr_i  , clinfo1=' lim_sbc: fr_i   : ', tab2d_2=tn_ice, clinfo2=' tn_ice  : ')
       ENDIF