- Timestamp:
- 2014-11-13T16:43:38+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_4728_CNRS04_coupled_interface/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r4734 r4855 177 177 ! --- Energy received in the lead, zqld is defined everywhere (J.m-2) --- ! 178 178 zqld = tms(ji,jj) * rdt_ice * & 179 & ( pfrld(ji,jj)* ( zqsr(ji,jj) * oatte(ji,jj) + zqns(ji,jj) ) &179 & ( pfrld(ji,jj) * ( zqsr(ji,jj) * oatte(ji,jj) + zqns(ji,jj) ) & 180 180 & + ( pfrld(ji,jj)**betas - pfrld(ji,jj) ) * sprecip(ji,jj) * & ! heat content of precip 181 181 & ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus ) & … … 212 212 hfx_in(ji,jj) = hfx_in(ji,jj) & 213 213 ! heat flux above the ocean 214 & + pfrld(ji,jj) * ( qns(ji,jj) + qsr(ji,jj) )&214 & + pfrld(ji,jj) * ( zqns(ji,jj) + zqsr(ji,jj) ) & 215 215 ! latent heat of precip (note that precip is included in qns but not in qns_ice) 216 216 & + ( 1._wp - pfrld(ji,jj) ) * sprecip(ji,jj) * ( cpic * ( MIN( tatm_ice(ji,jj), rt0_snow ) - rtt ) - lfus ) & … … 584 584 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in configuration namelist', lwp ) 585 585 IF(lwm) WRITE ( numoni, namicethd ) 586 587 IF( lk_cpl .AND. parsub /= 0.0 ) CALL ctl_stop( 'In coupled mode, use parsub = 0. or send dqla' ) 586 588 ! 587 589 IF(lwp) THEN ! control print
Note: See TracChangeset
for help on using the changeset viewer.