Changeset 13442 for NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/icedyn_adv_pra.F90

Timestamp:

2020-08-27T17:13:15+02:00 (4 years ago)

Author:

dancopsey

Message:

Put majority of code in to get melt pond freezing to feed back into ice thickness.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/icedyn_adv_pra.F90 (modified) (12 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_fix_cpl_retain_melt/src/ICE/icedyn_adv_pra.F90

@@ -45 +45 @@
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   sxvp , syvp , sxxvp , syyvp , sxyvp    ! melt pond volume
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   sxlhp, sylhp, sxxlhp, syylhp, sxylhp   ! melt pond lid thickness
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   sxsalp, sysalp, sxxsalp, syysalp, sxysalp   ! melt pond salinity
 
    !! * Substitutions
@@ -56 +57 @@
 
    SUBROUTINE ice_dyn_adv_pra( kt, pu_ice, pv_ice,  &
-      &                        pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, plh_ip, pe_s, pe_i )
+      &                        pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, plh_ip, psv_ip, pe_s, pe_i )
       !!----------------------------------------------------------------------
       !!                **  routine ice_dyn_adv_pra  **
@@ -80 +81 @@
       REAL(wp), DIMENSION(:,:,:)  , INTENT(inout) ::   pv_ip      ! melt pond volume
       REAL(wp), DIMENSION(:,:,:)  , INTENT(inout) ::   plh_ip     ! melt pond lid thickness
+      REAL(wp), DIMENSION(:,:,:)  , INTENT(inout) ::   psv_ip     ! melt pond salt content
       REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) ::   pe_s       ! snw heat content
       REAL(wp), DIMENSION(:,:,:,:), INTENT(inout) ::   pe_i       ! ice heat content
@@ -132 +134 @@
             z0vp(:,:,jl)  = pv_ip(:,:,jl) * e1e2t(:,:)     ! Melt pond volume
             z0lhp(:,:,jl)  = plh_ip(:,:,jl) * e1e2t(:,:)   ! Melt pond lid thickness
+            z0smp(:,:,jl)  = psv_ip(:,:,jl) * e1e2t(:,:)   ! Melt pond salt content
          ENDIF
       END DO
@@ -172 +175 @@
                CALL adv_x( zdt , pu_ice , 1._wp , zarea , z0lhp, sxlhp, sxxlhp, sylhp, syylhp, sxylhp)    !--- melt pond lid thickness
                CALL adv_y( zdt , pv_ice , 0._wp , zarea , z0lhp, sxlhp, sxxlhp, sylhp, syylhp, sxylhp)
+               CALL adv_x( zdt , pu_ice , 1._wp , zarea , z0smp, sxsalp, sxxsalp, sysalp, syysalp, sxysalp)    !--- melt pond salinity
+               CALL adv_y( zdt , pv_ice , 0._wp , zarea , z0smp, sxsalp, sxxsalp, sysalp, syysalp, sxysalp)
             ENDIF
          END DO
@@ -209 +214 @@
                CALL adv_y( zdt , pv_ice , 1._wp , zarea , z0lhp, sxlhp, sxxlhp, sylhp, syylhp, sxylhp)    !--- melt pond lid thickness
                CALL adv_x( zdt , pu_ice , 0._wp , zarea , z0lhp, sxlhp, sxxlhp, sylhp, syylhp, sxylhp) 
+               CALL adv_x( zdt , pu_ice , 1._wp , zarea , z0smp, sxsalp, sxxsalp, sysalp, syysalp, sxysalp)    !--- melt pond salinity
+               CALL adv_y( zdt , pv_ice , 0._wp , zarea , z0smp, sxsalp, sxxsalp, sysalp, syysalp, sxysalp)
             ENDIF
          END DO
@@ -233 +240 @@
             pv_ip(:,:,jl) = z0vp(:,:,jl) * r1_e1e2t(:,:) * tmask(:,:,1)
             plh_ip(:,:,jl) = z0lhp (:,:,jl) * r1_e1e2t(:,:) * tmask(:,:,1)
+            psv_ip(:,:,jl) = z0smp(:,:,jl) * r1_e1e2t(:,:) * tmask(:,:,1)
          ENDIF
       END DO
@@ -239 +247 @@
       !     Remove negative values (conservation is ensured)
       !     (because advected fields are not perfectly bounded and tiny negative values can occur, e.g. -1.e-20)
-      CALL ice_var_zapneg( zdt, pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, plh_ip, pe_s, pe_i )
+      CALL ice_var_zapneg( zdt, pato_i, pv_i, pv_s, psv_i, poa_i, pa_i, pa_ip, pv_ip, plh_ip, psv_ip, pe_s, pe_i )
       !
       ! --- Ensure snow load is not too big --- !
@@ -662 +670 @@
          &      sxvp(jpi,jpj,jpl)  , syvp (jpi,jpj,jpl) , sxxvp (jpi,jpj,jpl) , syyvp (jpi,jpj,jpl) , sxyvp (jpi,jpj,jpl) ,   &
          &      sxlhp(jpi,jpj,jpl) , sylhp (jpi,jpj,jpl), sxxlhp (jpi,jpj,jpl), syylhp (jpi,jpj,jpl), sxylhp (jpi,jpj,jpl),   &
+         &      sxsalp(jpi,jpj,jpl) , sysalp (jpi,jpj,jpl), sxxsalp (jpi,jpj,jpl), syysalp (jpi,jpj,jpl), sxysalp (jpi,jpj,jpl),   &
          !
          &      sxc0 (jpi,jpj,nlay_s,jpl) , syc0 (jpi,jpj,nlay_s,jpl) , sxxc0(jpi,jpj,nlay_s,jpl) , &
@@ -780 +789 @@
                CALL iom_get( numrir, jpdom_autoglo, 'syylhp', syylhp )
                CALL iom_get( numrir, jpdom_autoglo, 'sxylhp', sxylhp )
+               !                                                     ! melt pond salinity
+               CALL iom_get( numrir, jpdom_autoglo, 'sxsalp' , sxsalp  )
+               CALL iom_get( numrir, jpdom_autoglo, 'sysalp' , sysalp  )
+               CALL iom_get( numrir, jpdom_autoglo, 'sxxsalp', sxxsalp )
+               CALL iom_get( numrir, jpdom_autoglo, 'syysalp', syysalp )
+               CALL iom_get( numrir, jpdom_autoglo, 'sxysalp', sxysalp )
             ENDIF
             !
@@ -798 +813 @@
                sxvp  = 0._wp   ;   syvp  = 0._wp   ;   sxxvp  = 0._wp   ;   syyvp  = 0._wp   ;   sxyvp  = 0._wp   ! melt pond volume
                sxlhp  = 0._wp  ;   sylhp  = 0._wp  ;   sxxlhp  = 0._wp  ;   syylhp  = 0._wp  ;   sxylhp  = 0._wp  ! melt pond lid thickness
+               sxsalp  = 0._wp  ;   sysalp  = 0._wp  ;   sxxsalp  = 0._wp  ;   syysalp  = 0._wp  ;   sxysalp  = 0._wp  ! melt pond salinity
             ENDIF
          ENDIF
@@ -884 +900 @@
             CALL iom_rstput( iter, nitrst, numriw, 'syylhp', syylhp )
             CALL iom_rstput( iter, nitrst, numriw, 'sxylhp', sxylhp )
+            !                                                        ! melt pond salinity
+            CALL iom_rstput( iter, nitrst, numriw, 'sxsalp' , sxsalp  )
+            CALL iom_rstput( iter, nitrst, numriw, 'sysalp' , sysalp  )
+            CALL iom_rstput( iter, nitrst, numriw, 'sxxsalp', sxxsalp )
+            CALL iom_rstput( iter, nitrst, numriw, 'syysalp', syysalp )
+            CALL iom_rstput( iter, nitrst, numriw, 'sxysalp', sxysalp )
          ENDIF
          !