Changeset 7077


Ignore:
Timestamp:
2016-10-24T17:07:43+02:00 (4 years ago)
Author:
clem
Message:

enable analytic forcings with LIM3. It requires a change in namelist_ref.

Location:
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6861 r7077  
    278278&namsbc_ana    !   analytical surface boundary condition 
    279279!----------------------------------------------------------------------- 
     280   ! --- oce variables --- ! 
    280281   nn_tau000   =   0       !  gently increase the stress over the first ntau_rst time-steps 
    281282   rn_utau0    =   0.5     !  uniform value for the i-stress 
     
    284285   rn_qsr0     =   0.e0    !  uniform value for the solar radiation 
    285286   rn_emp0     =   0.e0    !  uniform value for the freswater budget (E-P) 
     287   ! --- ice variables --- ! 
     288   rn_iutau0   =   0.e0    !  uniform value for the i-stress over ice 
     289   rn_ivtau0   =   0.e0    !  uniform value for the j-stress over ice 
     290   rn_iqns0    =   0.e0    !  uniform value for the total heat flux over ice 
     291   rn_iqsr0    =   0.e0    !  uniform value for the solar radiation over ice 
     292   rn_sprec0   =   0.e0    !  uniform value for snow precip 
     293   rn_ievap0   =   0.e0    !  uniform value for sublimation 
    286294/ 
    287295!----------------------------------------------------------------------- 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90

    r4624 r7077  
    1515   USE dom_oce         ! ocean space and time domain 
    1616   USE sbc_oce         ! Surface boundary condition: ocean fields 
     17   USE sbc_ice         ! Surface boundary condition: ice   fields 
    1718   USE phycst          ! physical constants 
    1819   USE in_out_manager  ! I/O manager 
     
    2021   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    2122   USE lib_fortran 
    22  
     23   USE wrk_nemo 
     24#if defined key_lim3 
     25   USE ice, ONLY       : pfrld, a_i_b 
     26   USE limthd_dh       ! for CALL lim_thd_snwblow 
     27#endif 
     28    
    2329   IMPLICIT NONE 
    2430   PRIVATE 
    2531 
    26    PUBLIC   sbc_ana    ! routine called in sbcmod module 
    27    PUBLIC   sbc_gyre   ! routine called in sbcmod module 
     32   PUBLIC   sbc_ana         ! routine called in sbcmod module 
     33   PUBLIC   sbc_gyre        ! routine called in sbcmod module 
     34#if defined key_lim3 
     35   PUBLIC   ana_ice_tau     ! routine called in sbc_ice_lim module 
     36   PUBLIC   ana_ice_flx     ! routine called in sbc_ice_lim module 
     37#endif 
    2838 
    2939   !                       !!* Namelist namsbc_ana * 
    30    INTEGER  ::   nn_tau000  ! nb of time-step during which the surface stress 
    31    !                        ! increase from 0 to its nominal value  
    32    REAL(wp) ::   rn_utau0   ! constant wind stress value in i-direction 
    33    REAL(wp) ::   rn_vtau0   ! constant wind stress value in j-direction 
    34    REAL(wp) ::   rn_qns0    ! non solar heat flux 
    35    REAL(wp) ::   rn_qsr0    !     solar heat flux 
    36    REAL(wp) ::   rn_emp0    ! net freshwater flux 
     40   ! --- oce variables --- ! 
     41   INTEGER  ::   nn_tau000 ! nb of time-step during which the surface stress 
     42   !                       ! increase from 0 to its nominal value  
     43   REAL(wp) ::   rn_utau0  ! constant wind stress value in i-direction 
     44   REAL(wp) ::   rn_vtau0  ! constant wind stress value in j-direction 
     45   REAL(wp) ::   rn_qns0   ! non solar heat flux 
     46   REAL(wp) ::   rn_qsr0   !     solar heat flux 
     47   REAL(wp) ::   rn_emp0   ! net freshwater flux 
     48   ! --- ice variables --- ! 
     49   REAL(wp) ::   rn_iutau0 ! constant wind stress value in i-direction over ice 
     50   REAL(wp) ::   rn_ivtau0 ! constant wind stress value in j-direction over ice 
     51   REAL(wp) ::   rn_iqns0  ! non solar heat flux over ice 
     52   REAL(wp) ::   rn_iqsr0  !     solar heat flux over ice 
     53   REAL(wp) ::   rn_sprec0 ! snow precip 
     54   REAL(wp) ::   rn_ievap0 ! sublimation 
    3755    
    3856   !! * Substitutions 
     
    6987      REAL(wp) ::   zcoef, zty, zmod      !   -      - 
    7088      !! 
    71       NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0 
     89      NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0,  & 
     90         &                 rn_iutau0, rn_ivtau0, rn_iqsr0, rn_iqns0, rn_sprec0, rn_ievap0 
    7291      !!--------------------------------------------------------------------- 
    7392      ! 
     
    86105         IF(lwp) WRITE(numout,*)' sbc_ana : Constant surface fluxes read in namsbc_ana namelist' 
    87106         IF(lwp) WRITE(numout,*)' ~~~~~~~ ' 
    88          IF(lwp) WRITE(numout,*)'              spin up of the stress  nn_tau000 = ', nn_tau000, ' time-steps' 
    89          IF(lwp) WRITE(numout,*)'              constant i-stress      rn_utau0  = ', rn_utau0 , ' N/m2' 
    90          IF(lwp) WRITE(numout,*)'              constant j-stress      rn_vtau0  = ', rn_vtau0 , ' N/m2' 
    91          IF(lwp) WRITE(numout,*)'              non solar heat flux    rn_qns0   = ', rn_qns0  , ' W/m2' 
    92          IF(lwp) WRITE(numout,*)'              solar heat flux        rn_qsr0   = ', rn_qsr0  , ' W/m2' 
    93          IF(lwp) WRITE(numout,*)'              net heat flux          rn_emp0   = ', rn_emp0  , ' Kg/m2/s' 
     107         IF(lwp) WRITE(numout,*)'              spin up of the stress         nn_tau000 = ', nn_tau000 , ' time-steps' 
     108         IF(lwp) WRITE(numout,*)'              constant i-stress             rn_utau0  = ', rn_utau0  , ' N/m2' 
     109         IF(lwp) WRITE(numout,*)'              constant j-stress             rn_vtau0  = ', rn_vtau0  , ' N/m2' 
     110         IF(lwp) WRITE(numout,*)'              non solar heat flux           rn_qns0   = ', rn_qns0   , ' W/m2' 
     111         IF(lwp) WRITE(numout,*)'              solar heat flux               rn_qsr0   = ', rn_qsr0   , ' W/m2' 
     112         IF(lwp) WRITE(numout,*)'              net freshwater flux           rn_emp0   = ', rn_emp0   , ' Kg/m2/s' 
     113         IF(lwp) WRITE(numout,*)'              constant ice-atm stress       rn_iutau0 = ', rn_iutau0 , ' N/m2' 
     114         IF(lwp) WRITE(numout,*)'              constant ice-atm stress       rn_ivtau0 = ', rn_ivtau0 , ' N/m2' 
     115         IF(lwp) WRITE(numout,*)'              solar heat flux over ice      rn_iqsr0  = ', rn_iqsr0  , ' W/m2' 
     116         IF(lwp) WRITE(numout,*)'              non solar heat flux over ice  rn_iqns0  = ', rn_iqns0  , ' W/m2' 
     117         IF(lwp) WRITE(numout,*)'              snow precip                   rn_sprec0 = ', rn_sprec0 , ' Kg/m2/s' 
     118         IF(lwp) WRITE(numout,*)'              sublimation                   rn_ievap0 = ', rn_ievap0 , ' Kg/m2/s' 
    94119         ! 
    95120         nn_tau000 = MAX( nn_tau000, 1 )     ! must be >= 1 
     
    133158   END SUBROUTINE sbc_ana 
    134159 
    135  
     160#if defined key_lim3 
     161   SUBROUTINE ana_ice_tau 
     162      !!--------------------------------------------------------------------- 
     163      !!                     ***  ROUTINE ana_ice_tau  *** 
     164      !! 
     165      !! ** Purpose :   provide the surface boundary (momentum) condition over sea-ice 
     166      !!--------------------------------------------------------------------- 
     167      utau_ice(:,:) = rn_iutau0 
     168      vtau_ice(:,:) = rn_ivtau0 
     169      
     170   END SUBROUTINE ana_ice_tau 
     171    
     172   SUBROUTINE ana_ice_flx 
     173      !!--------------------------------------------------------------------- 
     174      !!                     ***  ROUTINE ana_ice_flx  *** 
     175      !! 
     176      !! ** Purpose :   provide the surface boundary (flux) condition over sea-ice 
     177      !!--------------------------------------------------------------------- 
     178      REAL(wp), DIMENSION(:,:), POINTER ::   zsnw       ! snw distribution after wind blowing 
     179      !!--------------------------------------------------------------------- 
     180      CALL wrk_alloc( jpi,jpj, zsnw )  
     181 
     182      ! ocean variables (renaming) 
     183      emp_oce (:,:)   = rn_emp0 
     184      qsr_oce (:,:)   = rn_qsr0 
     185      qns_oce (:,:)   = rn_qns0 
     186       
     187      ! ice variables 
     188      alb_ice (:,:,:) = 0.7_wp ! useless 
     189      qsr_ice (:,:,:) = rn_iqsr0 
     190      qns_ice (:,:,:) = rn_iqns0 
     191      sprecip (:,:)   = rn_sprec0 
     192      evap_ice(:,:,:) = rn_ievap0 
     193 
     194      ! ice variables deduced from above 
     195      zsnw(:,:) = 0._wp 
     196      CALL lim_thd_snwblow( pfrld, zsnw )  ! snow distribution over ice after wind blowing  
     197      emp_ice  (:,:)   = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw 
     198      emp_oce  (:,:)   = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw ) 
     199      qevap_ice(:,:,:) =   0._wp 
     200      qprec_ice(:,:)   =   rhosn * ( sst_m(:,:) * cpic - lfus ) * tmask(:,:,1) ! in J/m3 
     201      qemp_oce (:,:)   = - emp_oce(:,:) * sst_m(:,:) * rcp 
     202      qemp_ice (:,:)   =   sprecip(:,:) * zsnw * ( sst_m(:,:) * cpic - lfus ) * tmask(:,:,1) ! solid precip (only) 
     203 
     204      ! total fluxes 
     205      emp_tot (:,:) = emp_ice  + emp_oce  
     206      qns_tot (:,:) = pfrld(:,:) * qns_oce(:,:) + SUM( a_i_b(:,:,:) * qns_ice(:,:,:), dim=3 ) + qemp_ice(:,:) + qemp_oce(:,:) 
     207      qsr_tot (:,:) = pfrld(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 
     208 
     209      !-------------------------------------------------------------------- 
     210      ! FRACTIONs of net shortwave radiation which is not absorbed in the 
     211      ! thin surface layer and penetrates inside the ice cover 
     212      ! ( Maykut and Untersteiner, 1971 ; Ebert and Curry, 1993 ) 
     213      fr1_i0(:,:) = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) 
     214      fr2_i0(:,:) = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) 
     215 
     216      CALL wrk_dealloc( jpi,jpj, zsnw )  
     217       
     218   END SUBROUTINE ana_ice_flx 
     219#endif 
     220 
     221    
    136222   SUBROUTINE sbc_gyre( kt ) 
    137223      !!--------------------------------------------------------------------- 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r7060 r7077  
    3030   USE sbcblk_clio     ! Surface boundary condition: CLIO bulk 
    3131   USE sbccpl          ! Surface boundary condition: coupled interface 
     32   USE sbcana          ! Surface boundary condition: analytic formulation 
    3233   USE albedo          ! ocean & ice albedo 
    3334 
     
    109110      !!--------------------------------------------------------------------- 
    110111      INTEGER, INTENT(in) ::   kt      ! ocean time step 
    111       INTEGER, INTENT(in) ::   kblk    ! type of bulk (=3 CLIO, =4 CORE, =5 COUPLED) 
     112      INTEGER, INTENT(in) ::   kblk    ! type of bulk (=1 ANALYTIC, =3 CLIO, =4 CORE, =5 COUPLED) 
    112113      !! 
    113114      INTEGER  ::   jl                 ! dummy loop index 
     
    156157         !----------------------------------------------------------------- 
    157158         SELECT CASE( kblk ) 
     159            CASE( jp_ana     )   ;    CALL ana_ice_tau                              ! analytic formulation             
    158160            CASE( jp_clio    )   ;    CALL blk_ice_clio_tau                         ! CLIO bulk formulation             
    159161            CASE( jp_core    )   ;    CALL blk_ice_core_tau                         ! CORE bulk formulation 
     
    228230         SELECT CASE( kblk ) 
    229231 
     232            CASE( jp_ana )                                        ! analytic formulation 
     233                                      CALL ana_ice_flx 
     234                
    230235            CASE( jp_clio )                                       ! CLIO bulk formulation 
    231236               ! In CLIO the cloud fraction is read in the climatology and the all-sky albedo  
Note: See TracChangeset for help on using the changeset viewer.