Changeset 991 for branches/dev_003_CPL/NEMO/LIM_SRC_3/limsbc.F90
- Timestamp:
- 2008-05-23T17:55:55+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_003_CPL/NEMO/LIM_SRC_3/limsbc.F90
r990 r991 242 242 REAL(wp) :: zpme ! freshwater exchanges at the ice/ocean interface 243 243 REAL(wp), DIMENSION(jpi,jpj) :: zfcm1 , zfcm2 ! solar/non solar heat fluxes 244 #if defined key_coupled245 244 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalb ! albedo of ice under overcast sky 246 245 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalbp ! albedo of ice under clear sky 247 #endif248 246 !!--------------------------------------------------------------------- 249 247 … … 422 420 !-----------------------------------------------! 423 421 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 439 432 IF(ln_ctl) THEN 440 433 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.