Changeset 11355 for NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icethd.F90

Timestamp:

2019-07-26T11:32:13+02:00 (5 years ago)

Author:

dancopsey

Message:

Add lots of print statements printing out values of files throughout the ocean_ice timestep.. New output will be in files crash_stats_419.out and crash_stats_68.out.

File:

• NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icethd.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icethd.F90

@@ -329 +329 @@
       !
       INTEGER ::   jk   ! dummy loop indices
+      INTEGER ::   jn , jid, jjd
       !!-----------------------------------------------------------------------
       !
@@ -335 +336 @@
       CASE( 1 )            !==  from 2D to 1D  ==!
          !                 !---------------------!
+
+         DO jn = 1, npti
+           jid        = MOD( nptidx(jn) - 1 , jpi ) + 1
+           jjd        =    ( nptidx(jn) - 1 ) / jpi + 1
+           IF( jn == 54 ) THEN
+             IF( narea == 68) THEN
+               WRITE(968,*) 'in ice_thd_1d2d: ice_cat=',kl
+               WRITE(968,*) 'in ice_thd_1d2d: jid, jjd for point 54 is ',jid,jjd
+               WRITE(968,*) 'a_i for this point = ',a_i(jid,jjd)
+             END IF
+           END IF
+         END DO
+
          CALL tab_2d_1d( npti, nptidx(1:npti), at_i_1d(1:npti), at_i             )
          CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d (1:npti), a_i (:,:,kl)     )