Changeset 15456 for NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate

Timestamp:

2021-10-28T12:15:52+02:00 (3 years ago)

Author:

hadjt

Message:

Preparation for change in SBC/tide_mod.F90
comments added to SBC/tide_mod.F90 (and BDY/bdydyn2d.F90).

Location:

NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE

Files:

• BDY/bdydyn2d.F90 (modified) (1 diff)
• SBC/tide_mod.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/BDY/bdydyn2d.F90

@@ -174 +174 @@
       TYPE(OBC_DATA),               INTENT(in) ::   dta  ! OBC external data
       INTEGER,                      INTENT(in) ::   ib_bdy  ! BDY set index
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pua2d, pva2d
-      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   pssh, phur, phvr
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pua2d, pva2d     ! JT After barotropic velocities
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   pssh, phur, phvr  ! JT SSH after, and phur, phvr = inverse of depth of water 
       LOGICAL                     , INTENT(in) ::   llrim0   ! indicate if rim 0 is treated
       INTEGER  ::   ibeg, iend                       ! length of rim to be treated (rim 0 or rim 1)

NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/SBC/tide_mod.F90

@@ -71 +71 @@
       !!----------------------------------------------------------------------
       !
-      zqy = AINT( (nyear-1901.)/4. )
-      zsy = nyear - 1900.
-      !
-      zdj  = dayjul( nyear, nmonth, nday )
-      zday = zdj + zqy - 1.
-      !
-      zhfrac = nsec_day / 3600.
+      zqy = AINT( (nyear-1901.)/4. )        ! leap years since 1901
+      zsy = nyear - 1900.                   ! years since 1900
+      !
+      zdj  = dayjul( nyear, nmonth, nday )  ! day number of year
+      zday = zdj + zqy - 1.                 ! day number of year + No of leap yrs
+                                            ! i.e. what would doy if every year = 365 day??
+      !
+      zhfrac = nsec_day / 3600.             ! The seconds of the day/3600
       !
       !----------------------------------------------------------------------