Changeset 4270
- Timestamp:
- 2013-11-20T10:01:12+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90
r4161 r4270 51 51 !-------------------------------------------------------------------- 52 52 ijpij = SIZE(zhti,1) 53 zht_i(1:ijpij,1:jpl) = 0. d054 zht_s(1:ijpij,1:jpl) = 0. d055 za_i (1:ijpij,1:jpl) = 0. d053 zht_i(1:ijpij,1:jpl) = 0._wp 54 zht_s(1:ijpij,1:jpl) = 0._wp 55 za_i (1:ijpij,1:jpl) = 0._wp 56 56 57 57 !------------------------------------------------------------------------------------ … … 87 87 88 88 ! initialisation of ice variables for each try 89 zht_i(ji,1:jpl) = 0. d090 za_i (ji,1:jpl) = 0. d089 zht_i(ji,1:jpl) = 0._wp 90 za_i (ji,1:jpl) = 0._wp 91 91 92 92 ! *** case very thin ice: fill only category 1 … … 100 100 ! Fill ice thicknesses except the last one (i_fill) by (hmax-hmin)/2 101 101 DO jl = 1, i_fill - 1 102 zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) / 2.102 zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) * 0.5_wp 103 103 END DO 104 104 … … 116 116 DO jl = 1, i_fill - 1 117 117 IF ( jl == jl0 ) CYCLE 118 zarg = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) / 2.)118 zarg = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) * 0.5_wp ) 119 119 za_i(ji,jl) = za_i (ji,jl0) * EXP(-zarg**2) 120 120 END DO … … 146 146 ztest_4 = 1 147 147 DO jl = 1, i_fill 148 IF ( za_i(ji,jl) < 0. 0d0) ztest_4 = 0148 IF ( za_i(ji,jl) < 0._wp ) ztest_4 = 0 149 149 END DO 150 150
Note: See TracChangeset
for help on using the changeset viewer.