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 991 for branches/dev_003_CPL/NEMO/LIM_SRC_3/limsbc.F90 – NEMO

Ignore:
Timestamp:
2008-05-23T17:55:55+02:00 (16 years ago)
Author:
smasson
Message:

dev_003_CPL: preliminary draft (not working), see ticket #155

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/dev_003_CPL/NEMO/LIM_SRC_3/limsbc.F90

    r990 r991  
    242242      REAL(wp) ::   zpme             ! freshwater exchanges at the ice/ocean interface 
    243243      REAL(wp), DIMENSION(jpi,jpj) ::   zfcm1 , zfcm2    ! solar/non solar heat fluxes 
    244 #if defined key_coupled     
    245244      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb     ! albedo of ice under overcast sky 
    246245      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalbp    ! albedo of ice under clear sky 
    247 #endif 
    248246      !!--------------------------------------------------------------------- 
    249247 
     
    422420      !-----------------------------------------------! 
    423421 
    424       freeze(:,:)   = at_i(:,:)             ! Sea ice cover             
    425       tn_ice(:,:,:) = t_su(:,:,:)           ! Ice surface temperature                       
    426  
    427 #if defined key_coupled             
    428       !------------------------------------------------! 
    429       !    Computation of snow/ice and ocean albedo    ! 
    430       !------------------------------------------------! 
    431       zalb  (:,:,:) = 0.e0 
    432       zalbp (:,:,:) = 0.e0 
    433  
    434       CALL albedo_ice( t_su, ht_i, ht_s, zalbp, zalb ) 
    435  
    436       alb_ice(:,:,:) =  0.5 * zalbp(:,:,:) + 0.5 * zalb (:,:,:)   ! Ice albedo (mean clear and overcast skys) 
    437 #endif 
    438  
     422      freeze(:,:)   = at_i(:,:)             ! Sea ice cover   
     423 
     424      IF ( lk_cpl ) THEN                 
     425         ! Ice surface temperature 
     426         tn_ice(:,:,:) = t_su(:,:,:)            
     427         ! Computation of snow/ice and ocean albedo 
     428         CALL albedo_ice( t_su, ht_i, ht_s, zalbp, zalb ) 
     429         alb_ice(:,:,:) =  0.5 * ( zalbp(:,:,:) + zalb (:,:,:) )   ! Ice albedo (mean clear and overcast skys) 
     430      ENDIF 
     431          
    439432      IF(ln_ctl) THEN 
    440433         CALL prt_ctl( tab2d_1=qsr   , clinfo1=' lim_sbc: qsr    : ', tab2d_2=qns , clinfo2=' qns     : ' ) 
Note: See TracChangeset for help on using the changeset viewer.