Changeset 9119 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
- Timestamp:
- 2017-12-18T17:33:07+01:00 (7 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icethd_zdf.F90
r9019 r9119 119 119 ENDIF 120 120 ! 121 IF ( ( ln_cndi_U64 .AND. ln_cndi_P07 ) .OR. ( .NOT. ln_cndi_U64 .AND. .NOT.ln_cndi_P07 ) ) THEN122 CALL ctl_stop( 'ice_thd_zdf_init: choose one and only oneformulation for thermal conduction (ln_cndi_U64 or ln_cndi_P07)' )121 IF ( ( ln_cndi_U64 .AND. ln_cndi_P07 ) .OR. ( .NOT. ln_cndi_U64 .AND. .NOT. ln_cndi_P07 ) ) THEN 122 CALL ctl_stop( 'ice_thd_zdf_init: choose 1 and only 1 formulation for thermal conduction (ln_cndi_U64 or ln_cndi_P07)' ) 123 123 ENDIF 124 124 ! !== set the choice of ice vertical thermodynamic formulation ==! -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90
r9118 r9119 236 236 DO jk = 1, nlay_s 237 237 WHERE( v_s(:,:,:) > epsi20 ) !--- icy area 238 t_s(:,:,jk,:) = MAX( -100._wp , MIN( r1_cpic * ( -r1_rhosn * (e_s(:,:,jk,:)/v_s(:,:,:)*zlay_s) + lfus ) , 0._wp ) ) + rt0 238 t_s(:,:,jk,:) = MAX( -100._wp , & 239 & MIN( r1_cpic * ( -r1_rhosn * ( e_s(:,:,jk,:) / v_s(:,:,:) * zlay_s ) + lfus ) , & 240 & 0._wp ) ) & 241 & + rt0 239 242 ELSEWHERE !--- no ice 240 243 t_s(:,:,jk,:) = rt0
Note: See TracChangeset
for help on using the changeset viewer.