Changeset 7036
- Timestamp:
- 2016-10-18T10:49:34+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r6972 r7036 340 340 & Cd, Ch, Ce, zt_zu, zq_zu ) 341 341 342 ! Make ocean-atm. drag dependent on ice concentration (see Lupkes et al. 2012) (clem) 343 #if defined key_lim3 344 IF( ln_Cd_L12 ) THEN 345 346 Cd_oce(:,:) = Cd(:,:) ! record value of pure ocean-atm. drag 347 348 CALL Cdn10_Lupkes2012( Cd ) ! calculate new drag from Lupkes(2012) equations 349 350 ENDIF 351 #endif 352 342 Cd_oce(:,:) = Cd(:,:) ! record value of pure ocean-atm. drag (clem) 343 353 344 ! ... tau module, i and j component 354 345 DO jj = 1, jpj … … 483 474 #if defined key_lim3 484 475 IF( ln_Cd_L12 ) THEN 485 486 476 CALL Cdn10_Lupkes2012( Cd ) ! calculate new drag from Lupkes(2012) equations 487 488 477 ENDIF 489 478 #endif … … 964 953 !! *** ROUTINE Cdn10_Lupkes2012 *** 965 954 !! 966 !! ** Purpose : Recompute the ice-atm and ocean-atm dragsat 10m height to make967 !! themdependent on edges at leads, melt ponds and flows.955 !! ** Purpose : Recompute the ice-atm drag at 10m height to make 956 !! it dependent on edges at leads, melt ponds and flows. 968 957 !! After some approximations, this can be resumed to a dependency 969 958 !! on ice concentration. … … 971 960 !! ** Method : The parameterization is taken from Lupkes et al. (2012) eq.(50) 972 961 !! with the highest level of approximation: level4, eq.(59) 973 !! The dragcan be re-written as follows:962 !! The generic drag over a cell partly covered by ice can be re-written as follows: 974 963 !! 975 964 !! Cd = Cdw * (1-A) + Cdi * A + Ce * (1-A)**(nu+1/(10*beta)) * A**mu … … 999 988 zcoef = znu + 1._wp / ( 10._wp * zbeta ) 1000 989 1001 Cd(:,:) = Cd_oce(:,:) * ( 1._wp - at_i_b(:,:) ) + & ! pure ocean drag 1002 & Cd_ice * at_i_b(:,:) + & ! pure ice drag 1003 & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**zmu ! change due to sea-ice morphology 990 ! generic drag over a cell partly covered by ice 991 !!Cd(:,:) = Cd_oce(:,:) * ( 1._wp - at_i_b(:,:) ) + & ! pure ocean drag 992 !! & Cd_ice * at_i_b(:,:) + & ! pure ice drag 993 !! & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**zmu ! change due to sea-ice morphology 994 995 ! ice-atm drag 996 Cd(:,:) = Cd_ice + & ! pure ice drag 997 & zCe * ( 1._wp - at_i_b(:,:) )**zcoef * at_i_b(:,:)**(zmu-1._wp) ! change due to sea-ice morphology 1004 998 1005 999 END SUBROUTINE Cdn10_Lupkes2012
Note: See TracChangeset
for help on using the changeset viewer.