Changeset 12475 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn.F90

Timestamp:

2020-02-27T12:27:50+01:00 (4 years ago)

Author:

dancopsey

Message:

• Add more print statements.
• Move away from using snow to ice diagnostics and use a new snow to pond one instead.

File:

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

NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn.F90

@@ -105 +105 @@
       !
 
-      write(numout,*)'ice_dyn 0: e_i = ',e_i (42,26,1,1)
+      write(numout,*)'ice_dyn 0:sv_i, e_i, t_i = ',sv_i(3,4,1), ' ', e_i(3,4,1,1), ' ', t_i(3,4,1,1)
+      write(numout,*)'ice_dyn 0:a_i, h_i, v_i = ',a_i(3,4,1), ' ', h_i(3,4,1), ' ', v_i(3,4,1)
+      write(numout,*)'ice_dyn 0:s_i = ',s_i(3,4,1)
+      write(numout,*)'ice_dyn 0: MEAN(sv_i) = ',SUM(sv_i)/SIZE(sv_i)
+      write(numout,*)'ice_dyn 0: MEAN(e_i) = ',SUM(e_i)/SIZE(e_i)
+      write(numout,*)'ice_dyn 0: MEAN(t_i) = ',SUM(t_i)/SIZE(t_i)
+      write(numout,*)'ice_dyn 0: MEAN(a_i) = ',SUM(a_i)/SIZE(a_i)
+      write(numout,*)'ice_dyn 0: MEAN(h_i) = ',SUM(h_i)/SIZE(h_i)
+      write(numout,*)'ice_dyn 0: MEAN(v_i) = ',SUM(v_i)/SIZE(v_i)
+      write(numout,*)'ice_dyn 0: MEAN(s_i) = ',SUM(s_i)/SIZE(s_i)
+      write(numout,*)'ice_dyn 0: MEAN(u_ice) = ',SUM(u_ice)/SIZE(u_ice)
+      write(numout,*)'ice_dyn 0: MEAN(e2u) = ',SUM(e2u)/SIZE(e2u)
+      write(numout,*)'ice_dyn 0: u_ice, e2u = ',u_ice(3,4), e2u(3,4)
+      
       !
       SELECT CASE( nice_dyn )          !-- Set which dynamics is running
@@ -111 +124 @@
       CASE ( np_dynALL )           !==  all dynamical processes  ==!
          !
-         CALL ice_dyn_rhg   ( kt )                                          ! -- rheology  
+         CALL ice_dyn_rhg   ( kt )                                          ! -- rheology 
+          write(numout,*)'ice_dyn 1:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1)
+          write(numout,*)'ice_dyn 1: u_ice, e2u = ',u_ice(3,4), e2u(3,4)
          CALL ice_dyn_adv   ( kt )                                          ! -- advection of ice
+          write(numout,*)'ice_dyn 2:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1)
          CALL ice_dyn_rdgrft( kt )                                          ! -- ridging/rafting 
+          write(numout,*)'ice_dyn 3:sv_i, e_s = ',sv_i(3,4,1), ' ', e_i (3,4,1,1)
          CALL ice_cor       ( kt , 1 )                                      ! -- Corrections
          !
@@ -150 +167 @@
       END SELECT
 
-      write(numout,*)'ice_dyn 1: e_i = ',e_i (42,26,1,1)
+      write(numout,*)'ice_dyn 4: e_i = ',e_i (42,26,1,1)
       !
       !
@@ -174 +191 @@
          !
       ENDIF
-      write(numout,*)'ice_dyn 2: e_i = ',e_i (42,26,1,1)
+      write(numout,*)'ice_dyn 5: e_i = ',e_i (42,26,1,1)
       !
       ! controls