Changeset 5612 for branches/UKMO/dev_r5518_CICE_coupling_GSI7/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2015-07-20T11:11:11+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_CICE_coupling_GSI7/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Property svn:keywords set to Id
r5522 r5612 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,nt_Tsfc 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, & … … 106 108 # include "domzgr_substitute.h90" 107 109 110 !! $Id$ 108 111 CONTAINS 109 112 … … 147 150 IF ( ksbc == jp_flx ) THEN 148 151 CALL cice_sbc_force(kt) 149 ELSE IF ( ksbc == jp_ cpl ) THEN152 ELSE IF ( ksbc == jp_purecpl ) THEN 150 153 CALL sbc_cpl_ice_flx( 1.0-fr_i ) 151 154 ENDIF … … 155 158 CALL cice_sbc_out ( kt, ksbc ) 156 159 157 IF ( ksbc == jp_ cpl ) CALL cice_sbc_hadgam(kt+1)160 IF ( ksbc == jp_purecpl ) CALL cice_sbc_hadgam(kt+1) 158 161 159 162 ENDIF ! End sea-ice time step only … … 209 212 210 213 ! Do some CICE consistency checks 211 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN214 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 212 215 IF ( calc_strair .OR. calc_Tsfc ) THEN 213 216 CALL ctl_stop( 'STOP', 'cice_sbc_init : Forcing option requires calc_strair=F and calc_Tsfc=F in ice_in' ) … … 231 234 232 235 CALL cice2nemo(aice,fr_i, 'T', 1. ) 233 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN236 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 234 237 DO jl=1,ncat 235 238 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 339 342 ! forced and coupled case 340 343 341 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN344 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 342 345 343 346 ztmpn(:,:,:)=0.0 … … 521 524 CALL nemo2cice(ztmp,ss_tlty,'F', -1. ) 522 525 523 CALL wrk_dealloc( jpi,jpj, ztmp )526 CALL wrk_dealloc( jpi,jpj, ztmp, zpice ) 524 527 CALL wrk_dealloc( jpi,jpj,ncat, ztmpn ) 525 528 ! … … 575 578 ! Combine wind stress and ocean-ice stress 576 579 ! [Note that fr_iu hasn't yet been updated, so still from start of CICE timestep] 580 ! strocnx and strocny already weighted by ice fraction in CICE so not done here 577 581 578 582 utau(:,:)=(1.0-fr_iu(:,:))*utau(:,:)-ss_iou(:,:) 579 583 vtau(:,:)=(1.0-fr_iv(:,:))*vtau(:,:)-ss_iov(:,:) 584 585 ! Also need ice/ocean stress on T points so that taum can be updated 586 ! This interpolation is already done in CICE so best to use those values 587 CALL cice2nemo(strocnxT,ztmp1,'T',-1.) 588 CALL cice2nemo(strocnyT,ztmp2,'T',-1.) 589 590 ! Update taum with modulus of ice-ocean stress 591 ! strocnxT and strocnyT are not weighted by ice fraction in CICE so must be done here 592 taum(:,:)=(1.0-fr_i(:,:))*taum(:,:)+fr_i(:,:)*SQRT(ztmp1**2. + ztmp2**2.) 580 593 581 594 ! Freshwater fluxes … … 589 602 ELSE IF (ksbc == jp_core) THEN 590 603 emp(:,:) = (1.0-fr_i(:,:))*emp(:,:) 591 ELSE IF (ksbc == jp_ cpl) THEN604 ELSE IF (ksbc == jp_purecpl) THEN 592 605 ! emp_tot is set in sbc_cpl_ice_flx (called from cice_sbc_in above) 593 606 ! This is currently as required with the coupling fields from the UM atmosphere … … 625 638 ENDIF 626 639 ! Take into account snow melting except for fully coupled when already in qns_tot 627 IF (ksbc == jp_ cpl) THEN640 IF (ksbc == jp_purecpl) THEN 628 641 qsr(:,:)= qsr_tot(:,:) 629 642 qns(:,:)= qns_tot(:,:) … … 660 673 661 674 CALL cice2nemo(aice,fr_i,'T', 1. ) 662 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN675 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 663 676 DO jl=1,ncat 664 677 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 1127 1140 !! Default option Dummy module NO CICE sea-ice model 1128 1141 !!---------------------------------------------------------------------- 1142 !! $Id$ 1129 1143 CONTAINS 1130 1144
Note: See TracChangeset
for help on using the changeset viewer.