Changeset 9929 for NEMO/trunk/src/ICE

Timestamp:

2018-07-11T15:59:30+02:00 (6 years ago)

Author:

clem

Message:

remove unused variables

Location:

NEMO/trunk/src/ICE

Files:

• icealb.F90 (modified) (1 diff)
• icedyn_adv_umx.F90 (modified) (3 diffs)
• icethd_zdf_bl99.F90 (modified) (1 diff)

NEMO/trunk/src/ICE/icealb.F90

@@ -149 +149 @@
                !
                !                       !--- Snow-covered ice albedo (freezing, melting cases)
-               IF( pt_su(ji,jj,jl) < rt0_snow ) THEN
+               IF( pt_su(ji,jj,jl) < rt0 ) THEN
                   zalb_snw = rn_alb_sdry - ( rn_alb_sdry - zalb_ice ) * EXP( - ph_snw(ji,jj,jl) * z1_c3 )
                ELSE

NEMO/trunk/src/ICE/icedyn_adv_umx.F90

@@ -136 +136 @@
       !
    END SUBROUTINE ice_dyn_adv_umx
+
    
    SUBROUTINE adv_umx( k_order, kt, pdt, puc, pvc, pubox, pvbox, ptc )
@@ -159 +160 @@
       INTEGER  ::   ji, jj           ! dummy loop indices  
       REAL(wp) ::   ztra             ! local scalar
-      REAL(wp) ::   zfp_ui, zfp_vj   !   -      -
-      REAL(wp) ::   zfm_ui, zfm_vj   !   -      -
       REAL(wp), DIMENSION(jpi,jpj) ::   zfu_ups, zfu_ho, zt_u, zt_ups
       REAL(wp), DIMENSION(jpi,jpj) ::   zfv_ups, zfv_ho, zt_v, ztrd
@@ -169 +168 @@
       DO jj = 1, jpjm1         ! upstream tracer flux in the i and j direction
          DO ji = 1, fs_jpim1   ! vector opt.
-            zfp_ui = puc(ji,jj) + ABS( puc(ji,jj) )
-            zfm_ui = puc(ji,jj) - ABS( puc(ji,jj) )
-            zfp_vj = pvc(ji,jj) + ABS( pvc(ji,jj) )
-            zfm_vj = pvc(ji,jj) - ABS( pvc(ji,jj) )
-            zfu_ups(ji,jj) = 0.5_wp * ( zfp_ui * ptc(ji,jj) + zfm_ui * ptc(ji+1,jj  ) )
-            zfv_ups(ji,jj) = 0.5_wp * ( zfp_vj * ptc(ji,jj) + zfm_vj * ptc(ji  ,jj+1) )
+            zfu_ups(ji,jj) = MAX( puc(ji,jj), 0._wp ) * ptc(ji,jj) + MIN( puc(ji,jj), 0._wp ) * ptc(ji+1,jj)
+            zfv_ups(ji,jj) = MAX( pvc(ji,jj), 0._wp ) * ptc(ji,jj) + MIN( pvc(ji,jj), 0._wp ) * ptc(ji,jj+1)
          END DO
       END DO

NEMO/trunk/src/ICE/icethd_zdf_bl99.F90

@@ -805 +805 @@
                cnd_ice_1d(ji) = 2._wp * zkappa_i(ji,0)
             ELSE
-               cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) / 0.1_wp
+               cnd_ice_1d(ji) = 2._wp * ztcond_i(ji,0) * 10._wp
             ENDIF
          ENDIF