Changeset 4322

Timestamp:

2013-11-28T13:26:52+01:00 (10 years ago)

Author:

cetlod

Message:

v3.6_alpha : defines the surface cells in offline, see ticket #1191

Location:

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OFF_SRC

Files:

• domain.F90 (modified) (1 diff)
• domrea.F90 (modified) (1 diff)

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OFF_SRC/domain.F90

@@ -71 +71 @@
       CALL dom_zgr      ! Vertical mesh and bathymetry option
       CALL dom_rea      ! Create a domain file
+
+     !
+      ! - ML - Used in dom_vvl_sf_nxt and lateral diffusion routines
+      !        but could be usefull in many other routines
+      e12t    (:,:) = e1t(:,:) * e2t(:,:)
+      e1e2t   (:,:) = e1t(:,:) * e2t(:,:)
+      e12u    (:,:) = e1u(:,:) * e2u(:,:)
+      e12v    (:,:) = e1v(:,:) * e2v(:,:)
+      e12f    (:,:) = e1f(:,:) * e2f(:,:)
+      r1_e12t (:,:) = 1._wp    / e12t(:,:)
+      r1_e12u (:,:) = 1._wp    / e12u(:,:)
+      r1_e12v (:,:) = 1._wp    / e12v(:,:)
+      r1_e12f (:,:) = 1._wp    / e12f(:,:)
+      re2u_e1u(:,:) = e2u(:,:) / e1u(:,:)
+      re1v_e2v(:,:) = e1v(:,:) / e2v(:,:)
+      !
+      hu(:,:) = 0._wp                          ! Ocean depth at U- and V-points
+      hv(:,:) = 0._wp
+      DO jk = 1, jpk
+         hu(:,:) = hu(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk)
+         hv(:,:) = hv(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk)
+      END DO
+      !                                        ! Inverse of the local depth
+      hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask(:,:,1) ) * umask(:,:,1)
+      hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask(:,:,1) ) * vmask(:,:,1)
+
       CALL dom_stp      ! Time step
       CALL dom_msk      ! Masks

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OFF_SRC/domrea.F90

@@ -153 +153 @@
          CALL iom_get( inum3, jpdom_data, 'e2u', e2u )
          CALL iom_get( inum3, jpdom_data, 'e2v', e2v )
-
-         e1e2t(:,:) = e1t(:,:) * e2t(:,:)                              ! surface at T-points
 
          CALL iom_get( inum3, jpdom_data, 'ff', ff )