Changeset 13480 for NEMO/releases/r4.0/r4.0-HEAD/tests
- Timestamp:
- 2020-09-16T19:03:32+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/tests/BENCH/MY_SRC/usrdef_sbc.F90
r10179 r13480 126 126 REAL(wp), DIMENSION(:,:,:), INTENT(in) :: phi ! ice thickness 127 127 !! 128 REAL(wp) :: zfr1, zfr2 ! local variables129 128 REAL(wp), DIMENSION(jpi,jpj) :: zsnw ! snw distribution after wind blowing 130 129 !!--------------------------------------------------------------------- … … 147 146 ! ice fields deduced from above 148 147 zsnw(:,:) = 1._wp 149 !!CALL lim_thd_snwblow( at_i_b, zsnw ) ! snow distribution over ice after150 !wind blowing151 148 emp_ice (:,:) = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw(:,:) 152 149 emp_oce (:,:) = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) … … 161 158 qsr_tot (:,:) = at_i_b(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 162 159 163 ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 164 zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice ) ! transmission when hi>10cm 165 zfr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) ! zfr2 such that zfr1 + zfr2 to equal 1 166 ! 167 WHERE ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) < 0.1_wp ) ! linear decrease from hi=0 to 10cm 168 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 169 ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp ) ! constant (zfr1) when hi>10cm 170 qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 171 ELSEWHERE ! zero when hs>0 172 qtr_ice_top(:,:,:) = 0._wp 173 END WHERE 160 ! --- shortwave radiation transmitted below the surface (W/m2) 161 qtr_ice_top(:,:,:) = 0._wp 174 162 #endif 175 163
Note: See TracChangeset
for help on using the changeset viewer.