Changeset 1469
- Timestamp:
- 2009-06-11T14:20:08+02:00 (15 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 13 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/limrhg_2.F90
r1465 r1469 176 176 177 177 ! Wind stress. 178 ! always provide stress at I-point (ocean F-point)178 ! always provide stress at I-point 179 179 ztagnx = ( zi1(ji,jj ) * wght(ji,jj,2,2) + zi1(ji-1,jj ) * wght(ji,jj,1,2) & 180 & + zi1(ji,jj-1) * wght(ji,jj,2,1) + zi1(ji-1,jj-1) * wght(ji,jj,1,1) ) * zusw * utau i_ice(ji,jj)180 & + zi1(ji,jj-1) * wght(ji,jj,2,1) + zi1(ji-1,jj-1) * wght(ji,jj,1,1) ) * zusw * utau_ice(ji,jj) 181 181 ztagny = ( zi2(ji,jj ) * wght(ji,jj,2,2) + zi2(ji-1,jj ) * wght(ji,jj,1,2) & 182 & + zi2(ji,jj-1) * wght(ji,jj,2,1) + zi2(ji-1,jj-1) * wght(ji,jj,1,1) ) * zusw * vtau i_ice(ji,jj)182 & + zi2(ji,jj-1) * wght(ji,jj,2,1) + zi2(ji-1,jj-1) * wght(ji,jj,1,1) ) * zusw * vtau_ice(ji,jj) 183 183 184 184 ! Gradient of ice strength -
trunk/NEMO/LIM_SRC_2/limwri_2.F90
r1465 r1469 155 155 ! See thersf for the coefficient 156 156 zcmo(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce !!gm ??? 157 zcmo(ji,jj,15) = utau i_ice(ji,jj)158 zcmo(ji,jj,16) = vtau i_ice(ji,jj)157 zcmo(ji,jj,15) = utau_ice(ji,jj) 158 zcmo(ji,jj,16) = vtau_ice(ji,jj) 159 159 zcmo(ji,jj,17) = qsr_ice(ji,jj,1) 160 160 zcmo(ji,jj,18) = qns_ice(ji,jj,1) -
trunk/NEMO/LIM_SRC_2/limwri_dimg_2.h90
r1463 r1469 121 121 ! See thersf for the coefficient 122 122 zcmo(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 123 zcmo(ji,jj,15) = utau i_ice(ji,jj)124 zcmo(ji,jj,16) = vtau i_ice(ji,jj)123 zcmo(ji,jj,15) = utau_ice(ji,jj) 124 zcmo(ji,jj,16) = vtau_ice(ji,jj) 125 125 zcmo(ji,jj,17) = qsr_ice(ji,jj,1) 126 126 zcmo(ji,jj,18) = qns_ice(ji,jj,1) … … 164 164 ! See thersf for the coefficient 165 165 rcmoy(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 166 rcmoy(ji,jj,15) = utau i_ice(ji,jj)167 rcmoy(ji,jj,16) = vtau i_ice(ji,jj)166 rcmoy(ji,jj,15) = utau_ice(ji,jj) 167 rcmoy(ji,jj,16) = vtau_ice(ji,jj) 168 168 rcmoy(ji,jj,17) = qsr_ice(ji,jj,1) 169 169 rcmoy(ji,jj,18) = qns_ice(ji,jj,1) -
trunk/NEMO/LIM_SRC_3/limrhg.F90
r1465 r1469 281 281 & / (e2t(ji,jj+1)+e2t(ji,jj)) * tmv(ji,jj) 282 282 283 ! Wind stress .284 ztagnx = ( 1. - zfrld1(ji,jj) ) * utau i_ice(ji,jj)285 ztagny = ( 1. - zfrld2(ji,jj) ) * vtau i_ice(ji,jj)283 ! Wind stress at U,V-point 284 ztagnx = ( 1. - zfrld1(ji,jj) ) * utau_ice(ji,jj) 285 ztagny = ( 1. - zfrld2(ji,jj) ) * vtau_ice(ji,jj) 286 286 287 287 ! Computation of the velocity field taking into account the ice internal interaction. -
trunk/NEMO/LIM_SRC_3/limrst.F90
r1465 r1469 159 159 CALL iom_rstput( iter, nitrst, numriw, 'u_ice' , u_ice ) 160 160 CALL iom_rstput( iter, nitrst, numriw, 'v_ice' , v_ice ) 161 CALL iom_rstput( iter, nitrst, numriw, 'utaui_ice' , utaui_ice )162 CALL iom_rstput( iter, nitrst, numriw, 'vtaui_ice' , vtaui_ice )163 161 CALL iom_rstput( iter, nitrst, numriw, 'fsbbq' , fsbbq ) 164 162 CALL iom_rstput( iter, nitrst, numriw, 'stress1_i' , stress1_i ) … … 519 517 CALL iom_get( numrir, jpdom_autoglo, 'u_ice' , u_ice ) 520 518 CALL iom_get( numrir, jpdom_autoglo, 'v_ice' , v_ice ) 521 CALL iom_get( numrir, jpdom_autoglo, 'utaui_ice' , utaui_ice )522 CALL iom_get( numrir, jpdom_autoglo, 'vtaui_ice' , vtaui_ice )523 519 CALL iom_get( numrir, jpdom_autoglo, 'fsbbq' , fsbbq ) 524 520 CALL iom_get( numrir, jpdom_autoglo, 'stress1_i' , stress1_i ) -
trunk/NEMO/LIM_SRC_3/limthd_lac.F90
r1465 r1469 265 265 !------------- 266 266 ! C-grid wind stress components 267 ztaux = ( utau i_ice(ji-1,jj ) * tmu(ji-1,jj ) &268 + utaui_ice(ji ,jj ) * tmu(ji ,jj ) ) / 2.0269 ztauy = ( vtau i_ice(ji ,jj-1) * tmv(ji ,jj-1) &270 + vtaui_ice(ji ,jj ) * tmv(ji ,jj ) ) / 2.0267 ztaux = ( utau_ice(ji-1,jj ) * tmu(ji-1,jj ) & 268 & + utau_ice(ji ,jj ) * tmu(ji ,jj ) ) / 2.0 269 ztauy = ( vtau_ice(ji ,jj-1) * tmv(ji ,jj-1) & 270 & + vtau_ice(ji ,jj ) * tmv(ji ,jj ) ) / 2.0 271 271 ! Square root of wind stress 272 272 ztenagm = SQRT( SQRT( ztaux * ztaux + ztauy * ztauy ) ) -
trunk/NEMO/LIM_SRC_3/limupdate.F90
r1465 r1469 1087 1087 CALL prt_ctl_info(' ~~~~~~~~~~ ') 1088 1088 CALL prt_ctl(tab2d_1=utau , clinfo1= ' lim_update : utau : ', tab2d_2=vtau , clinfo2= ' vtau : ') 1089 CALL prt_ctl(tab2d_1=utau i_ice , clinfo1= ' lim_update : utaui_ice : ', tab2d_2=vtaui_ice , clinfo2= ' vtaui_ice: ')1089 CALL prt_ctl(tab2d_1=utau_ice , clinfo1= ' lim_update : utau_ice : ', tab2d_2=vtau_ice , clinfo2= ' vtau_ice : ') 1090 1090 CALL prt_ctl(tab2d_1=u_oce , clinfo1= ' lim_update : u_oce : ', tab2d_2=v_oce , clinfo2= ' v_oce : ') 1091 1091 ENDIF -
trunk/NEMO/LIM_SRC_3/limwri.F90
r1465 r1469 263 263 zcmo(ji,jj,13) = qns(ji,jj) 264 264 zcmo(ji,jj,14) = fhbri(ji,jj) 265 zcmo(ji,jj,15) = utau i_ice(ji,jj)266 zcmo(ji,jj,16) = vtau i_ice(ji,jj)265 zcmo(ji,jj,15) = utau_ice(ji,jj) 266 zcmo(ji,jj,16) = vtau_ice(ji,jj) 267 267 zcmo(ji,jj,17) = zcmo(ji,jj,17) + (1.0-at_i(ji,jj))*qsr(ji,jj) 268 268 zcmo(ji,jj,18) = zcmo(ji,jj,18) + (1.0-at_i(ji,jj))*qns(ji,jj) -
trunk/NEMO/LIM_SRC_3/limwri_dimg.h90
r1156 r1469 119 119 ! See thersf for the coefficient 120 120 zcmo(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 121 zcmo(ji,jj,15) = utau i_ice(ji,jj)122 zcmo(ji,jj,16) = vtau i_ice(ji,jj)121 zcmo(ji,jj,15) = utau_ice(ji,jj) 122 zcmo(ji,jj,16) = vtau_ice(ji,jj) 123 123 zcmo(ji,jj,17) = qsr (ji,jj) 124 124 zcmo(ji,jj,18) = qns(ji,jj) … … 162 162 ! See thersf for the coefficient 163 163 rcmoy(ji,jj,14) = - emps(ji,jj) * rday * ( sss_m(ji,jj) + epsi16 ) / soce 164 rcmoy(ji,jj,15) = utau i_ice(ji,jj)165 rcmoy(ji,jj,16) = vtau i_ice(ji,jj)164 rcmoy(ji,jj,15) = utau_ice(ji,jj) 165 rcmoy(ji,jj,16) = vtau_ice(ji,jj) 166 166 rcmoy(ji,jj,17) = qsr(ji,jj) 167 167 rcmoy(ji,jj,18) = qns(ji,jj) -
trunk/NEMO/OPA_SRC/SBC/sbc_ice.F90
r1467 r1469 4 4 !! parameter and variables defined in memory in forced mode 5 5 !!====================================================================== 6 !! History : 9.0 ! 06-08 (G. Modec) Surface module 6 !! History : 3.0 ! 2006-08 (G. Madec) Surface module 7 !! 3.2 ! 2009-06 (S. Masson) merge with ice_oce 7 8 !!---------------------------------------------------------------------- 8 9 #if defined key_lim3 || defined key_lim2 … … 40 41 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) :: alb_ice !: albedo of ice 41 42 42 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: utau i_ice !: u-stress over ice (I-point) [N/m2]43 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: vtau i_ice !: v-stress over ice (I-point) [N/m2]43 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: utau_ice !: u-stress over ice (I-point for LIM2 or U,V-point for LIM3) [N/m2] 44 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: vtau_ice !: v-stress over ice (I-point for LIM2 or U,V-point for LIM3) [N/m2] 44 45 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fr1_i0 !: 1st fraction of sol. rad. which penetrate inside the ice cover 45 46 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fr2_i0 !: 2nd fraction of sol. rad. which penetrate inside the ice cover … … 51 52 52 53 #else 53 54 54 !!---------------------------------------------------------------------- 55 55 !! Default option NO LIM 2.0 or 3.0 sea-ice model … … 58 58 LOGICAL , PUBLIC, PARAMETER :: lk_lim3 = .FALSE. !: no LIM-3 ice model 59 59 CHARACTER(len=1), PUBLIC :: cigr_type = '-' !: no grid ice-velocity 60 61 60 #endif 62 61 63 62 !!---------------------------------------------------------------------- 64 !! OPA 9.0 , LOCEAN-IPSL (2005)63 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 65 64 !! $Id$ 66 65 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 67 66 !!---------------------------------------------------------------------- 67 68 !!====================================================================== 68 69 END MODULE sbc_ice -
trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
r1468 r1469 6 6 !! History : 2.0 ! 06-2007 (R. Redler, N. Keenlyside, W. Park) Original code split into flxmod & taumod 7 7 !! 3.0 ! 02-2008 (G. Madec, C Talandier) surface module 8 !! - ! 08-2008 (S. Masson, E. ....) generic coupled interface8 !! 3.1 ! 02-2009 (S. Masson, E. Maisonave, A. Caubel) generic coupled interface 9 9 !!---------------------------------------------------------------------- 10 10 #if defined key_oasis3 || defined key_oasis4 … … 1094 1094 DO jj = 2, jpjm1 1095 1095 DO ji = fs_2, fs_jpim1 ! vector opt. 1096 zotx1(ji,jj) = 0.5 * ( un (ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj)1097 zoty1(ji,jj) = 0.5 * ( vn (ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj)1098 zitx1(ji,jj) = 0.5 * ( utau i_ice(ji,jj ) + utaui_ice(ji-1,jj ) ) * fr_i(ji,jj)1099 zity1(ji,jj) = 0.5 * ( vtau i_ice(ji,jj ) + vtaui_ice(ji ,jj-1 ) ) * fr_i(ji,jj)1096 zotx1(ji,jj) = 0.5 * ( un (ji,jj,1) + un (ji-1,jj ,1) ) * zfr_l(ji,jj) 1097 zoty1(ji,jj) = 0.5 * ( vn (ji,jj,1) + vn (ji ,jj-1,1) ) * zfr_l(ji,jj) 1098 zitx1(ji,jj) = 0.5 * ( utau_ice(ji,jj ) + utau_ice(ji-1,jj ) ) * fr_i(ji,jj) 1099 zity1(ji,jj) = 0.5 * ( vtau_ice(ji,jj ) + vtau_ice(ji ,jj-1 ) ) * fr_i(ji,jj) 1100 1100 END DO 1101 1101 END DO … … 1103 1103 DO jj = 2, jpjm1 1104 1104 DO ji = fs_2, fs_jpim1 ! vector opt. 1105 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1)) * zfr_l(ji,jj)1106 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1)) * zfr_l(ji,jj)1107 zitx1(ji,jj) = 0.25 * ( utau i_ice(ji+1,jj+1) + utaui_ice(ji,jj+1) &1108 & + utau i_ice(ji+1,jj ) + utaui_ice(ji,jj ) ) * fr_i(ji,jj)1109 zity1(ji,jj) = 0.25 * ( vtau i_ice(ji+1,jj+1) + vtaui_ice(ji,jj+1) &1110 & + vtau i_ice(ji+1,jj ) + vtaui_ice(ji,jj ) ) * fr_i(ji,jj)1105 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1) ) * zfr_l(ji,jj) 1106 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) 1107 zitx1(ji,jj) = 0.25 * ( utau_ice(ji+1,jj+1) + utau_ice(ji,jj+1) & 1108 & + utau_ice(ji+1,jj ) + utau_ice(ji,jj ) ) * fr_i(ji,jj) 1109 zity1(ji,jj) = 0.25 * ( vtau_ice(ji+1,jj+1) + vtau_ice(ji,jj+1) & 1110 & + vtau_ice(ji+1,jj ) + vtau_ice(ji,jj ) ) * fr_i(ji,jj) 1111 1111 END DO 1112 1112 END DO … … 1114 1114 DO jj = 2, jpjm1 1115 1115 DO ji = fs_2, fs_jpim1 ! vector opt. 1116 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1)) * zfr_l(ji,jj)1117 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1)) * zfr_l(ji,jj)1118 zitx1(ji,jj) = 0.25 * ( utau i_ice(ji-1,jj-1) + utaui_ice(ji,jj-1) &1119 & + utau i_ice(ji-1,jj ) + utaui_ice(ji,jj ) ) * fr_i(ji,jj)1120 zity1(ji,jj) = 0.25 * ( vtau i_ice(ji-1,jj-1) + vtaui_ice(ji,jj-1) &1121 & + vtau i_ice(ji-1,jj ) + vtaui_ice(ji,jj ) ) * fr_i(ji,jj)1116 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1) ) * zfr_l(ji,jj) 1117 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) 1118 zitx1(ji,jj) = 0.25 * ( utau_ice(ji-1,jj-1) + utau_ice(ji,jj-1) & 1119 & + utau_ice(ji-1,jj ) + utau_ice(ji,jj ) ) * fr_i(ji,jj) 1120 zity1(ji,jj) = 0.25 * ( vtau_ice(ji-1,jj-1) + vtau_ice(ji,jj-1) & 1121 & + vtau_ice(ji-1,jj ) + vtau_ice(ji,jj ) ) * fr_i(ji,jj) 1122 1122 END DO 1123 1123 END DO … … 1129 1129 DO jj = 2, jpjm1 1130 1130 DO ji = fs_2, fs_jpim1 ! vector opt. 1131 zotx1(ji,jj) = 0.5 * ( un (ji,jj,1) + un(ji-1,jj ,1) ) * zfr_l(ji,jj) &1132 & + 0.5 * ( utau i_ice(ji,jj ) + utaui_ice(ji-1,jj ) ) * fr_i(ji,jj)1133 zoty1(ji,jj) = 0.5 * ( vn (ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) &1134 & + 0.5 * ( vtau i_ice(ji,jj ) + vtaui_ice(ji ,jj-1 ) ) * fr_i(ji,jj)1131 zotx1(ji,jj) = 0.5 * ( un (ji,jj,1) + un (ji-1,jj ,1) ) * zfr_l(ji,jj) & 1132 & + 0.5 * ( utau_ice(ji,jj ) + utau_ice(ji-1,jj ) ) * fr_i(ji,jj) 1133 zoty1(ji,jj) = 0.5 * ( vn (ji,jj,1) + vn (ji ,jj-1,1) ) * zfr_l(ji,jj) & 1134 & + 0.5 * ( vtau_ice(ji,jj ) + vtau_ice(ji ,jj-1 ) ) * fr_i(ji,jj) 1135 1135 END DO 1136 1136 END DO … … 1138 1138 DO jj = 2, jpjm1 1139 1139 DO ji = fs_2, fs_jpim1 ! vector opt. 1140 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1)) * zfr_l(ji,jj) &1141 & + 0.25 * ( utau i_ice(ji+1,jj+1) + utaui_ice(ji,jj+1) &1142 & + utau i_ice(ji+1,jj ) + utaui_ice(ji,jj ) ) * fr_i(ji,jj)1143 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1)) * zfr_l(ji,jj) &1144 & + 0.25 * ( vtau i_ice(ji+1,jj+1) + vtaui_ice(ji,jj+1) &1145 & + vtau i_ice(ji+1,jj ) + vtaui_ice(ji,jj ) ) * fr_i(ji,jj)1140 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1) ) * zfr_l(ji,jj) & 1141 & + 0.25 * ( utau_ice(ji+1,jj+1) + utau_ice(ji,jj+1) & 1142 & + utau_ice(ji+1,jj ) + utau_ice(ji,jj ) ) * fr_i(ji,jj) 1143 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) & 1144 & + 0.25 * ( vtau_ice(ji+1,jj+1) + vtau_ice(ji,jj+1) & 1145 & + vtau_ice(ji+1,jj ) + vtau_ice(ji,jj ) ) * fr_i(ji,jj) 1146 1146 END DO 1147 1147 END DO … … 1149 1149 DO jj = 2, jpjm1 1150 1150 DO ji = fs_2, fs_jpim1 ! vector opt. 1151 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1)) * zfr_l(ji,jj) &1152 & + 0.25 * ( utau i_ice(ji-1,jj-1) + utaui_ice(ji,jj-1) &1153 & + utau i_ice(ji-1,jj ) + utaui_ice(ji,jj ) ) * fr_i(ji,jj)1154 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1)) * zfr_l(ji,jj) &1155 & + 0.25 * ( vtau i_ice(ji-1,jj-1) + vtaui_ice(ji,jj-1) &1156 & + vtau i_ice(ji-1,jj ) + vtaui_ice(ji,jj ) ) * fr_i(ji,jj)1151 zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj-1,1) ) * zfr_l(ji,jj) & 1152 & + 0.25 * ( utau_ice(ji-1,jj-1) + utau_ice(ji,jj-1) & 1153 & + utau_ice(ji-1,jj ) + utau_ice(ji,jj ) ) * fr_i(ji,jj) 1154 zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji ,jj-1,1) ) * zfr_l(ji,jj) & 1155 & + 0.25 * ( vtau_ice(ji-1,jj-1) + vtau_ice(ji,jj-1) & 1156 & + vtau_ice(ji-1,jj ) + vtau_ice(ji,jj ) ) * fr_i(ji,jj) 1157 1157 END DO 1158 1158 END DO -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r1465 r1469 131 131 END DO 132 132 ! Bulk formulea - provides the following fields: 133 ! utau i_ice, vtaui_ice : surface ice stress (U- & V-points) [N/m2]134 ! qsr_ice , qns_ice: solar & non solar heat flux over ice (T-point) [W/m2]135 ! qla_ice 136 ! dqns_ice , dqla_ice: non solar & latent heat sensistivity (T-point) [W/m2]137 ! tprecip , sprecip: total & solid precipitation (T-point) [Kg/m2/s]138 ! fr1_i0 , fr2_i0: 1sr & 2nd fraction of qsr penetration in ice [%]133 ! utau_ice, vtau_ice : surface ice stress (U- & V-points) [N/m2] 134 ! qsr_ice , qns_ice : solar & non solar heat flux over ice (T-point) [W/m2] 135 ! qla_ice : latent heat flux over ice (T-point) [W/m2] 136 ! dqns_ice, dqla_ice : non solar & latent heat sensistivity (T-point) [W/m2] 137 ! tprecip , sprecip : total & solid precipitation (T-point) [Kg/m2/s] 138 ! fr1_i0 , fr2_i0 : 1sr & 2nd fraction of qsr penetration in ice [%] 139 139 ! 140 140 SELECT CASE( kblk ) 141 141 CASE( 3 ) ! CLIO bulk formulation 142 142 CALL blk_ice_clio( t_su , alb_ice_cs, alb_ice_os, & 143 & utau i_ice , vtaui_ice, qns_ice , qsr_ice , &143 & utau_ice , vtau_ice , qns_ice , qsr_ice , & 144 144 & qla_ice , dqns_ice , dqla_ice , & 145 145 & tprecip , sprecip , & … … 147 147 ! 148 148 CASE( 4 ) ! CORE bulk formulation 149 CALL blk_ice_core( t_su , u_ice , v_ice , alb_ice_cs, &150 & utau i_ice, vtaui_ice, qns_ice , qsr_ice , &151 & qla_ice , dqns_ice , dqla_ice , &152 & tprecip , sprecip , &153 & fr1_i0 , fr2_i0 , cl_grid, jpl )149 CALL blk_ice_core( t_su , u_ice , v_ice , alb_ice_cs, & 150 & utau_ice , vtau_ice , qns_ice , qsr_ice , & 151 & qla_ice , dqns_ice , dqla_ice , & 152 & tprecip , sprecip , & 153 & fr1_i0 , fr2_i0 , cl_grid, jpl ) 154 154 END SELECT 155 155 … … 327 327 WRITE(numout,*) ' oceanic stress utau : ', utau(ji,jj) 328 328 WRITE(numout,*) ' oceanic stress vtau : ', vtau(ji,jj) 329 WRITE(numout,*) ' sea-ice stress utau i_ice : ', utaui_ice(ji,jj)330 WRITE(numout,*) ' sea-ice stress vtau i_ice : ', vtaui_ice(ji,jj)329 WRITE(numout,*) ' sea-ice stress utau_ice : ', utau_ice(ji,jj) 330 WRITE(numout,*) ' sea-ice stress vtau_ice : ', vtau_ice(ji,jj) 331 331 WRITE(numout,*) ' oceanic speed u : ', u_oce(ji,jj) 332 332 WRITE(numout,*) ' oceanic speed v : ', v_oce(ji,jj) … … 631 631 WRITE(numout,*) ' - Stresses ' 632 632 WRITE(numout,*) ' ~~~~~~~~ ' 633 WRITE(numout,*) ' utau i_ice : ', utaui_ice(ki,kj)634 WRITE(numout,*) ' vtau i_ice : ', vtaui_ice(ki,kj)633 WRITE(numout,*) ' utau_ice : ', utau_ice(ki,kj) 634 WRITE(numout,*) ' vtau_ice : ', vtau_ice(ki,kj) 635 635 WRITE(numout,*) ' utau : ', utau(ki,kj) 636 636 WRITE(numout,*) ' vtau : ', vtau(ki,kj) -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r1468 r1469 127 127 128 128 ! ... Sea-ice surface boundary conditions output from bulk formulae : 129 ! - utau i_ice! surface ice stress i-component (I-point) [N/m2]130 ! - vtau i_ice! surface ice stress j-component (I-point) [N/m2]129 ! - utau_ice ! surface ice stress i-component (I-point) [N/m2] 130 ! - vtau_ice ! surface ice stress j-component (I-point) [N/m2] 131 131 ! - qns_ice ! non solar heat flux over ice (T-point) [W/m2] 132 132 ! - qsr_ice ! solar heat flux over ice (T-point) [W/m2] … … 141 141 SELECT CASE( ksbc ) 142 142 CASE( 3 ) ! CLIO bulk formulation 143 CALL blk_ice_clio( zsist, zalb_ice_cs, zalb_ice_os ,&144 & utau i_ice , vtaui_ice , qns_ice , qsr_ice,&145 & qla_ice , dqns_ice , dqla_ice ,&146 & tprecip , sprecip ,&147 & fr1_i0 , fr2_i0 , cl_grid, jpl )143 CALL blk_ice_clio( zsist, zalb_ice_cs, zalb_ice_os, & 144 & utau_ice , vtau_ice , qns_ice , qsr_ice, & 145 & qla_ice , dqns_ice , dqla_ice , & 146 & tprecip , sprecip , & 147 & fr1_i0 , fr2_i0 , cl_grid , jpl ) 148 148 149 149 CASE( 4 ) ! CORE bulk formulation 150 CALL blk_ice_core( zsist, ui_ice , vi_ice , zalb_ice_cs, &151 & utau i_ice , vtaui_ice , qns_ice, qsr_ice, &152 & qla_ice , dqns_ice , dqla_ice, &153 & tprecip , sprecip ,&154 & fr1_i0 , fr2_i0 , cl_grid, jpl )150 CALL blk_ice_core( zsist, ui_ice , vi_ice , zalb_ice_cs, & 151 & utau_ice , vtau_ice , qns_ice , qsr_ice, & 152 & qla_ice , dqns_ice , dqla_ice , & 153 & tprecip , sprecip , & 154 & fr1_i0 , fr2_i0 , cl_grid , jpl ) 155 155 CASE( 5 ) ! Coupled formulation : atmosphere-ice stress only (fluxes provided after ice dynamics) 156 CALL sbc_cpl_ice_tau( utau i_ice , vtaui_ice )156 CALL sbc_cpl_ice_tau( utau_ice , vtau_ice ) 157 157 END SELECT 158 158 159 159 IF(ln_ctl) THEN ! print mean trends (used for debugging) 160 160 CALL prt_ctl_info( 'Ice Forcings ' ) 161 CALL prt_ctl( tab2d_1=tprecip ,clinfo1=' sbc_ice_lim: precip : ', tab2d_2=sprecip , clinfo2=' Snow: ' )162 CALL prt_ctl( tab2d_1=utau i_ice,clinfo1=' sbc_ice_lim: utaui_ice: ', tab2d_2=vtaui_ice, clinfo2=' vtaui_ice: ' )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=ui_oce ,clinfo1=' sbc_ice_lim: u_io : ', tab2d_2=vi_oce , clinfo2=' v_io: ' )165 CALL prt_ctl( tab2d_1=hsnif ,clinfo1=' sbc_ice_lim: hsnif 1 : ', tab2d_2=hicif , clinfo2=' hicif: ' )166 CALL prt_ctl( tab2d_1=frld ,clinfo1=' sbc_ice_lim: frld 1 : ', tab2d_2=sist , clinfo2=' sist: ' )161 CALL prt_ctl( tab2d_1=tprecip ,clinfo1=' sbc_ice_lim: precip : ', tab2d_2=sprecip , clinfo2=' Snow : ' ) 162 CALL prt_ctl( tab2d_1=utau_ice,clinfo1=' sbc_ice_lim: utau_ice: ', tab2d_2=vtau_ice, clinfo2=' vtau_ice: ' ) 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=ui_oce ,clinfo1=' sbc_ice_lim: u_io : ', tab2d_2=vi_oce , clinfo2=' v_io : ' ) 165 CALL prt_ctl( tab2d_1=hsnif ,clinfo1=' sbc_ice_lim: hsnif 1: ', tab2d_2=hicif , clinfo2=' hicif : ' ) 166 CALL prt_ctl( tab2d_1=frld ,clinfo1=' sbc_ice_lim: frld 1: ', tab2d_2=sist , clinfo2=' sist : ' ) 167 167 ENDIF 168 168
Note: See TracChangeset
for help on using the changeset viewer.