- Timestamp:
- 2017-12-13T15:58:53+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r7646 r9019 137 137 CALL cice_sbc_force(kt) 138 138 ELSE IF ( ksbc == jp_purecpl ) THEN 139 CALL sbc_cpl_ice_flx( 1.0-fr_i)139 CALL sbc_cpl_ice_flx( fr_i ) 140 140 ENDIF 141 141 … … 230 230 CALL lbc_lnk ( fr_iv , 'V', 1. ) 231 231 232 ! ! embedded sea ice 233 IF( nn_ice_embd /= 0 ) THEN ! mass exchanges between ice and ocean (case 1 or 2) set the snow+ice mass 234 CALL cice2nemo(vsno(:,:,:),ztmp1,'T', 1. ) 235 CALL cice2nemo(vice(:,:,:),ztmp2,'T', 1. ) 236 snwice_mass (:,:) = ( rhosn * ztmp1(:,:) + rhoic * ztmp2(:,:) ) 237 snwice_mass_b(:,:) = snwice_mass(:,:) 238 ELSE 239 snwice_mass (:,:) = 0.0_wp ! no mass exchanges 240 snwice_mass_b(:,:) = 0.0_wp ! no mass exchanges 241 ENDIF 232 ! set the snow+ice mass 233 CALL cice2nemo(vsno(:,:,:),ztmp1,'T', 1. ) 234 CALL cice2nemo(vice(:,:,:),ztmp2,'T', 1. ) 235 snwice_mass (:,:) = ( rhosn * ztmp1(:,:) + rhoic * ztmp2(:,:) ) 236 snwice_mass_b(:,:) = snwice_mass(:,:) 237 242 238 IF( .NOT.ln_rstart ) THEN 243 IF( nn_ice_embd == 2 ) THEN ! full embedment (case 2)deplete the initial ssh below sea-ice area239 IF( ln_ice_embd ) THEN ! embedded sea-ice: deplete the initial ssh below sea-ice area 244 240 sshn(:,:) = sshn(:,:) - snwice_mass(:,:) * r1_rau0 245 241 sshb(:,:) = sshb(:,:) - snwice_mass(:,:) * r1_rau0 … … 473 469 CALL nemo2cice(ztmp,vocn,'F', -1. ) 474 470 475 IF( nn_ice_embd == 2) THEN !== embedded sea ice: compute representative ice top surface ==!471 IF( ln_ice_embd ) THEN !== embedded sea ice: compute representative ice top surface ==! 476 472 ! 477 473 ! average interpolation coeff as used in dynspg = (1/nn_fsbc) * {SUM[n/nn_fsbc], n=0,nn_fsbc-1} … … 676 672 CALL lbc_lnk ( fr_iv , 'V', 1. ) 677 673 678 ! ! embedded sea ice 679 IF( nn_ice_embd /= 0 ) THEN ! mass exchanges between ice and ocean (case 1 or 2) set the snow+ice mass 680 CALL cice2nemo(vsno(:,:,:),ztmp1,'T', 1. ) 681 CALL cice2nemo(vice(:,:,:),ztmp2,'T', 1. ) 682 snwice_mass (:,:) = ( rhosn * ztmp1(:,:) + rhoic * ztmp2(:,:) ) 683 snwice_mass_b(:,:) = snwice_mass(:,:) 684 snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt 685 ENDIF 674 ! set the snow+ice mass 675 CALL cice2nemo(vsno(:,:,:),ztmp1,'T', 1. ) 676 CALL cice2nemo(vice(:,:,:),ztmp2,'T', 1. ) 677 snwice_mass (:,:) = ( rhosn * ztmp1(:,:) + rhoic * ztmp2(:,:) ) 678 snwice_mass_b(:,:) = snwice_mass(:,:) 679 snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt 686 680 687 681 ! Release work space … … 716 710 ! ! Prepare Coupling fields ! 717 711 ! ! =========================== ! 718 719 ! x and y comp of ice velocity720 712 ! 713 ! x and y comp of ice velocity 714 ! 721 715 CALL cice2nemo(uvel,u_ice,'F', -1. ) 722 716 CALL cice2nemo(vvel,v_ice,'F', -1. ) 723 724 ! Ice concentration (CO_1) = a_i calculated at end of cice_sbc_out725 726 ! Snow and ice thicknesses (CO_2 and CO_3)727 728 DO jl = 1, ncat729 CALL cice2nemo( vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. )730 CALL cice2nemo( vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1. )731 END DO717 ! 718 ! Ice concentration (CO_1) = a_i calculated at end of cice_sbc_out 719 ! 720 ! Snow and ice thicknesses (CO_2 and CO_3) 721 ! 722 DO jl = 1, ncat 723 CALL cice2nemo( vsnon(:,:,jl,:), h_s(:,:,jl),'T', 1. ) 724 CALL cice2nemo( vicen(:,:,jl,:), h_i(:,:,jl),'T', 1. ) 725 END DO 732 726 ! 733 727 IF( nn_timing == 1 ) CALL timing_stop('cice_sbc_hadgam')
Note: See TracChangeset
for help on using the changeset viewer.