Changeset 10531 for NEMO/trunk/src/ICE
- Timestamp:
- 2019-01-16T14:17:50+01:00 (5 years ago)
- Location:
- NEMO/trunk/src/ICE
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/ice.F90
r10425 r10531 114 114 INTEGER , PUBLIC :: nlay_i !: number of ice layers 115 115 INTEGER , PUBLIC :: nlay_s !: number of snow layers 116 INTEGER , PUBLIC :: nn_virtual_itd !: virtual ITD mono-category parameterizations (1-3) or not (0)116 LOGICAL , PUBLIC :: ln_virtual_itd !: virtual ITD mono-category parameterization (T) or not (F) 117 117 LOGICAL , PUBLIC :: ln_icedyn !: flag for ice dynamics (T) or not (F) 118 118 LOGICAL , PUBLIC :: ln_icethd !: flag for ice thermo (T) or not (F) -
NEMO/trunk/src/ICE/icedyn_rdgrft.F90
r10511 r10531 586 586 & + rhos * vsrft(ji) * ( 1._wp - rn_fsnwrft ) ) * r1_rdtice 587 587 588 ! Put the melt pond water into the ocean589 ! clem: I think the following lines must be commented since there590 ! is no net mass flux between melt ponds and the ocean (see icethd_pnd.F90 for ex.)591 !IF ( ln_pnd_fwb ) THEN592 ! wfx_pnd_1d(ji) = wfx_pnd_1d(ji) + ( rhow * vprdg(ji) * ( 1._wp - rn_fpndrdg ) & ! fresh water source for ocean593 ! & + rhow * vprft(ji) * ( 1._wp - rn_fpndrft ) ) * r1_rdtice594 !ENDIF595 596 588 ! virtual salt flux to keep salinity constant 597 589 IF( nn_icesal /= 2 ) THEN -
NEMO/trunk/src/ICE/iceforcing.F90
r10425 r10531 145 145 SELECT CASE( nice_jules ) 146 146 CASE( np_jules_ACTIVE ) ! compute conduction flux and surface temperature (as in Jules surface module) 147 CALL blk_ice_qcn ( nn_virtual_itd, t_su, t_bo, h_s, h_i )147 CALL blk_ice_qcn ( ln_virtual_itd, t_su, t_bo, h_s, h_i ) 148 148 END SELECT 149 149 CASE ( jp_purecpl ) !--- coupled formulation -
NEMO/trunk/src/ICE/icestp.F90
r10519 r10531 297 297 INTEGER :: ios ! Local integer 298 298 !! 299 NAMELIST/nampar/ jpl, nlay_i, nlay_s, nn_virtual_itd, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s, &299 NAMELIST/nampar/ jpl, nlay_i, nlay_s, ln_virtual_itd, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s, & 300 300 & cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir 301 301 !!------------------------------------------------------------------- … … 317 317 WRITE(numout,*) ' number of ice layers nlay_i = ', nlay_i 318 318 WRITE(numout,*) ' number of snow layers nlay_s = ', nlay_s 319 WRITE(numout,*) ' virtual ITD param for jpl=1 ( 1-3) or not (0) nn_virtual_itd = ', nn_virtual_itd319 WRITE(numout,*) ' virtual ITD param for jpl=1 (T) or not (F) ln_virtual_itd = ', ln_virtual_itd 320 320 WRITE(numout,*) ' Ice dynamics (T) or not (F) ln_icedyn = ', ln_icedyn 321 321 WRITE(numout,*) ' Ice thermodynamics (T) or not (F) ln_icethd = ', ln_icethd … … 324 324 ENDIF 325 325 ! !--- check consistency 326 IF ( jpl > 1 .AND. nn_virtual_itd == 1) THEN327 nn_virtual_itd = 0326 IF ( jpl > 1 .AND. ln_virtual_itd ) THEN 327 ln_virtual_itd = .FALSE. 328 328 IF(lwp) WRITE(numout,*) 329 IF(lwp) WRITE(numout,*) ' nn_virtual_itd forced to 0as jpl>1, no need with multiple categories to emulate them'329 IF(lwp) WRITE(numout,*) ' ln_virtual_itd forced to false as jpl>1, no need with multiple categories to emulate them' 330 330 ENDIF 331 331 ! -
NEMO/trunk/src/ICE/icethd.F90
r10511 r10531 233 233 CALL ice_thd_temp ! --- temperature update --- ! 234 234 ! 235 !!gm please create a new logical (l_thd_mono or a better explicit name) set one for all in icestp.F90 module 236 !!gm l_thd_mono = ln_icedH .AND. ( ( nn_virtual_itd == 1 .OR. nn_virtual_itd == 4 ) .AND. jpl == 1 ) 237 !!gm by the way, the different options associated with nn_virtual_itd =1 to 4 are quite impossible to identify 238 !!gm more comment to add when ready the namelist, with an explicit print in the ocean.output 239 IF( ln_icedH ) THEN 240 IF ( ( nn_virtual_itd == 1 .OR. nn_virtual_itd == 3 ) .AND. jpl == 1 ) THEN 241 CALL ice_thd_mono ! --- extra lateral melting if virtual_itd --- ! 242 END IF 243 END IF 235 IF( ln_icedH .AND. ln_virtual_itd ) & 236 & CALL ice_thd_mono ! --- extra lateral melting if virtual_itd --- ! 244 237 ! 245 238 IF( ln_icedA ) CALL ice_thd_da ! --- lateral melting --- ! -
NEMO/trunk/src/ICE/icethd_zdf_bl99.F90
r10425 r10531 265 265 zghe(1:npti) = 1._wp 266 266 ! 267 SELECT CASE ( nn_virtual_itd ) 268 ! 269 CASE ( 1 , 2 ) 267 IF( ln_virtual_itd ) THEN 270 268 ! 271 269 zepsilon = 0.1_wp … … 277 275 END DO 278 276 ! 279 END SELECT277 ENDIF 280 278 ! 281 279 !-----------------
Note: See TracChangeset
for help on using the changeset viewer.