Changeset 2412
- Timestamp:
- 2010-11-19T18:36:22+01:00 (14 years ago)
- Location:
- branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r2409 r2412 235 235 !------------------------------------! 236 236 !CDIR COLLAPSE 237 DO jj = 1 , jpj 238 DO ji = 1, jpi 239 utau(ji,jj) = sf(jp_utau)%fnow(ji,jj,1) 240 vtau(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1) 241 END DO 242 END DO 237 utau(:,:) = sf(jp_utau)%fnow(:,:,1) 238 !CDIR COLLAPSE 239 vtau(:,:) = sf(jp_vtau)%fnow(:,:,1) 243 240 244 241 !------------------------------------! … … 260 257 !------------------------------------! 261 258 !CDIR COLLAPSE 262 DO jj = 1 , jpj 263 DO ji = 1, jpi 264 wndm(ji,jj) = sf(jp_wndm)%fnow(ji,jj,1) 265 END DO 266 END DO 259 wndm(:,:) = sf(jp_wndm)%fnow(:,:,1) 267 260 268 261 !------------------------------------------------! … … 368 361 369 362 !CDIR COLLAPSE 370 !CDIR NOVERRCHK 371 DO jj = 1, jpj 372 !CDIR NOVERRCHK 373 DO ji = 1, jpi 374 qns (ji,jj) = zqlw(ji,jj) - zqsb(ji,jj) - zqla(ji,jj) ! Downward Non Solar flux 375 emp (ji,jj) = zqla(ji,jj) / cevap - sf(jp_prec)%fnow(ji,jj,1) / rday * tmask(ji,jj,1) 376 END DO 377 END DO 363 emp (:,:) = zqla(:,:) / cevap - sf(jp_prec)%fnow(:,:,1) / rday * tmask(:,:,1) 364 qns (:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) ! Downward Non Solar flux 378 365 emps(:,:) = emp(:,:) 379 366 ! … … 464 451 SELECT CASE( cd_grid ) 465 452 CASE( 'C' ) ! C-grid ice dynamics 466 zcoef = cai / cao ! Change from air-sea stress to air-ice stress 467 !CDIR COLLAPSE 468 DO jj = 1 , jpj 469 DO ji = 1, jpi 470 p_taui(ji,jj) = zcoef * utau(ji,jj) 471 p_tauj(ji,jj) = zcoef * vtau(ji,jj) 472 END DO 473 END DO 453 zcoef = cai / cao ! Change from air-sea stress to air-ice stress 454 p_taui(:,:) = zcoef * utau(:,:) 455 p_tauj(:,:) = zcoef * vtau(:,:) 474 456 CASE( 'I' ) ! I-grid ice dynamics: I-point (i.e. F-point lower-left corner) 475 zcoef = 0.5_wp * cai / cao ! Change from air-sea stress to air-ice stress!CDIR COLLAPSE457 zcoef = 0.5_wp * cai / cao ! Change from air-sea stress to air-ice stress 476 458 DO jj = 2, jpj ! stress from ocean U- and V-points to ice U,V point 477 459 DO ji = 2, jpi ! I-grid : no vector opt. -
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r2370 r2412 359 359 qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:) ! Downward Non Solar flux 360 360 !CDIR COLLAPSE 361 emp 361 emp(:,:) = zevap(:,:) - sf(jp_prec)%fnow(:,:,1) * rn_pfac * tmask(:,:,1) 362 362 !CDIR COLLAPSE 363 363 emps(:,:) = emp(:,:) … … 396 396 !! caution : the net upward water flux has with mm/day unit 397 397 !!--------------------------------------------------------------------- 398 REAL(wp), INTENT(in ), DIMENSION(:,:,:):: pst ! ice surface temperature (>0, =rt0 over land) [Kelvin]399 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pui ! ice surface velocity (i- and i- components [m/s]400 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj) :: pvi ! at I-point (B-grid) or U & V-point (C-grid)401 REAL(wp), INTENT(in ), DIMENSION(:,:,:):: palb ! ice albedo (clear sky) (alb_ice_cs) [%]402 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_taui ! i- & j-components of surface ice stress [N/m2]403 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tauj ! at I-point (B-grid) or U & V-point (C-grid)404 REAL(wp), INTENT( out), DIMENSION(:,:,:):: p_qns ! non solar heat flux over ice (T-point) [W/m2]405 REAL(wp), INTENT( out), DIMENSION(:,:,:):: p_qsr ! solar heat flux over ice (T-point) [W/m2]406 REAL(wp), INTENT( out), DIMENSION(:,:,:):: p_qla ! latent heat flux over ice (T-point) [W/m2]407 REAL(wp), INTENT( out), DIMENSION(:,:,:):: p_dqns ! non solar heat sensistivity (T-point) [W/m2]408 REAL(wp), INTENT( out), DIMENSION(:,:,:):: p_dqla ! latent heat sensistivity (T-point) [W/m2]409 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_tpr ! total precipitation (T-point) [Kg/m2/s]410 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_spr ! solid precipitation (T-point) [Kg/m2/s]411 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr1 ! 1sr fraction of qsr penetration in ice (T-point) [%]412 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: p_fr2 ! 2nd fraction of qsr penetration in ice (T-point) [%]413 CHARACTER(len=1) , INTENT(in ):: cd_grid ! ice grid ( C or B-grid)414 INTEGER , INTENT(in ):: pdim ! number of ice categories398 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: pst ! ice surface temperature (>0, =rt0 over land) [Kelvin] 399 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: pui ! ice surface velocity (i- and i- components [m/s] 400 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: pvi ! at I-point (B-grid) or U & V-point (C-grid) 401 REAL(wp), DIMENSION(:,:,:) , INTENT(in ) :: palb ! ice albedo (clear sky) (alb_ice_cs) [%] 402 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_taui ! i- & j-components of surface ice stress [N/m2] 403 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_tauj ! at I-point (B-grid) or U & V-point (C-grid) 404 REAL(wp), DIMENSION(:,:,:) , INTENT( out) :: p_qns ! non solar heat flux over ice (T-point) [W/m2] 405 REAL(wp), DIMENSION(:,:,:) , INTENT( out) :: p_qsr ! solar heat flux over ice (T-point) [W/m2] 406 REAL(wp), DIMENSION(:,:,:) , INTENT( out) :: p_qla ! latent heat flux over ice (T-point) [W/m2] 407 REAL(wp), DIMENSION(:,:,:) , INTENT( out) :: p_dqns ! non solar heat sensistivity (T-point) [W/m2] 408 REAL(wp), DIMENSION(:,:,:) , INTENT( out) :: p_dqla ! latent heat sensistivity (T-point) [W/m2] 409 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_tpr ! total precipitation (T-point) [Kg/m2/s] 410 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_spr ! solid precipitation (T-point) [Kg/m2/s] 411 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_fr1 ! 1sr fraction of qsr penetration in ice (T-point) [%] 412 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: p_fr2 ! 2nd fraction of qsr penetration in ice (T-point) [%] 413 CHARACTER(len=1) , INTENT(in ) :: cd_grid ! ice grid ( C or B-grid) 414 INTEGER , INTENT(in ) :: pdim ! number of ice categories 415 415 !! 416 416 INTEGER :: ji, jj, jl ! dummy loop indices … … 432 432 433 433 ! local scalars ( place there for vector optimisation purposes) 434 zcoef_wnorm = rhoa * Cice434 zcoef_wnorm = rhoa * Cice 435 435 zcoef_wnorm2 = rhoa * Cice * 0.5 436 zcoef_dqlw = 4.0 * 0.95 * Stef437 zcoef_dqla = -Ls * Cice * 11637800. * (-5897.8)438 zcoef_dqsb = rhoa * cpa * Cice439 zcoef_frca = 1.0 - 0.3436 zcoef_dqlw = 4.0 * 0.95 * Stef 437 zcoef_dqla = -Ls * Cice * 11637800. * (-5897.8) 438 zcoef_dqsb = rhoa * cpa * Cice 439 zcoef_frca = 1.0 - 0.3 440 440 441 441 !!gm brutal.... … … 451 451 CASE( 'I' ) ! B-grid ice dynamics : I-point (i.e. F-point with sea-ice indexation) 452 452 ! and scalar wind at T-point ( = | U10m - U_ice | ) (masked) 453 #if defined key_vectopt_loop454 !CDIR COLLAPSE455 #endif456 453 !CDIR NOVERRCHK 457 454 DO jj = 2, jpjm1 458 DO ji = 2, jpim1 ! B grid : novector opt455 DO ji = 2, jpim1 ! B grid : NO vector opt 459 456 ! ... scalar wind at I-point (fld being at T-point) 460 457 zwndi_f = 0.25 * ( sf(jp_wndi)%fnow(ji-1,jj ,1) + sf(jp_wndi)%fnow(ji ,jj ,1) &
Note: See TracChangeset
for help on using the changeset viewer.