- Timestamp:
- 2017-08-08T15:58:05+02:00 (7 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90
r8321 r8422 41 41 #if defined key_lim3 42 42 USE ice , ONLY : u_ice, v_ice, jpl, a_i_b, at_i_b 43 USE limthd_dh ! for CALL lim_thd_snwblow43 USE icethd_dh ! for CALL ice_thd_snwblow 44 44 #endif 45 45 USE sbcblk_algo_ncar ! => turb_ncar : NCAR - CORE (Large & Yeager, 2009) … … 789 789 ! --- evaporation minus precipitation --- ! 790 790 zsnw(:,:) = 0._wp 791 CALL lim_thd_snwblow( (1.-at_i_b(:,:)), zsnw ) ! snow distribution over ice after wind blowing791 CALL ice_thd_snwblow( (1.-at_i_b(:,:)), zsnw ) ! snow distribution over ice after wind blowing 792 792 emp_oce(:,:) = ( 1._wp - at_i_b(:,:) ) * zevap(:,:) - ( tprecip(:,:) - sprecip(:,:) ) - sprecip(:,:) * (1._wp - zsnw ) 793 793 emp_ice(:,:) = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8321 r8422 40 40 #endif 41 41 #if defined key_lim3 42 USE limthd_dh ! for CALL lim_thd_snwblow42 USE icethd_dh ! for CALL ice_thd_snwblow 43 43 #endif 44 44 ! … … 1621 1621 #if defined key_lim3 1622 1622 ! zsnw = snow fraction over ice after wind blowing (=picefr if no blowing) 1623 zsnw(:,:) = 0._wp ; CALL lim_thd_snwblow( ziceld, zsnw )1623 zsnw(:,:) = 0._wp ; CALL ice_thd_snwblow( ziceld, zsnw ) 1624 1624 1625 1625 ! --- evaporation minus precipitation corrected (because of wind blowing on snow) --- !
Note: See TracChangeset
for help on using the changeset viewer.