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

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

@@ -145 +145 @@
 #endif
                                         CALL store_fields             ! Store now ice values
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_sbc_tau = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+         WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31)
+      ENDIF
          !
          !------------------------------------------------!
@@ -159 +165 @@
                                         CALL ice_rst_opn( kt )        ! Open Ice restart file (if necessary) 
          !
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_dyn = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+         WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31)
+      ENDIF
+
          IF( ln_icedyn .AND. .NOT.lk_c1d )   &
             &                           CALL ice_dyn( kt )            ! -- Ice dynamics
@@ -165 +178 @@
          IF( ln_icethd .AND. ln_bdy )   CALL bdy_ice( kt )            ! -- bdy ice thermo
          !
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+         WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31)
+      ENDIF
+
          !                          !==  previous lead fraction and ice volume for flux calculations
                                         CALL ice_var_glo2eqv          ! h_i and h_s for ice albedo calculation
@@ -183 +203 @@
          !    qprec_ice, qevap_ice
          !------------------------------------------------------!
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_sbc_flx = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+         WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31)
+      ENDIF
+
                                         CALL ice_sbc_flx( kt, ksbc )
          !----------------------------!
          ! --- ice thermodynamics --- !
          !----------------------------!
-         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics      
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_thd = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+         WRITE(9419,*) 'qlead(48,31) = ',qlead(48,31)
+      ENDIF
+
+         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics  
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_cor = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+      ENDIF
+    
          !
                                         CALL ice_cor( kt , 2 )        ! -- Corrections
          !
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_var_glo2eqv = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+      ENDIF
+
                                         CALL ice_var_glo2eqv          ! necessary calls (at least for coupling)
                                         CALL ice_var_agg( 2 )         ! necessary calls (at least for coupling)
          !
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_update_flx = ',emp(48,31)
+         WRITE(9419,*) 'wfx_opw(48,31) = ',wfx_opw(48,31)
+      ENDIF
+
                                         CALL ice_update_flx( kt )     ! -- Update ocean surface mass, heat and salt fluxes
          !
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp before ice_dia = ',emp(48,31)
+      ENDIF
+
          IF( ln_icediahsb )             CALL ice_dia( kt )            ! -- Diagnostics outputs 
          !
@@ -203 +260 @@
          !
          IF( ln_icectl )                CALL ice_ctl( kt )            ! -- alerts in case of model crash
+
+      IF(narea == 419) THEN
+         WRITE(9419,*) 'max emp ice_stp at end = ',emp(48,31)
+      ENDIF
+
          !
       ENDIF   ! End sea-ice time step only