- Timestamp:
- 2017-09-08T17:09:25+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_lac.F90
r8486 r8514 26 26 USE icethd_ent ! sea-ice: thermodynamics, enthalpy 27 27 USE icevar ! sea-ice: operations 28 USE icethd_sal ! sea-ice: salinity profiles 28 29 ! 29 30 USE in_out_manager ! I/O manager … … 112 113 !!-----------------------------------------------------------------------! 113 114 114 IF( ln_ limdiachk ) CALL ice_cons_hsm(0, 'icethd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)115 IF( ln_icediachk ) CALL ice_cons_hsm(0, 'icethd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 115 116 116 117 CALL ice_var_agg(1) … … 133 134 134 135 ! Default new ice thickness 135 WHERE( qlead(:,:) < 0._wp ) ; hicol(:,:) = rn_h newice136 WHERE( qlead(:,:) < 0._wp ) ; hicol(:,:) = rn_hinew 136 137 ELSEWHERE ; hicol(:,:) = 0._wp 137 138 END WHERE … … 348 349 DO ji = 1, nidx 349 350 rswitch = 1._wp - MAX( 0._wp, SIGN( 1._wp , - zat_i_1d(ji) ) ) 350 zfrazb = rswitch * ( TANH( rn_Cfraz b * ( zvrel_1d(ji) - rn_vfrazb ) ) + 1.0 ) * 0.5 * rn_maxfrazb351 zfrazb = rswitch * ( TANH( rn_Cfraz * ( zvrel_1d(ji) - rn_vfraz ) ) + 1.0 ) * 0.5 * rn_maxfraz 351 352 zv_frazb(ji) = zfrazb * zv_newice(ji) 352 353 zv_newice(ji) = ( 1.0 - zfrazb ) * zv_newice(ji) … … 478 479 ENDIF ! nidx > 0 479 480 ! 480 IF( ln_ limdiachk ) CALL ice_cons_hsm(1, 'icethd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)481 IF( ln_icediachk ) CALL ice_cons_hsm(1, 'icethd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 481 482 ! 482 483 END SUBROUTINE ice_thd_lac
Note: See TracChangeset
for help on using the changeset viewer.