Changeset 878 for branches

Timestamp:

2008-04-04T14:53:06+02:00 (16 years ago)

Author:

ctlod

Message:

dev_001_SBC: clean and add freeze(:,:) & tn_ice(:,:) update, see ticket:#103

File:

• branches/dev_001_SBC/NEMO/LIM_SRC/limsbc.F90 (modified) (4 diffs)

branches/dev_001_SBC/NEMO/LIM_SRC/limsbc.F90

@@ -152 +152 @@
                !                                                   !  ice-covered fraction:
 
-!!gm old    zemp = + evap(ji,jj) * frld(ji,jj)                 &   !  evaporation over oceanic fraction
-!!gm old       &   - tprecip(ji,jj)                            &   !  total precipitation
-!!gm old       &   + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) )    &   !  remov. snow precip over ice (store over sea-ice)
-!!gm old       &   + rdmsnif(ji,jj) / rdt_ice                      !  freshwaterflux due to snow melting 
-            
             !  computing salt exchanges at the ice/ocean interface
             zfons =  ( soce - sice ) * ( rdmicif(ji,jj) / rdt_ice ) 
             
             !  converting the salt flux from ice to a freshwater flux from ocean
-!!gm old    zfm  = zfons / ( sss_io(ji,jj) + epsi16 )
             zfm  = zfons / ( sss_m(ji,jj) + epsi16 )
             
@@ -185 +179 @@
             DO ji = 1, jpi
                ! ... change the cosinus angle sign in the south hemisphere
-!!              zsang  = SIGN(1.e0, gphif(ji-1,jj-1) ) * sangvg    ! do the full loop and avoid lbc_lnk
                zsang  = SIGN(1.e0, gphif(ji,jj) ) * sangvg
                ! ... ice velocity relative to the ocean
@@ -219 +212 @@
       ENDIF
 
+      !-----------------------------------------------!
+      !   Storing the transmitted variables           !
+      !-----------------------------------------------!
+
+      freeze(:,:) = 1.0 - frld(:,:)       ! Sea ice cover            
+      tn_ice(:,:) = sist(:,:)             ! Ice surface temperature                      
+
 #if defined key_coupled            
       !------------------------------------------------!
@@ -234 +234 @@
 
       IF(ln_ctl) THEN
-         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=utau, clinfo1=' lim_sbc: utau : ', tab2d_2=vtau, clinfo2=' vtau : ')
-!!gm tn_ice ???  what is it now?
-!!gm     CALL prt_ctl(tab2d_1=freeze, clinfo1=' lim_sbc: freeze : ', tab2d_2=tn_ice , clinfo2=' tn_ice  : ')
+         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=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  : ')
       ENDIF