Changeset 874 for branches

Timestamp:

2008-04-01T17:40:36+02:00 (16 years ago)

Author:

ctlod

Message:

dev_001_SBC: Bound the latent heat flux (over ocean & sea-ice) in CORE bulks and mask properly the wind stress over the ocean, see ticket:#99

File:

• branches/dev_001_SBC/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (4 diffs)

branches/dev_001_SBC/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -276 +276 @@
          &                   Cd, Ch, Ce )
    
-      ! ...  umasked Momentum : utau, vtau at U- and V_points, resp.
-      !      Note the use of 2-tmask in order to umask the stress along coastlines
+      ! ... utau, vtau at U- and V_points, resp.
+      !     Note the use of 0.5*(2-umask) in order to unmask the stress along coastlines
       zwnd_i(:,:) = rhoa * zwind_speed_t(:,:) * Cd(:,:) * zwnd_i(:,:)
       zwnd_j(:,:) = rhoa * zwind_speed_t(:,:) * Cd(:,:) * zwnd_j(:,:)
       DO jj = 1, jpjm1
          DO ji = 1, fs_jpim1
-            utau(ji,jj) = 0.5 * ( zwnd_i(ji,jj) + zwnd_i(ji+1,jj  ) ) * ( 2. - tmask(ji,jj,1) )
-            vtau(ji,jj) = 0.5 * ( zwnd_j(ji,jj) + zwnd_j(ji  ,jj+1) ) * ( 2. - tmask(ji,jj,1) )
+            utau(ji,jj) = 0.5 * ( 2. - umask(ji,jj,1) ) * ( zwnd_i(ji,jj) + zwnd_i(ji+1,jj  ) )
+            vtau(ji,jj) = 0.5 * ( 2. - vmask(ji,jj,1) ) * ( zwnd_j(ji,jj) + zwnd_j(ji  ,jj+1) )
          END DO
       END DO
@@ -292 +292 @@
       ! -----------------------------
 !CDIR COLLAPSE
-      zevap(:,:) = rhoa    *Ce(:,:)*( zqsatw(:,:) - sf(jp_humi)%fnow(:,:) ) * zwind_speed_t(:,:)   ! Evaporation
-!CDIR COLLAPSE
-      zqsb (:,:) = rhoa*cpa*Ch(:,:)*( zst   (:,:) - sf(jp_tair)%fnow(:,:) ) * zwind_speed_t(:,:)   ! Sensible Heat
-!CDIR COLLAPSE
-      zqla (:,:) = Lv * zevap(:,:)                                                                 ! Latent Heat
-
+      zevap(:,:) = MAX( 0.e0, rhoa    *Ce(:,:)*( zqsatw(:,:) - sf(jp_humi)%fnow(:,:) ) * zwind_speed_t(:,:) )   ! Evaporation
+!CDIR COLLAPSE
+      zqsb (:,:) =            rhoa*cpa*Ch(:,:)*( zst   (:,:) - sf(jp_tair)%fnow(:,:) ) * zwind_speed_t(:,:)     ! Sensible Heat
+!CDIR COLLAPSE
+      zqla (:,:) = Lv * zevap(:,:)                                                              ! Latent Heat
        
       ! ----------------------------------------------------------------------------- !
@@ -305 +304 @@
 !CDIR COLLAPSE
       qns(:,:) = zqlw(:,:) - zqsb(:,:) - zqla(:,:)      ! Downward Non Solar flux
-       
 !CDIR COLLAPSE
       emp (:,:) = zevap(:,:) - sf(jp_prec)%fnow(:,:) * tmask(:,:,1)
@@ -428 +426 @@
 !CDIR NOVERRCHK
 !CDIR COLLAPSE
-      p_qla(:,:) = rhoa * Ls  * Cice * z_wnds_t(:,:)   &                                    ! Latent Heat
-         &       * (  11637800. * EXP( -5897.8 / pst(:,:) ) / rhoa - sf(jp_humi)%fnow(:,:)  )
+      p_qla(:,:) = MAX( 0.e0, rhoa * Ls  * Cice * z_wnds_t(:,:)   &                           ! Latent Heat
+         &                    * (  11637800. * EXP( -5897.8 / pst(:,:) ) / rhoa - sf(jp_humi)%fnow(:,:)  ) )
       
       ! Latent heat sensitivity for ice (Dqla/Dt)