- Timestamp:
- 2016-12-19T12:06:58+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90
r7504 r7507 641 641 zwndj_t = sf(jp_wndj)%fnow(ji,jj,1) - rn_vfac * 0.25 * ( v_ice(ji,jj+1) + v_ice(ji+1,jj+1) & 642 642 & + v_ice(ji,jj ) + v_ice(ji+1,jj ) ) 643 wndm_ice(ji,jj) 643 wndm_ice(ji,jj) = SQRT( zwndi_t * zwndi_t + zwndj_t * zwndj_t ) * tmask(ji,jj,1) 644 644 END DO 645 645 END DO … … 760 760 ! Latent heat sensitivity for ice (Dqla/Dt) 761 761 IF( qla_ice(ji,jj,jl) > 0._wp ) THEN 762 dqla_ice(ji,jj,jl) = rn_efac * zcoef_dqla * wndm_ice(ji,jj) / ( zst2 ) * EXP( -5897.8 / ptsu(ji,jj,jl) )762 dqla_ice(ji,jj,jl) = rn_efac * zcoef_dqla * Cd(ji,jj) * wndm_ice(ji,jj) / ( zst2 ) * EXP( -5897.8 / ptsu(ji,jj,jl) ) 763 763 ELSE 764 764 dqla_ice(ji,jj,jl) = 0._wp … … 818 818 819 819 ! --- heat content of evap over ice in W/m2 (to be used in 1D-thermo) --- 820 ! !821 820 DO jl = 1, jpl 822 qevap_ice(:,:,jl) = 0._wp ! should be -evap_ice(:,:,jl)*( ( Tice - rt0) * cpic * tmask(:,:,1) ) 823 ! But we do not have Tice => consider it at 824 ! 0??C => evap=0 821 qevap_ice(:,:,jl) = 0._wp ! should be -evap_ice(:,:,jl)*( ( Tice - rt0 ) * cpic * tmask(:,:,1) ) 822 ! But we do not have Tice => consider it at 0degC => evap=0 825 823 END DO 826 824 … … 998 996 !! and going down to Cdi(say 1.4e-3) for A=1 999 997 !! 1000 !! It is theoretically applicable to all ice conditions !(not only MIZ)998 !! It is theoretically applicable to all ice conditions (not only MIZ) 1001 999 !! => see Lupkes et al (2013) 1002 1000 !! … … 1015 1013 1016 1014 ! generic drag over a cell partly covered by ice 1017 !!Cd(:,:) = Cd_oce(:,:) * ( 1._wp - at_i_b(:,:) ) + & !! pure ocean drag1018 !! & Cd_ice * at_i_b(:,:) + & !! pure ice drag1019 !! & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * !at_i_b(:,:)**zmu ! change due to sea-ice morphology1015 !!Cd(:,:) = Cd_oce(:,:) * ( 1._wp - at_i_b(:,:) ) + & ! pure ocean drag 1016 !! & Cd_ice * at_i_b(:,:) + & ! pure ice drag 1017 !! & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**zmu ! change due to sea-ice morphology 1020 1018 1021 1019 ! ice-atm drag 1022 Cd(:,:) = Cd_ice + & ! pure ice drag1023 & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**(zmu-1._wp) 1020 Cd(:,:) = Cd_ice + & ! pure ice drag 1021 & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**(zmu-1._wp) ! change due to sea-ice morphology 1024 1022 1025 1023 END SUBROUTINE Cdn10_Lupkes2012
Note: See TracChangeset
for help on using the changeset viewer.