Changeset 8971

Timestamp:

2017-12-11T11:43:16+01:00 (6 years ago)

Author:

clem

Message:

debug for conservation when Jules option=2 (not the reference). This commit is on behalf of Martin

File:

• branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90 (modified) (3 diffs)

branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90

@@ -1003 +1003 @@
       INTEGER  ::   iter                 ! local integer
       REAL(wp) ::   zfac, zfac2, zfac3   ! local scalars
-      REAL(wp) ::   zkeff_h, ztsu        !
+      REAL(wp) ::   zkeff_h, ztsu, ztsu0 !
       REAL(wp) ::   zqc, zqnet           !
       REAL(wp) ::   zhe, zqa0            !
@@ -1050 +1050 @@
                zkeff_h = zfac * zgfac(ji,jj,jl) / &                                    ! Effective conductivity of the snow-ice system divided by thickness
                   &      ( rcdic * phs(ji,jj,jl) + rn_cnd_s * MAX( 0.01, phi(ji,jj,jl) ) )
-               ztsu    = ptsu(ji,jj,jl)                                                ! Store initial surface temperature
+               ztsu    = ptsu(ji,jj,jl)                                                ! Store current iteration temperature
+               ztsu0   = ptsu(ji,jj,jl)                                                ! Store initial surface temperature
                zqa0    = qsr_ice(ji,jj,jl) - qsr_ice_tr(ji,jj,jl) + qns_ice(ji,jj,jl)  ! Net initial atmospheric heat flux
                !
@@ -1061 +1062 @@
                ptsu   (ji,jj,jl) = MIN( rt0, ztsu )
                qcn_ice(ji,jj,jl) = zkeff_h * ( ptsu(ji,jj,jl) - ptb(ji,jj) )
-               !
+               qns_ice(ji,jj,jl) = qns_ice(ji,jj,jl) + dqns_ice(ji,jj,jl) * ( ptsu(ji,jj,jl) - ztsu0 )
+               qml_ice(ji,jj,jl) = ( qsr_ice(ji,jj,jl) - qsr_ice_tr(ji,jj,jl) + qns_ice(ji,jj,jl) - qcn_ice(ji,jj,jl) )   &
+                             &   * MAX( 0._wp , SIGN( 1._wp, ptsu(ji,jj,jl) - rt0 ) )
+
             END DO
          END DO