Changeset 2387 for trunk

Timestamp:

2010-11-15T16:18:00+01:00 (13 years ago)

Author:

smasson

Message:

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

File:

• trunk/NEMO/LIM_SRC_3/limthd.F90 (modified) (2 diffs)

trunk/NEMO/LIM_SRC_3/limthd.F90

@@ -86 +86 @@
       REAL(wp) ::   zfric_umax = 2e-02    ! upper bound for the friction velocity
       REAL(wp) ::   zinda, zindb, zthsnice, zfric_u    ! temporary scalar
-      REAL(wp) ::   zfnsol, zfontn, zfntlat, zpareff   !    -         -
+      REAL(wp) ::   zfntlat, zpareff                   !    -         -
       REAL(wp) ::   zeps, zareamin, zcoef
       REAL(wp), DIMENSION(jpi,jpj) ::   zqlbsbq   ! link with lead energy budget qldif
@@ -190 +190 @@
             !           !-- Lead heat budget (part 1, next one is in limthd_dh
             !           !-- qldif -- (or qldif_1d in 1d routines)
-            zfontn         = sprecip(ji,jj) * lfus              ! energy of melting
-            zfnsol         = qns(ji,jj)                         ! total non solar flux
-            qldif(ji,jj)   = tms(ji,jj) * ( qsr(ji,jj)                               &
-               &                               + zfnsol + fdtcn(ji,jj) - zfontn      &
-               &                               + ( 1.0 - zindb ) * fsbbq(ji,jj)  )   &
-               &                               * ( 1.0 - at_i(ji,jj) ) * rdt_ice    
-
+            qldif(ji,jj) =  tms(ji,jj) * rdt_ice  * (                             & 
+               &   pfrld(ji,jj)        * (  qsr(ji,jj)                            &   ! solar heat 
+               &                            + qns(ji,jj)                          &   ! non solar heat 
+               &                            + fdtcn(ji,jj)                        &   ! turbulent ice-ocean heat 
+               &                            + fsbbq(ji,jj) * ( 1.0 - zindb )  )   &   ! residual heat from previous step 
+               & - pfrld(ji,jj)**betas * sprecip(ji,jj) * lfus                    )   ! latent heat of sprecip melting
+            
             ! Positive heat budget is used for bottom ablation
             zfntlat        = 1.0 - MAX( zzero , SIGN( zone ,  - qldif(ji,jj) ) )