Changeset 5443 for branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2015-06-19T17:18:00+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Property svn:keywords set to Id
r5162 r5443 41 41 USE ice_state, only: aice,aicen,uvel,vvel,vsno,vsnon,vice,vicen 42 42 USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow, & 43 strocnxT,strocnyT, & 43 44 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_gbm, & 44 45 fresh_gbm,fhocn_gbm,fswthru_gbm,frzmlt, & … … 51 52 vsnon,vice,vicen 52 53 USE ice_flux, only: strax,stray,strocnx,strocny,frain,fsnow, & 54 strocnxT,strocnyT, & 53 55 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai, & 54 56 fresh_ai,fhocn_ai,fswthru_ai,frzmlt, & … … 98 100 # include "domzgr_substitute.h90" 99 101 102 !! $Id$ 100 103 CONTAINS 101 104 … … 139 142 IF ( ksbc == jp_flx ) THEN 140 143 CALL cice_sbc_force(kt) 141 ELSE IF ( ksbc == jp_ cpl ) THEN144 ELSE IF ( ksbc == jp_purecpl ) THEN 142 145 CALL sbc_cpl_ice_flx( 1.0-fr_i ) 143 146 ENDIF … … 147 150 CALL cice_sbc_out ( kt, ksbc ) 148 151 149 IF ( ksbc == jp_ cpl ) CALL cice_sbc_hadgam(kt+1)152 IF ( ksbc == jp_purecpl ) CALL cice_sbc_hadgam(kt+1) 150 153 151 154 ENDIF ! End sea-ice time step only … … 202 205 203 206 ! Do some CICE consistency checks 204 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN207 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 205 208 IF ( calc_strair .OR. calc_Tsfc ) THEN 206 209 CALL ctl_stop( 'STOP', 'cice_sbc_init : Forcing option requires calc_strair=F and calc_Tsfc=F in ice_in' ) … … 223 226 224 227 CALL cice2nemo(aice,fr_i, 'T', 1. ) 225 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN228 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 226 229 DO jl=1,ncat 227 230 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 331 334 ! forced and coupled case 332 335 333 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN336 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 334 337 335 338 ztmpn(:,:,:)=0.0 … … 562 565 ! Combine wind stress and ocean-ice stress 563 566 ! [Note that fr_iu hasn't yet been updated, so still from start of CICE timestep] 567 ! strocnx and strocny already weighted by ice fraction in CICE so not done here 564 568 565 569 utau(:,:)=(1.0-fr_iu(:,:))*utau(:,:)-ss_iou(:,:) 566 570 vtau(:,:)=(1.0-fr_iv(:,:))*vtau(:,:)-ss_iov(:,:) 571 572 ! Also need ice/ocean stress on T points so that taum can be updated 573 ! This interpolation is already done in CICE so best to use those values 574 CALL cice2nemo(strocnxT,ztmp1,'T',-1.) 575 CALL cice2nemo(strocnyT,ztmp2,'T',-1.) 576 577 ! Update taum with modulus of ice-ocean stress 578 ! strocnxT and strocnyT are not weighted by ice fraction in CICE so must be done here 579 taum(:,:)=(1.0-fr_i(:,:))*taum(:,:)+fr_i(:,:)*SQRT(ztmp1**2. + ztmp2**2.) 567 580 568 581 ! Freshwater fluxes … … 576 589 ELSE IF (ksbc == jp_core) THEN 577 590 emp(:,:) = (1.0-fr_i(:,:))*emp(:,:) 578 ELSE IF (ksbc == jp_ cpl) THEN591 ELSE IF (ksbc == jp_purecpl) THEN 579 592 ! emp_tot is set in sbc_cpl_ice_flx (called from cice_sbc_in above) 580 593 ! This is currently as required with the coupling fields from the UM atmosphere … … 612 625 ENDIF 613 626 ! Take into account snow melting except for fully coupled when already in qns_tot 614 IF (ksbc == jp_ cpl) THEN627 IF (ksbc == jp_purecpl) THEN 615 628 qsr(:,:)= qsr_tot(:,:) 616 629 qns(:,:)= qns_tot(:,:) … … 647 660 648 661 CALL cice2nemo(aice,fr_i,'T', 1. ) 649 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN662 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 650 663 DO jl=1,ncat 651 664 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 1094 1107 !! Default option Dummy module NO CICE sea-ice model 1095 1108 !!---------------------------------------------------------------------- 1109 !! $Id$ 1096 1110 CONTAINS 1097 1111
Note: See TracChangeset
for help on using the changeset viewer.