- Timestamp:
- 2017-10-18T19:14:32+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90
r8586 r8637 94 94 REAL(wp) :: ztmelts, zdh 95 95 INTEGER :: i_hemis, i_fill, jl0 96 REAL(wp) :: zarg, zV, zconv, zdv 96 REAL(wp) :: zarg, zV, zconv, zdv, zfac 97 97 INTEGER , DIMENSION(4) :: itest 98 98 REAL(wp), DIMENSION(jpi,jpj) :: z2d … … 314 314 315 315 ! for constant salinity in time 316 IF( nn_icesal == 1 .OR. nn_icesal == 3) THEN316 IF( nn_icesal /= 2 ) THEN 317 317 CALL ice_var_salprof 318 318 sv_i = s_i * v_i … … 358 358 tn_ice (:,:,:) = t_su (:,:,:) 359 359 360 ! MV MP 2016361 360 ! Melt pond volume and fraction 362 IF ( ln_pnd ) THEN 363 DO jl = 1, jpl 364 a_ip_frac(:,:,jl) = 0.2 * zswitch(:,:) 365 h_ip (:,:,jl) = 0.05 * zswitch(:,:) 366 a_ip(:,:,jl) = a_ip_frac(:,:,jl) * a_i (:,:,jl) 367 v_ip(:,:,jl) = h_ip (:,:,jl) * a_ip(:,:,jl) 368 END DO 369 ELSE 370 a_ip(:,:,:) = 0._wp 371 v_ip(:,:,:) = 0._wp 372 a_ip_frac(:,:,:) = 0._wp 373 h_ip (:,:,:) = 0._wp 374 ENDIF 375 ! END MV MP 2016 361 IF ( ln_pnd_CST .OR. ln_pnd_H12 ) THEN ; zfac = 1._wp 362 ELSE ; zfac = 0._wp 363 ENDIF 364 DO jl = 1, jpl 365 a_ip_frac(:,:,jl) = rn_apnd * zswitch(:,:) * zfac 366 h_ip (:,:,jl) = rn_hpnd * zswitch(:,:) * zfac 367 END DO 368 a_ip(:,:,:) = a_ip_frac(:,:,:) * a_i (:,:,:) 369 v_ip(:,:,:) = h_ip (:,:,:) * a_ip(:,:,:) 376 370 377 371 ELSE ! if ln_iceini=false
Note: See TracChangeset
for help on using the changeset viewer.