Changeset 14259 for NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbutl.F90

Timestamp:

2021-01-05T12:11:16+01:00 (3 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.3_icb_speed_limit : Bug fix and extra write-outs for error trap in icbutl_bilin_e.

File:

• NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbutl.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.3_icb_speed_limit/src/OCE/ICB/icbutl.F90

@@ -422 +422 @@
       IF    ( ij < mjg( 1 ) ) THEN   ;   ij = 1       ; ierr = ierr + 1
       ELSEIF( ij > mjg(jpj) ) THEN   ;   ij = jpj     ; ierr = ierr + 1
-      ELSE                           ;   ij  = mj1(ij)
+      ELSE                           ;   ij = mj1(ij)
       ENDIF
       !
@@ -428 +428 @@
       IF( ij == jpj ) THEN ; ij = ij-1 ; ierr = ierr + 1 ; END IF
       !
-      IF ( ierr > 0 ) CALL ctl_stop('STOP','icb_utl_bilin_e: an icebergs coordinates is out of valid range (out of bound error)')
+      IF ( ierr > 0 ) THEN
+         WRITE(numicb,*) '>>>> Out of valid range crash <<<<'
+         WRITE(numicb,*) 'new global position estimate (pi,pj) = ',pi,pj 
+         WRITE(numicb,*) 'mig(1), mig(jpi) : ',mig(1),mig(jpi)
+         WRITE(numicb,*) 'mjg(1), mjg(jpj) : ',mjg(1),mjg(jpj)
+         WRITE(numicb,*) 'new local position estimate (ii,ij) = ',ii,ij 
+         WRITE(numicb,*) 'jpi, jpj : ',jpi,jpj
+         WRITE(numicb,*) 'nlei, nlej : ',nlei,nlej
+         CALL ctl_stop('STOP','icb_utl_bilin_e: an icebergs coordinates is out of valid range (out of bound error)')
+      ENDIF
       !
       IF(    0.0_wp <= zi .AND. zi < 0.5_wp   ) THEN