diff --git a/columnphysics/icepack_flux.F90 b/columnphysics/icepack_flux.F90 index fa880545..d5578076 100644 --- a/columnphysics/icepack_flux.F90 +++ b/columnphysics/icepack_flux.F90 @@ -117,10 +117,10 @@ subroutine merge_fluxes (aicen, & fswthrun_vdf, & ! vis dif sw radiation through ice bot (W/m**2) fswthrun_idr, & ! nir dir sw radiation through ice bot (W/m**2) fswthrun_idf, & ! nir dif sw radiation through ice bot (W/m**2) - fswthrun_uvrdr, & ! > 700nm vis uvr dir sw radiation through ice bot (W/m**2) - fswthrun_uvrdf, & ! > 700nm vis uvr dif sw radiation through ice bot (W/m**2) - fswthrun_pardr, & ! 400-700nm vis par dir sw radiation through ice bot (W/m**2) - fswthrun_pardf, & ! 400-700nm vis par dif sw radiation through ice bot (W/m**2) + fswthrun_uvrdr, & ! < 400nm uv dir sw radiation through ice bot (W/m**2) + fswthrun_uvrdf, & ! < 400nm uv dif sw radiation through ice bot (W/m**2) + fswthrun_pardr, & ! 400-700nm par dir sw radiation through ice bot (W/m**2) + fswthrun_pardf, & ! 400-700nm par dif sw radiation through ice bot (W/m**2) Urefn ! air speed reference level (m/s) ! cumulative fluxes @@ -159,10 +159,10 @@ subroutine merge_fluxes (aicen, & fswthru_vdf, & ! vis dif sw radiation through ice bot (W/m**2) fswthru_idr, & ! nir dir sw radiation through ice bot (W/m**2) fswthru_idf, & ! nir dif sw radiation through ice bot (W/m**2) - fswthru_uvrdr, & ! > 700nm vis uvr dir sw radiation through ice bot (W/m**2) - fswthru_uvrdf, & ! > 700nm vis uvr dif sw radiation through ice bot (W/m**2) - fswthru_pardr, & ! 400-700nm vis par dir sw radiation through ice bot (W/m**2) - fswthru_pardf, & ! 400-700nm vis par dif sw radiation through ice bot (W/m**2) + fswthru_uvrdr, & ! < 400nm uv dir sw radiation through ice bot (W/m**2) + fswthru_uvrdf, & ! < 400nm uv dif sw radiation through ice bot (W/m**2) + fswthru_pardr, & ! 400-700nm par dir sw radiation through ice bot (W/m**2) + fswthru_pardf, & ! 400-700nm par dif sw radiation through ice bot (W/m**2) dsnow, & ! change in snow depth (m) Uref ! air speed reference level (m/s) diff --git a/columnphysics/icepack_shortwave.F90 b/columnphysics/icepack_shortwave.F90 index 6ecdb692..af3f0f65 100644 --- a/columnphysics/icepack_shortwave.F90 +++ b/columnphysics/icepack_shortwave.F90 @@ -250,10 +250,10 @@ subroutine shortwave_ccsm3 (aicen, vicen, & swidf ! sw down, near IR, diffuse (W/m^2) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr , & ! sw down, vis uvr dir (W/m^2) - swuvrdf , & ! sw down, vis uvr dif (W/m^2) - swpardr , & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr , & ! sw down, ultraviolet, direct (W/m^2) + swuvrdf , & ! sw down, ultraviolet, diffuse (W/m^2) + swpardr , & ! sw down, PAR, direct (W/m^2) + swpardf ! sw down, PAR, diffuse (W/m^2) ! baseline albedos for ccsm3 shortwave, set in namelist real (kind=dbl_kind), intent(in) :: & @@ -282,10 +282,10 @@ subroutine shortwave_ccsm3 (aicen, vicen, & fswthrun_vdf, & ! vis dif SW through ice to ocean (W m-2) fswthrun_idr, & ! nir dir SW through ice to ocean (W m-2) fswthrun_idf, & ! nir dif SW through ice to ocean (W m-2) - fswthrun_uvrdr,&! vis dir uvr SW through ice to ocean (W m-2) - fswthrun_uvrdf,&! vis dif uvr SW through ice to ocean (W m-2) - fswthrun_pardr,&! vis dir par SW through ice to ocean (W m-2) - fswthrun_pardf ! vis dif par SW through ice to ocean (W m-2) + fswthrun_uvrdr,&! uv dir SW through ice to ocean (W m-2) + fswthrun_uvrdf,&! uv dif SW through ice to ocean (W m-2) + fswthrun_pardr,&! par dir SW through ice to ocean (W m-2) + fswthrun_pardf ! par dif SW through ice to ocean (W m-2) real (kind=dbl_kind), intent(inout) :: & coszen ! cosine(zenith angle) @@ -318,10 +318,10 @@ subroutine shortwave_ccsm3 (aicen, vicen, & l_fswthru_vdf, & ! vis dif SW through ice to ocean (W m-2) l_fswthru_idr, & ! nir dir SW through ice to ocean (W m-2) l_fswthru_idf, & ! nir dif SW through ice to ocean (W m-2) - l_fswthru_uvrdr,&! vis uvr dir SW through ice to ocean (W m-2) - l_fswthru_uvrdf,&! vis uvr dif SW through ice to ocean (W m-2) - l_fswthru_pardr,&! vis par dir SW through ice to ocean (W m-2) - l_fswthru_pardf ! vis par dif SW through ice to ocean (W m-2) + l_fswthru_uvrdr,&! uv dir SW through ice to ocean (W m-2) + l_fswthru_uvrdf,&! uv dif SW through ice to ocean (W m-2) + l_fswthru_pardr,&! par dir SW through ice to ocean (W m-2) + l_fswthru_pardf ! par dif SW through ice to ocean (W m-2) character(len=*),parameter :: subname='(shortwave_ccsm3)' !----------------------------------------------------------------- @@ -739,10 +739,10 @@ subroutine absorbed_solar (aicen, & alidfns ! near-ir, diffuse albedo, snow real (kind=dbl_kind), intent(in), optional :: & - swuvrdr , & ! sw down, vis uvr dir (W/m^2) - swuvrdf , & ! sw down, vis uvr dif (W/m^2) - swpardr , & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr , & ! sw down, uv dir (W/m^2) + swuvrdf , & ! sw down, uv dif (W/m^2) + swpardr , & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) real (kind=dbl_kind), intent(out):: & fswsfc , & ! SW absorbed at ice/snow surface (W m-2) @@ -754,10 +754,10 @@ subroutine absorbed_solar (aicen, & fswthru_vdf , & ! vis dif SW through ice to ocean (W m-2) fswthru_idr , & ! nir dir SW through ice to ocean (W m-2) fswthru_idf , & ! nir dif SW through ice to ocean (W m-2) - fswthru_uvrdr, & ! vis dir uvr SW through ice to ocean (W m-2) - fswthru_uvrdf, & ! vis dif uvr SW through ice to ocean (W m-2) - fswthru_pardr, & ! vis dir par SW through ice to ocean (W m-2) - fswthru_pardf ! vis dif par SW through ice to ocean (W m-2) + fswthru_uvrdr, & ! uv dir SW through ice to ocean (W m-2) + fswthru_uvrdf, & ! uv dif SW through ice to ocean (W m-2) + fswthru_pardr, & ! par dir SW through ice to ocean (W m-2) + fswthru_pardf ! par dif SW through ice to ocean (W m-2) real (kind=dbl_kind), dimension (:), intent(out) :: & Iswabs , & ! SW absorbed in particular layer (W m-2) @@ -788,10 +788,10 @@ subroutine absorbed_solar (aicen, & asnow ! fractional area of snow cover real (kind=dbl_kind) :: & - swuvrdrpen , & ! penetrating SW, vis uvr dir (W/m^2) - swuvrdfpen , & ! penetrating SW, vis uvr dif (W/m^2) - swpardrpen , & ! penetrating SW, vis par dir (W/m^2) - swpardfpen ! penetrating SW, vis par dif (W/m^2) + swuvrdrpen , & ! penetrating SW, uv dir (W/m^2) + swuvrdfpen , & ! penetrating SW, uv dif (W/m^2) + swpardrpen , & ! penetrating SW, par dir (W/m^2) + swpardfpen ! penetrating SW, par dif (W/m^2) character(len=*),parameter :: subname='(absorbed_solar)' @@ -972,10 +972,10 @@ subroutine run_dEdd(dt, & fsnow ! snowfall rate (kg/m^2 s) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr, & ! sw down, vis uvr dir (W/m^2) - swuvrdf, & ! sw down, vis uvr dif (W/m^2) - swpardr, & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr, & ! sw down, uv dir (W/m^2) + swuvrdf, & ! sw down, uv dif (W/m^2) + swpardr, & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) real(kind=dbl_kind), dimension(:), intent(in) :: & aicen, & ! concentration of ice @@ -1017,8 +1017,8 @@ subroutine run_dEdd(dt, & fswthrun_vdf, & ! vis dif SW through ice to ocean (W/m^2) fswthrun_idr, & ! nir dir SW through ice to ocean (W/m^2) fswthrun_idf, & ! nir dif SW through ice to ocean (W/m^2) - fswthrun_uvrdr,&! uvr dir SW through ice to ocean (W/m^2) - fswthrun_uvrdf,&! uvr dif SW through ice to ocean (W/m^2) + fswthrun_uvrdr,&! uv dir SW through ice to ocean (W/m^2) + fswthrun_uvrdf,&! uv dif SW through ice to ocean (W/m^2) fswthrun_pardr,&! par dir SW through ice to ocean (W/m^2) fswthrun_pardf ! par dif SW through ice to ocean (W/m^2) @@ -1062,10 +1062,10 @@ subroutine run_dEdd(dt, & l_fswthru_vdf , & ! vis dif SW through ice to ocean (W m-2) l_fswthru_idr , & ! nir dir SW through ice to ocean (W m-2) l_fswthru_idf , & ! nir dif SW through ice to ocean (W m-2) - l_fswthru_uvrdr, & ! vis uvr dir SW through ice to ocean (W m-2) - l_fswthru_uvrdf, & ! vis uvr dif SW through ice to ocean (W m-2) - l_fswthru_pardr, & ! vis par dir SW through ice to ocean (W m-2) - l_fswthru_pardf ! vis par dif SW through ice to ocean (W m-2) + l_fswthru_uvrdr, & ! uv dir SW through ice to ocean (W m-2) + l_fswthru_uvrdf, & ! uv dif SW through ice to ocean (W m-2) + l_fswthru_pardr, & ! par dir SW through ice to ocean (W m-2) + l_fswthru_pardf ! par dif SW through ice to ocean (W m-2) logical (kind=log_kind) :: & l_initonly ! local initonly value @@ -1087,7 +1087,7 @@ subroutine run_dEdd(dt, & call compute_coszen (TLAT, TLON, yday, sec, coszen, & days_per_year, nextsw_cday, calendar_type) #else - if (.not.semi_implicit_Tsfc) then ! geos sets solar angles in driver level + if (.not.semi_implicit_Tsfc) then ! GEOS sets solar angles in driver level call compute_coszen (TLAT, TLON, yday, sec, coszen) endif #endif @@ -1318,10 +1318,10 @@ subroutine shortwave_dEdd (coszen, & swidf ! sw down, near IR, diffuse (W/m^2) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr , & ! sw down, vis uvr dir (W/m^2) - swuvrdf , & ! sw down, vis uvr dif (W/m^2) - swpardr , & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr , & ! sw down, uv dir (W/m^2) + swuvrdf , & ! sw down, uv dif (W/m^2) + swpardr , & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) real (kind=dbl_kind), intent(inout) :: & coszen , & ! cosine of solar zenith angle @@ -1338,10 +1338,10 @@ subroutine shortwave_dEdd (coszen, & fswthru_vdf , & ! vis dif SW through snow/bare ice/ponded ice into ocean (W m-2) fswthru_idr , & ! nir dir SW through snow/bare ice/ponded ice into ocean (W m-2) fswthru_idf , & ! nir dif SW through snow/bare ice/ponded ice into ocean (W m-2) - fswthru_uvrdr,& ! vis uvr dir sw radiation through ice bot (GEOS) (W/m**2) - fswthru_uvrdf,& ! vis uvr dif sw radiation through ice bot (GEOS) (W/m**2) - fswthru_pardr,& ! vis par dir sw radiation through ice bot (GEOS) (W/m**2) - fswthru_pardf ! vis par dif sw radiation through ice bot (GEOS) (W/m**2) + fswthru_uvrdr,& ! uv dir SW through ice into ocean (W m-2) + fswthru_uvrdf,& ! uv dif SW through ice into ocean (W m-2) + fswthru_pardr,& ! par dir SW through ice into ocean (W m-2) + fswthru_pardf ! par dif SW through ice into ocean (W m-2) real (kind=dbl_kind), dimension (:), intent(inout) :: & fswpenl , & ! visible SW entering ice layers (W m-2) @@ -1731,10 +1731,10 @@ subroutine compute_dEdd_3bd( & fswthru_vdf, & ! vis dif SW through snow/bare ice/ponded ice into ocean (W m-2) fswthru_idr, & ! nir dir SW through snow/bare ice/ponded ice into ocean (W m-2) fswthru_idf, & ! nir dif SW through snow/bare ice/ponded ice into ocean (W m-2) - fswthru_uvrdr,&! vis uvr dir sw radiation through ice bot (GEOS) (W/m**2) - fswthru_uvrdf,&! vis uvr dif sw radiation through ice bot (GEOS) (W/m**2) - fswthru_pardr,&! vis par dir sw radiation through ice bot (GEOS) (W/m**2) - fswthru_pardf ! vis par dif sw radiation through ice bot (GEOS) (W/m**2) + fswthru_uvrdr,&! uv dir SW through ice into ocean (W m-2) + fswthru_uvrdf,&! uv dif SW through ice into ocean (W m-2) + fswthru_pardr,&! par dir SW through ice into ocean (W m-2) + fswthru_pardf ! par dif SW through ice into ocean (W m-2) real (kind=dbl_kind), dimension (:), intent(inout) :: & fswpenl, & ! visible SW entering ice layers (W m-2) @@ -1742,10 +1742,10 @@ subroutine compute_dEdd_3bd( & Iswabs ! SW absorbed in ice layer (W m-2) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr, & ! sw down, vis uvr dir (W/m^2) - swuvrdf, & ! sw down, vis uvr dif (W/m^2) - swpardr, & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr, & ! sw down, uv dir (W/m^2) + swuvrdf, & ! sw down, uv dif (W/m^2) + swpardr, & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) !----------------------------------------------------------------------- ! @@ -1876,10 +1876,10 @@ subroutine compute_dEdd_3bd( & fthruvdf, & ! vis dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) fthruidr, & ! nir dir shortwave through snow/bare ice/ponded ice to ocean (W/m^2) fthruidf, & ! nir dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) - fthruuvrdr,&! vis uvr dir shortwave through snow/bare ice/ponded ice to ocean (W/m^2) - fthruuvrdf,&! vis uvr dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) - fthrupardr,&! vis par dir shortwave through snow/bare ice/ponded ice to ocean (W/m^2) - fthrupardf ! vis par dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) + fthruuvrdr,&! uv dir shortwave through snow/bare ice/ponded ice to ocean (W/m^2) + fthruuvrdf,&! uv dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) + fthrupardr,&! par dir shortwave through snow/bare ice/ponded ice to ocean (W/m^2) + fthrupardf ! par dif shortwave through snow/bare ice/ponded ice to ocean (W/m^2) real (kind=dbl_kind), dimension(nslyr) :: & Sabs ! shortwave absorbed in snow layer (W m-2) @@ -3998,10 +3998,10 @@ subroutine icepack_step_radiation (dt, & TLAT, TLON ! latitude and longitude (radian) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr , & ! sw down, vis uvr dir (W/m^2) - swuvrdf , & ! sw down, vis uvr dif (W/m^2) - swpardr , & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr , & ! sw down, uv dir (W/m^2) + swuvrdf , & ! sw down, uv dif (W/m^2) + swpardr , & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) integer (kind=int_kind), intent(in) :: & sec ! elapsed seconds into date @@ -4062,10 +4062,10 @@ subroutine icepack_step_radiation (dt, & fswthrun_vdf , & ! vis dif SW through ice to ocean (W/m^2) fswthrun_idr , & ! nir dir SW through ice to ocean (W/m^2) fswthrun_idf , & ! nir dif SW through ice to ocean (W/m^2) - fswthrun_uvrdr,& ! vis uvr dir SW through ice to ocean (W/m^2) - fswthrun_uvrdf,& ! vis uvr dif SW through ice to ocean (W/m^2) - fswthrun_pardr,& ! vis par dir SW through ice to ocean (W/m^2) - fswthrun_pardf ! vis par dif SW through ice to ocean (W/m^2) + fswthrun_uvrdr,& ! uv dir SW through ice to ocean (W/m^2) + fswthrun_uvrdf,& ! uv dif SW through ice to ocean (W/m^2) + fswthrun_pardr,& ! par dir SW through ice to ocean (W/m^2) + fswthrun_pardf ! par dif SW through ice to ocean (W/m^2) real (kind=dbl_kind), dimension(:,:), intent(inout) :: & fswpenln , & ! visible SW entering ice layers (W m-2) diff --git a/columnphysics/icepack_therm_vertical.F90 b/columnphysics/icepack_therm_vertical.F90 index 53f082f6..00b3756d 100644 --- a/columnphysics/icepack_therm_vertical.F90 +++ b/columnphysics/icepack_therm_vertical.F90 @@ -2370,10 +2370,10 @@ subroutine icepack_step_therm1(dt, & fswthru_vdf , & ! vis dif shortwave penetrating to ocean (W/m^2) fswthru_idr , & ! nir dir shortwave penetrating to ocean (W/m^2) fswthru_idf , & ! nir dif shortwave penetrating to ocean (W/m^2) - fswthru_uvrdr,& ! vis uvr dir shortwave penetrating to ocean (W/m^2) - fswthru_uvrdf,& ! vis uvr dif shortwave penetrating to ocean (W/m^2) - fswthru_pardr,& ! vis par dir shortwave penetrating to ocean (W/m^2) - fswthru_pardf,& ! vis par dif shortwave penetrating to ocean (W/m^2) + fswthru_uvrdr,& ! uv dir shortwave penetrating to ocean (W/m^2) + fswthru_uvrdf,& ! uv dif shortwave penetrating to ocean (W/m^2) + fswthru_pardr,& ! par dir shortwave penetrating to ocean (W/m^2) + fswthru_pardf,& ! par dif shortwave penetrating to ocean (W/m^2) dsnow , & ! change in snow depth (m/step-->cm/day) fsloss ! rate of snow loss to leads (kg/m^2/s) @@ -2461,10 +2461,10 @@ subroutine icepack_step_therm1(dt, & fswthrun_vdf , & ! vis dif SW through ice to ocean (W/m^2) fswthrun_idr , & ! nir dir SW through ice to ocean (W/m^2) fswthrun_idf , & ! nir dif SW through ice to ocean (W/m^2) - fswthrun_uvrdr,& ! vis uvr dir SW through ice to ocean (W/m^2) - fswthrun_uvrdf,& ! vis uvr dif SW through ice to ocean (W/m^2) - fswthrun_pardr,& ! vis par dir SW through ice to ocean (W/m^2) - fswthrun_pardf ! via par dif SW through ice to ocean (W/m^2) + fswthrun_uvrdr,& ! uv dir SW through ice to ocean (W/m^2) + fswthrun_uvrdf,& ! uv dif SW through ice to ocean (W/m^2) + fswthrun_pardr,& ! par dir SW through ice to ocean (W/m^2) + fswthrun_pardf ! par dif SW through ice to ocean (W/m^2) real (kind=dbl_kind), dimension(:,:), intent(inout) :: & zqsn , & ! snow layer enthalpy (J m-3) @@ -2532,17 +2532,17 @@ subroutine icepack_step_therm1(dt, & l_meltsliqn ! mass of snow melt local (kg/m^2) real (kind=dbl_kind) :: & - l_fswthrun_vdr, & ! vis dir SW local n ice to ocean (W/m^2) - l_fswthrun_vdf, & ! vis dif SW local n ice to ocean (W/m^2) - l_fswthrun_idr, & ! nir dir SW local n ice to ocean (W/m^2) - l_fswthrun_idf, & ! nir dif SW local n ice to ocean (W/m^2) - l_fswthrun_uvrdr,& ! vis uvr dir SW through ice to ocean (W/m^2) - l_fswthrun_uvrdf,& ! vis uvr dif SW through ice to ocean (W/m^2) - l_fswthrun_pardr,& ! vis par dir SW through ice to ocean (W/m^2) - l_fswthrun_pardf,& ! vis par dif SW through ice to ocean (W/m^2) - l_dsnow, & ! local snow change - l_dsnown, & ! local snow change category - l_meltsliq ! mass of snow melt local (kg/m^2) + l_fswthrun_vdr, & ! vis dir SW category ice to ocean (W/m^2) + l_fswthrun_vdf, & ! vis dif SW category ice to ocean (W/m^2) + l_fswthrun_idr, & ! nir dir SW category ice to ocean (W/m^2) + l_fswthrun_idf, & ! nir dif SW category ice to ocean (W/m^2) + l_fswthrun_uvrdr, & ! uv dir SW category ice to ocean (W/m^2) + l_fswthrun_uvrdf, & ! uv dif SW category ice to ocean (W/m^2) + l_fswthrun_pardr, & ! par dir SW category ice to ocean (W/m^2) + l_fswthrun_pardf, & ! par dif SW category ice to ocean (W/m^2) + l_dsnow, & ! local snow change + l_dsnown, & ! local snow change category + l_meltsliq ! mass of snow melt local (kg/m^2) real (kind=dbl_kind) :: & l_dpnd_flushn, & ! category pond flushing rate (m/step) diff --git a/doc/source/icepack_index.rst b/doc/source/icepack_index.rst index 6a46aeb9..38b3a4cf 100755 --- a/doc/source/icepack_index.rst +++ b/doc/source/icepack_index.rst @@ -194,10 +194,10 @@ section :ref:`tabnamelist`. "fswthru_ai", "grid-box-mean shortwave penetrating to ocean (fswthru)", "W/m\ :math:`^2`" "fswthru_idf", "near IR diffuse shortwave penetrating to ocean", "W/m\ :math:`^2`" "fswthru_idr", "near IR direct shortwave penetrating to ocean", "W/m\ :math:`^2`" - "fswthru_pardf", "visible photosynthetically active shortwave diffuse radiation 400-700nm penetrating to ocean", "W/m\ :math:`^2`" - "fswthru_pardr", "visible photosynthetically active shortwave direct radiation 400-700nm penetrating to ocean", "W/m\ :math:`^2`" - "fswthru_uvrdf", "visible ultraviolet shortwave diffuse radiation > 700nm penetrating to ocean", "W/m\ :math:`^2`" - "fswthru_uvrdr", "visible ultraviolet shortwave direct radiation > 700nm penetrating to ocean", "W/m\ :math:`^2`" + "fswthru_pardf", "photosynthetically active diffuse shortwave radiation 400-700nm penetrating to ocean", "W/m\ :math:`^2`" + "fswthru_pardr", "photosynthetically active direct shortwave radiation 400-700nm penetrating to ocean", "W/m\ :math:`^2`" + "fswthru_uvrdf", "ultraviolet diffuse shortwave radiation < 400nm penetrating to ocean", "W/m\ :math:`^2`" + "fswthru_uvrdr", "ultraviolet direct shortwave radiation < 400nm penetrating to ocean", "W/m\ :math:`^2`" "fswthru_vdr", "visible direct shortwave penetrating to ocean", "W/m\ :math:`^2`" "fswthru_vdf", "visible diffuse shortwave penetrating to ocean", "W/m\ :math:`^2`" "fyear", "current data year", "" @@ -254,7 +254,7 @@ section :ref:`tabnamelist`. "istep1", "total number of steps at current time step", "" "Iswabs", "shortwave radiation absorbed in ice layers", "W/m\ :math:`^2`" "itd_area_min", ":math:`\bullet` zap residual ice below a minimum area", "1.e-11" - "itd_mass_min", ":math:`\bullet` zap residual ice below a minimum mass", "1.e-10 + "itd_mass_min", ":math:`\bullet` zap residual ice below a minimum mass", "1.e-10" "**J**", "", "" "**K**", "", "" "kalg", ":math:`\bullet` absorption coefficient for algae", "" diff --git a/doc/source/science_guide/sg_thermo.rst b/doc/source/science_guide/sg_thermo.rst index 874b0653..bf3d3c25 100755 --- a/doc/source/science_guide/sg_thermo.rst +++ b/doc/source/science_guide/sg_thermo.rst @@ -889,19 +889,27 @@ heat flux, :math:`F_{L\downarrow}` is the incoming longwave flux, :math:`F_{L\uparrow}` is the outgoing longwave flux, :math:`F_{sw}` is the incoming shortwave flux, :math:`\alpha` is the shortwave albedo, and :math:`i_0` is the fraction of absorbed shortwave flux that penetrates -into the ice. The albedo may be altered by the presence of melt ponds. +into the ice. + +Two methods for computing albedo and shortwave fluxes are available, the +"CCSM3" method and a multiple scattering radiative transfer scheme that +uses a Delta-Eddington approach. Shortwave forcing input is generally +divided into two wavelength bands (:math:`< 700` \ nm, visible and below) and +(:math:`> 700` \ nm, infrared), with an additional capability to accept +shortwave divided into ultraviolet (:math:`< 400` \ nm) and visible / +photosynthetically active (:math:`400 - 700` \ nm) bands. +Radiation calculations may occur in the original two bands +or for up to five bands, as described below. + +The albedo may be altered by the presence of melt ponds. Each of the explicit melt pond parameterizations (topo, level-ice and sealvl ponds) should be used in conjunction with the Delta-Eddington -shortwave scheme, described below. +shortwave scheme. Shortwave radiation: Delta-Eddington ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ -Two methods for computing albedo and shortwave fluxes are available, the -"ccsm3" method, described in the next section, and a multiple scattering -radiative transfer scheme that uses a Delta-Eddington approach -(``shortwave`` = ``dEdd``). - +In the Delta-Eddington approach (``shortwave`` = ``dEdd``), "Inherent" optical properties (IOPs) for snow and sea ice, such as extinction coefficient and single scattering albedo, are prescribed based on physical measurements; reflected, absorbed and transmitted diff --git a/doc/source/user_guide/interfaces.include b/doc/source/user_guide/interfaces.include index 27e93917..3620e5c1 100644 --- a/doc/source/user_guide/interfaces.include +++ b/doc/source/user_guide/interfaces.include @@ -2061,10 +2061,10 @@ icepack_step_radiation TLAT, TLON ! latitude and longitude (radian) real (kind=dbl_kind), intent(in), optional :: & - swuvrdr , & ! sw down, vis uvr dir (W/m^2) - swuvrdf , & ! sw down, vis uvr dif (W/m^2) - swpardr , & ! sw down, vis par dir (W/m^2) - swpardf ! sw down, vis par dif (W/m^2) + swuvrdr , & ! sw down, uv dir (W/m^2) + swuvrdf , & ! sw down, uv dif (W/m^2) + swpardr , & ! sw down, par dir (W/m^2) + swpardf ! sw down, par dif (W/m^2) integer (kind=int_kind), intent(in) :: & sec ! elapsed seconds into date @@ -2125,10 +2125,10 @@ icepack_step_radiation fswthrun_vdf , & ! vis dif SW through ice to ocean (W/m^2) fswthrun_idr , & ! nir dir SW through ice to ocean (W/m^2) fswthrun_idf , & ! nir dif SW through ice to ocean (W/m^2) - fswthrun_uvrdr,& ! vis uvr dir SW through ice to ocean (W/m^2) - fswthrun_uvrdf,& ! vis uvr dif SW through ice to ocean (W/m^2) - fswthrun_pardr,& ! vis par dir SW through ice to ocean (W/m^2) - fswthrun_pardf ! vis par dif SW through ice to ocean (W/m^2) + fswthrun_uvrdr,& ! uv dir SW through ice to ocean (W/m^2) + fswthrun_uvrdf,& ! uv dif SW through ice to ocean (W/m^2) + fswthrun_pardr,& ! par dir SW through ice to ocean (W/m^2) + fswthrun_pardf ! par dif SW through ice to ocean (W/m^2) real (kind=dbl_kind), dimension(:,:), intent(inout) :: & fswpenln , & ! visible SW entering ice layers (W m-2) @@ -2712,10 +2712,10 @@ icepack_step_therm1 fswthru_vdf , & ! vis dif shortwave penetrating to ocean (W/m^2) fswthru_idr , & ! nir dir shortwave penetrating to ocean (W/m^2) fswthru_idf , & ! nir dif shortwave penetrating to ocean (W/m^2) - fswthru_uvrdr,& ! vis uvr dir shortwave penetrating to ocean (W/m^2) - fswthru_uvrdf,& ! vis uvr dif shortwave penetrating to ocean (W/m^2) - fswthru_pardr,& ! vis par dir shortwave penetrating to ocean (W/m^2) - fswthru_pardf,& ! vis par dif shortwave penetrating to ocean (W/m^2) + fswthru_uvrdr,& ! uv dir shortwave penetrating to ocean (W/m^2) + fswthru_uvrdf,& ! uv dif shortwave penetrating to ocean (W/m^2) + fswthru_pardr,& ! par dir shortwave penetrating to ocean (W/m^2) + fswthru_pardf,& ! par dif shortwave penetrating to ocean (W/m^2) dsnow , & ! change in snow depth (m/step-->cm/day) fsloss ! rate of snow loss to leads (kg/m^2/s) @@ -2803,10 +2803,10 @@ icepack_step_therm1 fswthrun_vdf , & ! vis dif SW through ice to ocean (W/m^2) fswthrun_idr , & ! nir dir SW through ice to ocean (W/m^2) fswthrun_idf , & ! nir dif SW through ice to ocean (W/m^2) - fswthrun_uvrdr,& ! vis uvr dir SW through ice to ocean (W/m^2) - fswthrun_uvrdf,& ! vis uvr dif SW through ice to ocean (W/m^2) - fswthrun_pardr,& ! vis par dir SW through ice to ocean (W/m^2) - fswthrun_pardf ! via par dif SW through ice to ocean (W/m^2) + fswthrun_uvrdr,& ! uv dir SW through ice to ocean (W/m^2) + fswthrun_uvrdf,& ! uv dif SW through ice to ocean (W/m^2) + fswthrun_pardr,& ! par dir SW through ice to ocean (W/m^2) + fswthrun_pardf ! par dif SW through ice to ocean (W/m^2) real (kind=dbl_kind), dimension(:,:), intent(inout) :: & zqsn , & ! snow layer enthalpy (J m-3)