[1529] | 1 | module radcommon_h |
---|
| 2 | use radinc_h, only: L_NSPECTI, L_NSPECTV, L_NGAUSS, NTstar, NTstop, & |
---|
| 3 | naerkind, nsizemax |
---|
[135] | 4 | implicit none |
---|
| 5 | |
---|
| 6 | !----------------------------------------------------------------------C |
---|
| 7 | ! |
---|
| 8 | ! radcommon.h |
---|
[1648] | 9 | !v |
---|
[135] | 10 | !----------------------------------------------------------------------C |
---|
| 11 | ! |
---|
| 12 | ! "Include" grid.h and radinc.h before this file in code that uses |
---|
| 13 | ! some or all of this common data set |
---|
| 14 | ! |
---|
| 15 | ! WNOI - Array of wavenumbers at the spectral interval |
---|
| 16 | ! centers for the infrared. Array is NSPECTI |
---|
| 17 | ! elements long. |
---|
| 18 | ! DWNI - Array of "delta wavenumber", i.e., the width, |
---|
| 19 | ! in wavenumbers (cm^-1) of each IR spectral |
---|
| 20 | ! interval. NSPECTI elements long. |
---|
| 21 | ! WAVEI - Array (NSPECTI elements long) of the wavelenght |
---|
| 22 | ! (in microns) at the center of each IR spectral |
---|
| 23 | ! interval. |
---|
| 24 | ! WNOV - Array of wavenumbers at the spectral interval |
---|
| 25 | ! center for the VISUAL. Array is NSPECTV |
---|
| 26 | ! elements long. |
---|
| 27 | ! DWNV - Array of "delta wavenumber", i.e., the width, |
---|
| 28 | ! in wavenumbers (cm^-1) of each VISUAL spectral |
---|
| 29 | ! interval. NSPECTV elements long. |
---|
| 30 | ! WAVEV - Array (NSPECTV elements long) of the wavelenght |
---|
| 31 | ! (in microns) at the center of each VISUAL spectral |
---|
| 32 | ! interval. |
---|
| 33 | ! STELLARF - Array (NSPECTV elements) of stellar flux (W/M^2) in |
---|
| 34 | ! each spectral interval. Values are for 1 AU, |
---|
| 35 | ! scaled to the planetary distance elsewhere. |
---|
| 36 | ! TAURAY - Array (NSPECTV elements) of the pressure-independent |
---|
| 37 | ! part of Rayleigh scattering optical depth. |
---|
| 38 | ! FZEROI - Fraction of zeros in the IR CO2 k-coefficients, for |
---|
| 39 | ! each temperature, pressure, and spectral interval |
---|
| 40 | ! FZEROV - Fraction of zeros in the VISUAL CO2 k-coefficients, for |
---|
| 41 | ! each temperature, pressure, and spectral interval |
---|
| 42 | ! |
---|
| 43 | ! AEROSOL RADIATIVE OPTICAL CONSTANTS |
---|
| 44 | ! |
---|
| 45 | ! Shortwave |
---|
| 46 | ! ~~~~~~~~~ |
---|
| 47 | ! |
---|
| 48 | ! For the "naerkind" kind of aerosol radiative properties : |
---|
| 49 | ! QVISsQREF : Qext / Qext("longrefvis") |
---|
| 50 | ! omegavis : single scattering albedo |
---|
| 51 | ! gvis : assymetry factor |
---|
| 52 | ! |
---|
| 53 | ! Longwave |
---|
| 54 | ! ~~~~~~~~ |
---|
| 55 | ! |
---|
| 56 | ! For the "naerkind" kind of aerosol radiative properties : |
---|
| 57 | ! QIRsQREF : Qext / Qext("longrefvis") |
---|
| 58 | ! omegaIR : mean single scattering albedo |
---|
| 59 | ! gIR : mean assymetry factor |
---|
| 60 | |
---|
[1315] | 61 | REAL*8 BWNI(L_NSPECTI+1), WNOI(L_NSPECTI), DWNI(L_NSPECTI), WAVEI(L_NSPECTI) !BWNI read by master in setspi |
---|
| 62 | REAL*8 BWNV(L_NSPECTV+1), WNOV(L_NSPECTV), DWNV(L_NSPECTV), WAVEV(L_NSPECTV) !BWNV read by master in setspv |
---|
[1648] | 63 | REAL*8 STELLARF(L_NSPECTV), TAURAY(L_NSPECTV) |
---|
[1315] | 64 | !$OMP THREADPRIVATE(WNOI,DWNI,WAVEI,& |
---|
| 65 | !$OMP WNOV,DWNV,WAVEV,& |
---|
[1648] | 66 | !$OMP STELLARF,TAURAY) |
---|
[135] | 67 | |
---|
| 68 | REAL*8 blami(L_NSPECTI+1) |
---|
| 69 | REAL*8 blamv(L_NSPECTV+1) ! these are needed by suaer.F90 |
---|
[1315] | 70 | !$OMP THREADPRIVATE(blami,blamv) |
---|
[135] | 71 | |
---|
[873] | 72 | !! AS: introduced to avoid doing same computations again for continuum |
---|
[878] | 73 | INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: indi |
---|
| 74 | INTEGER, DIMENSION(:,:,:), ALLOCATABLE :: indv |
---|
[1315] | 75 | !$OMP THREADPRIVATE(indi,indv) |
---|
[873] | 76 | |
---|
[470] | 77 | !!! ALLOCATABLE STUFF SO THAT DIMENSIONS ARE READ in *.dat FILES -- AS 12/2011 |
---|
| 78 | REAL*8, DIMENSION(:,:,:,:,:), ALLOCATABLE :: gasi, gasv |
---|
[1648] | 79 | REAL*8, DIMENSION(:), ALLOCATABLE :: PGASREF, TGASREF, PFGASREF |
---|
[135] | 80 | real*8 FZEROI(L_NSPECTI) |
---|
| 81 | real*8 FZEROV(L_NSPECTV) |
---|
| 82 | real*8 pgasmin, pgasmax |
---|
| 83 | real*8 tgasmin, tgasmax |
---|
[1648] | 84 | !$OMP THREADPRIVATE(gasi,gasv,& !pgasref,tgasref,pfgasref read by master in sugas_corrk |
---|
[1315] | 85 | !$OMP FZEROI,FZEROV) !pgasmin,pgasmax,tgasmin,tgasmax read by master in sugas_corrk |
---|
[135] | 86 | |
---|
| 87 | real QVISsQREF(L_NSPECTV,naerkind,nsizemax) |
---|
| 88 | real omegavis(L_NSPECTV,naerkind,nsizemax) |
---|
| 89 | real gvis(L_NSPECTV,naerkind,nsizemax) |
---|
| 90 | real QIRsQREF(L_NSPECTI,naerkind,nsizemax) |
---|
| 91 | real omegair(L_NSPECTI,naerkind,nsizemax) |
---|
| 92 | real gir(L_NSPECTI,naerkind,nsizemax) |
---|
[1315] | 93 | !$OMP THREADPRIVATE(QVISsQREF,omegavis,gvis,QIRsQREF,omegair,gir) |
---|
[135] | 94 | |
---|
| 95 | |
---|
| 96 | ! Reference wavelengths used to compute reference optical depth (m) |
---|
| 97 | ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
---|
| 98 | |
---|
| 99 | REAL lamrefir(naerkind),lamrefvis(naerkind) |
---|
| 100 | |
---|
| 101 | ! Actual number of grain size classes in each domain for a |
---|
| 102 | ! given aerosol: |
---|
| 103 | |
---|
| 104 | INTEGER :: nsize(naerkind,2) |
---|
| 105 | |
---|
| 106 | ! Particle size axis (depend on the kind of aerosol and the |
---|
| 107 | ! radiation domain) |
---|
| 108 | |
---|
| 109 | DOUBLE PRECISION :: radiustab(naerkind,2,nsizemax) |
---|
[1315] | 110 | !$OMP THREADPRIVATE(lamrefir,lamrefvis,radiustab) !nsize read by suaer_corrk |
---|
[135] | 111 | |
---|
| 112 | ! Extinction coefficient at reference wavelengths; |
---|
| 113 | ! These wavelengths are defined in aeroptproperties, and called |
---|
| 114 | ! longrefvis and longrefir. |
---|
| 115 | |
---|
| 116 | REAL :: QREFvis(naerkind,nsizemax) |
---|
| 117 | REAL :: QREFir(naerkind,nsizemax) |
---|
| 118 | REAL :: omegaREFvis(naerkind,nsizemax) |
---|
| 119 | REAL :: omegaREFir(naerkind,nsizemax) |
---|
| 120 | |
---|
[1529] | 121 | REAL,SAVE :: tstellar ! Stellar brightness temperature (SW) |
---|
[135] | 122 | |
---|
[1529] | 123 | real*8,save :: planckir(L_NSPECTI,NTstop-NTstar+1) |
---|
[135] | 124 | |
---|
[1529] | 125 | real*8,save :: PTOP |
---|
| 126 | real*8,save,allocatable :: TAUREF(:) |
---|
[135] | 127 | |
---|
[1529] | 128 | real*8,parameter :: UBARI = 0.5D0 |
---|
[135] | 129 | |
---|
[1529] | 130 | real*8,save :: gweight(L_NGAUSS) |
---|
[1315] | 131 | !$OMP THREADPRIVATE(QREFvis,QREFir,omegaREFvis,omegaREFir,& ! gweight read by master in sugas_corrk |
---|
| 132 | !$OMP tstellar,planckir,PTOP,TAUREF) |
---|
[135] | 133 | |
---|
| 134 | ! If the gas optical depth (top to the surface) is less than |
---|
| 135 | ! this value, we place that Gauss-point into the "zeros" |
---|
| 136 | ! channel. |
---|
| 137 | real*8, parameter :: TLIMIT = 1.0D-30 |
---|
| 138 | |
---|
| 139 | ! Factor to convert pressures from millibars to Pascals |
---|
| 140 | real*8, parameter :: SCALEP = 1.00D+2 |
---|
| 141 | |
---|
[959] | 142 | real*8, parameter :: sigma = 5.67032D-8 |
---|
[1194] | 143 | real*8, parameter :: grav = 6.672E-11 |
---|
[135] | 144 | |
---|
[1529] | 145 | real*8,save :: Cmk |
---|
| 146 | real*8,save :: glat_ig |
---|
[1315] | 147 | !$OMP THREADPRIVATE(Cmk,glat_ig) |
---|
[135] | 148 | |
---|
[1133] | 149 | ! extinction of incoming sunlight (Saturn's rings, eclipses, etc...) |
---|
[1529] | 150 | REAL, DIMENSION(:), ALLOCATABLE ,SAVE :: eclipse |
---|
[135] | 151 | |
---|
[1194] | 152 | !Latitude-dependent gravity |
---|
[1529] | 153 | REAL, DIMENSION(:), ALLOCATABLE , SAVE :: glat |
---|
[1315] | 154 | !$OMP THREADPRIVATE(glat,eclipse) |
---|
[1194] | 155 | |
---|
[1529] | 156 | contains |
---|
| 157 | |
---|
| 158 | subroutine ini_radcommon_h |
---|
| 159 | use radinc_h, only: L_LEVELS |
---|
| 160 | implicit none |
---|
| 161 | |
---|
| 162 | allocate(TAUREF(L_LEVELS+1)) |
---|
| 163 | |
---|
| 164 | end subroutine ini_radcommon_h |
---|
| 165 | |
---|
| 166 | end module radcommon_h |
---|