Changeset 1470 for trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
- Timestamp:
- 2009-06-11T15:01:25+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r1469 r1470 108 108 ! 109 109 ! ... mean surface ocean current at ice dynamics point 110 ! B-grid dynamics : I-point (F-point with sea-ice indexation)110 ! B-grid dynamics : I-point 111 111 DO jj = 2, jpj 112 112 DO ji = 2, jpi ! B grid : no vector opt. 113 u i_oce(ji,jj) = 0.5 * ( ssu_m(ji-1,jj ) + ssu_m(ji-1,jj-1) ) * tmu(ji,jj)114 v i_oce(ji,jj) = 0.5 * ( ssv_m(ji ,jj-1) + ssv_m(ji-1,jj-1) ) * tmu(ji,jj)113 u_oce(ji,jj) = 0.5 * ( ssu_m(ji-1,jj ) + ssu_m(ji-1,jj-1) ) * tmu(ji,jj) 114 v_oce(ji,jj) = 0.5 * ( ssv_m(ji ,jj-1) + ssv_m(ji-1,jj-1) ) * tmu(ji,jj) 115 115 END DO 116 116 END DO 117 CALL lbc_lnk( u i_oce, 'I', -1. ) ! I-point (i.e. F-point with ice indices)118 CALL lbc_lnk( v i_oce, 'I', -1. ) ! I-point (i.e. F-point with ice indices)117 CALL lbc_lnk( u_oce, 'I', -1. ) ! I-point (i.e. F-point with ice indices) 118 CALL lbc_lnk( v_oce, 'I', -1. ) ! I-point (i.e. F-point with ice indices) 119 119 120 120 ! ... masked sea surface freezing temperature [Kelvin] (set to rt0 over land) … … 148 148 149 149 CASE( 4 ) ! CORE bulk formulation 150 CALL blk_ice_core( zsist, u i_ice , vi_ice, zalb_ice_cs, &150 CALL blk_ice_core( zsist, u_ice , v_ice , zalb_ice_cs, & 151 151 & utau_ice , vtau_ice , qns_ice , qsr_ice, & 152 152 & qla_ice , dqns_ice , dqla_ice , & … … 162 162 CALL prt_ctl( tab2d_1=utau_ice,clinfo1=' sbc_ice_lim: utau_ice: ', tab2d_2=vtau_ice, clinfo2=' vtau_ice: ' ) 163 163 CALL prt_ctl( tab2d_1=sst_m ,clinfo1=' sbc_ice_lim: sst : ', tab2d_2=sss_m , clinfo2=' sss : ' ) 164 CALL prt_ctl( tab2d_1=u i_oce ,clinfo1=' sbc_ice_lim: u_io : ', tab2d_2=vi_oce, clinfo2=' v_io : ' )164 CALL prt_ctl( tab2d_1=u_oce ,clinfo1=' sbc_ice_lim: u_io : ', tab2d_2=v_oce , clinfo2=' v_io : ' ) 165 165 CALL prt_ctl( tab2d_1=hsnif ,clinfo1=' sbc_ice_lim: hsnif 1: ', tab2d_2=hicif , clinfo2=' hicif : ' ) 166 166 CALL prt_ctl( tab2d_1=frld ,clinfo1=' sbc_ice_lim: frld 1: ', tab2d_2=sist , clinfo2=' sist : ' )
Note: See TracChangeset
for help on using the changeset viewer.