- Timestamp:
- 2011-03-09T16:02:46+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r2612 r2678 79 79 !! profile of the ice/snow layers : z_i, z_s 80 80 !! total ice/snow thickness : ht_i_b, ht_s_b 81 !! 82 !! ** External : 83 !! 84 !! ** References : 85 !! 86 !! ** History : 87 !! (02-2003) Martin Vancoppenolle, Louvain-la-Neuve, Belgium 88 !! (06-2005) Martin Vancoppenolle, 3d version 89 !! (11-2006) Vectorized by Xavier Fettweis (UCL-ASTR) 90 !! (04-2007) Energy conservation tested by M. Vancoppenolle 81 91 !!------------------------------------------------------------------ 82 92 INTEGER , INTENT (in) :: & … … 333 343 END DO 334 344 END DO 345 ENDIF 346 347 IF ( thcon_i_swi .EQ. 1 ) THEN ! Pringle (0.011/2=0.0055) 348 DO layer = 1, nlay_i-1 349 DO ji = kideb , kiut 350 ztcond_i(ji,layer) = rcdic + 0.09*( s_i_b(ji,layer) & 351 + s_i_b(ji,layer+1) ) / MIN(-2.0*zeps, & 352 t_i_b(ji,layer)+t_i_b(ji,layer+1)-2.0*rtt) - & 353 0.0055* ( t_i_b(ji,layer) + t_i_b(ji,layer+1) - 2.0*rtt ) 354 ztcond_i(ji,layer) = MAX(ztcond_i(ji,layer),zkimin) 355 END DO 356 END DO 357 ENDIF 358 359 IF ( thcon_i_swi .EQ. 0 ) THEN ! Untersteiner 335 360 DO ji = kideb , kiut 336 361 ztcond_i(ji,nlay_i) = rcdic + zbeta*s_i_b(ji,nlay_i) / &
Note: See TracChangeset
for help on using the changeset viewer.