Changeset 1756 for trunk/NEMO/LIM_SRC_2/limsbc_2.F90
- Timestamp:
- 2009-11-25T15:15:20+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_2/limsbc_2.F90
r1695 r1756 24 24 USE lbclnk ! ocean lateral boundary condition 25 25 USE in_out_manager ! I/O manager 26 USE diaar5, ONLY : lk_diaar5 26 27 USE iom ! 27 28 USE albedo ! albedo parameters … … 80 81 INTEGER :: iflt, ial, iadv, ifral, ifrdv 81 82 INTEGER :: ii0, ii1, ij0, ij1 ! temporary integers 83 REAL(wp) :: zrdtir ! 1. / rdt_ice 82 84 REAL(wp) :: zqsr , zqns ! solar & non solar heat flux 83 85 REAL(wp) :: zinda ! switch for testing the values of ice concentration … … 97 99 !!--------------------------------------------------------------------- 98 100 101 zrdtir = 1. / rdt_ice 102 99 103 IF( kt == nit000 ) THEN 100 104 IF(lwp) WRITE(numout,*) … … 180 184 zqns = - ( 1. - thcm(ji,jj) ) * zqsr & ! part of the solar energy used in leads 181 185 & + iflt * ( fscmbq(ji,jj) + ffltbif(ji,jj) ) & 182 & + ifral * ( ial * qcmif(ji,jj) + (1 - ial) * qldif(ji,jj) ) / rdt_ice&183 & + ifrdv * ( qfvbq(ji,jj) + qdtcn(ji,jj) ) / rdt_ice186 & + ifral * ( ial * qcmif(ji,jj) + (1 - ial) * qldif(ji,jj) ) * zrdtir & 187 & + ifrdv * ( qfvbq(ji,jj) + qdtcn(ji,jj) ) * zrdtir 184 188 185 189 fsbbq(ji,jj) = ( 1.0 - ( ifvt + iflt ) ) * fscmbq(ji,jj) ! ??? … … 190 194 END DO 191 195 196 CALL iom_put( 'hflx_ice_cea', - fdtcn(:,:) ) 192 197 CALL iom_put( 'qns_io_cea', qns(:,:) - zqnsoce(:,:) * pfrld(:,:) ) 193 198 CALL iom_put( 'qsr_io_cea', fstric(:,:) * (1. - pfrld(:,:)) ) … … 206 211 #if defined key_coupled 207 212 zemp = emp_tot(ji,jj) - emp_ice(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! 208 & + rdmsnif(ji,jj) / rdt_ice! freshwaterflux due to snow melting213 & + rdmsnif(ji,jj) * zrdtir ! freshwaterflux due to snow melting 209 214 #else 210 215 !!$ ! computing freshwater exchanges at the ice/ocean interface … … 217 222 & - tprecip(ji,jj) * ( 1. - frld(ji,jj) ) & ! liquid precipitation reaches directly the ocean 218 223 & + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! taking into account change in ice cover within the time step 219 & + rdmsnif(ji,jj) / rdt_ice! freshwaterflux due to snow melting224 & + rdmsnif(ji,jj) * zrdtir ! freshwaterflux due to snow melting 220 225 ! ! ice-covered fraction: 221 226 #endif 222 227 223 228 ! computing salt exchanges at the ice/ocean interface 224 zfons = ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) / rdt_ice)229 zfons = ( soce_r(ji,jj) - sice_r(ji,jj) ) * ( rdmicif(ji,jj) * zrdtir ) 225 230 226 231 ! converting the salt flux from ice to a freshwater flux from ocean … … 232 237 END DO 233 238 END DO 239 240 IF( lk_diaar5 ) THEN 241 CALL iom_put( 'isnwmlt_cea' , rdmsnif(:,:) * zrdtir ) 242 CALL iom_put( 'fsal_virt_cea', soce_r(:,:) * rdmicif(:,:) * zrdtir ) 243 CALL iom_put( 'fsal_real_cea', - sice_r(:,:) * rdmicif(:,:) * zrdtir ) 244 ENDIF 234 245 235 246 !------------------------------------------!
Note: See TracChangeset
for help on using the changeset viewer.