Add new 'thermo' namelist variable to disable thermodynamics. Add cod…

…e to set coriolis term to zero. Update box2001 namelist option file

cicecore/cicedynB/dynamics/ice_dyn_shared.F90

@@ -28,6 +28,7 @@ module ice_dyn_shared
 
       integer (kind=int_kind), public :: &
          kdyn     , & ! type of dynamics ( 1 = evp, 2 = eap )
+         thermo   , & ! thermodynamics off (0) or on (1)
          ndte         ! number of subcycles:  ndte=dt/dte
 
          character (len=char_len), public :: &
@@ -135,6 +136,8 @@ subroutine init_evp (dt)
          ! Coriolis parameter
          if (trim(coriolis) == 'constant') then
          fcor_blk(i,j,iblk) = 1.46e-4_dbl_kind ! Hibler 1979, N. Hem; 1/s
+         else if (trim(coriolis) == 'zero') then
+         fcor_blk(i,j,iblk) = 0.0
          else
          fcor_blk(i,j,iblk) = c2*omega*sin(ULAT(i,j,iblk)) ! 1/s
          endif

cicecore/cicedynB/general/ice_init.F90

@@ -79,7 +79,7 @@ subroutine input_data
       use ice_grid, only: grid_file, gridcpl_file, kmt_file, grid_type, grid_format, &
                           dxrect, dyrect
       use ice_dyn_shared, only: ndte, kdyn, revised_evp, yield_curve, &
-                                basalstress, Ktens, e_ratio, coriolis
+                                basalstress, Ktens, e_ratio, coriolis, thermo
       use ice_transport_driver, only: advection
       use ice_restoring, only: restore_ice
 #ifdef CESMCOUPLED
@@ -151,7 +151,7 @@ subroutine input_data
 
       namelist /dynamics_nml/ &
         kdyn,           ndte,           revised_evp,    yield_curve,    &
-        advection,      coriolis,                                       &
+        advection,      coriolis,       thermo,                         &
         kstrength,      krdg_partic,    krdg_redist,    mu_rdg,         &
         e_ratio,        Ktens,          Cf,             basalstress
 
@@ -246,6 +246,7 @@ subroutine input_data
       kitd = 1           ! type of itd conversions (0 = delta, 1 = linear)
       kcatbound = 1      ! category boundary formula (0 = old, 1 = new, etc)
       kdyn = 1           ! type of dynamics (1 = evp, 2 = eap)
+      thermo = 1         ! set to 1 for thoermodynamics on, and thermo off (0=thermo off)
       ndtd = 1           ! dynamic time steps per thermodynamic time step
       ndte = 120         ! subcycles per dynamics timestep:  ndte=dt_dyn/dte
       revised_evp = .false.  ! if true, use revised procedure for evp dynamics
@@ -484,6 +485,7 @@ subroutine input_data
       call broadcast_scalar(kitd,               master_task)
       call broadcast_scalar(kcatbound,          master_task)
       call broadcast_scalar(kdyn,               master_task)
+      call broadcast_scalar(thermo,             master_task)
       call broadcast_scalar(ndtd,               master_task)
       call broadcast_scalar(ndte,               master_task)
       call broadcast_scalar(revised_evp,        master_task)
@@ -863,6 +865,7 @@ subroutine input_data
          write(nu_diag,1020) ' kcatbound                 = ', &
                                kcatbound
          write(nu_diag,1020) ' kdyn                      = ', kdyn
+         write(nu_diag,1020) ' thermo                    = ', thermo
          write(nu_diag,1020) ' ndtd                      = ', ndtd
          write(nu_diag,1020) ' ndte                      = ', ndte
          write(nu_diag,1010) ' revised_evp               = ', &

cicecore/drivers/cice/CICE_RunMod.F90

@@ -25,6 +25,8 @@ module CICE_RunMod
       use icepack_intfc, only: icepack_max_aero
       use icepack_intfc, only: icepack_query_parameters
       use icepack_intfc, only: icepack_query_tracer_flags, icepack_query_tracer_numbers
+      use ice_blocks, only: nx_block, ny_block
+      use ice_flux, only: uatm, vatm, strairxT, strairyT, wind, rhoa
 
       implicit none
       private
@@ -136,7 +138,7 @@ subroutine ice_step
       use ice_domain, only: halo_info, nblocks
       use ice_domain_size, only: nslyr
       use ice_dyn_eap, only: write_restart_eap
-      use ice_dyn_shared, only: kdyn
+      use ice_dyn_shared, only: kdyn, thermo
       use ice_flux, only: scale_factor, init_history_therm, &
           daidtt, daidtd, dvidtt, dvidtd, dagedtt, dagedtd
       use ice_history, only: accum_hist
@@ -158,10 +160,12 @@ subroutine ice_step
 
       integer (kind=int_kind) :: &
          iblk        , & ! block index 
+         i,j         , & ! loop counter
          k               ! dynamics supercycling index
 
       real (kind=dbl_kind) :: &
-         offset          ! d(age)/dt time offset
+         offset      , & ! d(age)/dt time offset
+         tau             ! part of wind stress
 
       logical (kind=log_kind) :: &
           tr_iage, tr_FY, tr_lvl, &
@@ -202,7 +206,7 @@ subroutine ice_step
       !-----------------------------------------------------------------
       ! Scale radiation fields
       !-----------------------------------------------------------------
-
+            if(thermo.ne.0) then
             if (calc_Tsfc) call prep_radiation (dt, iblk)
 
       !-----------------------------------------------------------------
@@ -213,6 +217,20 @@ subroutine ice_step
             call biogeochemistry (dt, iblk) ! biogeochemistry
             call step_therm2     (dt, iblk) ! ice thickness distribution thermo
 
+            else  ! for box problem
+            !wind stress
+
+               do j = 1, ny_block  
+               do i = 1, nx_block  
+                   wind(i,j,iblk) = sqrt(uatm(i,j,iblk)**2 + vatm(i,j,iblk)**2)
+                   tau = rhoa(i,j,iblk) * 0.0012_dbl_kind * wind(i,j,iblk)
+                   strairxT(i,j,iblk) = tau * uatm(i,j,iblk)
+                   strairyT(i,j,iblk) = tau * vatm(i,j,iblk)
+            enddo
+            enddo
+
+            endif
+
          enddo ! iblk
          !$OMP END PARALLEL DO
 
@@ -255,7 +273,7 @@ subroutine ice_step
          !$OMP PARALLEL DO PRIVATE(iblk)
          do iblk = 1, nblocks
 
-            call step_radiation (dt, iblk)
+            if (thermo.ne.0) call step_radiation (dt, iblk)
 
       !-----------------------------------------------------------------
       ! get ready for coupling and the next time step

configuration/scripts/ice_in

@@ -95,6 +95,7 @@
 
 &dynamics_nml
     kdyn            = 1
+    thermo          = 1
     ndte            = 120
     revised_evp     = .false.
     advection       = 'remap'

configuration/scripts/options/set_nml.box2001

@@ -4,7 +4,7 @@ ice_ic = 'default'
 restart = .false.
 restart_ext = .true.
 restart_format = 'bin'
-dumpfreq_n = 5
+dumpfreq_n = 1
 histfreq       = 'd','x','x','x','x'
 grid_type    = 'rectangular'
 grid_file = 'grid'
@@ -15,8 +15,8 @@ tr_iage      = .false.
 tr_FY        = .false.
 tr_lvl       = .false.
 tr_pond_lvl  = .false.
-kitd              = 0
-ktherm            = 0
+kitd              = 1
+ktherm            = 2
 kstrength       = 0
 atmbndy         = 'constant'
 atm_data_type   = 'box'
@@ -30,7 +30,8 @@ f_uatm         = 'd'
 f_vatm         = 'd'
 f_uocn         = 'd'
 f_vocn         = 'd'
-coriolis       = 'constant'
+coriolis       = 'zero'
 dxrect         = 16.e5
 dyrect         = 16.e5
-kdyn           = 0
+thermo         = 0
+f_Tsfc         = 'd'