Changeset 5682 for branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Timestamp:
- 2015-08-12T17:46:45+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5072_UKMO2_OBS_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
- Property svn:keywords set to Id
r4990 r5682 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, & … … 94 96 # include "domzgr_substitute.h90" 95 97 98 !! $Id$ 96 99 CONTAINS 97 100 … … 135 138 IF ( ksbc == jp_flx ) THEN 136 139 CALL cice_sbc_force(kt) 137 ELSE IF ( ksbc == jp_ cpl ) THEN140 ELSE IF ( ksbc == jp_purecpl ) THEN 138 141 CALL sbc_cpl_ice_flx( 1.0-fr_i ) 139 142 ENDIF … … 143 146 CALL cice_sbc_out ( kt, ksbc ) 144 147 145 IF ( ksbc == jp_ cpl ) CALL cice_sbc_hadgam(kt+1)148 IF ( ksbc == jp_purecpl ) CALL cice_sbc_hadgam(kt+1) 146 149 147 150 ENDIF ! End sea-ice time step only … … 184 187 185 188 ! Do some CICE consistency checks 186 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN189 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 187 190 IF ( calc_strair .OR. calc_Tsfc ) THEN 188 191 CALL ctl_stop( 'STOP', 'cice_sbc_init : Forcing option requires calc_strair=F and calc_Tsfc=F in ice_in' ) … … 209 212 210 213 CALL cice2nemo(aice,fr_i, 'T', 1. ) 211 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_ cpl) ) THEN214 IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN 212 215 DO jl=1,ncat 213 216 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 316 319 ! forced and coupled case 317 320 318 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN321 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 319 322 320 323 ztmpn(:,:,:)=0.0 … … 506 509 CALL nemo2cice(ztmp,ss_tlty,'F', -1. ) 507 510 508 CALL wrk_dealloc( jpi,jpj, ztmp )511 CALL wrk_dealloc( jpi,jpj, ztmp, zpice ) 509 512 CALL wrk_dealloc( jpi,jpj,ncat, ztmpn ) 510 513 ! … … 560 563 ! Combine wind stress and ocean-ice stress 561 564 ! [Note that fr_iu hasn't yet been updated, so still from start of CICE timestep] 565 ! strocnx and strocny already weighted by ice fraction in CICE so not done here 562 566 563 567 utau(:,:)=(1.0-fr_iu(:,:))*utau(:,:)-ss_iou(:,:) 564 568 vtau(:,:)=(1.0-fr_iv(:,:))*vtau(:,:)-ss_iov(:,:) 569 570 ! Also need ice/ocean stress on T points so that taum can be updated 571 ! This interpolation is already done in CICE so best to use those values 572 CALL cice2nemo(strocnxT,ztmp1,'T',-1.) 573 CALL cice2nemo(strocnyT,ztmp2,'T',-1.) 574 575 ! Update taum with modulus of ice-ocean stress 576 ! strocnxT and strocnyT are not weighted by ice fraction in CICE so must be done here 577 taum(:,:)=(1.0-fr_i(:,:))*taum(:,:)+fr_i(:,:)*SQRT(ztmp1**2. + ztmp2**2.) 565 578 566 579 ! Freshwater fluxes … … 574 587 ELSE IF (ksbc == jp_core) THEN 575 588 emp(:,:) = (1.0-fr_i(:,:))*emp(:,:) 576 ELSE IF (ksbc == jp_ cpl) THEN589 ELSE IF (ksbc == jp_purecpl) THEN 577 590 ! emp_tot is set in sbc_cpl_ice_flx (called from cice_sbc_in above) 578 591 ! This is currently as required with the coupling fields from the UM atmosphere … … 610 623 ENDIF 611 624 ! Take into account snow melting except for fully coupled when already in qns_tot 612 IF (ksbc == jp_ cpl) THEN625 IF (ksbc == jp_purecpl) THEN 613 626 qsr(:,:)= qsr_tot(:,:) 614 627 qns(:,:)= qns_tot(:,:) … … 645 658 646 659 CALL cice2nemo(aice,fr_i,'T', 1. ) 647 IF ( (ksbc == jp_flx).OR.(ksbc == jp_ cpl) ) THEN660 IF ( (ksbc == jp_flx).OR.(ksbc == jp_purecpl) ) THEN 648 661 DO jl=1,ncat 649 662 CALL cice2nemo(aicen(:,:,jl,:),a_i(:,:,jl), 'T', 1. ) … … 1083 1096 !! Default option Dummy module NO CICE sea-ice model 1084 1097 !!---------------------------------------------------------------------- 1098 !! $Id$ 1085 1099 CONTAINS 1086 1100
Note: See TracChangeset
for help on using the changeset viewer.