Changeset 8092


Ignore:
Timestamp:
2017-05-30T10:10:54+02:00 (3 years ago)
Author:
cbricaud
Message:

use masked horizontal surface to coarsene surface fluxes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsfld.F90

    r7805 r8092  
    304304      nk_rnf_crs=INT(z2d_crs) 
    305305 
    306       CALL crs_dom_ope( qsr  , 'SUM', 'T', tmask, qsr_crs  , p_e12=e1e2t, p_surf_crs=e1e2t_crs, psgn=1.0 ) 
    307       CALL crs_dom_ope( wndm , 'SUM', 'T', tmask, wndm_crs , p_e12=e1e2t, p_surf_crs=e1e2t_crs, psgn=1.0 ) 
    308       CALL crs_dom_ope( emp   ,'SUM', 'T', tmask, emp_crs   , p_e12=e1e2t, p_surf_crs=e1e2t_crs, psgn=1.0 ) 
    309       CALL crs_dom_ope( fmmflx,'SUM', 'T', tmask, fmmflx_crs, p_e12=e1e2t, p_surf_crs=e1e2t_crs, psgn=1.0 ) 
    310  
    311       CALL crs_dom_ope( fr_i  ,'SUM', 'T', tmask, fr_i_crs  , p_e12=e1e2t, p_surf_crs=e1e2t_crs, psgn=1.0 ) 
     306      CALL crs_dom_ope( qsr  , 'SUM', 'T', tmask, qsr_crs   , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 
     307      CALL crs_dom_ope( wndm , 'SUM', 'T', tmask, wndm_crs  , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 
     308      CALL crs_dom_ope( emp   ,'SUM', 'T', tmask, emp_crs   , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 
     309      CALL crs_dom_ope( fmmflx,'SUM', 'T', tmask, fmmflx_crs, p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 
     310 
     311      CALL crs_dom_ope( fr_i  ,'SUM', 'T', tmask, fr_i_crs  , p_e12=e1e2t, p_surf_crs=e1e2w_msk(:,:,1), psgn=1.0 ) 
    312312      fr_i_crs=MAX( 0._wp, MIN( fr_i_crs , 1._wp ) ) 
    313313 
     
    323323      CALL iom_put( "qsr"      , qsr_crs  )   ! qsr output  
    324324      CALL iom_put( "empmr"    , emp_crs  )   ! water flux output  
    325       CALL iom_put( "saltflx"  , sfx_crs  )   ! salt flux output  
     325      CALL iom_put( "saltflx"  , fmmflx_crs  )   ! salt flux output  
    326326      CALL iom_put( "ice_cover", fr_i_crs )   ! ice cover output  
    327327 
Note: See TracChangeset for help on using the changeset viewer.