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

Timestamp:

2019-07-29T18:10:24+02:00 (5 years ago)

Author:

dancopsey

Message:

More print statements including fix

File:

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

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

@@ -115 +115 @@
 
       at_i(:,:) = SUM( a_i, dim=3 )
+
+      IF(narea == 68) THEN
+         WRITE(968,*) 'utau_ice in icethd_do = ',utau_ice(54,1)
+         WRITE(968,*) 'vtau_ice in icethd_do = ',vtau_ice(54,1)
+      ENDIF
+
       !------------------------------------------------------------------------------!
       ! 1) Collection thickness of ice formed in leads and polynyas
@@ -207 +213 @@
       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)
+      IF(narea == 68) THEN
+         WRITE(968,*) 'npti in icetd_do = ',npti
+         WRITE(968,*) 'max qlead = ',MAXVAL(  qlead(:,:) )
+         WRITE(968,*) 'min qlead = ',MINVAL(  qlead(:,:) )
+         WRITE(968,*) 'qlead(48,31) = ', qlead(48,31)
+         WRITE(968,*) 'tau_icebfr(48,31) = ',tau_icebfr(48,31)
+         WRITE(968,*) 'sum wfx_opw = ',SUM( wfx_opw )
       ENDIF
 
@@ -239 +246 @@
          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(:) )  )
+         IF(narea == 68) THEN
+            WRITE(968,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+            WRITE(968,*) 'max wfx_opw_1d icethd_do after to 1d = ',MAXVAL(  ABS( wfx_opw_1d(:) )  )
+            WRITE(968,*) 'sum wfx_opw_1d = ',SUM( wfx_opw_1d )
          ENDIF
 
@@ -309 +317 @@
          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(:) )  )
+         IF(narea == 68) THEN
+            WRITE(968,*) 'sum wfx_opw_1d icethd_do after calculation = ',SUM( wfx_opw_1d )
+            WRITE(968,*) 'sum zv_newice icethd_do after calculation = ',SUM(  zv_newice)
+            WRITE(968,*) 'rhoi = ',rhoi
+            WRITE(968,*) 'r1_rdtice = ',r1_rdtice
+            WRITE(968,*) 'sum qlead_1d icethd_do after calculation = ',SUM( qlead_1d)
+            WRITE(968,*) 'sum ze_newice icethd_do after calculation = ',SUM( ze_newice)
+            WRITE(968,*) 'sum t_bo_1d icethd_do after calculation = ',SUM( t_bo_1d)
          ENDIF