- Timestamp:
- 2014-05-27T11:28:12+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r4634 r4649 69 69 70 70 PUBLIC sbc_ice_lim ! routine called by sbcmod.F90 71 PUBLIC lim_prt_state 71 72 72 73 !! * Substitutions … … 170 171 ! 171 172 IF( ln_nicep ) THEN ! control print at a given point 172 jiindx = 3 ; jjindx = 49173 jiindx = 15 ; jjindx = 44 173 174 IF(lwp) WRITE(numout,*) ' The debugging point is : jiindx : ',jiindx, ' jjindx : ',jjindx 174 175 ENDIF … … 183 184 u_oce(:,:) = ssu_m(:,:) ! mean surface ocean current at ice velocity point 184 185 v_oce(:,:) = ssv_m(:,:) ! (C-grid dynamics : U- & V-points as the ocean) 185 ! 186 t_bo(:,:) = tfreez( sss_m ) + rt0 ! masked sea surface freezing temperature [Kelvin] 187 ! ! (set to rt0 over land) 186 187 ! masked sea surface freezing temperature [Kelvin] 188 t_bo(:,:) = ( tfreez( sss_m ) + rt0 ) * tmask(:,:,1) + rt0 * ( 1. - tmask(:,:,1) ) 189 188 190 CALL albedo_ice( t_su, ht_i, ht_s, zalb_ice_cs, zalb_ice_os ) ! ... ice albedo 189 191 … … 310 312 ! salt, heat and mass fluxes 311 313 sfx (:,:) = 0._wp ; 312 sfx_bri(:,:) = 0._wp ; sfx_dyn(:,:) = 0._wp ; sfx_res(:,:) = 0._wp314 sfx_bri(:,:) = 0._wp ; sfx_dyn(:,:) = 0._wp 313 315 sfx_sni(:,:) = 0._wp ; sfx_opw(:,:) = 0._wp 314 316 sfx_bog(:,:) = 0._wp ; sfx_dyn(:,:) = 0._wp 315 317 sfx_bom(:,:) = 0._wp ; sfx_sum(:,:) = 0._wp 316 317 hfx_thd(:,:) = 0._wp ; hfx_dyn(:,:) = 0._wp ; hfx_snw(:,:) = 0._wp 318 hfx_tot(:,:) = 0._wp ; hfx_spr(:,:) = 0._wp ; hfx_res(:,:) = 0._wp 319 hfx_sub(:,:) = 0._wp ; hfx_err(:,:) = 0._wp ; hfx_in (:,:) = 0._wp ; hfx_out(:,:) = 0._wp 320 hfx_err_rem(:,:) = 0._wp 321 322 wfx_snw(:,:) = 0._wp ; wfx_ice(:,:) = 0._wp ; wfx_sub(:,:) = 0._wp 318 sfx_res(:,:) = 0._wp 319 320 wfx_snw(:,:) = 0._wp ; wfx_ice(:,:) = 0._wp 323 321 wfx_sni(:,:) = 0._wp ; wfx_opw(:,:) = 0._wp 324 322 wfx_bog(:,:) = 0._wp ; wfx_dyn(:,:) = 0._wp 325 323 wfx_bom(:,:) = 0._wp ; wfx_sum(:,:) = 0._wp 326 wfx_res(:,:) = 0._wp ; 327 ! 328 fhld (:,:) = 0._wp 329 fmmflx(:,:) = 0._wp 330 ftr_ice(:,:,:) = 0._wp ! part of solar radiation transmitted through the ice 324 wfx_res(:,:) = 0._wp ; wfx_sub(:,:) = 0._wp 325 wfx_spr(:,:) = 0._wp ; 326 327 hfx_in (:,:) = 0._wp ; hfx_out(:,:) = 0._wp 328 hfx_thd(:,:) = 0._wp ; 329 hfx_snw(:,:) = 0._wp ; hfx_opw(:,:) = 0._wp 330 hfx_bog(:,:) = 0._wp ; hfx_dyn(:,:) = 0._wp 331 hfx_bom(:,:) = 0._wp ; hfx_sum(:,:) = 0._wp 332 hfx_res(:,:) = 0._wp ; hfx_sub(:,:) = 0._wp 333 hfx_spr(:,:) = 0._wp ; hfx_dif(:,:) = 0._wp 334 hfx_err(:,:) = 0._wp ; hfx_err_rem(:,:) = 0._wp 335 336 ! 337 fhld (:,:) = 0._wp 338 fmmflx(:,:) = 0._wp 339 ! part of solar radiation transmitted through the ice 340 ftr_ice(:,:,:) = 0._wp 331 341 332 342 ! diags 333 diag_trp_vi (:,:) = 0._wp ; diag_trp_vs(:,:) = 0._wp ; diag_trp_ei(:,:) = 0._wp ; diag_trp_es(:,:) = 0._wp ;334 diag_heat_dhc 1(:,:) = 0._wp ;343 diag_trp_vi (:,:) = 0._wp ; diag_trp_vs(:,:) = 0._wp ; diag_trp_ei(:,:) = 0._wp ; diag_trp_es(:,:) = 0._wp 344 diag_heat_dhc(:,:) = 0._wp 335 345 336 346 ! dynamical invariants … … 820 830 WRITE(numout,*) ' hfx_in : ', hfx_in(ji,jj) 821 831 WRITE(numout,*) ' hfx_out : ', hfx_out(ji,jj) 822 WRITE(numout,*) ' hfx_tot : ', hfx_tot(ji,jj) 823 WRITE(numout,*) ' dhc : ', diag_heat_dhc1(ji,jj) 832 WRITE(numout,*) ' dhc : ', diag_heat_dhc(ji,jj) 824 833 WRITE(numout,*) 825 834 WRITE(numout,*) ' hfx_dyn : ', hfx_dyn(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.