Index | Rundeck: E6TomaF40 | Created: Fri May 17 03:30:39 EDT 2024 |
SURF_ALBEDO | File: ALBEDO.f |
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/) )