Changeset 1270 for trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

Timestamp:

2009-01-14T19:35:02+01:00 (15 years ago)

Author:

rblod

Message:

Light bug in blk_ice, see ticket #289

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

@@ -30 +30 @@
    USE prtctl          ! Print control
 #if defined key_lim3
-   USE par_ice
    USE ice
    USE ice_oce         ! For ice surface temperature
 #elif defined key_lim2
-   USE par_ice_2
    USE ice_2
 #endif
@@ -414 +412 @@
       &                      p_qla , p_dqns, p_dqla,          &
       &                      p_tpr , p_spr ,                  &
-      &                      p_fr1 , p_fr2 , cd_grid  )
+      &                      p_fr1 , p_fr2 , cd_grid, pdim  )
       !!---------------------------------------------------------------------------
       !!                     ***  ROUTINE blk_ice_clio  ***
@@ -452 +450 @@
       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj) ::   p_fr2    ! 2nd fraction of qsr penetration in ice         [%]
       CHARACTER(len=1), INTENT(in   )             ::   cd_grid  ! type of sea-ice grid ("C" or "B" grid)
+      INTEGER, INTENT(in   )                      ::   pdim     ! number of ice categories
       !!
       INTEGER  ::   ji, jj, jl    ! dummy loop indices
@@ -467 +466 @@
       REAL(wp), DIMENSION(jpi,jpj) ::   zevsqr  ! vapour pressure square-root
       REAL(wp), DIMENSION(jpi,jpj) ::   zrhoa   ! air density
-      REAL(wp), DIMENSION(jpi,jpj,SIZE(pst,3)) ::   z_qlw, z_qsb
+      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_qlw, z_qsb
       !!---------------------------------------------------------------------
 
-      ijpl  = SIZE( pst, 3 )                 ! number of ice categories
+      ijpl  = pdim                           ! number of ice categories
       zpatm = 101000.                        ! atmospheric pressure  (assumed constant  here)