Changeset 1037 for trunk/NEMO/LIM_SRC_3/limsbc.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_3/limsbc.F90 (modified) (4 diffs)

trunk/NEMO/LIM_SRC_3/limsbc.F90

@@ -21 +21 @@
    USE sbc_oce          ! Surface boundary condition: ocean fields
    USE phycst           ! physical constants
-   USE ocfzpt           ! surface ocean freezing point
    USE ice_oce          ! sea-ice variable
    USE ice              ! LIM sea-ice variables
@@ -225 +224 @@
       !!              - Update the ocean sbc
       !!     
-      !! ** Outputs : - qsr    : sea heat flux:     solar 
-      !!              - qns    : sea heat flux: non solar
-      !!              - emp    : freshwater budget: volume flux 
-      !!              - emps   : freshwater budget: concentration/dillution 
+      !! ** Outputs : - qsr     : sea heat flux:     solar 
+      !!              - qns     : sea heat flux: non solar
+      !!              - emp     : freshwater budget: volume flux 
+      !!              - emps    : freshwater budget: concentration/dillution 
+      !!              - 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.
@@ -422 +424 @@
       !-----------------------------------------------!
 
-      freeze(:,:)   = at_i(:,:)             ! Sea ice cover            
+      fr_i  (:,:)   = at_i(:,:)             ! Sea-ice fraction            
       tn_ice(:,:,:) = t_su(:,:,:)           ! Ice surface temperature                      
 
@@ -440 +442 @@
          CALL prt_ctl( tab2d_1=qsr   , clinfo1=' lim_sbc: qsr    : ', tab2d_2=qns , clinfo2=' qns     : ' )
          CALL prt_ctl( tab2d_1=emp   , clinfo1=' lim_sbc: emp    : ', tab2d_2=emps, clinfo2=' emps    : ' )
-         CALL prt_ctl( tab2d_1=freeze, clinfo1=' lim_sbc: freeze : ' )
+         CALL prt_ctl( tab2d_1=fr_i  , clinfo1=' lim_sbc: fr_i   : ' )
          CALL prt_ctl( tab3d_1=tn_ice, clinfo1=' lim_sbc: tn_ice : ', kdim=jpl )
       ENDIF