Changeset 1234
- Timestamp:
- 2008-12-29T17:37:02+01:00 (15 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 2 edited
-
LIM_SRC_2/par_ice_2.F90 (modified) (2 diffs)
-
OPA_SRC/SBC/sbcblk_core.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/par_ice_2.F90
r1156 r1234 1 1 MODULE par_ice_2 2 #if defined key_lim23 2 !!====================================================================== 4 3 !! *** MODULE par_ice_2 *** 5 4 !! Sea-Ice model : definition of the parameters 6 5 !!====================================================================== 6 #if defined key_lim2 7 7 !!---------------------------------------------------------------------- 8 !! LIM 2.0, UCL-LOCEAN-IPSL (2005) 9 !! $Id$ 10 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 8 !! 'key_lim2' LIM-2 sea-ice model 11 9 !!---------------------------------------------------------------------- 12 !! * Modules used13 10 USE par_oce 14 11 … … 16 13 PUBLIC ! allows par_oce and par_kind to be known in ice modules 17 14 18 INTEGER, PUBLIC, PARAMETER :: & !:19 jpkmax = 1 , &!: ???20 jpsmax =2 !: ???15 INTEGER, PUBLIC, PARAMETER :: jpl = 1 !: number of ice categories (only 1 in LIM-2) 16 INTEGER, PUBLIC, PARAMETER :: jpkmax = 1 !: ??? 17 INTEGER, PUBLIC, PARAMETER :: jpsmax = 2 !: ??? 21 18 22 INTEGER, PUBLIC, PARAMETER :: & !: 23 jplayers = 2 , & !: number of vertical ice layers 24 jplayersp1 = jplayers + 1 !: ??? 19 INTEGER, PUBLIC, PARAMETER :: jplayers = 2 !: number of vertical ice layers 20 INTEGER, PUBLIC, PARAMETER :: jplayersp1 = jplayers + 1 !: ??? 25 21 22 #else 23 !!---------------------------------------------------------------------- 24 !! Default option No LIM-2 sea-ice model 25 !!---------------------------------------------------------------------- 26 #endif 27 !!---------------------------------------------------------------------- 28 !! NEMO/LIM 2.0, UCL-LOCEAN-IPSL (2008) 29 !! $Id$ 30 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 26 31 !!====================================================================== 27 #endif28 32 END MODULE par_ice_2 -
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.
