- Timestamp:
- 2015-02-11T11:50:34+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO7_STARTHOUR/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r4624 r5075 114 114 !! - utau, vtau i- and j-component of the wind stress 115 115 !! - taum wind stress module at T-point 116 !! - wndm 10m wind module at T-point 116 !! - wndm 10m wind module at T-point over free ocean or leads in presence of sea-ice 117 117 !! - qns non-solar heat flux including latent heat of solid 118 118 !! precip. melting and emp heat content … … 204 204 !! - utau, vtau i- and j-component of the wind stress 205 205 !! - taum wind stress module at T-point 206 !! - wndm 10m wind module at T-point 206 !! - wndm 10m wind module at T-point over free ocean or leads in presence of sea-ice 207 207 !! - qns non-solar heat flux including latent heat of solid 208 208 !! precip. melting and emp heat content … … 257 257 END DO 258 258 END DO 259 utau(:,:) = utau(:,:) * umask(:,:,1) 260 vtau(:,:) = vtau(:,:) * vmask(:,:,1) 261 taum(:,:) = taum(:,:) * tmask(:,:,1) 259 262 CALL lbc_lnk( taum, 'T', 1. ) 260 263 … … 264 267 !CDIR COLLAPSE 265 268 wndm(:,:) = sf(jp_wndm)%fnow(:,:,1) 269 wndm(:,:) = wndm(:,:) * tmask(:,:,1) 266 270 267 271 !------------------------------------------------! … … 270 274 271 275 CALL blk_clio_qsr_oce( qsr ) 272 276 qsr(:,:) = qsr(:,:) * tmask(:,:,1) ! no shortwave radiation into the ocean beneath ice shelf 273 277 !------------------------! 274 278 ! Other ocean fluxes ! … … 376 380 & - zqla(:,:) * pst(:,:) * zcevap & ! remove evap. heat content at SST in Celcius 377 381 & + sf(jp_prec)%fnow(:,:,1) * sf(jp_tair)%fnow(:,:,1) * zcprec ! add precip. heat content at Tair in Celcius 382 qns(:,:) = qns(:,:) * tmask(:,:,1) 378 383 ! NB: if sea-ice model, the snow precip are computed and the associated heat is added to qns (see blk_ice_clio) 379 384 … … 398 403 399 404 400 SUBROUTINE blk_ice_clio( pst , palb_cs, palb_os ,&405 SUBROUTINE blk_ice_clio( pst , palb_cs, palb_os, palb, & 401 406 & p_taui, p_tauj, p_qns , p_qsr, & 402 407 & p_qla , p_dqns, p_dqla, & … … 427 432 !!---------------------------------------------------------------------- 428 433 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: pst ! ice surface temperature [Kelvin] 429 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb_cs ! ice albedo (clear sky) (alb_ice_cs) [%] 430 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb_os ! ice albedo (overcast sky) (alb_ice_os) [%] 434 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb_cs ! ice albedo (clear sky) (alb_ice_cs) [-] 435 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb_os ! ice albedo (overcast sky) (alb_ice_os) [-] 436 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: palb ! ice albedo (actual value) [-] 431 437 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_taui ! surface ice stress at I-point (i-component) [N/m2] 432 438 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tauj ! surface ice stress at I-point (j-component) [N/m2] … … 438 444 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tpr ! total precipitation (T-point) [Kg/m2/s] 439 445 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_spr ! solid precipitation (T-point) [Kg/m2/s] 440 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice [ %]441 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice [ %]446 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice [-] 447 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice [-] 442 448 CHARACTER(len=1), INTENT(in ) :: cd_grid ! type of sea-ice grid ("C" or "B" grid) 443 449 INTEGER, INTENT(in ) :: pdim ! number of ice categories … … 542 548 !-----------------------------------------------------------! 543 549 CALL blk_clio_qsr_ice( palb_cs, palb_os, p_qsr ) 550 551 DO jl = 1, ijpl 552 palb(:,:,jl) = ( palb_cs(:,:,jl) * ( 1.e0 - sf(jp_ccov)%fnow(:,:,1) ) & 553 & + palb_os(:,:,jl) * sf(jp_ccov)%fnow(:,:,1) ) 554 END DO 544 555 545 556 ! ! ========================== !
Note: See TracChangeset
for help on using the changeset viewer.