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