Changeset 4270
 Timestamp:
 20131120T10:01:12+01:00 (9 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 (hmaxhmin)/2 101 101 DO jl = 1, i_fill  1 102 zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl1) ) / 2.102 zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl1) ) * 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.