Changeset 9567 for branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

Timestamp:

2018-05-10T16:26:43+02:00 (6 years ago)

Author:

mathiot

Message:

apply tmask(:,:,1) to erp, qrp, precip, snow, rnf, qns to fix #2083

File:

• branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90 (modified) (2 diffs)

branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

@@ -461 +461 @@
       !!----------------------------------------------------------------------
 
-            esre  = shms(ji,jj) - shnow(ji,jj)   ! --- calculates the term : qs(Ta)-qa(t2d)
-
-            cseep = ce(ji,jj) * wndm(ji,jj) * esre     ! --- calculates the term : Ce*|V|*[qs(Ts)-qa(t2d)]
-
-            evap(ji,jj) = (cseep * rhom(ji,jj))  ! in [kg/m2/sec] !! --- calculates the EVAPORATION RATE [m/yr]
-
-            elat(ji,jj) = rhom(ji,jj) * cseep * heatlat(sst(ji,jj))
+            esre  = shms(ji,jj) - shnow(ji,jj) * tmask(ji,jj,1)   ! --- calculates the term : qs(Ta)-qa(t2d)
+
+            cseep = ce(ji,jj) * wndm(ji,jj) * esre * tmask(ji,jj,1)     ! --- calculates the term : Ce*|V|*[qs(Ts)-qa(t2d)]
+
+            evap(ji,jj) = (cseep * rhom(ji,jj)) * tmask(ji,jj,1)  ! in [kg/m2/sec] !! --- calculates the EVAPORATION RATE [m/yr]
+
+            elat(ji,jj) = rhom(ji,jj) * cseep * heatlat(sst(ji,jj)) * tmask(ji,jj,1)
 
       !!----------------------------------------------------------------------
@@ -635 +635 @@
 
           IF(cldnow(ji,jj).LT.0.3) THEN
-             qsw(ji,jj) = qtot * (1.-albedo)
+             qsw(ji,jj) = qtot * (1.-albedo) * tmask(ji,jj,1)
           ELSE
              qsw(ji,jj) = qtot*(1.-0.62*cldnow(ji,jj)              &
-                                + .0019*sunbet)*(1.-albedo)
+                                + .0019*sunbet)*(1.-albedo) * tmask(ji,jj,1)
           ENDIF