Index

Rundeck: E6F40

Created: Fri May 17 03:30:27 EDT 2024

SURF_ALBEDO

File: ALBEDO.f

Summary: SURF_ALBEDO contains parameters/variables needed for albedo calc
Author : A. Lacis/V. Oinas (modifications by I. Aleinov/G. Schmidt)
Version:

---

Subroutines:

get_albedo_data

returns some internal SURF_ALBEDO data. Needed to provide corresponding data to WRITER in RADPAR.

getsur

GETSUR computes surface albedo for each grid box

updsur

UPDSUR updates variables for surface albedo once a day

---

Global Variables:

nkband : integer, parameter

number of K-bands
Initial Value = 33

nv : integer, parameter

total number of vegetation types
Initial Value = 12

srfoam : real*8, parameter, dimension(25)

look up table for ocean foam as a function of wind speed
Initial Value = (/ 0.000,0.000,0.000,0.000,0.001,0.002,0.003,0.005,0.007,0.010, 0.014,0.019,0.025,0.032,0.041,0.051,0.063,0.077,0.094,0.112, 0.138,0.164,0.191,0.218,0.246/)

season : real*8, parameter, dimension(4)

julian day for start of season (used for veg albedo calc)
Initial Value = (/ 15.00, 105.0, 196.0, 288.0/)

albvnd : real*8, parameter, dimension(NV, 4, 6)

veg alb by veg type, season and band
Initial Value = RESHAPE( (/ .500,.067,.089,.089,.078,.100,.067,.061,.089,.000,.200,.089, .500,.062,.100,.100,.073,.055,.067,.061,.100,.000,.200,.100, .500,.085,.091,.139,.085,.058,.083,.061,.091,.000,.200,.091, .500,.080,.090,.111,.064,.055,.061,.061,.090,.000,.200,.090, .500,.200,.267,.267,.233,.300,.200,.183,.267,.000,.200,.267, .500,.206,.350,.300,.241,.218,.200,.183,.350,.000,.200,.350, .500,.297,.364,.417,.297,.288,.250,.183,.364,.000,.200,.364, .500,.255,.315,.333,.204,.218,.183,.183,.315,.000,.200,.315, .500,.200,.267,.267,.233,.300,.200,.183,.267,.000,.200,.267, .500,.206,.350,.300,.241,.218,.200,.183,.350,.000,.200,.350, .500,.297,.364,.417,.297,.288,.250,.183,.364,.000,.200,.364, .500,.255,.315,.333,.204,.218,.183,.183,.315,.000,.200,.315, .500,.080,.107,.107,.093,.120,.080,.073,.107,.000,.200,.107, .500,.082,.140,.120,.096,.083,.080,.073,.140,.000,.200,.140, .500,.119,.145,.167,.119,.115,.100,.073,.145,.000,.200,.145, .500,.102,.126,.132,.081,.087,.073,.073,.126,.000,.200,.126, .500,.100,.133,.133,.116,.150,.100,.091,.133,.000,.200,.133, .500,.103,.175,.150,.120,.109,.100,.091,.175,.000,.200,.175, .500,.148,.182,.208,.148,.144,.125,.091,.182,.000,.200,.182, .500,.127,.157,.166,.102,.109,.091,.091,.157,.000,.200,.157, .500,.020,.027,.027,.023,.030,.020,.018,.027,.000,.200,.027, .500,.021,.035,.030,.024,.022,.020,.018,.035,.000,.200,.035, .500,.030,.036,.042,.030,.029,.025,.018,.036,.000,.200,.036, .500,.026,.032,.033,.020,.022,.018,.018,.032,.000,.200,.032 /),(/NV,4,6/) )

vtmask : real*8, parameter, dimension(NV)

vegetation depth mask by type (m)
Initial Value = (/ .1d0,.2d0,.2d0,.5d0,2d0, 5d0, 10d0,25d0,.2d0,.1d0,1d-5,.2d0 /)

anhzoi : real*8

hemisph.Ice Albedo half-max depth (m) (orig.version)
Initial Value = .1d0

ashzoi : real*8

hemisph.Ice Albedo half-max depth (m) (orig.version)
Initial Value = .1d0

dmoice : real*8

masking depth for snow on sea ice (orig.version)
Initial Value = 10.

dmlice : real*8

masking depth for snow on land ice
Initial Value = 10.

anfoam : real*8

for ocean albedo calc
Initial Value = .1514d0

avfoam : real*8

for ocean albedo calc
Initial Value = .2197d0

anscat : real*8

for ocean albedo calc
Initial Value = 0d0

avscat : real*8

for ocean albedo calc
Initial Value = .0156d0

asnalb : real*8, parameter, dimension(7)

snow albedo for old snow
Initial Value = (/.60d0,.55d0,.55d0,.30d0,.10d0,.05d0, .35d0/)

aoialb : real*8, parameter, dimension(7)

seaice albedo (original version)
Initial Value = (/.55d0,.50d0,.45d0,.25d0,.10d0,.05d0, .30d0/)

alialb : real*8, parameter, dimension(7)

land ice albedo
Initial Value = (/.60d0,.55d0,.50d0,.30d0,.10d0,.05d0, .35d0/)

agexpf : real*8, dimension(3, 2)

exponent in snowage calculation depends on hemi/surf type
Initial Value = RESHAPE( (/ 0.2d0, 0.2d0, 0.2d0, 0.2d0, 0.2d0, 0.2d0 /), (/3,2/) )

albdif : real*8, dimension(3, 2)

difference in albedo as function of snowage
Initial Value = RESHAPE( (/ 0.35d0, 0.35d0, 0.35d0, 0.35d0, 0.35d0, 0.35d0/), (/3,2/) )

aoimax : real*8, dimension(6)

range for seaice albedo
Initial Value = (/ .62d0, .42d0, .30d0, .120d0, .05d0, .03d0/)

aoimin : real*8, dimension(6)

range for seaice albedo
Initial Value = (/ .05d0, .05d0, .05d0, .050d0, .05d0, .03d0/)

asndry : real*8, dimension(6)

wet,dry snow albedo over sea ice
Initial Value = (/ .90d0, .85d0, .65d0, .450d0, .10d0, .10d0/)

asnwet : real*8, dimension(6)

wet,dry snow albedo over sea ice
Initial Value = (/ .85d0, .75d0, .50d0, .175d0, .03d0, .01d0/)

ampmin : real*8, dimension(6)

mininimal melt pond albedo
Initial Value = (/ .10d0, .05d0, .05d0, .050d0, .05d0, .03d0/)

aocean : real*8, parameter, dimension(NKBAND)

K-band dependent Thermal radiation characteristics for ocn
Initial Value = (/ 0.04000,0.09566,0.10273,0.10389,0.10464,0.10555,0.10637, 0.10666,0.10697,0.10665,0.10719,0.10728,0.11007,0.04009, 0.04553,0.05554,0.08178,0.09012,0.09464,0.09548,0.09532, 0.09558,0.09558,0.09568,0.09565,0.05771,0.04985,0.04670, 0.04630,0.04575,0.04474,0.04468,0.04500/)

agsidv : real*8, parameter, dimension(NKBAND, 4)

K-band dependent Thermal radiation for other types
Initial Value = RESHAPE( (/ 0.01400,0.09262,0.09170,0.07767,0.07130,0.06603,0.06540, 0.06397,0.06358,0.06361,0.06365,0.06386,0.06564,0.01354, 0.01537,0.02320,0.04156,0.03702,0.03633,0.03417,0.03346, 0.03342,0.03322,0.03350,0.03170,0.01967,0.01845,0.01977, 0.01986,0.01994,0.02013,0.02041,0.02100, 0.01400,0.09262,0.09170,0.07767,0.07130,0.06603,0.06540, 0.06397,0.06358,0.06361,0.06365,0.06386,0.06564,0.01354, 0.01537,0.02320,0.04156,0.03702,0.03633,0.03417,0.03346, 0.03342,0.03322,0.03350,0.03170,0.01967,0.01845,0.01977, 0.01986,0.01994,0.02013,0.02041,0.02100, 0.04500,0.10209,0.08806,0.05856,0.04835,0.04052,0.04001, 0.03775,0.03687,0.03740,0.03637,0.03692,0.03570,0.07001, 0.05665,0.05326,0.05349,0.04356,0.03845,0.03589,0.03615, 0.03610,0.03602,0.03613,0.03471,0.13687,0.14927,0.16484, 0.16649,0.16820,0.17199,0.17484,0.18000, 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0. /), (/ NKBAND,4 /) )

wetsra : real*8

adjustment factors for wet earth albedo calc
Initial Value = 1.0

wettra : real*8

adjustment factors for wet earth albedo calc
Initial Value = 1.0

zvgsra : real*8

adjustment factors for
Initial Value = 1.0

zdssra : real*8

adjustment factors for
Initial Value = 1.0

zicsra : real*8

adjustment factors for
Initial Value = 1.0

zsnsra : real*8

adjustment factors for
Initial Value = 1.0

zocsra : real*8

adjustment factors for
Initial Value = 1.0

evgtra : real*8

adjustment factors for
Initial Value = 1.0

edstra : real*8

adjustment factors for
Initial Value = 1.0

eictra : real*8

adjustment factors for
Initial Value = 1.0

esntra : real*8

adjustment factors for
Initial Value = 1.0

eoctra : real*8

adjustment factors for
Initial Value = 1.0

albvnh : real*8, dimension(NV, 6, 2)

hemispherically varying vegetation albedo
Used by: | SOIL_DRV:daily_earth |

gzsnow : real*8, parameter, dimension(7, 3, 2)

asymmetry parameter for snow over three types
Initial Value = RESHAPE( (/ 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0, 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0, 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0, 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0, 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0, 0.95d0, 0.94d0, 0.905d0, 0.896d0, 0.894d0, 0.89d0, 0.91d0 /), (/7,3,2/) )