Changeset 1103


Ignore:
Timestamp:
2008-06-10T15:11:32+02:00 (12 years ago)
Author:
rblod
Message:

Fixes in kil3, change results, see ticket #206

Location:
trunk/NEMO/LIM_SRC_3
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limdyn.F90

    r921 r1103  
    115115                  i_j1 = i_j1 + 1 
    116116               END DO 
    117                i_j1 = MAX( 1, i_j1-1 ) 
     117               i_j1 = MAX( 1, i_j1-2 ) 
    118118               IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : NH  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj ) 
    119119               CALL lim_rhg( i_j1, i_jpj ) 
     
    125125                  i_jpj = i_jpj - 1 
    126126               END DO 
    127                i_jpj = MIN( jpj, i_jpj+2 ) 
     127               i_jpj = MIN( jpj, i_jpj+1 ) 
    128128               IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : SH  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj ) 
    129129 
    130                CALL lim_rhg( i_j1, i_jpj ) 
    131  
    132             ELSE                                 ! local domain extends over one hemisphere only 
    133                !                                 ! Rheology is computed only over the ice cover 
    134                !                                 ! latitude strip 
    135                i_j1  = 1 
     130       CALL lim_rhg( i_j1, i_jpj ) 
     131 
     132    ELSE                                 ! local domain extends over one hemisphere only 
     133       !                                 ! Rheology is computed only over the ice cover 
     134       !                                 ! latitude strip 
     135       i_j1  = 1 
    136136               DO WHILE ( i_j1 <= jpj .AND. zind(i_j1) == FLOAT(jpi) .AND. zmsk(i_j1) /=0 ) 
    137137                  i_j1 = i_j1 + 1 
    138138               END DO 
    139                i_j1 = MAX( 1, i_j1-1 ) 
     139               i_j1 = MAX( 1, i_j1-2 ) 
    140140 
    141141               i_jpj  = jpj 
     
    143143                  i_jpj = i_jpj - 1 
    144144               END DO 
    145                i_jpj = MIN( jpj, i_jpj+2) 
     145               i_jpj = MIN( jpj, i_jpj+1) 
    146146 
    147147               IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : one hemisphere:  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj ) 
  • trunk/NEMO/LIM_SRC_3/limistate.F90

    r1055 r1103  
    290290                     t_s(ji,jj,jk,jl) = zidto(ji,jj) * 270.00 + ( 1.0 - zidto(ji,jj) ) * rtt 
    291291                     ! Snow energy of melting 
    292                      e_s(ji,jk,jk,jl) = zidto(ji,jj) * rhosn * ( cpic * ( rtt - t_s(ji,jj,jk,jl) ) + lfus ) 
     292                     e_s(ji,jj,jk,jl) = zidto(ji,jj) * rhosn * ( cpic * ( rtt - t_s(ji,jj,jk,jl) ) + lfus ) 
    293293                     ! Change dimensions 
    294294                     e_s(ji,jj,jk,jl) = e_s(ji,jj,jk,jl) / unit_fac 
  • trunk/NEMO/LIM_SRC_3/limthd_dif.F90

    r1055 r1103  
    112112         maxnumeqmax 
    113113 
    114       INTEGER , DIMENSION(jpij) :: & 
     114      INTEGER , DIMENSION(kiut) :: & 
    115115         numeqmin, &   ! reference number of top equation 
    116116         numeqmax, &   ! reference number of bottom equation 
     
    118118 
    119119      !! * New local variables        
    120       REAL(wp) , DIMENSION(jpij,0:nlay_i) ::    & 
     120      REAL(wp) , DIMENSION(kiut,0:nlay_i) ::    & 
    121121         ztcond_i,    & !Ice thermal conductivity 
    122122         zradtr_i,    & !Radiation transmitted through the ice 
     
    124124         zkappa_i       !Kappa factor in the ice 
    125125 
    126       REAL(wp) , DIMENSION(jpij,0:nlay_s) ::    & 
     126      REAL(wp) , DIMENSION(kiut,0:nlay_s) ::    & 
    127127         zradtr_s,    & !Radiation transmited through the snow 
    128128         zradab_s,    & !Radiation absorbed in the snow 
    129129         zkappa_s       !Kappa factor in the snow 
    130130 
    131       REAL(wp) , DIMENSION(jpij,0:nlay_i) :: & 
     131      REAL(wp) , DIMENSION(kiut,0:nlay_i) :: & 
    132132         ztiold,      & !Old temperature in the ice 
    133133         zeta_i,      & !Eta factor in the ice  
     
    136136         z_i            !Vertical cotes of the layers in the ice 
    137137 
    138       REAL(wp) , DIMENSION(jpij,0:nlay_s) :: & 
     138      REAL(wp) , DIMENSION(kiut,0:nlay_s) :: & 
    139139         zeta_s,      & !Eta factor in the snow 
    140140         ztstemp,     & !Temporary temperature in the snow to check the convergence 
     
    142142         z_s            !Vertical cotes of the layers in the snow 
    143143 
    144       REAL(wp) , DIMENSION(jpij,jkmax+2) ::    & 
     144      REAL(wp) , DIMENSION(kiut,jkmax+2) ::    & 
    145145         zindterm,    & ! Independent term 
    146146         zindtbis,    & ! temporary independent term 
    147147         zdiagbis 
    148148 
    149       REAL(wp) , DIMENSION(jpij,jkmax+2,3) ::  & 
     149      REAL(wp) , DIMENSION(kiut,jkmax+2,3) ::  & 
    150150         ztrid          ! tridiagonal system terms 
    151151 
    152       REAL(wp), DIMENSION(jpij) ::  & 
     152      REAL(wp), DIMENSION(kiut) ::  & 
    153153         ztfs     ,   & ! ice melting point 
    154154         ztsuold  ,   & ! old surface temperature (before the iterative 
     
    175175         zerritmax              ! current maximal error on temperature  
    176176 
    177       REAL(wp), DIMENSION(jpij)  :: & 
     177      REAL(wp), DIMENSION(kiut)  :: & 
    178178         zerrit,             &  ! current error on temperature  
    179179         zdifcase,           &  ! case of the equation resolution (1->4) 
Note: See TracChangeset for help on using the changeset viewer.