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

Timestamp:

2009-06-09T16:45:31+02:00 (15 years ago)

Author:

smasson

Message:

force 3rd dimension in sea-ice coupling fields, see ticket:444

File:

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

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -173 +173 @@
             !   computation the solar flux at ocean surface
 #if defined key_coupled 
-            zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj) ) * ( 1.0 - pfrld(ji,jj) )
+            zqsr = qsr_tot(ji,jj) + ( fstric(ji,jj) - qsr_ice(ji,jj,1) ) * ( 1.0 - pfrld(ji,jj) )
 #else
             zqsr = pfrld(ji,jj) * qsr(ji,jj)  + ( 1.  - pfrld(ji,jj) ) * fstric(ji,jj)
@@ -286 +286 @@
       IF ( lk_cpl ) THEN           
          ! Ice surface temperature 
-         tn_ice(:,:) = sist(:,:)          ! sea-ice surface temperature       
+         tn_ice(:,:,1) = sist(:,:)          ! sea-ice surface temperature       
          ! Computation of snow/ice and ocean albedo
          ! INTERFACE 3D versus 2D
@@ -292 +292 @@
          zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:)
          CALL albedo_ice( zsist, zhicif, zhsnif, zalbp, zalb )
-         alb_ice(:,:) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
+         alb_ice(:,:,1) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
       ENDIF
 
@@ -300 +300 @@
          CALL prt_ctl(tab2d_1=utau  , clinfo1=' lim_sbc: utau   : ', mask1=umask,   &
             &         tab2d_2=vtau  , clinfo2=' vtau    : '        , mask2=vmask )
-         CALL prt_ctl(tab2d_1=fr_i  , clinfo1=' lim_sbc: fr_i   : ', tab2d_2=tn_ice, clinfo2=' tn_ice  : ')
+         CALL prt_ctl(tab2d_1=fr_i  , clinfo1=' lim_sbc: fr_i   : ', tab2d_2=tn_ice(:,:,1), clinfo2=' tn_ice  : ')
       ENDIF