- Timestamp:
- 2010-05-04T10:39:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limsbc_2.F90
r1857 r1858 4 4 !! computation of the flux at the sea ice/ocean interface 5 5 !!====================================================================== 6 !! History : 00-01 (H. Goosse) Original code 7 !! 02-07 (C. Ethe, G. Madec) re-writing F90 8 !! 06-07 (G. Madec) surface module 6 !! History : LIM ! 2000-01 (H. Goosse) Original code 7 !! 2.0 ! 2002-07 (C. Ethe, G. Madec) re-writing F90 8 !! - ! 2006-07 (G. Madec) surface module 9 !! 2.1 ! 2010-05 (Y. Aksenov, M. Vancoppenolle, G. Madec) add heat content exchanges 9 10 !!---------------------------------------------------------------------- 10 11 #if defined key_lim2 11 12 !!---------------------------------------------------------------------- 12 13 !! 'key_lim2' LIM 2.0 sea-ice model 13 !!----------------------------------------------------------------------14 14 !!---------------------------------------------------------------------- 15 15 !! lim_sbc_2 : flux at the ice / ocean interface … … 211 211 #if defined key_coupled 212 212 zemp = emp_tot(ji,jj) - emp_ice(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! 213 & + rdm snif(ji,jj) * zrdtir ! freshwaterflux due to snow melting213 & + rdm_snw(ji,jj) * zrdtir ! freshwaterflux due to snow melting 214 214 #else 215 215 !!$ ! computing freshwater exchanges at the ice/ocean interface … … 217 217 !!$ & + tprecip(ji,jj) & ! total precipitation 218 218 !!$ & - sprecip(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! remov. snow precip over ice 219 !!$ & - rdm snif(ji,jj) / rdt_ice ! freshwaterflux due to snow melting219 !!$ & - rdm_snw(ji,jj) / rdt_ice ! freshwaterflux due to snow melting 220 220 ! computing freshwater exchanges at the ice/ocean interface 221 221 zemp = + emp(ji,jj) * frld(ji,jj) & ! e-p budget over open ocean fraction 222 222 & - tprecip(ji,jj) * ( 1. - frld(ji,jj) ) & ! liquid precipitation reaches directly the ocean 223 223 & + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! taking into account change in ice cover within the time step 224 & + rdm snif(ji,jj) * zrdtir ! freshwaterflux due to snow melting224 & + rdm_snw(ji,jj) * zrdtir ! freshwaterflux due to snow melting 225 225 ! ! ice-covered fraction: 226 226 #endif 227 227 228 228 ! computing salt exchanges at the ice/ocean interface 229 zfons = ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdm icif(ji,jj) * zrdtir )229 zfons = ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdm_ice(ji,jj) * zrdtir ) 230 230 231 231 ! converting the salt flux from ice to a freshwater flux from ocean … … 239 239 240 240 IF( lk_diaar5 ) THEN 241 CALL iom_put( 'isnwmlt_cea' , rdm snif(:,:) * zrdtir )242 CALL iom_put( 'fsal_virt_cea', soce_r(:,:) * rdm icif(:,:) * zrdtir )243 CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdm icif(:,:) * zrdtir )241 CALL iom_put( 'isnwmlt_cea' , rdm_snw(:,:) * zrdtir ) 242 CALL iom_put( 'fsal_virt_cea', soce_r(:,:) * rdm_ice(:,:) * zrdtir ) 243 CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdm_ice(:,:) * zrdtir ) 244 244 ENDIF 245 245
Note: See TracChangeset
for help on using the changeset viewer.