Changeset 13472 for NEMO/trunk/src/ICE/iceupdate.F90
- Timestamp:
- 2020-09-16T15:05:19+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/iceupdate.F90
r13295 r13472 25 25 USE icectl ! sea-ice: control prints 26 26 USE bdy_oce , ONLY : ln_bdy 27 USE zdfdrg , ONLY : ln_drgice_imp 27 28 ! 28 29 USE in_out_manager ! I/O manager … … 93 94 REAL(wp) :: zqmass ! Heat flux associated with mass exchange ice->ocean (W.m-2) 94 95 REAL(wp) :: zqsr ! New solar flux received by the ocean 95 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace 96 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalb_cs, zalb_os ! 3D workspace 96 REAL(wp), DIMENSION(jpi,jpj) :: z2d ! 2D workspace 97 97 !!--------------------------------------------------------------------- 98 98 IF( ln_timing ) CALL timing_start('ice_update') … … 182 182 ! Snow/ice albedo (only if sent to coupler, useless in forced mode) 183 183 !------------------------------------------------------------------ 184 CALL ice_alb( t_su, h_i, h_s, ln_pnd_alb, a_ip_frac, h_ip, zalb_cs, zalb_os ) ! cloud-sky and overcast-sky ice albedos 185 ! 186 alb_ice(:,:,:) = ( 1._wp - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 184 CALL ice_alb( t_su, h_i, h_s, ln_pnd_alb, a_ip_eff, h_ip, cloud_fra, alb_ice ) ! ice albedo 185 187 186 ! 188 187 IF( lrst_ice ) THEN !* write snwice_mass fields in the restart file … … 320 319 REAL(wp) :: zat_u, zutau_ice, zu_t, zmodt ! local scalar 321 320 REAL(wp) :: zat_v, zvtau_ice, zv_t, zrhoco ! - - 321 REAL(wp) :: zflagi ! - - 322 322 !!--------------------------------------------------------------------- 323 323 IF( ln_timing ) CALL timing_start('ice_update_tau') … … 350 350 ! 351 351 ! !== every ocean time-step ==! 352 IF ( ln_drgice_imp ) THEN 353 ! Save drag with right sign to update top drag in the ocean implicit friction 354 rCdU_ice(:,:) = -r1_rho0 * tmod_io(:,:) * at_i(:,:) * tmask(:,:,1) 355 zflagi = 0._wp 356 ELSE 357 zflagi = 1._wp 358 ENDIF 352 359 ! 353 360 DO_2D( 0, 0, 0, 0 )
Note: See TracChangeset
for help on using the changeset viewer.