- Timestamp:
- 2017-10-05T16:44:46+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90
r8564 r8597 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 … … 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 ; ENDIF 363 DO jl = 1, jpl 364 a_ip_frac(:,:,jl) = rn_apnd * zswitch(:,:) * zfac 365 h_ip (:,:,jl) = rn_hpnd * zswitch(:,:) * zfac 366 END DO 367 a_ip(:,:,:) = a_ip_frac(:,:,:) * a_i (:,:,:) 368 v_ip(:,:,:) = h_ip (:,:,:) * a_ip(:,:,:) 376 369 377 370 ELSE ! if ln_iceini=false
Note: See TracChangeset
for help on using the changeset viewer.