Changeset 5133 for trunk/NEMOGCM
- Timestamp:
- 2015-03-09T10:10:16+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r4990 r5133 40 40 # if defined key_cice4 41 41 USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow, & 42 strocnxT,strocnyT, & 42 43 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_gbm, & 43 44 fresh_gbm,fhocn_gbm,fswthru_gbm,frzmlt, & … … 48 49 #else 49 50 USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow, & 51 strocnxT,strocnyT, & 50 52 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai, & 51 53 fresh_ai,fhocn_ai,fswthru_ai,frzmlt, & … … 560 562 ! Combine wind stress and ocean-ice stress 561 563 ! [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 562 565 563 566 utau(:,:)=(1.0-fr_iu(:,:))*utau(:,:)-ss_iou(:,:) 564 567 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 576 taum(:,:)=(1.0-fr_i(:,:))*taum(:,:)+fr_i(:,:)*SQRT(ztmp1**2. + ztmp2**2.) 565 577 566 578 ! Freshwater fluxes
Note: See TracChangeset
for help on using the changeset viewer.