Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 8 additions & 8 deletions columnphysics/icepack_flux.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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)

Expand Down
126 changes: 63 additions & 63 deletions columnphysics/icepack_shortwave.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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) :: &
Expand Down Expand Up @@ -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)
Expand Down Expand Up @@ -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)'

!-----------------------------------------------------------------
Expand Down Expand Up @@ -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)
Expand All @@ -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)
Expand Down Expand Up @@ -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)'

Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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)

Expand Down Expand Up @@ -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
Expand All @@ -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
Expand Down Expand Up @@ -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
Expand All @@ -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)
Expand Down Expand Up @@ -1731,21 +1731,21 @@ 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)
Sswabs , & ! SW absorbed in snow layer (W m-2)
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)

!-----------------------------------------------------------------------
!
Expand Down Expand Up @@ -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)
Expand Down Expand Up @@ -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
Expand Down Expand Up @@ -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)
Expand Down
38 changes: 19 additions & 19 deletions columnphysics/icepack_therm_vertical.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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)

Expand Down Expand Up @@ -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)
Expand Down Expand Up @@ -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)
Expand Down
Loading