Updated Numerical Solution for Wave Attenuation in sea ice (#1294)

* This PR adds a write statement in ww3_grid.out and therefore will result in changes in all ww3_grid.out files in regtests. 
 
*add IC4_NUMERICS namelist to w3gdatmd.F90

* add IC4_NUMERICS namelist to w3gridmd.F90

* make ICNUMERICS namelist setting (to be put under MISC namelist) - Ic numerics fix should be applicable to ALL IC settings.

* IC_NUMERICS Fix for All sea ice source terms

* remove lines that modify timestep in sea ice

* remove pointer on line 704 w3gdatmd.F90

Author: erinethomas

model/src/w3gdatmd.F90

@@ -240,6 +240,7 @@ MODULE W3GDATMD
   !                               for individual grid points.
   !      IICEDISP   Log.  Public   Flag for use of the ice covered dispertion relation.
   !      IICESMOOTH Log.  Public   Flag to smooth the ice covered dispertion relation in broken ice.
+  !      IC_NUMERICS Log. Public  Turn on/off IC numerics fix
   !
   !
   !      GNAME     C*30  Public   Grid name.
@@ -700,6 +701,7 @@ MODULE W3GDATMD
 
     LOGICAL          :: GINIT, FLDRY, FLCX, FLCY, FLCTH, FLCK, FLSOU, IICEDISP,&
          IICESMOOTH
+    LOGICAL          :: IC_NUMERICS
     LOGICAL          :: FLAGLL
     LOGICAL          :: CMPRTRCK
     LOGICAL, POINTER :: FLAGST(:)
@@ -1216,6 +1218,7 @@ MODULE W3GDATMD
 
   LOGICAL, POINTER :: GINIT, FLDRY, FLCX, FLCY, FLCTH, FLCK, FLSOU, IICEDISP,&
        IICESMOOTH
+  LOGICAL, POINTER :: IC_NUMERICS
   LOGICAL, POINTER :: FLAGLL
   LOGICAL, POINTER :: CMPRTRCK
   LOGICAL, POINTER :: FLAGST(:)
@@ -2385,6 +2388,7 @@ SUBROUTINE W3SETG ( IMOD, NDSE, NDST )
     FLSOU  => GRIDS(IMOD)%FLSOU
     IICEDISP => GRIDS(IMOD)%IICEDISP
     IICESMOOTH => GRIDS(IMOD)%IICESMOOTH
+    IC_NUMERICS => GRIDS(IMOD)%IC_NUMERICS
     !
     GNAME  => GRIDS(IMOD)%GNAME
     FILEXT => GRIDS(IMOD)%FILEXT

model/src/w3gridmd.F90

@@ -809,6 +809,7 @@ MODULE W3GRIDMD
   !
   REAL(8)                 :: GSHIFT ! see notes in WMGHGH
   LOGICAL                 :: FLC, ICEDISP, TRCKCMPR
+  LOGICAL                 :: ICNUMERICS
   INTEGER                 :: PTM   ! Partitioning method
   REAL                    :: PTFC  ! Part. cut off freq (for method 5)
   REAL                    :: AIRCMIN, AIRGB
@@ -1113,7 +1114,7 @@ MODULE W3GRIDMD
        STDX, STDY, STDT, ICEHMIN, ICEHINIT, ICEDISP,   &
        ICESLN, ICEWIND, ICESNL, ICESDS, ICEHFAC,       &
        ICEHDISP, ICEDDISP, ICEFDISP, CALTYPE,          &
-       TRCKCMPR, PTM, PTFC, BTBET
+       TRCKCMPR, PTM, PTFC, BTBET, ICNUMERICS
   NAMELIST /OUTS/ P2SF, I1P2SF, I2P2SF,                      &
        US3D, I1US3D, I2US3D,                    &
        USSP, IUSSP, STK_WN,                     &
@@ -2765,6 +2766,7 @@ SUBROUTINE W3GRID()
     STDY = -1.
     STDT = -1.
     ICEDISP = .FALSE.
+    ICNUMERICS=.FALSE.
     CALTYPE = 'standard'
     ! Variables for 3D array output
     E3D=0
@@ -3035,6 +3037,7 @@ SUBROUTINE W3GRID()
     IICEHDISP  = ICEHDISP
     IICEDDISP  = ICEDDISP
     IICEFDISP  = ICEFDISP
+    IC_NUMERICS=ICNUMERICS
     PMOVE  = MAX ( 0. , PMOVE )
     PFMOVE = PMOVE
     !
@@ -3424,7 +3427,7 @@ SUBROUTINE W3GRID()
              ICEHINIT, ICEDISP, ICEHDISP,          &
              ICESLN, ICEWIND, ICESNL, ICESDS,      &
              ICEDDISP,ICEFDISP, CALTYPE, TRCKCMPR, &
-             BTBETA
+             BTBETA,ICNUMERICS
       ELSE
         WRITE (NDSO,2966) CICE0, CICEN, LICE, PMOVE, XSEED, FLAGTR, &
              XP, XR, XFILT, IHMAX, HSPMIN, WSMULT, &
@@ -3434,7 +3437,7 @@ SUBROUTINE W3GRID()
              ICEHINIT, ICEDISP, ICEHDISP,          &
              ICESLN, ICEWIND, ICESNL, ICESDS,      &
              ICEDDISP, ICEFDISP, CALTYPE, TRCKCMPR,&
-             BTBETA
+             BTBETA,ICNUMERICS
       END IF
       !
 #ifdef W3_FLD1
@@ -6833,7 +6836,7 @@ SUBROUTINE W3GRID()
          ', ICESNL = ',F6.2,', ICESDS = ',F5.2,','/       &
          '        ICEDDISP = ',F5.2,', ICEFDISP = ',F5.2,       &
          ', CALTYPE = ',A8,' , TRCKCMPR = ', L3,','/      &
-         '        BTBET  = ', F6.2, ' /')
+         '        BTBET  = ', F6.2, ', ICNUMERICS =',L3,' /')
     !
 2976 FORMAT ( '  &OUTS P2SF  =',I2,', I1P2SF =',I2,', I2P2SF =',I3,','/&
          '        US3D  =',I2,', I1US3D =',I3,', I2US3D =',I3,','/&