Changeset 4923 for branches/2014/dev_CNRS_2014
- Timestamp:
- 2014-11-28T16:07:52+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS_2014/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r4901 r4923 351 351 !! q_air and t_air are not given at 10m (wind reference height) 352 352 ! Values of temp. and hum. adjusted to height of wind during bulk algorithm iteration must be used!!! 353 zevap(:,:) = rn_efac*MAX( 0. e0, rhoa*Ce(:,:)*( zqsatw(:,:) - zq_zu(:,:) )*wndm(:,:) ) ! Evaporation354 zqsb (:,:) = cpa*rhoa*Ch(:,:)*( zst (:,:) - zt_zu(:,:) )*wndm(:,:) ! Sensible Heat353 zevap(:,:) = rn_efac*MAX( 0._wp, rhoa*Ce(:,:)*( zqsatw(:,:) - zq_zu(:,:) )*wndm(:,:) ) ! Evaporation 354 zqsb (:,:) = cpa*rhoa*Ch(:,:)*( zst (:,:) - zt_zu(:,:) )*wndm(:,:) ! Sensible Heat 355 355 ENDIF 356 356 zqla (:,:) = Lv * zevap(:,:) ! Latent Heat … … 541 541 p_qsr(ji,jj,jl) = zztmp * ( 1. - palb(ji,jj,jl) ) * qsr(ji,jj) 542 542 ! Long Wave (lw) 543 ! iovino 544 IF( ff(ji,jj) .GT. 0._wp ) THEN 545 z_qlw(ji,jj,jl) = ( 0.95 * sf(jp_qlw)%fnow(ji,jj,1) - Stef * pst(ji,jj,jl) * zst3 ) * tmask(ji,jj,1) 546 ELSE 547 z_qlw(ji,jj,jl) = 0.95 * ( sf(jp_qlw)%fnow(ji,jj,1) - Stef * pst(ji,jj,jl) * zst3 ) * tmask(ji,jj,1) 548 ENDIF 543 z_qlw(ji,jj,jl) = 0.95 * ( sf(jp_qlw)%fnow(ji,jj,1) - Stef * pst(ji,jj,jl) * zst3 ) * tmask(ji,jj,1) 549 544 ! lw sensitivity 550 545 z_dqlw(ji,jj,jl) = zcoef_dqlw * zst3 … … 561 556 & * ( 11637800. * EXP( -5897.8 / pst(ji,jj,jl) ) / rhoa - sf(jp_humi)%fnow(ji,jj,1) ) ) 562 557 ! Latent heat sensitivity for ice (Dqla/Dt) 563 p_dqla(ji,jj,jl) = rn_efac * zcoef_dqla * z_wnds_t(ji,jj) / ( zst2 ) * EXP( -5897.8 / pst(ji,jj,jl) ) 558 IF( p_qla(ji,jj,jl) > 0._wp ) THEN 559 p_dqla(ji,jj,jl) = rn_efac * zcoef_dqla * z_wnds_t(ji,jj) / ( zst2 ) * EXP( -5897.8 / pst(ji,jj,jl) ) 560 ELSE 561 p_dqla(ji,jj,jl) = 0._wp 562 ENDIF 563 564 564 ! Sensible heat sensitivity (Dqsb_ice/Dtn_ice) 565 565 z_dqsb(ji,jj,jl) = zcoef_dqsb * z_wnds_t(ji,jj) … … 582 582 ! thin surface layer and penetrates inside the ice cover 583 583 ! ( Maykut and Untersteiner, 1971 ; Ebert and Curry, 1993 ) 584 584 ! 585 585 p_fr1(:,:) = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) 586 586 p_fr2(:,:) = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) 587 587 ! 588 588 p_tpr(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac ! total precipitation [kg/m2/s] 589 589 p_spr(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac ! solid precipitation [kg/m2/s] … … 602 602 ENDIF 603 603 604 CALL wrk_dealloc( jpi,jpj, z_wnds_t )605 CALL wrk_dealloc( jpi,jpj, pdim, z_qlw, z_qsb, z_dqlw, z_dqsb )604 CALL wrk_dealloc( jpi,jpj, z_wnds_t ) 605 CALL wrk_dealloc( jpi,jpj, pdim, z_qlw, z_qsb, z_dqlw, z_dqsb ) 606 606 ! 607 607 IF( nn_timing == 1 ) CALL timing_stop('blk_ice_core')
Note: See TracChangeset
for help on using the changeset viewer.