Changeset 7507 for branches/2016/dev_merge_2016
- Timestamp:
- 2016-12-19T12:06:58+01:00 (8 years ago)
- Location:
- branches/2016/dev_merge_2016/NEMOGCM/NEMO
- Files:
-
- 2 deleted
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
r7421 r7507 218 218 zh_i_ini(ji,jj,jl) = hi_mean(jl) 219 219 END DO 220 220 ! 221 221 !--- Concentrations 222 222 za_i_ini(ji,jj,jl0) = zat_i_ini(ji,jj) / SQRT(REAL(jpl)) … … 227 227 ENDIF 228 228 END DO 229 229 ! 230 230 ! Concentration in the last (i_fill) category 231 231 za_i_ini(ji,jj,i_fill) = zat_i_ini(ji,jj) - SUM( za_i_ini(ji,jj,1:i_fill-1) ) … … 274 274 END DO ! end iteration on categories 275 275 ! !============================ 276 276 ! 277 277 IF( lwp .AND. SUM(itest) /= 4 ) THEN 278 278 WRITE(numout,*) -
branches/2016/dev_merge_2016/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r7421 r7507 34 34 USE traqsr ! add penetration of solar flux in the calculation of heat budget 35 35 USE domvvl ! Variable volume 36 USE limctl ! 37 USE limcons ! 36 USE limctl ! 37 USE limcons ! 38 38 USE bdy_oce , ONLY: ln_bdy 39 39 ! -
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.