- Timestamp:
- 2018-06-04T15:58:26+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r9499 r9734 1895 1895 1896 1896 ! Heat content per unit mass of snow (J/kg) 1897 WHERE( SUM( a_i, dim=3 ) > 1.e-10 ) ; zcptsnw(:,:) = cpic * SUM( (t n_ice - rt0) * a_i, dim=3 ) / SUM( a_i, dim=3 )1897 WHERE( SUM( a_i, dim=3 ) > 1.e-10 ) ; zcptsnw(:,:) = cpic * SUM( (t1_ice - rt0) * a_i, dim=3 ) / SUM( a_i, dim=3 ) 1898 1898 ELSEWHERE ; zcptsnw(:,:) = zcptn(:,:) 1899 1899 ENDWHERE 1900 1900 ! Heat content per unit mass of rain (J/kg) 1901 zcptrain(:,:) = rcp * ( SUM( (t n_ice(:,:,:) - rt0) * a_i(:,:,:), dim=3 ) + sst_m(:,:) * ziceld(:,:) )1901 zcptrain(:,:) = rcp * ( SUM( (t1_ice(:,:,:) - rt0) * a_i(:,:,:), dim=3 ) + sst_m(:,:) * ziceld(:,:) ) 1902 1902 1903 1903 ! --- enthalpy of snow precip over ice in J/m3 (to be used in 1D-thermo) --- ! … … 2140 2140 SELECT CASE( sn_snd_temp%clcat ) 2141 2141 CASE( 'yes' ) 2142 ztmp3(:,:,1:jpl) = t n_ice(:,:,1:jpl)2142 ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) 2143 2143 CASE( 'no' ) 2144 2144 WHERE( SUM( a_i, dim=3 ) /= 0. ) 2145 ztmp3(:,:,1) = SUM( t n_ice * a_i, dim=3 ) / SUM( a_i, dim=3 )2145 ztmp3(:,:,1) = SUM( t1_ice * a_i, dim=3 ) / SUM( a_i, dim=3 ) 2146 2146 ELSEWHERE 2147 2147 ztmp3(:,:,1) = rt0 … … 2152 2152 SELECT CASE( sn_snd_temp%clcat ) 2153 2153 CASE( 'yes' ) 2154 ztmp3(:,:,1:jpl) = t n_ice(:,:,1:jpl) * a_i(:,:,1:jpl)2154 ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) * a_i(:,:,1:jpl) 2155 2155 CASE( 'no' ) 2156 2156 ztmp3(:,:,:) = 0.0 2157 2157 DO jl=1,jpl 2158 ztmp3(:,:,1) = ztmp3(:,:,1) + t n_ice(:,:,jl) * a_i(:,:,jl)2158 ztmp3(:,:,1) = ztmp3(:,:,1) + t1_ice(:,:,jl) * a_i(:,:,jl) 2159 2159 ENDDO 2160 2160 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%clcat' ) … … 2163 2163 SELECT CASE( sn_snd_temp%clcat ) 2164 2164 CASE( 'yes' ) 2165 ztmp3(:,:,1:jpl) = t n_ice(:,:,1:jpl) * a_i(:,:,1:jpl)2165 ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) * a_i(:,:,1:jpl) 2166 2166 CASE( 'no' ) 2167 2167 ztmp3(:,:,:) = 0.0 2168 2168 DO jl=1,jpl 2169 ztmp3(:,:,1) = ztmp3(:,:,1) + t n_ice(:,:,jl) * a_i(:,:,jl)2169 ztmp3(:,:,1) = ztmp3(:,:,1) + t1_ice(:,:,jl) * a_i(:,:,jl) 2170 2170 ENDDO 2171 2171 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%clcat' ) … … 2174 2174 ztmp1(:,:) = ( ztmp1(:,:) + rt0 ) * zfr_l(:,:) 2175 2175 DO jl=1,jpl 2176 ztmp1(:,:) = ztmp1(:,:) + t n_ice(:,:,jl) * a_i(:,:,jl)2176 ztmp1(:,:) = ztmp1(:,:) + t1_ice(:,:,jl) * a_i(:,:,jl) 2177 2177 ENDDO 2178 2178 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%cldes' )
Note: See TracChangeset
for help on using the changeset viewer.