Changeset 7036


Ignore:
Timestamp:
2016-10-18T10:49:34+02:00 (4 years ago)
Author:
clem
Message:

correct Lupkes param on ice-atm drags

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  
    340340         &               Cd, Ch, Ce, zt_zu, zq_zu ) 
    341341     
    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      
    353344      ! ... tau module, i and j component 
    354345      DO jj = 1, jpj 
     
    483474#if defined key_lim3 
    484475      IF( ln_Cd_L12 ) THEN 
    485  
    486476         CALL Cdn10_Lupkes2012( Cd ) ! calculate new drag from Lupkes(2012) equations 
    487  
    488477      ENDIF 
    489478#endif 
     
    964953      !!                      ***  ROUTINE  Cdn10_Lupkes2012  *** 
    965954      !! 
    966       !! ** Purpose :    Recompute the ice-atm and ocean-atm drags at 10m height to make 
    967       !!                 them dependent 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. 
    968957      !!                 After some approximations, this can be resumed to a dependency 
    969958      !!                 on ice concentration. 
     
    971960      !! ** Method :     The parameterization is taken from Lupkes et al. (2012) eq.(50) 
    972961      !!                 with the highest level of approximation: level4, eq.(59) 
    973       !!                 The drag can be re-written as follows: 
     962      !!                 The generic drag over a cell partly covered by ice can be re-written as follows: 
    974963      !! 
    975964      !!                 Cd = Cdw * (1-A) + Cdi * A + Ce * (1-A)**(nu+1/(10*beta)) * A**mu 
     
    999988      zcoef = znu + 1._wp / ( 10._wp * zbeta ) 
    1000989 
    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 
    1004998       
    1005999   END SUBROUTINE Cdn10_Lupkes2012 
Note: See TracChangeset for help on using the changeset viewer.