Changeset 12369 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icecor.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/icecor.F90 (modified) (2 diffs)

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

@@ -85 +85 @@
       IF ( nn_icesal == 2 ) THEN    !  salinity must stay in bounds [Simin,Simax]        !
       !                             !-----------------------------------------------------
+
+      WRITE(numout,*) 'icecor: sv_i before bounding = ',sv_i(3,4,1)
+
          zzc = rhoi * r1_rdtice
          DO jl = 1, jpl
@@ -92 +95 @@
                   sv_i(ji,jj,jl) = MIN(  MAX( rn_simin*v_i(ji,jj,jl) , sv_i(ji,jj,jl) ) , rn_simax*v_i(ji,jj,jl)  )
                   sfx_res(ji,jj) = sfx_res(ji,jj) - ( sv_i(ji,jj,jl) - zsal ) * zzc   ! associated salt flux
+
+                  IF ( ji == 3 .AND. jj == 4 .AND. jl == 1 ) THEN
+                     WRITE(numout,*) 'icecor: zsal, v_i(ji,jj,jl), sv_i(ji,jj,jl) = ',zsal, v_i(ji,jj,jl), sv_i(ji,jj,jl)
+                  ENDIF
                END DO
             END DO
          END DO
       ENDIF
+
+
+      WRITE(numout,*) 'icecor before redist: v_i(42,26,1), a_i(42,26,1), h_i = ',v_i(42,26,1), ' ', a_i(42,26,1), ' ', v_i(42,26,1)/a_i(42,26,1)
+
       !                             !-----------------------------------------------------
       !                             !  Rebin categories with thickness out of bounds     !
       !                             !-----------------------------------------------------
       IF ( jpl > 1 )   CALL ice_itd_reb( kt )
+
+      WRITE(numout,*) 'icecor after redist: v_i(42,26,1), a_i(42,26,1), h_i = ',v_i(42,26,1), ' ', a_i(42,26,1), ' ', v_i(42,26,1)/a_i(42,26,1)
 
       !                             !-----------------------------------------------------