Changeset 12369 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icestp.F90

Timestamp:

2020-02-12T13:06:25+01:00 (4 years ago)

Author:

dancopsey

Message:

Add print statements

File:

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

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

@@ -80 +80 @@
    USE prtctl         ! Print control
 
+   USE lbclnk         ! lateral boundary conditions (or mpp links)
+
    IMPLICIT NONE
    PRIVATE
@@ -118 +120 @@
       INTEGER, INTENT(in) ::   ksbc    ! flux formulation (user defined, bulk, or Pure Coupled)
       !
-      INTEGER ::   jl   ! dummy loop index
+      INTEGER  :: ji, jj, jk, jl, jn, gi, gj   ! dummy loop indices
+      REAL(wp) :: delta_lon, delta_lat, lon_west, lon_east, lat_south, lat_north, lon_print, lat_print
       !!----------------------------------------------------------------------
       !
       IF( ln_timing )   CALL timing_start('ice_stp')
+
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+
+            gi = nimpp + ji - 1
+            gj = njmpp + jj - 1
+            delta_lon = (glamt(2,jj) - glamt(1,jj))/2.0
+            delta_lat = (gphit(ji,2) - gphit(ji,1))/2.0
+
+            lon_west = glamt(ji,jj) - delta_lon
+            lon_east = glamt(ji,jj) + delta_lon
+            lat_south = gphit(ji,jj) - delta_lat
+            lat_north = gphit(ji,jj) + delta_lat
+            to_print_2d(ji,jj) = 0
+
+            lon_print = 95.0
+            lat_print = -63.0  ! For bottom melt investigation
+
+            IF ( lon_west <= lon_print .AND. lon_print <= lon_east .AND. lat_south <= lat_print .AND. lat_print <= lat_north ) THEN
+              to_print_2d(ji,jj) = 1
+              write(numout,*)'icethd: location to print: ji, jj, gi,gj, to_print_2d(ji,jj) = ', ji, jj, gi,gj, to_print_2d(ji,jj)
+            ENDIF
+
+            lon_print = 10.0
+            lat_print = -62.0   ! For top melt investigation
+
+            IF ( lon_west <= lon_print .AND. lon_print <= lon_east .AND. lat_south <= lat_print .AND. lat_print <= lat_north ) THEN
+              to_print_2d(ji,jj) = 2
+              write(numout,*)'icethd: location to print: ji, jj, gi,gj, to_print_2d(ji,jj) = ', ji, jj, gi,gj, to_print_2d(ji,jj)
+            ENDIF
+
+            lon_print = 137.81
+            lat_print = -66   ! For conductivity investigation
+
+            IF ( lon_west <= lon_print .AND. lon_print <= lon_east .AND. lat_south <= lat_print .AND. lat_print <= lat_north ) THEN
+              to_print_2d(ji,jj) = 3
+              write(numout,*)'icethd: location to print: ji, jj, gi,gj, to_print_2d(ji,jj) = ', ji, jj, gi,gj, to_print_2d(ji,jj)
+            ENDIF
+
+            lon_print = 137.81 - 360.0
+
+            IF ( lon_west <= lon_print .AND. lon_print <= lon_east .AND. lat_south <= lat_print .AND. lat_print <= lat_north ) THEN
+              to_print_2d(ji,jj) = 3
+              write(numout,*)'icethd: location to print: ji, jj, gi,gj, to_print_2d(ji,jj) = ', ji, jj, gi,gj, to_print_2d(ji,jj)
+            ENDIF
+
+            lon_print = 8
+            lat_print = -67   ! For top melt investigation
+
+            IF ( lon_west <= lon_print .AND. lon_print <= lon_east .AND. lat_south <= lat_print .AND. lat_print <= lat_north ) THEN
+              to_print_2d(ji,jj) = 10
+              write(numout,*)'icethd: location to print: ji, jj, gi,gj, to_print_2d(ji,jj) = ', ji, jj, gi,gj, to_print_2d(ji,jj)
+            ENDIF
+          END DO
+        END DO
+
       !
       !                                      !-----------------------!
@@ -128 +187 @@
          !
          kt_ice = kt                              ! -- Ice model time step
+
+         write(numout,*)'icestp 0: e_s = ',e_s (42,26,1,1)
          !
          u_oce(:,:) = ssu_m(:,:)                  ! -- mean surface ocean current
@@ -152 +213 @@
          !    utau_ice, vtau_ice = surface ice stress [N/m2]
          !------------------------------------------------!
-                                        CALL ice_sbc_tau( kt, ksbc, utau_ice, vtau_ice )          
+                                        CALL ice_sbc_tau( kt, ksbc, utau_ice, vtau_ice )      
+
+         write(numout,*)'icestp 1: e_s = ',e_s (42,26,1,1)
+    
          !-------------------------------------!
          ! --- ice dynamics and advection  --- !
@@ -161 +225 @@
          IF( ln_icedyn .AND. .NOT.lk_c1d )   &
             &                           CALL ice_dyn( kt )            ! -- Ice dynamics
+
+         write(numout,*)'icestp 2: e_s = ',e_s (42,26,1,1)
          !
          !                          !==  lateral boundary conditions  ==!
          IF( ln_icethd .AND. ln_bdy )   CALL bdy_ice( kt )            ! -- bdy ice thermo
+
+         write(numout,*)'icestp 3: e_s = ',e_s (42,26,1,1)
          !
          !                          !==  previous lead fraction and ice volume for flux calculations
                                         CALL ice_var_glo2eqv          ! h_i and h_s for ice albedo calculation
+
+         write(numout,*)'icestp 4: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
+
                                         CALL ice_var_agg(1)           ! at_i for coupling 
+
+         write(numout,*)'icestp 5: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
+
                                         CALL store_fields             ! Store now ice values
+
+         write(numout,*)'icestp 6: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !
          !------------------------------------------------------!
@@ -184 +260 @@
          !------------------------------------------------------!
                                         CALL ice_sbc_flx( kt, ksbc )
+         write(numout,*)'icestp 7: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !----------------------------!
          ! --- ice thermodynamics --- !
          !----------------------------!
-         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics      
+         IF( ln_icethd )                CALL ice_thd( kt )            ! -- Ice thermodynamics   
+         write(numout,*)'icestp 8: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)   
          !
                                         CALL ice_cor( kt , 2 )        ! -- Corrections
+         write(numout,*)'icestp 9: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !
                                         CALL ice_var_glo2eqv          ! necessary calls (at least for coupling)
+         write(numout,*)'icestp 10: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
                                         CALL ice_var_agg( 2 )         ! necessary calls (at least for coupling)
+         write(numout,*)'icestp 11: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !
                                         CALL ice_update_flx( kt )     ! -- Update ocean surface mass, heat and salt fluxes
+         write(numout,*)'icestp 12: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !
          IF( ln_icediahsb )             CALL ice_dia( kt )            ! -- Diagnostics outputs 
@@ -203 +285 @@
          !
          IF( ln_icectl )                CALL ice_ctl( kt )            ! -- alerts in case of model crash
+         write(numout,*)'icestp 13: t_s, e_s = ',t_s(ji,jj,jk,jl) - rt0, ' ', e_s (42,26,1,1)
          !
       ENDIF   ! End sea-ice time step only