Changeset 1270
- Timestamp:
- 2009-01-14T19:35:02+01:00 (15 years ago)
- Location:
- trunk/NEMO/OPA_SRC/SBC
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r1269 r1270 30 30 USE prtctl ! Print control 31 31 #if defined key_lim3 32 USE par_ice33 32 USE ice 34 33 USE ice_oce ! For ice surface temperature 35 34 #elif defined key_lim2 36 USE par_ice_237 35 USE ice_2 38 36 #endif … … 414 412 & p_qla , p_dqns, p_dqla, & 415 413 & p_tpr , p_spr , & 416 & p_fr1 , p_fr2 , cd_grid )414 & p_fr1 , p_fr2 , cd_grid, pdim ) 417 415 !!--------------------------------------------------------------------------- 418 416 !! *** ROUTINE blk_ice_clio *** … … 452 450 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice [%] 453 451 CHARACTER(len=1), INTENT(in ) :: cd_grid ! type of sea-ice grid ("C" or "B" grid) 452 INTEGER, INTENT(in ) :: pdim ! number of ice categories 454 453 !! 455 454 INTEGER :: ji, jj, jl ! dummy loop indices … … 467 466 REAL(wp), DIMENSION(jpi,jpj) :: zevsqr ! vapour pressure square-root 468 467 REAL(wp), DIMENSION(jpi,jpj) :: zrhoa ! air density 469 REAL(wp), DIMENSION(jpi,jpj, SIZE(pst,3)) :: z_qlw, z_qsb468 REAL(wp), DIMENSION(jpi,jpj,pdim) :: z_qlw, z_qsb 470 469 !!--------------------------------------------------------------------- 471 470 472 ijpl = SIZE( pst, 3 )! number of ice categories471 ijpl = pdim ! number of ice categories 473 472 zpatm = 101000. ! atmospheric pressure (assumed constant here) 474 473 -
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) -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r1240 r1270 145 145 & qla_ice , dqns_ice , dqla_ice , & 146 146 & tprecip , sprecip , & 147 & fr1_i0 , fr2_i0 , cl_grid )147 & fr1_i0 , fr2_i0 , cl_grid, jpl ) 148 148 ! 149 149 CASE( 4 ) ! CORE bulk formulation … … 152 152 & qla_ice , dqns_ice , dqla_ice , & 153 153 & tprecip , sprecip , & 154 & fr1_i0 , fr2_i0 , cl_grid )154 & fr1_i0 , fr2_i0 , cl_grid, jpl ) 155 155 END SELECT 156 156 -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r1239 r1270 21 21 USE lib_mpp 22 22 USE ice_2 23 USE par_ice_2 23 24 USE iceini_2 24 25 USE ice_oce ! ice variables … … 154 155 & zqla_ice , zdqns_ice , zdqla_ice , & 155 156 & tprecip , sprecip , & 156 & fr1_i0 , fr2_i0 , cl_grid )157 & fr1_i0 , fr2_i0 , cl_grid, jpl ) 157 158 158 159 CASE( 4 ) ! CORE bulk formulation … … 161 162 & zqla_ice , zdqns_ice , zdqla_ice , & 162 163 & tprecip , sprecip , & 163 & fr1_i0 , fr2_i0 , cl_grid )164 & fr1_i0 , fr2_i0 , cl_grid, jpl ) 164 165 CASE( 5 ) ! Coupled formulation : atmosphere-ice stress only (fluxes provided after ice dynamics) 165 166 CALL sbc_cpl_ice_tau( utaui_ice , vtaui_ice )
Note: See TracChangeset
for help on using the changeset viewer.