Changeset 2291 for branches/nemo_v3_3_beta/NEMOGCM/NEMO
- Timestamp:
- 2010-10-19T09:50:32+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r2287 r2291 62 62 REAL(wp), PARAMETER :: Stef = 5.67e-8 ! Stefan Boltzmann constant 63 63 REAL(wp), PARAMETER :: Cice = 1.63e-3 ! transfer coefficient over ice 64 REAL(wp), PARAMETER :: albo = 0.066 ! ocean albedo assumed to be contant 64 65 65 66 ! !!* Namelist namsbc_core : CORE bulk parameters … … 262 263 ! ----------------------------------------------------------------------------- ! 263 264 264 IF( ln_dm2dc ) THEN ; qsr(:,:) = sbc_dcy( sf(jp_qsr)%fnow(:,:,1) ) ! modify now Qsr to include the diurnal cycle 265 ELSE ; qsr(:,:) = sf(jp_qsr)%fnow(:,:,1) 265 ! ocean albedo assumed to be constant + modify now Qsr to include the diurnal cycle ! Short Wave 266 zztmp = 1. - albo 267 IF( ln_dm2dc ) THEN ; qsr(:,:) = zztmp * sbc_dcy( sf(jp_qsr)%fnow(:,:,1) ) * tmask(:,:,1) 268 ELSE ; qsr(:,:) = zztmp * sf(jp_qsr)%fnow(:,:,1) * tmask(:,:,1) 266 269 ENDIF 267 ! ocean albedo assumed to be 0.066268 qsr (:,:) = ( 1. - 0.066 ) * qsr(:,:) * tmask(:,:,1) ! Short Wave269 270 !CDIR COLLAPSE 270 271 zqlw(:,:) = ( sf(jp_qlw)%fnow(:,:,1) - Stef * zst(:,:)*zst(:,:)*zst(:,:)*zst(:,:) ) * tmask(:,:,1) ! Long Wave … … 417 418 REAL(wp) :: zst2, zst3 418 419 REAL(wp) :: zcoef_wnorm, zcoef_wnorm2, zcoef_dqlw, zcoef_dqla, zcoef_dqsb 420 REAL(wp) :: zztmp ! temporary variable 419 421 REAL(wp) :: zcoef_frca ! fractional cloud amount 420 422 REAL(wp) :: zwnorm_f, zwndi_f , zwndj_f ! relative wind module and components at F-point … … 504 506 END SELECT 505 507 508 zztmp = 1. / ( 1. - albo ) 506 509 ! ! ========================== ! 507 510 DO jl = 1, ijpl ! Loop over ice categories ! … … 518 521 zst3 = pst(ji,jj,jl) * zst2 519 522 ! Short Wave (sw) 520 p_qsr(ji,jj,jl) = ( 1. - palb(ji,jj,jl) ) * sf(jp_qsr)%fnow(ji,jj,1) * tmask(ji,jj,1)523 p_qsr(ji,jj,jl) = zztmp * ( 1. - palb(ji,jj,jl) ) * qsr(ji,jj) 521 524 ! Long Wave (lw) 522 525 z_qlw(ji,jj,jl) = 0.95 * ( sf(jp_qlw)%fnow(ji,jj,1) - Stef * pst(ji,jj,jl) * zst3 ) * tmask(ji,jj,1)
Note: See TracChangeset
for help on using the changeset viewer.