Changeset 8498 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceupdate.F90

Timestamp:

2017-09-05T19:53:41+02:00 (7 years ago)

Author:

clem

Message:

changes in style - part2 -

File:

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

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

@@ -86 +86 @@
       !!  
       !! ** Purpose :   Update the surface ocean boundary condition for heat 
-      !!              salt and mass over areas where sea-ice is non-zero
+      !!                salt and mass over areas where sea-ice is non-zero
       !!         
       !! ** Action  : - computes the heat and freshwater/salt fluxes
-      !!              at the ice-ocean interface.
+      !!                at the ice-ocean interface.
       !!              - Update the ocean sbc
       !!     
@@ -133 +133 @@
          DO ji = 1, jpi
 
-            !------------------------------------------!
-            !      heat flux at the ocean surface      !
-            !------------------------------------------!
             ! Solar heat flux reaching the ocean = zqsr (W.m-2) 
             !---------------------------------------------------
-            zqsr = qsr_tot(ji,jj)
-            DO jl = 1, jpl
-               zqsr = zqsr - a_i_b(ji,jj,jl) * (  qsr_ice(ji,jj,jl) - ftr_ice(ji,jj,jl) ) 
-            END DO
-!!gm  why not like almost everywhere else :
-!!gm        zqsr = qsr_tot(ji,jj) - SUM( a_i_b(ji,jj,:) * (  qsr_ice(ji,jj,:) - ftr_ice(ji,jj,:) )
+            zqsr = qsr_tot(ji,jj) - SUM( a_i_b(ji,jj,:) * ( qsr_ice(ji,jj,:) - ftr_ice(ji,jj,:) ) )
 
             ! Total heat flux reaching the ocean = hfx_out (W.m-2) 
@@ -204 +196 @@
       END DO
 
-      !-----------------------------------------------!
-      !   Storing the transmitted variables           !
-      !-----------------------------------------------!
+      ! Storing the transmitted variables
+      !----------------------------------
       fr_i  (:,:)   = at_i(:,:)             ! Sea-ice fraction            
       tn_ice(:,:,:) = t_su(:,:,:)           ! Ice surface temperature                      
 
-      !------------------------------------------------------------------------!
-      !    Snow/ice albedo (only if sent to coupler, useless in forced mode)   !
-      !------------------------------------------------------------------------!
+      ! Snow/ice albedo (only if sent to coupler, useless in forced mode)
+      !------------------------------------------------------------------
       CALL ice_alb( t_su, ht_i, ht_s, a_ip_frac, h_ip, ln_pnd_rad, zalb_cs, zalb_os ) ! cloud-sky and overcast-sky ice albedos
       !