New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1474 – NEMO

Changeset 1474


Ignore:
Timestamp:
2009-06-12T16:59:31+02:00 (15 years ago)
Author:
smasson
Message:

continue changeset:1463, see ticket:444

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_2/limsbc_2.F90

    r1470 r1474  
    9090      REAL(wp), DIMENSION(jpi,jpj,1) ::   zalb     ! albedo of ice under overcast sky 
    9191      REAL(wp), DIMENSION(jpi,jpj,1) ::   zalbp    ! albedo of ice under clear sky 
    92       REAL(wp), DIMENSION(jpi,jpj,1) ::   zsist    ! surface ice temperature (K) 
    93       REAL(wp), DIMENSION(jpi,jpj,1) ::   zhicif   ! ice thickness 
    94       REAL(wp), DIMENSION(jpi,jpj,1) ::   zhsnif   ! snow thickness 
    9592      REAL(wp) ::   zsang, zmod, zfm 
    9693      REAL(wp), DIMENSION(jpi,jpj) ::   ztio_u, ztio_v   ! ocean stress below sea-ice 
     
    281278      !-----------------------------------------------! 
    282279 
    283       fr_i  (:,:) = 1.0 - frld(:,:)       ! sea-ice fraction 
     280      fr_i(:,:) = 1.0 - frld(:,:)       ! sea-ice fraction 
    284281 
    285282      IF ( lk_cpl ) THEN            
     
    287284         tn_ice(:,:,1) = sist(:,:)          ! sea-ice surface temperature        
    288285         ! Computation of snow/ice and ocean albedo 
    289          ! INTERFACE 3D versus 2D 
    290          zsist (:,:,1) = sist (:,:) 
    291          zhicif(:,:,1) = hicif(:,:)   ;   zhsnif(:,:,1) = hsnif(:,:) 
    292          CALL albedo_ice( zsist, zhicif, zhsnif, zalbp, zalb ) 
     286         CALL albedo_ice( tn_ice(:,:,1), reshape( hicif, (/jpi,jpj,1/) ), reshape( hsnif, (/jpi,jpj,1/) ), zalbp, zalb ) 
    293287         alb_ice(:,:,1) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys) 
    294288      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.