Changeset 2387


Ignore:
Timestamp:
2010-11-15T16:18:00+01:00 (10 years ago)
Author:
smasson
Message:

ticket #670 lim3: small bug in latent heat of solid precipitation melting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limthd.F90

    r1758 r2387  
    8686      REAL(wp) ::   zfric_umax = 2e-02    ! upper bound for the friction velocity 
    8787      REAL(wp) ::   zinda, zindb, zthsnice, zfric_u    ! temporary scalar 
    88       REAL(wp) ::   zfnsol, zfontn, zfntlat, zpareff   !    -         - 
     88      REAL(wp) ::   zfntlat, zpareff                   !    -         - 
    8989      REAL(wp) ::   zeps, zareamin, zcoef 
    9090      REAL(wp), DIMENSION(jpi,jpj) ::   zqlbsbq   ! link with lead energy budget qldif 
     
    190190            !           !-- Lead heat budget (part 1, next one is in limthd_dh 
    191191            !           !-- qldif -- (or qldif_1d in 1d routines) 
    192             zfontn         = sprecip(ji,jj) * lfus              ! energy of melting 
    193             zfnsol         = qns(ji,jj)                         ! total non solar flux 
    194             qldif(ji,jj)   = tms(ji,jj) * ( qsr(ji,jj)                               & 
    195                &                               + zfnsol + fdtcn(ji,jj) - zfontn      & 
    196                &                               + ( 1.0 - zindb ) * fsbbq(ji,jj)  )   & 
    197                &                               * ( 1.0 - at_i(ji,jj) ) * rdt_ice     
    198  
     192            qldif(ji,jj) =  tms(ji,jj) * rdt_ice  * (                             &  
     193               &   pfrld(ji,jj)        * (  qsr(ji,jj)                            &   ! solar heat  
     194               &                            + qns(ji,jj)                          &   ! non solar heat  
     195               &                            + fdtcn(ji,jj)                        &   ! turbulent ice-ocean heat  
     196               &                            + fsbbq(ji,jj) * ( 1.0 - zindb )  )   &   ! residual heat from previous step  
     197               & - pfrld(ji,jj)**betas * sprecip(ji,jj) * lfus                    )   ! latent heat of sprecip melting 
     198             
    199199            ! Positive heat budget is used for bottom ablation 
    200200            zfntlat        = 1.0 - MAX( zzero , SIGN( zone ,  - qldif(ji,jj) ) ) 
Note: See TracChangeset for help on using the changeset viewer.