New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#1060 (minor alterations to iceberg trajectory component) – NEMO

Opened 11 years ago

Closed 11 years ago

#1060 closed Defect (fixed)

minor alterations to iceberg trajectory component

Reported by: acc Owned by: acc
Priority: low Milestone:
Component: OCE Version: v3.6
Severity: Keywords:
Cc:

Description

A couple of small corrections to the iceberg trajectory modules courtesy of Vladimir Ivchenko.

  1. A parameter correction in icbclv.F90. Line 59 corrected to use the iceberg density instead of the pure water density: 
         zfact = ( (1000._wp)**3 / ( NINT(rday) * nyear_len(1) ) ) * 850._wp
  1. Limiters placed on index calculations to prevent array bounds errors near cyclic boundaries. For example: 
    --- NEMOGCM/NEMO/OPA_SRC/ICB/icbutl.F90 (revision 3820)
+++ NEMOGCM/NEMO/OPA_SRC/ICB/icbutl.F90 (working copy)
@@ -195,23 +195,23 @@
             ! note that here there is no +0.5 added
             ! since we're looking for four T points containing quadrant we're in of 
             ! current T cell
-            ii = INT( pi     )
-            ij = INT( pj      )    ! T-point
+            ii = MAX(1, INT( pi     ))
+            ij = MAX(1, INT( pj     ))    ! T-point
             zi = pi - REAL(ii,wp)
             zj = pj - REAL(ij,wp)

Code now reflects the version in active use.

Commit History (1)

ChangesetAuthorTimeChangeLog
3821acc2013-02-25T15:13:39+01:00

Branch dev_MERGE_2012. #1060. Minor alterations to icbutl.F90 and icbclv.F90 (ICeBerg trajectory component)

Change History (1)

comment:1 Changed 11 years ago by acc

  • Resolution set to fixed
  • Status changed from new to closed

Changes committed at Changeset:3821
Ticket closed.

Note: See TracTickets for help on using tickets.