Changeset 5836 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90
- Timestamp:
- 2015-10-26T15:49:40+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90
r5183 r5836 185 185 zfs_b = glob_sum( ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) + & 186 186 & sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) & 187 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )187 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 188 188 189 189 ! water flux 190 190 zfw_b = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) + & 191 191 & wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:) & 192 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )192 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 193 193 194 194 ! heat flux 195 195 zft_b = glob_sum( ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:) & 196 196 & - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:) & 197 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv )198 199 zvi_b = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e1 2t * tmask(:,:,1) * zconv )200 201 zsmv_b = glob_sum( SUM( smv_i * rhoic , dim=3 ) * e1 2t * tmask(:,:,1) * zconv )197 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) 198 199 zvi_b = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 200 201 zsmv_b = glob_sum( SUM( smv_i * rhoic , dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 202 202 203 203 zei_b = glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) + & 204 204 & SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 ) & 205 ) * e1 2t * tmask(:,:,1) * zconv )205 ) * e1e2t * tmask(:,:,1) * zconv ) 206 206 207 207 ELSEIF( icount == 1 ) THEN … … 210 210 zfs = glob_sum( ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) + & 211 211 & sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) & 212 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zfs_b212 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfs_b 213 213 214 214 ! water flux 215 215 zfw = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) + & 216 216 & wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:) & 217 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zfw_b217 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfw_b 218 218 219 219 ! heat flux 220 220 zft = glob_sum( ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:) & 221 221 & - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:) & 222 & ) * e1 2t(:,:) * tmask(:,:,1) * zconv ) - zft_b222 & ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zft_b 223 223 224 224 ! outputs 225 225 zvi = ( ( glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) & 226 & * e1 2t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday226 & * e1e2t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday 227 227 228 228 zsmv = ( ( glob_sum( SUM( smv_i * rhoic , dim=3 ) & 229 & * e1 2t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday229 & * e1e2t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday 230 230 231 231 zei = glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) + & 232 232 & SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 ) & 233 & ) * e1 2t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft233 & ) * e1e2t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft 234 234 235 235 ! zvtrp and zetrp must be close to 0 if the advection scheme is conservative 236 zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e1 2t * tmask(:,:,1) * zconv ) * rday237 zetrp = glob_sum( ( diag_trp_ei + diag_trp_es ) * e1 2t * tmask(:,:,1) * zconv )236 zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e1e2t * tmask(:,:,1) * zconv ) * rday 237 zetrp = glob_sum( ( diag_trp_ei + diag_trp_es ) * e1e2t * tmask(:,:,1) * zconv ) 238 238 239 239 zvmin = glob_min( v_i ) … … 242 242 243 243 ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice) 244 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1 2t * zconv ) ! in 1.e9 m2244 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 245 245 zv_sill = zarea * 2.5e-5 246 246 zs_sill = zarea * 25.e-5 … … 286 286 #if ! defined key_bdy 287 287 ! heat flux 288 zhfx = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e1 2t * tmask(:,:,1) * zconv )288 zhfx = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e1e2t * tmask(:,:,1) * zconv ) 289 289 ! salt flux 290 zsfx = glob_sum( ( sfx + diag_smvi ) * e1 2t * tmask(:,:,1) * zconv ) * rday290 zsfx = glob_sum( ( sfx + diag_smvi ) * e1e2t * tmask(:,:,1) * zconv ) * rday 291 291 ! water flux 292 zvfx = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e1 2t * tmask(:,:,1) * zconv ) * rday292 zvfx = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e1e2t * tmask(:,:,1) * zconv ) * rday 293 293 294 294 ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice) 295 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1 2t * zconv ) ! in 1.e9 m2295 zarea = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 296 296 zv_sill = zarea * 2.5e-5 297 297 zs_sill = zarea * 25.e-5
Note: See TracChangeset
for help on using the changeset viewer.