Changeset 12326 for NEMO/branches/UKMO/NEMO_4.0.1_MIRROR_WAD_ZENV/src/OCE/DOM/domain.F90

Timestamp:

2020-01-16T07:56:47+01:00 (4 years ago)

Author:

deazer

Message:

Adjusted to account for flux form advection
should now work with both types of advection and enveloping bathymetry
moved one off caluculations of depth scales to domain.F90

File:

• NEMO/branches/UKMO/NEMO_4.0.1_MIRROR_WAD_ZENV/src/OCE/DOM/domain.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_MIRROR_WAD_ZENV/src/OCE/DOM/domain.F90

@@ -79 +79 @@
       INTEGER , DIMENSION(jpi,jpj) ::   ik_top , ik_bot       ! top and bottom ocean level
       REAL(wp), DIMENSION(jpi,jpj) ::   z1_hu_0, z1_hv_0
+!CEOD
+      REAL(wp), DIMENSION(jpi,jpj) ::  k_bot_i_min, k_bot_j_min
+      REAL(wp), DIMENSION(jpi,jpj) :: sum_e3t_0_min_ik
+      REAL(wp), DIMENSION(jpi,jpj) :: sum_e3t_0_min_jk
+      REAL(wp), DIMENSION(jpi,jpj) :: sum_e3t_0_min_ip1k
+      REAL(wp), DIMENSION(jpi,jpj) :: sum_e3t_0_min_jp1k
       !!----------------------------------------------------------------------
       !
@@ -155 +161 @@
          hv_0(:,:) = hv_0(:,:) + e3v_0(:,:,jk) * vmask(:,:,jk)
       END DO
+!CEOD Here we need to work out whihc is the the shallowest point in terms of k
+!levels at neighbouring points in i and je directions.
+! then we need to wkr out the proprtion of the water column down to ht_0 that
+! that level is for both i points.
+! e.g point i could say go down to level 10 but point i+1 only level 5
+! then work out depth of bottom of level 5 at point i over depth at point i to
+! level 10
+! This is needed later to work out where a u point depth should be when using
+! enveloping coordinates, as it wil be 1/2 depth at t point at level 5 for i,i+1
+! points Because level 5 at i can change in time, thus u bed can change in time!
+! What about under ice shelves? Need rto return to that later
+
+! 1) get the shallowest k level at nighbrougin i and j points store in ! k_bot_i_min
+
+         DO jj = 1, jpjm1
+            DO ji = 1, jpim1                ! SPG with the application of W/D gravity filters
+                 k_bot_i_min(ji,jj) = MIN( ik_bot(ji,jj), ik_bot(ji+1,jj  ))
+                 k_bot_j_min(ji,jj) = MIN( ik_bot(ji,jj), ik_bot(ji,  jj+1))
+            ENDDO
+        ENDDO
+! 2) Work out the depth down to the shallowest k levle both at point i and i+1
+! (likewise for j) 
+
+        sum_e3t_0_min_ik(:,:) = 0
+        sum_e3t_0_min_jk(:,:) = 0
+        sum_e3t_0_min_ip1k(:,:) = 0
+        sum_e3t_0_min_jp1k(:,:) = 0
+
+        !Get sum of e3t_0s down to local min
+        DO jj = 1, jpjm1
+           DO ji = 1, jpim1               
+             DO jk = 1, k_bot_i_min(ji,jj)
+ 
+                sum_e3t_0_min_ik  (ji,jj) =  sum_e3t_0_min_ik  (ji,jj) + e3t_0(ji  ,jj,jk)*tmask(ji,jj,jk)
+                sum_e3t_0_min_ip1k(ji,jj) =  sum_e3t_0_min_ip1k(ji,jj) + e3t_0(ji+1,jj,jk)*tmask(ji+1,jj,jk)
+             ENDDO
+           
+             DO jk = 1, k_bot_j_min(ji,jj)
+                sum_e3t_0_min_jk  (ji,jj) =  sum_e3t_0_min_jk  (ji,jj) + e3t_0(ji,jj  ,jk)*tmask(ji,jj,jk)
+                sum_e3t_0_min_jp1k(ji,jj) =  sum_e3t_0_min_jp1k(ji,jj) + e3t_0(ji,jj+1,jk)*tmask(ji,jj+1,jk)
+             ENDDO
+! 3)  Then work out what fraction of the at rest water column that is, we later
+! multiply the now water depth by this scale to work out the bottom of the kth
+! level at time now in dynspg_ts 
+             scaled_e3t_0_ik  (ji,jj) = sum_e3t_0_min_ik  (ji,jj)/( ht_0(ji  ,jj  ) + 1._wp - ssmask(ji,jj))
+             scaled_e3t_0_jk  (ji,jj) = sum_e3t_0_min_jk  (ji,jj)/( ht_0(ji  ,jj  ) + 1._wp - ssmask(ji,jj))
+             scaled_e3t_0_ip1k(ji,jj) = sum_e3t_0_min_ip1k(ji,jj)/( ht_0(ji+1,jj  ) + 1._wp - ssmask(ji+1,jj))
+             scaled_e3t_0_jp1k(ji,jj) = sum_e3t_0_min_jp1k(ji,jj)/( ht_0(ji  ,jj+1) + 1._wp - ssmask(ji,jj+1))
+          ENDDO
+       ENDDO
       !
       !           !==  time varying part of coordinate system  ==!