Changeset 810

Timestamp:

2008-02-06T17:36:03+01:00 (16 years ago)

Author:

ctlod

Message:

dev_001_SBC: problem in the computation of p_dqla(:,:), prf1(:,:) p_fr2(:,:) fields in CORE bulk formulae, see ticket:#63

File:

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

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

@@ -349 +349 @@
       INTEGER  ::   ji, jj                    ! dummy loop indices
       REAL(wp) ::   zst3
-      REAL(wp) ::   zcoef_wnorm, zcoef_dqlw, zcoef_dqla, zcoef_dqsb, zcoef_fr12
+      REAL(wp) ::   zcoef_wnorm, zcoef_dqlw, zcoef_dqla, zcoef_dqsb
+      REAL(wp) ::   zcoef_frca                       ! fractional cloud amount
       REAL(wp) ::   zwnorm_f, zwndi_f , zwndj_f      ! relative wind module and components at F-point
       REAL(wp) ::             zwndi_t , zwndj_t      ! relative wind components at T-point
@@ -362 +363 @@
       zcoef_wnorm = rhoa * Cice
       zcoef_dqlw = 4.0 * 0.95 * Stef
-      zcoef_dqla = -Ls * Cice * 0.98 * 11637800. / (rhoa*rhoa) * (-5897.8)
+      zcoef_dqla = -Ls * Cice * 11637800. * (-5897.8)
       zcoef_dqsb = rhoa * cpa * Cice
-      zcoef_fr12 = 1.0  - 0.3       !!!gm ??? ! sf(jp_snow)%fnow(:,:)
+      zcoef_frca = 1.0  - 0.3
 
 !!gm brutal....
@@ -452 +453 @@
       ! FRACTIONs of net shortwave radiation which is not absorbed in the
       ! thin surface layer and penetrates inside the ice cover
-      ! ( Maykut and Untersteiner, 1971 ; Elbert and Curry, 1993 )
+      ! ( Maykut and Untersteiner, 1971 ; Ebert and Curry, 1993 )
     
 !CDIR COLLAPSE
-      p_fr1(:,:) = ( 0.18 * zcoef_fr12 + 0.35 * sf(jp_snow)%fnow(:,:) )
-!CDIR COLLAPSE
-      p_fr2(:,:) = ( 0.82 * zcoef_fr12 + 0.65 * sf(jp_snow)%fnow(:,:) )
+      p_fr1(:,:) = ( 0.18 * ( 1.0 - zcoef_frca ) + 0.35 * zcoef_frca )
+!CDIR COLLAPSE
+      p_fr2(:,:) = ( 0.82 * ( 1.0 - zcoef_frca ) + 0.65 * zcoef_frca )
        
 !CDIR COLLAPSE