Changeset 8500 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

Timestamp:

2017-09-06T12:01:34+02:00 (7 years ago)

Author:

clem

Message:

changes in style - part3 - move output into proper files and correct a bug in debug mode

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90 (modified) (3 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

@@ -116 +116 @@
          !
          ALLOCATE( z1_at_i(jpi,jpj) , z1_vt_i(jpi,jpj) )
-         WHERE( at_i(:,:) > epsi10 )   ;   z1_at_i(:,:) = 1._wp / at_i(:,:)
+         WHERE( at_i(:,:) > epsi20 )   ;   z1_at_i(:,:) = 1._wp / at_i(:,:)
          ELSEWHERE                     ;   z1_at_i(:,:) = 0._wp
          END WHERE
-         WHERE( vt_i(:,:) > epsi10 )   ;   z1_vt_i(:,:) = 1._wp / vt_i(:,:)
+         WHERE( vt_i(:,:) > epsi20 )   ;   z1_vt_i(:,:) = 1._wp / vt_i(:,:)
          ELSEWHERE                     ;   z1_vt_i(:,:) = 0._wp
          END WHERE
@@ -163 +163 @@
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   z1_a_i, z1_v_i
       !!------------------------------------------------------------------
-
-!!gm Question 1:  here use epsi20 , in ice_var_agg it is epsi10 which is used...  why ??
 
 !!gm Question 2:  It is possible to define existence of sea-ice in a common way between 
@@ -549 +547 @@
       ! to be sure that at_i is the sum of a_i(jl)
       at_i (:,:) = a_i(:,:,1)
-     vt_i (:,:) = v_i(:,:,1)
+      vt_i (:,:) = v_i(:,:,1)
       DO jl = 2, jpl
          at_i(:,:) = at_i(:,:) + a_i(:,:,jl)
-        vt_i(:,:) = vt_i(:,:) + v_i(:,:,jl)
+         vt_i(:,:) = vt_i(:,:) + v_i(:,:,jl)
       END DO