Changeset 1037 for trunk/NEMO/LIM_SRC_2

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

Location:

trunk/NEMO/LIM_SRC_2

Files:

• iceini_2.F90 (modified) (2 diffs)
• limistate_2.F90 (modified) (3 diffs)
• limsbc_2.F90 (modified) (4 diffs)

trunk/NEMO/LIM_SRC_2/iceini_2.F90

@@ -21 +21 @@
    USE sbc_ice         ! surface boundary condition: ice
    USE phycst          ! Define parameters for the routines
-   USE ocfzpt
    USE ice_2
    USE limmsh_2
@@ -75 +74 @@
       
       tn_ice(:,:) = sist(:,:)         ! initialisation of ice temperature   
-      freeze(:,:) = 1.0 - frld(:,:)   ! initialisation of sea/ice cover    
+      fr_i  (:,:) = 1.0 - frld(:,:)   ! initialisation of sea-ice fraction    
 # if defined key_coupled
       alb_ice(:,:) = albege(:,:)      ! sea-ice albedo

trunk/NEMO/LIM_SRC_2/limistate_2.F90

@@ -19 +19 @@
    !!----------------------------------------------------------------------
    USE phycst
-   USE ocfzpt
    USE par_ice_2       ! ice parameters
    USE ice_oce         ! ice variables
    USE dom_ice_2
+   USE eosbn2          ! equation of state
    USE lbclnk
    USE oce
@@ -72 +72 @@
       IF( .NOT. ln_limini ) THEN  
          
-         ! Initialisation at tn or -2 if ice
+         tfu(:,:) = tfreez( sn(:,:,1) ) * tmask(:,:,1)       ! freezing/melting point of sea water [Celcius]
+
          DO jj = 1, jpj
             DO ji = 1, jpi
-               zbin = MAX( 0., SIGN( 1., fzptn(ji,jj) - tn(ji,jj,1) ) )
-               ztn(ji,jj) = ( (1.-zbin) * tn(ji,jj,1) - 2. * zbin + rt0 ) * tmask(ji,jj,1)
-            END DO
-         END DO
-                  
-         !  tfu: Melting point of sea water [Kelvin]
-         zs0 = 34.e0
-         ztf = rt0 + ( - 0.0575 + 1.710523e-3 * SQRT( zs0 ) - 2.154996e-4 * zs0 ) * zs0
-         tfu(:,:) = ztf
-         
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-               !--- Criterion for presence (zidto=1) or absence (zidto=0) of ice
-               zidto  = tms(ji,jj) * ( 1.0 - MAX(zzero, SIGN( zone, ztn(ji,jj) - tfu(ji,jj) - ttest) ) )
-               
+               !                     ! ice if sst <= t-freez + ttest
+               IF( tn(ji,jj,1)  - tfu(ji,jj) >= ttest ) THEN   ;   zidto = 0.e0      ! no ice
+               ELSE                                            ;   zidto = 1.e0      !    ice
+               ENDIF
+               !
                IF( fcor(ji,jj) >= 0.e0 ) THEN     !--  Northern hemisphere.
                   hicif(ji,jj)   = zidto * hginn
@@ -101 +92 @@
             END DO
          END DO
+
+         tfu(:,:) = tfu(:,:) + rt0       ! ftu converted from Celsius to Kelvin (rt0 over land)
          
          sist  (:,:)   = tfu(:,:)

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