Dummy and unused variables.

cicecore/cicedynB/analysis/ice_diagnostics.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_diagnostics.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 
 ! Diagnostic information output during run
@@ -14,7 +13,7 @@ module ice_diagnostics
       use ice_kinds_mod
       use ice_communicate, only: my_task, master_task
       use ice_constants, only: c0
-      use ice_calendar, only: diagfreq, istep1, istep
+      use ice_calendar, only: istep1
       use ice_fileunits, only: nu_diag
       use ice_fileunits, only: flush_fileunit
       use ice_exit, only: abort_ice
@@ -86,7 +85,7 @@ module ice_diagnostics
 
       ! printing info for routine print_state
       ! iblkp, ip, jp, mtask identify the grid cell to print
-      character (char_len) :: plabel
+!     character (char_len) :: plabel
       integer (kind=int_kind), parameter, public :: &
          check_step = 999999999, & ! begin printing at istep1=check_step
          iblkp = 1, &      ! block number 
@@ -116,14 +115,14 @@ subroutine runtime_diags (dt)
       use ice_domain_size, only: ncat, n_aero, max_blocks
       use ice_flux, only: alvdr, alidr, alvdf, alidf, evap, fsnow, frazil, &
           fswabs, fswthru, flw, flwout, fsens, fsurf, flat, frzmlt_init, frain, fpond, &
-          coszen, fhocn_ai, fsalt_ai, fresh_ai, frazil_diag, &
+          fhocn_ai, fsalt_ai, fresh_ai, frazil_diag, &
           update_ocn_f, Tair, Qa, fsw, fcondtop, meltt, meltb, meltl, snoice, &
           dsnow, congel, sst, sss, Tf, fhocn, &
           swvdr, swvdf, swidr, swidf, &
           alvdr_init, alvdf_init, alidr_init, alidf_init
       use ice_flux_bgc, only: faero_atm, faero_ocn
       use ice_global_reductions, only: global_sum, global_sum_prod, global_maxval
-      use ice_grid, only: lmask_n, lmask_s, tarean, tareas, grid_type
+      use ice_grid, only: lmask_n, lmask_s, tarean, tareas
       use ice_state   ! everything
 #ifdef CESMCOUPLED
       use ice_prescribed_mod, only: prescribed_ice
@@ -1216,7 +1215,7 @@ subroutine total_salt (work)
 
       use ice_blocks, only: nx_block, ny_block
       use ice_domain, only: nblocks
-      use ice_domain_size, only: ncat, nilyr, nslyr, max_blocks
+      use ice_domain_size, only: ncat, nilyr, max_blocks
       use ice_grid, only: tmask
       use ice_state, only: vicen, trcrn
 
@@ -1595,7 +1594,7 @@ subroutine print_points_state(plabel,ilabel)
 
       real (kind=dbl_kind) :: &
            eidebug, esdebug, &
-           qi, qs, Tsnow, &
+           qi, qs, &
            puny
 
       integer (kind=int_kind) :: m, n, k, i, j, iblk, nt_Tsfc, nt_qice, nt_qsno
@@ -1675,41 +1674,41 @@ subroutine print_points_state(plabel,ilabel)
       write(nu_diag,*) trim(llabel),'uvel=',uvel(i,j,iblk)
       write(nu_diag,*) trim(llabel),'vvel=',vvel(i,j,iblk)
 
-      !  write(nu_diag,*) ' '
-      !  write(nu_diag,*) 'atm states and fluxes'
-      !  write(nu_diag,*) '            uatm    = ',uatm (i,j,iblk)
-      !  write(nu_diag,*) '            vatm    = ',vatm (i,j,iblk)
-      !  write(nu_diag,*) '            potT    = ',potT (i,j,iblk)
-      !  write(nu_diag,*) '            Tair    = ',Tair (i,j,iblk)
-      !  write(nu_diag,*) '            Qa      = ',Qa   (i,j,iblk)
-      !  write(nu_diag,*) '            rhoa    = ',rhoa (i,j,iblk)
-      !  write(nu_diag,*) '            swvdr   = ',swvdr(i,j,iblk)
-      !  write(nu_diag,*) '            swvdf   = ',swvdf(i,j,iblk)
-      !  write(nu_diag,*) '            swidr   = ',swidr(i,j,iblk)
-      !  write(nu_diag,*) '            swidf   = ',swidf(i,j,iblk)
-      !  write(nu_diag,*) '            flw     = ',flw  (i,j,iblk)
-      !  write(nu_diag,*) '            frain   = ',frain(i,j,iblk)
-      !  write(nu_diag,*) '            fsnow   = ',fsnow(i,j,iblk)
-      !  write(nu_diag,*) ' '
-      !  write(nu_diag,*) 'ocn states and fluxes'
-      !  write(nu_diag,*) '            frzmlt  = ',frzmlt (i,j,iblk)
-      !  write(nu_diag,*) '            sst     = ',sst    (i,j,iblk)
-      !  write(nu_diag,*) '            sss     = ',sss    (i,j,iblk)
-      !  write(nu_diag,*) '            Tf      = ',Tf     (i,j,iblk)
-      !  write(nu_diag,*) '            uocn    = ',uocn   (i,j,iblk)
-      !  write(nu_diag,*) '            vocn    = ',vocn   (i,j,iblk)
-      !  write(nu_diag,*) '            strtltx = ',strtltx(i,j,iblk)
-      !  write(nu_diag,*) '            strtlty = ',strtlty(i,j,iblk)
-      !  write(nu_diag,*) ' '
-      !  write(nu_diag,*) 'srf states and fluxes'
-      !  write(nu_diag,*) '            Tref    = ',Tref  (i,j,iblk)
-      !  write(nu_diag,*) '            Qref    = ',Qref  (i,j,iblk)
-      !  write(nu_diag,*) '            Uref    = ',Uref  (i,j,iblk)
-      !  write(nu_diag,*) '            fsens   = ',fsens (i,j,iblk)
-      !  write(nu_diag,*) '            flat    = ',flat  (i,j,iblk)
-      !  write(nu_diag,*) '            evap    = ',evap  (i,j,iblk)
-      !  write(nu_diag,*) '            flwout  = ',flwout(i,j,iblk)
-      !  write(nu_diag,*) ' '
+      write(nu_diag,*) ' '
+      write(nu_diag,*) 'atm states and fluxes'
+      write(nu_diag,*) '            uatm    = ',uatm (i,j,iblk)
+      write(nu_diag,*) '            vatm    = ',vatm (i,j,iblk)
+      write(nu_diag,*) '            potT    = ',potT (i,j,iblk)
+      write(nu_diag,*) '            Tair    = ',Tair (i,j,iblk)
+      write(nu_diag,*) '            Qa      = ',Qa   (i,j,iblk)
+      write(nu_diag,*) '            rhoa    = ',rhoa (i,j,iblk)
+      write(nu_diag,*) '            swvdr   = ',swvdr(i,j,iblk)
+      write(nu_diag,*) '            swvdf   = ',swvdf(i,j,iblk)
+      write(nu_diag,*) '            swidr   = ',swidr(i,j,iblk)
+      write(nu_diag,*) '            swidf   = ',swidf(i,j,iblk)
+      write(nu_diag,*) '            flw     = ',flw  (i,j,iblk)
+      write(nu_diag,*) '            frain   = ',frain(i,j,iblk)
+      write(nu_diag,*) '            fsnow   = ',fsnow(i,j,iblk)
+      write(nu_diag,*) ' '
+      write(nu_diag,*) 'ocn states and fluxes'
+      write(nu_diag,*) '            frzmlt  = ',frzmlt (i,j,iblk)
+      write(nu_diag,*) '            sst     = ',sst    (i,j,iblk)
+      write(nu_diag,*) '            sss     = ',sss    (i,j,iblk)
+      write(nu_diag,*) '            Tf      = ',Tf     (i,j,iblk)
+      write(nu_diag,*) '            uocn    = ',uocn   (i,j,iblk)
+      write(nu_diag,*) '            vocn    = ',vocn   (i,j,iblk)
+      write(nu_diag,*) '            strtltx = ',strtltx(i,j,iblk)
+      write(nu_diag,*) '            strtlty = ',strtlty(i,j,iblk)
+      write(nu_diag,*) ' '
+      write(nu_diag,*) 'srf states and fluxes'
+      write(nu_diag,*) '            Tref    = ',Tref  (i,j,iblk)
+      write(nu_diag,*) '            Qref    = ',Qref  (i,j,iblk)
+      write(nu_diag,*) '            Uref    = ',Uref  (i,j,iblk)
+      write(nu_diag,*) '            fsens   = ',fsens (i,j,iblk)
+      write(nu_diag,*) '            flat    = ',flat  (i,j,iblk)
+      write(nu_diag,*) '            evap    = ',evap  (i,j,iblk)
+      write(nu_diag,*) '            flwout  = ',flwout(i,j,iblk)
+      write(nu_diag,*) ' '
 
       endif   ! my_task
       enddo   ! ncnt

cicecore/cicedynB/analysis/ice_diagnostics_bgc.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_diagnostics_bgc.F90 1447 2016-04-28 18:09:53Z afrobert@nps.edu $
 !=======================================================================
 
 ! Diagnostic information output during run
@@ -14,12 +13,11 @@ module ice_diagnostics_bgc
       use ice_kinds_mod
       use ice_communicate, only: my_task, master_task
       use ice_constants, only: c0, mps_to_cmpdy, c100, p5, c1
-      use ice_calendar, only: diagfreq, istep1, istep
       use ice_fileunits, only: nu_diag
       use ice_fileunits, only: flush_fileunit
       use ice_exit, only: abort_ice
       use icepack_intfc, only: icepack_warnings_flush, icepack_warnings_aborted
-      use icepack_intfc, only: icepack_max_algae, icepack_max_aero, icepack_max_dic
+      use icepack_intfc, only: icepack_max_algae, icepack_max_aero
       use icepack_intfc, only: icepack_max_doc, icepack_max_don, icepack_max_fe
       use icepack_intfc, only: icepack_query_parameters, icepack_query_tracer_flags
       use icepack_intfc, only: icepack_query_tracer_indices
@@ -42,18 +40,14 @@ module ice_diagnostics_bgc
 !          Cecilia M. Bitz, UW
 !          Nicole Jeffery, LANL
 
-      subroutine hbrine_diags (dt)
+      subroutine hbrine_diags
 
       use ice_arrays_column, only: darcy_V
       use ice_broadcast, only: broadcast_scalar, broadcast_array
-      use ice_diagnostics, only: npnt, print_points, pmloc, piloc, pjloc, pbloc, &
-                                plat, plon
-      use ice_domain_size, only: ncat, nltrcr, nilyr
+      use ice_diagnostics, only: npnt, print_points, pmloc, piloc, pjloc, pbloc
+      use ice_domain_size, only: nilyr
       use ice_state, only: aice, aicen, vicen, vice, trcr, trcrn
 
-      real (kind=dbl_kind), intent(in) :: &
-         dt      ! time step
-
       ! local variables
 
       integer (kind=int_kind) :: &
@@ -161,25 +155,21 @@ end subroutine hbrine_diags
 ! authors: Elizabeth C. Hunke, LANL
 !          Nicole Jeffery, LANL
 
-      subroutine bgc_diags (dt)
+      subroutine bgc_diags
 
       use ice_arrays_column, only: ocean_bio, zfswin, fbio_atmice, fbio_snoice, &
           Zoo, grow_net, ice_bio_net, trcrn_sw
       use ice_broadcast, only: broadcast_scalar, broadcast_array
       use ice_diagnostics, only: npnt, print_points, pmloc, piloc, pjloc, pbloc
-      use ice_domain_size, only: ncat, nltrcr, nblyr, n_algae, n_zaero, &
-          n_dic, n_doc, n_don, n_fed, n_fep, nilyr, nslyr
+      use ice_domain_size, only: ncat, nblyr, n_algae, n_zaero, &
+          n_doc, n_don, n_fed, n_fep, nilyr, nslyr
       use ice_flux_bgc, only: flux_bio, flux_bio_atm
-      use ice_state, only:aice, vicen, vice, trcr
-      use ice_timers, only: timer_bgc, ice_timer_start, ice_timer_stop
-
-      real (kind=dbl_kind), intent(in) :: &
-         dt      ! time step
+      use ice_state, only: vicen, vice, trcr
 
       ! local variables
 
       integer (kind=int_kind) :: &
-         i, j, k, n, nn, ii,jj, iblk,kk, klev
+         i, j, k, n, nn, iblk,kk, klev
       ! fields at diagnostic points
       real (kind=dbl_kind), dimension(npnt) :: &
          pNit_sk, pAm_sk, pSil_sk, phum_sk, &
@@ -870,24 +860,20 @@ end subroutine bgc_diags
 !          Cecilia M. Bitz, UW
 !          Nicole Jeffery, LANL
 
-      subroutine zsal_diags (dt)
+      subroutine zsal_diags
 
       use ice_arrays_column, only: fzsal, fzsal_g, sice_rho, bTiz, &
           iDi, bphi, dhbr_top, dhbr_bot, darcy_V
-      use ice_blocks, only: nx_block, ny_block
       use ice_broadcast, only: broadcast_scalar, broadcast_array
       use ice_diagnostics, only: npnt, print_points, pmloc, piloc, pjloc, &
-          pbloc, plat, plon
-      use ice_domain_size, only: max_blocks, nblyr, ncat, nilyr
+          pbloc
+      use ice_domain_size, only: nblyr, ncat, nilyr
       use ice_state, only: aicen, aice, vice, trcr, trcrn, vicen, vsno
 
-      real (kind=dbl_kind), intent(in) :: &
-         dt      ! time step
-
       ! local variables
 
       integer (kind=int_kind) :: &
-         i, j, k, n, nn, ii,jj, iblk
+         i, j, k, n, nn, iblk
 
       ! fields at diagnostic points
       real (kind=dbl_kind), dimension(npnt) :: &
@@ -897,16 +883,13 @@ subroutine zsal_diags (dt)
 
       ! vertical  fields of category 1 at diagnostic points for bgc layer model
       real (kind=dbl_kind), dimension(npnt,nblyr+2) :: &
-         pphin, pgrid, pphin1
+         pphin, pphin1
       real (kind=dbl_kind), dimension(npnt,nblyr) :: &
          pSin, pSice, pSin1
 
       real (kind=dbl_kind), dimension(npnt,nblyr+1) :: &
          pbTiz, piDin
 
-      real (kind=dbl_kind), dimension (nx_block,ny_block,max_blocks) :: &
-         work1, work2
-
       real (kind=dbl_kind) :: &
          rhosi, rhow, rhos
 

cicecore/cicedynB/analysis/ice_history.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_history.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 
 ! Driver for core history output
@@ -69,7 +68,7 @@ subroutine init_hist (dt)
           histfreq_n, nstreams
       use ice_domain_size, only: max_blocks, max_nstrm
       use ice_dyn_shared, only: kdyn
-      use ice_flux, only: mlt_onset, frz_onset, albcnt, taubx, tauby
+      use ice_flux, only: mlt_onset, frz_onset, albcnt
       use ice_history_shared ! everything
       use ice_history_mechred, only: init_hist_mechred_2D, init_hist_mechred_3Dc
       use ice_history_pond, only: init_hist_pond_2D, init_hist_pond_3Dc
@@ -1208,9 +1207,9 @@ subroutine accum_hist (dt)
           taubx, tauby, strocnx, strocny, fm, daidtt, dvidtt, daidtd, dvidtd, fsurf, &
           fcondtop, fsurfn, fcondtopn, flatn, fsensn, albcnt, &
           stressp_1, stressm_1, stress12_1, &
-          stressp_2, stressm_2, stress12_2, &
-          stressp_3, stressm_3, stress12_3, &
-          stressp_4, stressm_4, stress12_4, sig1, sig2, sigP, &
+          stressp_2, &
+          stressp_3, &
+          stressp_4, sig1, sig2, sigP, &
           mlt_onset, frz_onset, dagedtt, dagedtd, fswint_ai, keffn_top, &
           snowfrac, alvdr_ai, alvdf_ai, alidr_ai, alidf_ai
       use ice_arrays_column, only: snowfracn
@@ -1240,9 +1239,7 @@ subroutine accum_hist (dt)
 
       real (kind=dbl_kind) :: & 
            qn                , & ! temporary variable for enthalpy
-           sn                , & ! temporary variable for salinity
-           Tmlts             , & !  temporary variable for melting temperature
-           Tn                    !  temporary variable for ice temperature
+           sn                    ! temporary variable for salinity
 
       real (kind=dbl_kind), dimension (nx_block,ny_block) :: &
          worka, workb
@@ -1331,7 +1328,7 @@ subroutine accum_hist (dt)
       !---------------------------------------------------------------
 
       !$OMP PARALLEL DO PRIVATE(iblk,i,j,ilo,ihi,jlo,jhi,this_block, &
-      !$OMP             k,n,qn,ns,worka,workb,Tinz4d,Sinz4d,Tsnz4d)
+      !$OMP             k,n,qn,ns,sn,worka,workb,Tinz4d,Sinz4d,Tsnz4d)
       do iblk = 1, nblocks
          this_block = get_block(blocks_ice(iblk),iblk)         
          ilo = this_block%ilo

cicecore/cicedynB/analysis/ice_history_bgc.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_history_bgc.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 ! Biogeochemistry history output
 !
@@ -21,7 +20,7 @@ module ice_history_bgc
       use icepack_intfc, only: icepack_query_tracer_flags, &
           icepack_query_tracer_indices, icepack_query_parameters, &
           icepack_query_parameters
-      use ice_domain_size, only: max_nstrm, n_aero, nblyr, &
+      use ice_domain_size, only: max_nstrm, n_aero, &
           n_algae, n_dic, n_doc, n_don, n_zaero, n_fed, n_fep 
 
       implicit none
@@ -1730,8 +1729,6 @@ end subroutine init_hist_bgc_2D
 
       subroutine init_hist_bgc_3Dc
 
-      use ice_broadcast, only: broadcast_scalar
-      use ice_broadcast, only: broadcast_scalar
       use ice_calendar, only: nstreams
       use ice_history_shared, only: tstr3Dc, tcstr, define_hist_field
 
@@ -1816,14 +1813,14 @@ subroutine accum_hist_bgc (iblk)
           R_C2N, R_chl2N
       use ice_blocks, only: block, get_block, nx_block, ny_block
       use ice_domain, only: blocks_ice
-      use ice_domain_size, only: ncat, nblyr
+      use ice_domain_size, only: nblyr
       use ice_flux, only: sss
       use ice_flux_bgc, only: faero_atm, faero_ocn, flux_bio, flux_bio_ai, &
           fzsal_ai, fzsal_g_ai
-      use ice_history_shared, only: n2D, a2D, a3Dc, n3Dccum, &         
-          n3Dzcum, n3Dbcum, n3Dacum, a3Db, a3Da, &    
+      use ice_history_shared, only: n2D, a2D, a3Dc, &         
+          n3Dzcum, n3Dbcum, a3Db, a3Da, &    
           ncat_hist, accum_hist_field, nzblyr, nzalyr
-      use ice_state, only: trcrn, trcr, aicen, aice, vice, vicen
+      use ice_state, only: trcrn, trcr, aicen, aice, vicen
 
       integer (kind=int_kind), intent(in) :: &
            iblk                 ! block index

cicecore/cicedynB/analysis/ice_history_drag.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_history_drag.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 
 ! 2013 module for form drag parameters
@@ -8,10 +7,9 @@ module ice_history_drag
 
       use ice_kinds_mod
       use ice_domain_size, only: max_nstrm
-      use ice_constants, only: c0, c1, c100, mps_to_cmpdy
+      use ice_constants, only: c0, c1
       use ice_fileunits, only: nu_nml, nml_filename, &
           get_fileunit, release_fileunit
-      use ice_fileunits, only: nu_diag
       use ice_exit, only: abort_ice
       use icepack_intfc, only: icepack_warnings_flush, icepack_warnings_aborted
 
@@ -230,8 +228,7 @@ end subroutine init_hist_drag_2D
 
       subroutine accum_hist_drag (iblk)
 
-      use ice_history_shared, only: n2D, a2D, a3Dc, ncat_hist, &
-          accum_hist_field
+      use ice_history_shared, only: a2D, accum_hist_field
       use ice_arrays_column, only: hfreebd, hdraft, hridge, distrdg, hkeel, &
           dkeel, lfloe, dfloe, Cdn_atm, Cdn_atm_skin, Cdn_atm_floe, &
           Cdn_atm_pond, Cdn_atm_rdg, Cdn_atm_ratio, Cdn_ocn_skin, &

cicecore/cicedynB/analysis/ice_history_mechred.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_history_mechred.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 
 ! Mechanical redistribution history output

cicecore/cicedynB/analysis/ice_history_pond.F90

@@ -1,4 +1,3 @@
-!  SVN:$Id: ice_history_pond.F90 1228 2017-05-23 21:33:34Z tcraig $
 !=======================================================================
 
 ! Melt pond history output