| 1 | |
|---|
| 2 | module surfdat_h |
|---|
| 3 | |
|---|
| 4 | implicit none |
|---|
| 5 | |
|---|
| 6 | real,allocatable,dimension(:) :: n2frac ! n2 fraction of ice at the surfacestocked in startfi.nc file. |
|---|
| 7 | !$OMP THREADPRIVATE(n2frac) |
|---|
| 8 | |
|---|
| 9 | real,allocatable,dimension(:) :: albedodat ! albedo of bare ground stocked in startfi.nc file. |
|---|
| 10 | !$OMP THREADPRIVATE(albedodat) |
|---|
| 11 | ! Ehouarn: moved inertiedat to comsoil.h |
|---|
| 12 | ! real inertiedat, ! thermal inertia |
|---|
| 13 | real,allocatable,dimension(:) :: phisfi ! geopotential at ground level |
|---|
| 14 | !$OMP THREADPRIVATE(phisfi) |
|---|
| 15 | real,allocatable,dimension(:) :: phisfibed ! (paleo) geopotential of the bedrock (= phisfi-qsurf/1000*g) |
|---|
| 16 | !$OMP THREADPRIVATE(phisfibed) |
|---|
| 17 | real,allocatable,dimension(:,:) :: qsurfyear ! (paleo) kg.m-2 averaged mass of ice lost/gained in the last year |
|---|
| 18 | !$OMP THREADPRIVATE(qsurfyear) |
|---|
| 19 | real,dimension(2) :: emisice ! ice emissivity; 1:Northern hemisphere 2:Southern hemisphere |
|---|
| 20 | real emissiv |
|---|
| 21 | real,dimension(2) :: iceradius, dtemisice |
|---|
| 22 | !$OMP THREADPRIVATE(emisice,emissiv,iceradius,dtemisice) |
|---|
| 23 | real,allocatable,dimension(:) :: zmea,zstd,zsig,zgam,zthe |
|---|
| 24 | !$OMP THREADPRIVATE(zmea,zstd,zsig,zgam,zthe) |
|---|
| 25 | real phitop ! heatflux at top of atmosphere? for Triton !AF24 |
|---|
| 26 | real ttop |
|---|
| 27 | real,allocatable,dimension(:) :: kp ! TB ref pressure |
|---|
| 28 | real p00 |
|---|
| 29 | !$OMP THREADPRIVATE(phitop,ttop,kp,p00) |
|---|
| 30 | ! surface properties ! TB16 |
|---|
| 31 | real alb_n2b,alb_n2a,alb_ch4,alb_co,alb_tho,emis_n2b,emis_n2a |
|---|
| 32 | !$OMP THREADPRIVATE(alb_n2b,alb_n2a,alb_ch4,alb_co,alb_tho,emis_n2b,emis_n2a) |
|---|
| 33 | real emis_ch4,emis_co,emis_tho,emis_tho_eq,alb_tho_eq,alb_ch4_eq |
|---|
| 34 | !$OMP THREADPRIVATE(emis_ch4,emis_co,emis_tho,emis_tho_eq,alb_tho_eq,alb_ch4_eq) |
|---|
| 35 | real ITN2,ITCH4,ITH2O,ITN2d,ITCH4d,ITH2Od,alb_ch4_s,albspe,emispe |
|---|
| 36 | !$OMP THREADPRIVATE(ITN2,ITCH4,ITH2O,ITN2d,ITCH4d,ITH2Od,alb_ch4_s,albspe,emispe) |
|---|
| 37 | real alb_tho_spe,emis_tho_spe |
|---|
| 38 | !$OMP THREADPRIVATE(alb_tho_spe,emis_tho_spe) |
|---|
| 39 | |
|---|
| 40 | real,allocatable,dimension(:) :: dryness !"Dryness coefficient" for grnd water ice sublimation |
|---|
| 41 | ! AS: previously in tracer.h. it is more logical here. |
|---|
| 42 | !$OMP THREADPRIVATE(dryness) |
|---|
| 43 | |
|---|
| 44 | end module surfdat_h |
|---|
| 45 | |
|---|