Ignore:
Timestamp:
2016-03-22T18:17:23+01:00 (5 years ago)
Author:
clem
Message:

solve issues for 1)wrong fresh water budget sent to passive tracers and 2)ice-air coupling. See tickets #1694 #1696 #1697 and #1698

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r6316 r6399  
    125125      IF( iom_use('emp_ice' ) )  CALL iom_put( "emp_ice"  , emp_ice(:,:) )   ! emp over ice   (taking into account the snow blown away from the ice) 
    126126 
    127       ! clem 2016: albedo output 
     127      ! albedo output 
    128128      CALL wrk_alloc( jpi,jpj, zalb )     
    129129 
     
    158158            hfx_out(ji,jj) = hfx_out(ji,jj) + zqmass + zqsr 
    159159 
    160             ! Add the residual from heat diffusion equation (W.m-2) 
    161             !------------------------------------------------------- 
    162             hfx_out(ji,jj) = hfx_out(ji,jj) + hfx_err_dif(ji,jj) 
     160            ! Add the residual from heat diffusion equation and sublimation (W.m-2) 
     161            !---------------------------------------------------------------------- 
     162            hfx_out(ji,jj) = hfx_out(ji,jj) + hfx_err_dif(ji,jj) +   & 
     163               &           ( hfx_sub(ji,jj) - SUM( qevap_ice(ji,jj,:) * a_i_b(ji,jj,:) ) ) 
    163164 
    164165            ! New qsr and qns used to compute the oceanic heat flux at the next time step 
    165             !--------------------------------------------------- 
     166            !---------------------------------------------------------------------------- 
    166167            qsr(ji,jj) = zqsr                                       
    167168            qns(ji,jj) = hfx_out(ji,jj) - zqsr               
     
    183184 
    184185            ! mass flux at the ocean/ice interface 
    185             fmmflx(ji,jj) = - ( wfx_ice(ji,jj) + wfx_snw(ji,jj) ) * r1_rdtice  ! F/M mass flux save at least for biogeochemical model 
    186             emp(ji,jj)    = emp_oce(ji,jj) - wfx_ice(ji,jj) - wfx_snw(ji,jj)   ! mass flux + F/M mass flux (always ice/ocean mass exchange) 
    187              
     186            fmmflx(ji,jj) = - ( wfx_ice(ji,jj) + wfx_snw(ji,jj) + wfx_err_sub(ji,jj) )              ! F/M mass flux save at least for biogeochemical model 
     187            emp(ji,jj)    = emp_oce(ji,jj) - wfx_ice(ji,jj) - wfx_snw(ji,jj) - wfx_err_sub(ji,jj)   ! mass flux + F/M mass flux (always ice/ocean mass exchange) 
    188188         END DO 
    189189      END DO 
     
    193193      !------------------------------------------! 
    194194      sfx(:,:) = sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) + sfx_opw(:,:)   & 
    195          &     + sfx_res(:,:) + sfx_dyn(:,:) + sfx_bri(:,:) 
     195         &     + sfx_res(:,:) + sfx_dyn(:,:) + sfx_bri(:,:) + sfx_sub(:,:) 
    196196 
    197197      !-------------------------------------------------------------! 
Note: See TracChangeset for help on using the changeset viewer.