Index

Rundeck: P1SoM40

Created: Wed May 8 03:30:18 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

number of K-bands (33)
Used by: | RADPAR:rcomp1 |

ntalb : integer

number of albedo types in AGSIDV (4)
Used by: | RADPAR:rcomp1 |

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 (kg/m^2)
Initial Value = (/ 1d1, 2d1, 2d1, 5d1, 2d2, 5d2, 1d3, 25d2,2d1, 1d1,.001d0,2d1 /)

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, allocatable, dimension(:)

K-band dependent Thermal radiation characteristics for ocn
Used by: | RADPAR:rcomp1 |

agsidv : real*8, allocatable, dimension(:, :)

K-band dependent Thermal radiation for other types
Used by: | RADPAR:rcomp1 |

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/) )