- Timestamp:
- 2015-02-06T19:12:57+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r5064 r5067 26 26 USE phycst ! physical constants 27 27 USE dom_oce ! ocean domain 28 USE dom_ice, ONLY : tms29 28 USE ice ! LIM sea-ice variables 30 29 USE sbc_ice ! Surface boundary condition: sea-ice fields … … 170 169 zemp = emp(ji,jj) * pfrld(ji,jj) & ! evaporation over oceanic fraction 171 170 & - tprecip(ji,jj) * ( 1._wp - pfrld(ji,jj) ) & ! all precipitation reach the ocean 172 & + sprecip(ji,jj) * ( 1._wp - pfrld(ji,jj)** betas ) ! except solid precip intercepted by sea-ice171 & + sprecip(ji,jj) * ( 1._wp - pfrld(ji,jj)**rn_betas ) ! except solid precip intercepted by sea-ice 173 172 ENDIF 174 173 … … 197 196 snwice_mass_b(:,:) = snwice_mass(:,:) 198 197 ! new mass per unit area 199 snwice_mass (:,:) = tm s(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:) )198 snwice_mass (:,:) = tmask(:,:,1) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:) ) 200 199 ! time evolution of snow+ice mass 201 200 snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) * r1_rdtice … … 346 345 ! ! embedded sea ice 347 346 IF( nn_ice_embd /= 0 ) THEN ! mass exchanges between ice and ocean (case 1 or 2) set the snow+ice mass 348 snwice_mass (:,:) = tm s(:,:) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:) )347 snwice_mass (:,:) = tmask(:,:,1) * ( rhosn * vt_s(:,:) + rhoic * vt_i(:,:) ) 349 348 snwice_mass_b(:,:) = snwice_mass(:,:) 350 349 ELSE
Note: See TracChangeset
for help on using the changeset viewer.