Ignore:
Timestamp:
2015-02-02T11:28:50+01:00 (6 years ago)
Author:
vancop
Message:

new itd boundaries, namelist changes, mono-category and comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r4990 r5048  
    191191            CALL sbc_cpl_ice_tau( utau_ice , vtau_ice ) 
    192192 
    193             ! MV -> seb  
    194 !           CALL sbc_cpl_ice_flx( p_frld=ato_i, palbi=zalb_ice, psst=sst_m, pist=t_su    ) 
    195  
    196 !           IF( nn_limflx == 2 )   CALL ice_lim_flx( t_su, zalb_ice, qns_ice, qsr_ice ,   & 
    197 !              &                                           dqns_ice, qla_ice, dqla_ice, nn_limflx ) 
    198 !           ! Latent heat flux is forced to 0 in coupled : 
    199 !           !  it is included in qns (non-solar heat flux) 
    200 !           qla_ice  (:,:,:) = 0._wp 
    201 !           dqla_ice (:,:,:) = 0._wp 
    202             ! END MV -> seb 
    203             ! 
    204193         END SELECT 
    205194          
     
    236225         wfx_spr(:,:) = 0._wp   ;    
    237226 
     227         afx_tot(:,:) = at_i(:,:) ;  afx_dyn(:,:) = 0._wp 
     228         afx_thd(:,:) = 0._wp 
     229 
    238230         hfx_in (:,:) = 0._wp   ;   hfx_out(:,:) = 0._wp 
    239231         hfx_thd(:,:) = 0._wp   ;    
     
    252244         ! ---------------------------------------------- 
    253245         IF( .NOT. lk_c1d ) THEN 
     246 
     247         IF ( ln_limdyn ) afx_dyn(:,:)     = at_i(:,:) 
     248 
    254249                          CALL lim_dyn( kt )              ! Ice dynamics    ( rheology/dynamics ) 
    255250                          CALL lim_trp( kt )              ! Ice transport   ( Advection/diffusion ) 
    256251                          CALL lim_var_glo2eqv            ! equivalent variables, requested for rafting 
    257252         IF( ln_nicep )   CALL lim_prt_state( kt, jiindx, jjindx,-1, ' - ice dyn & trp - ' )   ! control print 
    258                          CALL lim_itd_me                 ! Mechanical redistribution ! (ridging/rafting) 
     253         IF( nn_monocat /= 2 ) CALL lim_itd_me                 ! Mechanical redistribution ! (ridging/rafting) 
    259254                          CALL lim_var_agg( 1 )  
    260255#if defined key_bdy 
     
    278273                          oa_i_b (:,:,:)   = oa_i (:,:,:) 
    279274                          smv_i_b(:,:,:)   = smv_i(:,:,:) 
     275 
     276         IF ( ln_limdyn ) afx_dyn(:,:)     = ( at_i(:,:) - afx_dyn(:,:) ) * r1_rdtice 
     277                          afx_thd(:,:)     = at_i(:,:) 
    280278  
    281279         ! ---------------------------------------------- 
    282          ! ice thermodynamic 
     280         ! ice thermodynamics 
    283281         ! ---------------------------------------------- 
    284282                          CALL lim_var_glo2eqv            ! equivalent variables 
     
    288286                          phicif(:,:)  = vt_i(:,:) 
    289287 
    290                           ! MV -> seb 
    291288                          SELECT CASE( kblk ) 
    292289                             CASE ( jp_cpl ) 
     
    299296                             dqla_ice (:,:,:) = 0._wp 
    300297                          END SELECT 
    301                           ! END MV -> seb 
    302298                          ! 
    303299                          CALL lim_var_bv                 ! bulk brine volume (diag) 
     
    320316         !                                           ! Diagnostics and outputs  
    321317         IF (ln_limdiaout) CALL lim_diahsb 
     318 
     319                          afx_thd(:,:)     = ( at_i(:,:) - afx_thd(:,:) ) * r1_rdtice 
     320                          afx_tot(:,:)     = ( at_i(:,:) - afx_tot(:,:) ) * r1_rdtice 
    322321 
    323322                          CALL lim_wri( 1  )              ! Ice outputs  
Note: See TracChangeset for help on using the changeset viewer.