Changeset 3383 for branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

Timestamp:

2012-05-04T16:06:37+02:00 (12 years ago)

Author:

sga

Message:

NEMO branch dev_r3337_NOCS10_ICB: rationalise date calculation (which also helps with key_agrif problems)

File:

• branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90 (modified) (4 diffs)

branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbclv.F90

@@ -117 +117 @@
       TYPE(iceberg) ::   newberg
       TYPE(point)   ::   newpt
-      REAL(wp)      ::   zddt, zcalved_to_berg, zheat_to_berg
+      REAL(wp)      ::   zday, zcalved_to_berg, zheat_to_berg
       !!----------------------------------------------------------------------
       !
       icntmax = 0
+      zday    = REAL(nday_year,wp) + REAL(nsec_day,wp)/86400.0_wp
       !
       DO jn = 1, nclasses
@@ -126 +127 @@
             DO ji = nicbdi, nicbei
                !
-               zddt  = 0._wp
                icnt = 0
                !
@@ -146 +146 @@
                   newpt%mass_of_bits   = 0._wp                          ! no bergy
                   !
-                  newpt%year   = current_year
-                  newpt%day    = current_yearday + zddt / rday
+                  newpt%year   = nyear
+                  newpt%day    = zday
                   newpt%heat_density = berg_grid%stored_heat(ji,jj) / berg_grid%stored_ice(ji,jj,jn)   ! This is in J/kg
                   !
@@ -162 +162 @@
                   berg_grid%stored_ice(ji,jj,jn) = berg_grid%stored_ice(ji,jj,jn) - zcalved_to_berg
                   !
-                  zddt  = zddt + berg_dt * 2._wp / 17._wp    ! Minor offset to start day    !!gm why???
                   icnt = icnt + 1
                   !