Changeset 809

Timestamp:

2008-02-06T17:19:08+01:00 (16 years ago)

Author:

ctlod

Message:

problem in the computation of dqla_ice(:,:), fr1_i0(:,:) fr2_i0(:,:) fields in CORE bulk formulae, see ticket: #62

File:

• trunk/NEMO/OPA_SRC/SBC/flx_core.h90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/SBC/flx_core.h90

@@ -564 +564 @@
 
          ! d qla_ice/ d zsss
-         dqla_ice(:,:) = -Ls*Cice*0.98*11637800/(rhoa*rhoa) &
+         dqla_ice(:,:) = -Ls*Cice*11637800
               * (-5897.8)/(zsss(:,:)*zsss(:,:)) &
               * exp(-5897.8/zsss(:,:)) * dUnormt(:,:)
@@ -576 +576 @@
          ! FRACTION 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 )
-
-         catm1(:,:) = 1.0  - 0.3  ! flxnow(:,:,8)
-         fr1_i0(:,:) = & 
-              (0.18 * catm1(:,:) + 0.35 * flxnow(:,:,8))
-         fr2_i0(:,:) = &
-              (0.82 * catm1(:,:) + 0.65 * flxnow(:,:,8))
-
+         ! ( Maykut and Untersteiner, 1971 ; Ebert and Curry, 1993 )
+
+         catm1(:,:) = 1.0 - 0.3
+         fr1_i0(:,:) = 0.18 * ( 1.0 - catm1(:,:) ) + 0.35 * catm1(:,:)
+         fr2_i0(:,:) = 0.82 * ( 1.0 - catm1(:,:) ) + 0.65 * catm1(:,:)
 
          !  Total  PRECIPITATION (kg/m2/s)