Changeset 4270

Timestamp:

2013-11-20T10:01:12+01:00 (10 years ago)

Author:

clem

Message:

few rewritings limcat_1D

File:

• branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90 (modified) (5 diffs)

branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/LIM_SRC_3/limcat_1D.F90

@@ -51 +51 @@
       !--------------------------------------------------------------------
       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
 
       !------------------------------------------------------------------------------------
@@ -87 +87 @@
             
             ! 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
@@ -100 +100 @@
                ! 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
                
@@ -116 +116 @@
                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
@@ -146 +146 @@
             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