Changeset 8158 for branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

Timestamp:

2017-06-09T09:09:35+02:00 (7 years ago)

Author:

vancop

Message:

SIMIP outputs, phase 2, commit#5

File:

• branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90 (modified) (3 diffs)

branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

@@ -282 +282 @@
             &                              ) * a_i_1d(ji) * r1_rdtice
          ! Mass flux by sublimation
-         wfx_sub_1d(ji) =  wfx_sub_1d(ji) - rhosn * a_i_1d(ji) * zdh_s_sub(ji) * r1_rdtice
+         wfx_snw_sub_1d(ji) = wfx_snw_sub_1d(ji) - rhosn * a_i_1d(ji) * zdh_s_sub(ji) * r1_rdtice
+
          ! new snow thickness
          ht_s_1d(ji)    =  MAX( 0._wp , ht_s_1d(ji) + zdh_s_sub(ji) )
@@ -383 +384 @@
             hfx_sub_1d(ji) = hfx_sub_1d(ji) + zdum * q_i_1d(ji,jk) * a_i_1d(ji) * r1_rdtice
             ! Mass flux > 0
-            wfx_sub_1d(ji) =  wfx_sub_1d(ji) - rhoic * a_i_1d(ji) * zdum * r1_rdtice
+            wfx_ice_sub_1d(ji) =  wfx_ice_sub_1d(ji) - rhoic * a_i_1d(ji) * zdum * r1_rdtice
+
             ! update remaining mass flux
             zevap_rema(ji)  = zevap_rema(ji) + zdum * rhoic
@@ -621 +623 @@
       END DO
 
-      ! add snow melt to mass snow mass flux to the ocean
-      DO ji = kideb, kiut
-         wfx_snw_1d(ji) = wfx_snw_1d(ji) + wfx_snw_sum_1d(ji)
+      ! Water fluxes
+      DO ji = kideb, kiut
+         wfx_snw_1d(ji) = wfx_snw_1d(ji)     + wfx_snw_sum_1d(ji) ! add snow melt to snow mass flux to the ocean
+         wfx_sub_1d(ji) = wfx_snw_sub_1d(ji) + wfx_ice_sub_1d(ji) ! sum ice and snow sublimation contributions
       END DO