Changeset 1270 for trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90
- Timestamp:
- 2009-01-14T19:35:02+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r1242 r1270 31 31 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 32 32 USE prtctl ! Print control 33 #if defined key_lim234 USE par_ice_2 ! For number of ice categories35 #endif36 33 #if defined key_lim3 37 USE par_ice ! For number of ice categories38 34 USE ice_oce ! For ice surface temperature 39 35 #endif … … 342 338 & p_qla , p_dqns, p_dqla, & 343 339 & p_tpr , p_spr , & 344 & p_fr1 , p_fr2 , cd_grid )340 & p_fr1 , p_fr2 , cd_grid, pdim ) 345 341 !!--------------------------------------------------------------------- 346 342 !! *** ROUTINE blk_ice_core *** … … 355 351 !! caution : the net upward water flux has with mm/day unit 356 352 !!--------------------------------------------------------------------- 357 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj,jpl) :: pst ! ice surface temperature (>0, =rt0 over land) [Kelvin] 358 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pui ! ice surface velocity (i- and i- components [m/s] 359 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pvi ! at I-point (B-grid) or U & V-point (C-grid) 360 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj,jpl) :: palb ! ice albedo (clear sky) (alb_ice_cs) [%] 361 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_taui ! i- & j-components of surface ice stress [N/m2] 362 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tauj ! at I-point (B-grid) or U & V-point (C-grid) 363 REAL(wp), INTENT( out), DIMENSION(jpi,jpj,jpl) :: p_qns ! non solar heat flux over ice (T-point) [W/m2] 364 REAL(wp), INTENT( out), DIMENSION(jpi,jpj,jpl) :: p_qsr ! solar heat flux over ice (T-point) [W/m2] 365 REAL(wp), INTENT( out), DIMENSION(jpi,jpj,jpl) :: p_qla ! latent heat flux over ice (T-point) [W/m2] 366 REAL(wp), INTENT( out), DIMENSION(jpi,jpj,jpl) :: p_dqns ! non solar heat sensistivity (T-point) [W/m2] 367 REAL(wp), INTENT( out), DIMENSION(jpi,jpj,jpl) :: p_dqla ! latent heat sensistivity (T-point) [W/m2] 368 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tpr ! total precipitation (T-point) [Kg/m2/s] 369 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_spr ! solid precipitation (T-point) [Kg/m2/s] 370 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice (T-point) [%] 371 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice (T-point) [%] 372 CHARACTER(len=1), INTENT(in ) :: cd_grid ! ice grid ( C or B-grid) 353 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: pst ! ice surface temperature (>0, =rt0 over land) [Kelvin] 354 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pui ! ice surface velocity (i- and i- components [m/s] 355 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pvi ! at I-point (B-grid) or U & V-point (C-grid) 356 REAL(wp), INTENT(in ), DIMENSION(:,:,:) :: palb ! ice albedo (clear sky) (alb_ice_cs) [%] 357 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_taui ! i- & j-components of surface ice stress [N/m2] 358 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tauj ! at I-point (B-grid) or U & V-point (C-grid) 359 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qns ! non solar heat flux over ice (T-point) [W/m2] 360 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qsr ! solar heat flux over ice (T-point) [W/m2] 361 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_qla ! latent heat flux over ice (T-point) [W/m2] 362 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_dqns ! non solar heat sensistivity (T-point) [W/m2] 363 REAL(wp), INTENT( out), DIMENSION(:,:,:) :: p_dqla ! latent heat sensistivity (T-point) [W/m2] 364 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tpr ! total precipitation (T-point) [Kg/m2/s] 365 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_spr ! solid precipitation (T-point) [Kg/m2/s] 366 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice (T-point) [%] 367 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice (T-point) [%] 368 CHARACTER(len=1), INTENT(in ) :: cd_grid ! ice grid ( C or B-grid) 369 INTEGER, INTENT(in ) :: pdim ! number of ice categories 373 370 !! 374 371 INTEGER :: ji, jj, jl ! dummy loop indices … … 380 377 REAL(wp) :: zwndi_t , zwndj_t ! relative wind components at T-point 381 378 REAL(wp), DIMENSION(jpi,jpj) :: z_wnds_t ! wind speed ( = | U10m - U_ice | ) at T-point 382 REAL(wp), DIMENSION(jpi,jpj, jpl) :: z_qlw ! long wave heat flux over ice383 REAL(wp), DIMENSION(jpi,jpj, jpl) :: z_qsb ! sensible heat flux over ice384 REAL(wp), DIMENSION(jpi,jpj, jpl) :: z_dqlw ! sensible heat flux over ice385 REAL(wp), DIMENSION(jpi,jpj, jpl) :: z_dqsb ! sensible heat flux over ice379 REAL(wp), DIMENSION(jpi,jpj,pdim) :: z_qlw ! long wave heat flux over ice 380 REAL(wp), DIMENSION(jpi,jpj,pdim) :: z_qsb ! sensible heat flux over ice 381 REAL(wp), DIMENSION(jpi,jpj,pdim) :: z_dqlw ! sensible heat flux over ice 382 REAL(wp), DIMENSION(jpi,jpj,pdim) :: z_dqsb ! sensible heat flux over ice 386 383 !!--------------------------------------------------------------------- 387 384 388 ijpl = jpl! number of ice categories385 ijpl = pdim ! number of ice categories 389 386 390 387 ! local scalars ( place there for vector optimisation purposes)
Note: See TracChangeset
for help on using the changeset viewer.