Changeset 11355 for NEMO/branches/UKMO/NEMO_4.0_GO8_coupled_iodef/src/ICE/icethd_do.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_do.F90 (modified) (3 diffs)

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

@@ -207 +207 @@
       END DO
 
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'npti in icetd_do = ',npti
+         WRITE(9419,*) 'max qlead = ',MAXVAL(  qlead(:,:) )
+         WRITE(9419,*) 'min qlead = ',MINVAL(  qlead(:,:) )
+         WRITE(9419,*) 'qlead(48,31) = ', qlead(48,31)
+         WRITE(9419,*) 'tau_icebfr(48,31) = ',tau_icebfr(48,31)
+      ENDIF
+
       ! Move from 2-D to 1-D vectors
       IF ( npti > 0 ) THEN
@@ -230 +238 @@
          CALL tab_2d_1d( npti, nptidx(1:npti), rn_amax_1d(1:npti) , rn_amax_2d )
          CALL tab_2d_1d( npti, nptidx(1:npti), sss_1d    (1:npti) , sss_m      )
+
+         IF(narea == 419) THEN
+            WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+            WRITE(9419,*) 'max wfx_opw_1d icethd_do after to 1d = ',MAXVAL(  ABS( wfx_opw_1d(:) )  )
+         ENDIF
 
          ! Convert units for ice internal energy
@@ -295 +308 @@
             sfx_opw_1d(ji) = sfx_opw_1d(ji) - zv_newice(ji) * rhoi * zs_newice(ji) * r1_rdtice
          END DO
+
+         IF(narea == 419) THEN
+            WRITE(9419,*) 'max wfx_opw_1d icethd_do after calculation = ',MAXVAL(  ABS( wfx_opw_1d(:) )  )
+            WRITE(9419,*) 'max zv_newice icethd_do after calculation = ',MAXVAL(  ABS( zv_newice(:) )  )
+            WRITE(9419,*) 'rhoi = ',rhoi
+            WRITE(9419,*) 'r1_rdtice = ',r1_rdtice
+            WRITE(9419,*) 'max qlead_1d icethd_do after calculation = ',MAXVAL(  ABS( qlead_1d(:) )  )
+            WRITE(9419,*) 'max ze_newice icethd_do after calculation = ',MAXVAL(  ABS( ze_newice(:) )  )
+            WRITE(9419,*) 'max t_bo_1d icethd_do after calculation = ',MAXVAL(  ABS( t_bo_1d(:) )  )
+         ENDIF
          
          zv_frazb(1:npti) = 0._wp