Changeset 8172 for branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

Timestamp:

2017-06-14T15:44:31+02:00 (7 years ago)

Author:

vancop

Message:

SIMIP branch, generic 2D transport fluxes

File:

• branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90 (modified) (2 diffs)

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

@@ -129 +129 @@
       REAL(wp) ::   zresm                                                    ! Maximal error on ice velocity
       REAL(wp) ::   zintb, zintn                                             ! dummy argument
-      REAL(wp) ::   zswi
+      REAL(wp) ::   zswi, zfac_x, zfac_y
       
       REAL(wp), POINTER, DIMENSION(:,:) ::   zpresh                          ! temporary array for ice strength
@@ -672 +672 @@
              diag_vtau_oi(ji,jj) = ztauy_oi(ji,jj) * zswi
 
+             ! 2D ice mass, snow mass, area transport arrays (X, Y)
+             zfac_x = 0.5 * u_ice(ji,jj) * e2u(ji,jj) * zswi
+             zfac_y = 0.5 * v_ice(ji,jj) * e1v(ji,jj) * zswi
+
+             diag_xmtrp_ice(ji,jj) = rhoic * zfac_x * ( vt_i(ji+1,jj) + vt_i(ji,jj) ) ! ice mass transport, X-component
+             diag_ymtrp_ice(ji,jj) = rhoic * zfac_y * ( vt_i(ji,jj+1) + vt_i(ji,jj) ) !        ''           Y-   ''
+
+             diag_xmtrp_snw(ji,jj) = rhosn * zfac_x * ( vt_s(ji+1,jj) + vt_s(ji,jj) ) ! snow mass transport, X-component
+             diag_ymtrp_snw(ji,jj) = rhosn * zfac_y * ( vt_s(ji,jj+1) + vt_s(ji,jj) ) !          ''          Y-   ''
+
+             diag_xatrp(ji,jj)     = zfac_x * ( at_i(ji+1,jj) + at_i(ji,jj) )         ! area transport,      X-component
+             diag_yatrp(ji,jj)     = zfac_y * ( at_i(ji,jj+1) + at_i(ji,jj) )         !        ''            Y-   ''
+
          END DO
       END DO