- Timestamp:
- 2012-11-02T07:13:40+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/LIM_SRC_3/limthd_sal.F90
r3523 r3524 12 12 !! 'key_lim3' LIM-3 sea-ice model 13 13 !!---------------------------------------------------------------------- 14 !! lim_thd_sal : salinity variations in the ice 15 !!---------------------------------------------------------------------- 16 USE par_oce ! ocean parameters 17 USE phycst ! physical constants (ocean directory) 18 USE sbc_oce ! Surface boundary condition: ocean fields 19 USE ice ! LIM variables 20 USE par_ice ! LIM parameters 21 USE thd_ice ! LIM thermodynamics 22 USE limvar ! LIM variables 23 USE in_out_manager ! I/O manager 24 USE lib_mpp ! MPP library 25 USE wrk_nemo ! work arrays 14 !! lim_thd_sal : salinity variations in the ice 15 !!---------------------------------------------------------------------- 16 USE par_oce ! ocean parameters 17 USE phycst ! physical constants (ocean directory) 18 USE sbc_oce ! Surface boundary condition: ocean fields 19 USE ice ! LIM variables 20 USE par_ice ! LIM parameters 21 USE thd_ice ! LIM thermodynamics 22 USE limvar ! LIM variables 23 USE in_out_manager ! I/O manager 24 USE lib_mpp ! MPP library 25 USE wrk_nemo ! work arrays 26 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 26 27 27 28 IMPLICIT NONE … … 32 33 33 34 !!---------------------------------------------------------------------- 34 !! NEMO/LIM3 4.0, UCL - NEMO Consortium (2011)35 !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2011) 35 36 !! $Id$ 36 37 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 149 150 DO ji = kideb, kiut 150 151 i_ice_switch = 1._wp - MAX( 0._wp, SIGN( 1._wp , - ht_i_b(ji) ) ) 151 fsbri_1d(ji) = fsbri_1d(ji) - i_ice_switch * rhoic * a_i_b(ji) * ht_i_b(ji) &152 & * ( MAX( dsm_i_gd_1d(ji) + dsm_i_fl_1d(ji) , sm_i_b(ji) - zsiold(ji) ) ) * r1_rdtice153 END DO ! ji152 sfx_bri_1d(ji) = sfx_bri_1d(ji) - i_ice_switch * rhoic * a_i_b(ji) * ht_i_b(ji) & 153 & * ( MAX( dsm_i_gd_1d(ji) + dsm_i_fl_1d(ji) , sm_i_b(ji) - zsiold(ji) ) ) * r1_rdtice 154 END DO 154 155 155 156 ! Only necessary for conservation check since salinity is modified … … 183 184 ! note: s_i_new = bulk_sal in constant salinity case 184 185 DO ji = kideb, kiut 185 fseqv_1d(ji) = fseqv_1d(ji) - s_i_new(ji) * rhoic * a_i_b(ji) * MAX( dh_i_bott(ji) , 0._wp ) * r1_rdtice186 sfx_thd_1d(ji) = sfx_thd_1d(ji) - s_i_new(ji) * rhoic * a_i_b(ji) * MAX( dh_i_bott(ji) , 0._wp ) * r1_rdtice 186 187 END DO 187 188 !
Note: See TracChangeset
for help on using the changeset viewer.