Changeset 5133


Ignore:
Timestamp:
2015-03-09T10:10:16+01:00 (6 years ago)
Author:
timgraham
Message:

Fix for ticket #1446. Taum now updated with ice-ocean stress under sea-ice when using cice.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r4990 r5133  
    4040# if defined key_cice4 
    4141   USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow,  & 
     42                strocnxT,strocnyT,                               &  
    4243                sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_gbm,     & 
    4344                fresh_gbm,fhocn_gbm,fswthru_gbm,frzmlt,          & 
     
    4849#else 
    4950   USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow,  & 
     51                strocnxT,strocnyT,                               &  
    5052                sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai,     & 
    5153                fresh_ai,fhocn_ai,fswthru_ai,frzmlt,          & 
     
    560562! Combine wind stress and ocean-ice stress 
    561563! [Note that fr_iu hasn't yet been updated, so still from start of CICE timestep] 
     564! strocnx and strocny already weighted by ice fraction in CICE so not done here  
    562565 
    563566      utau(:,:)=(1.0-fr_iu(:,:))*utau(:,:)-ss_iou(:,:) 
    564567      vtau(:,:)=(1.0-fr_iv(:,:))*vtau(:,:)-ss_iov(:,:)      
     568  
     569! Also need ice/ocean stress on T points so that taum can be updated  
     570! This interpolation is already done in CICE so best to use those values  
     571      CALL cice2nemo(strocnxT,ztmp1,'T',-1.)  
     572      CALL cice2nemo(strocnyT,ztmp2,'T',-1.)  
     573  
     574! Update taum with modulus of ice-ocean stress  
     575! strocnxT and strocnyT are not weighted by ice fraction in CICE so must be done here  
     576taum(:,:)=(1.0-fr_i(:,:))*taum(:,:)+fr_i(:,:)*SQRT(ztmp1**2. + ztmp2**2.)  
    565577 
    566578! Freshwater fluxes  
Note: See TracChangeset for help on using the changeset viewer.