Changeset 1103 for trunk/NEMO
- Timestamp:
- 2008-06-10T15:11:32+02:00 (16 years ago)
- Location:
- trunk/NEMO/LIM_SRC_3
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limdyn.F90
r921 r1103 115 115 i_j1 = i_j1 + 1 116 116 END DO 117 i_j1 = MAX( 1, i_j1- 1)117 i_j1 = MAX( 1, i_j1-2 ) 118 118 IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn : NH i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj ) 119 119 CALL lim_rhg( i_j1, i_jpj ) … … 125 125 i_jpj = i_jpj - 1 126 126 END DO 127 i_jpj = MIN( jpj, i_jpj+ 2)127 i_jpj = MIN( jpj, i_jpj+1 ) 128 128 IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn : SH i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj ) 129 129 130 131 132 133 134 135 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 136 136 DO WHILE ( i_j1 <= jpj .AND. zind(i_j1) == FLOAT(jpi) .AND. zmsk(i_j1) /=0 ) 137 137 i_j1 = i_j1 + 1 138 138 END DO 139 i_j1 = MAX( 1, i_j1- 1)139 i_j1 = MAX( 1, i_j1-2 ) 140 140 141 141 i_jpj = jpj … … 143 143 i_jpj = i_jpj - 1 144 144 END DO 145 i_jpj = MIN( jpj, i_jpj+ 2)145 i_jpj = MIN( jpj, i_jpj+1) 146 146 147 147 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 290 290 t_s(ji,jj,jk,jl) = zidto(ji,jj) * 270.00 + ( 1.0 - zidto(ji,jj) ) * rtt 291 291 ! Snow energy of melting 292 e_s(ji,j k,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 ) 293 293 ! Change dimensions 294 294 e_s(ji,jj,jk,jl) = e_s(ji,jj,jk,jl) / unit_fac -
trunk/NEMO/LIM_SRC_3/limthd_dif.F90
r1055 r1103 112 112 maxnumeqmax 113 113 114 INTEGER , DIMENSION( jpij) :: &114 INTEGER , DIMENSION(kiut) :: & 115 115 numeqmin, & ! reference number of top equation 116 116 numeqmax, & ! reference number of bottom equation … … 118 118 119 119 !! * New local variables 120 REAL(wp) , DIMENSION( jpij,0:nlay_i) :: &120 REAL(wp) , DIMENSION(kiut,0:nlay_i) :: & 121 121 ztcond_i, & !Ice thermal conductivity 122 122 zradtr_i, & !Radiation transmitted through the ice … … 124 124 zkappa_i !Kappa factor in the ice 125 125 126 REAL(wp) , DIMENSION( jpij,0:nlay_s) :: &126 REAL(wp) , DIMENSION(kiut,0:nlay_s) :: & 127 127 zradtr_s, & !Radiation transmited through the snow 128 128 zradab_s, & !Radiation absorbed in the snow 129 129 zkappa_s !Kappa factor in the snow 130 130 131 REAL(wp) , DIMENSION( jpij,0:nlay_i) :: &131 REAL(wp) , DIMENSION(kiut,0:nlay_i) :: & 132 132 ztiold, & !Old temperature in the ice 133 133 zeta_i, & !Eta factor in the ice … … 136 136 z_i !Vertical cotes of the layers in the ice 137 137 138 REAL(wp) , DIMENSION( jpij,0:nlay_s) :: &138 REAL(wp) , DIMENSION(kiut,0:nlay_s) :: & 139 139 zeta_s, & !Eta factor in the snow 140 140 ztstemp, & !Temporary temperature in the snow to check the convergence … … 142 142 z_s !Vertical cotes of the layers in the snow 143 143 144 REAL(wp) , DIMENSION( jpij,jkmax+2) :: &144 REAL(wp) , DIMENSION(kiut,jkmax+2) :: & 145 145 zindterm, & ! Independent term 146 146 zindtbis, & ! temporary independent term 147 147 zdiagbis 148 148 149 REAL(wp) , DIMENSION( jpij,jkmax+2,3) :: &149 REAL(wp) , DIMENSION(kiut,jkmax+2,3) :: & 150 150 ztrid ! tridiagonal system terms 151 151 152 REAL(wp), DIMENSION( jpij) :: &152 REAL(wp), DIMENSION(kiut) :: & 153 153 ztfs , & ! ice melting point 154 154 ztsuold , & ! old surface temperature (before the iterative … … 175 175 zerritmax ! current maximal error on temperature 176 176 177 REAL(wp), DIMENSION( jpij) :: &177 REAL(wp), DIMENSION(kiut) :: & 178 178 zerrit, & ! current error on temperature 179 179 zdifcase, & ! case of the equation resolution (1->4)
Note: See TracChangeset
for help on using the changeset viewer.