# Changeset 4270 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3

Ignore:
Timestamp:
2013-11-20T10:01:12+01:00 (8 years ago)
Message:

few rewritings limcat_1D

File:
1 edited

Unmodified
Added
Removed
• ## branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90

 r4161 !-------------------------------------------------------------------- ijpij = SIZE(zhti,1) zht_i(1:ijpij,1:jpl) = 0.d0 zht_s(1:ijpij,1:jpl) = 0.d0 za_i (1:ijpij,1:jpl) = 0.d0 zht_i(1:ijpij,1:jpl) = 0._wp zht_s(1:ijpij,1:jpl) = 0._wp za_i (1:ijpij,1:jpl) = 0._wp !------------------------------------------------------------------------------------ ! initialisation of ice variables for each try zht_i(ji,1:jpl) = 0.d0 za_i (ji,1:jpl) = 0.d0 zht_i(ji,1:jpl) = 0._wp za_i (ji,1:jpl) = 0._wp ! *** case very thin ice: fill only category 1 ! Fill ice thicknesses except the last one (i_fill) by (hmax-hmin)/2 DO jl = 1, i_fill - 1 zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) / 2. zht_i(ji,jl) = ( hi_max(jl) + hi_max(jl-1) ) * 0.5_wp END DO DO jl = 1, i_fill - 1 IF ( jl == jl0 ) CYCLE zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) / 2. ) zarg           = ( zht_i(ji,jl) - zhti(ji) ) / ( zhti(ji) * 0.5_wp ) za_i(ji,jl) =   za_i (ji,jl0) * EXP(-zarg**2) END DO ztest_4 = 1 DO jl = 1, i_fill IF ( za_i(ji,jl) < 0.0d0 ) ztest_4 = 0 IF ( za_i(ji,jl) < 0._wp ) ztest_4 = 0 END DO
Note: See TracChangeset for help on using the changeset viewer.