Skip to content
Snippets Groups Projects
Commit d76dea77 authored by Debolskiy Andrey's avatar Debolskiy Andrey :bicyclist_tone5:
Browse files

initial commit

parent 3c8e8e19
No related branches found
No related tags found
No related merge requests found
Expand all Show whitespace changes
Inline Side-by-side
Showing
with 3790 additions and 79 deletions
.idea/.gitignore 0 → 100644
+ 8
− 0
View file @ d76dea77
# Default ignored files
/shelf/
/workspace.xml
# Editor-based HTTP Client requests
/httpRequests/
# Datasource local storage ignored files
/dataSources/
/dataSources.local.xml
.idea/misc.xml 0 → 100644
+ 8
− 0
View file @ d76dea77
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="Black">
<option name="sdkName" value="Python 3.8 (untitled)" />
</component>
<component name="JsBuildToolPackageJson" sorting="DEFINITION_ORDER" />
<component name="ProjectRootManager" version="2" project-jdk-name="Python 3.8 (untitled)" project-jdk-type="Python SDK" />
</project>
\ No newline at end of file
.idea/modules.xml 0 → 100644
+ 8
− 0
View file @ d76dea77
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="ProjectModuleManager">
<modules>
<module fileurl="file://$PROJECT_DIR$/.idea/scm-abl.iml" filepath="$PROJECT_DIR$/.idea/scm-abl.iml" />
</modules>
</component>
</project>
\ No newline at end of file
.idea/scm-abl.iml 0 → 100644
+ 8
− 0
View file @ d76dea77
<?xml version="1.0" encoding="UTF-8"?>
<module type="CPP_MODULE" version="4">
<component name="NewModuleRootManager">
<content url="file://$MODULE_DIR$" />
<orderEntry type="inheritedJdk" />
<orderEntry type="sourceFolder" forTests="false" />
</component>
</module>
\ No newline at end of file
.idea/vcs.xml 0 → 100644
+ 6
− 0
View file @ d76dea77
<?xml version="1.0" encoding="UTF-8"?>
<project version="4">
<component name="VcsDirectoryMappings">
<mapping directory="" vcs="Git" />
</component>
</project>
\ No newline at end of file
CMakeLists.txt 0 → 100644
+ 27
− 0
View file @ d76dea77
cmake_minimum_required(VERSION 3.23)
option(BUILD_DOC "Build documentation" OFF)
option(USE_CONFIG_PARSER "Build config parser" ON)
option(USE_NETCDF "Enable netcdf library" ON)
option(WITH_RRTMG "with(out) RRTMG radiation transfer" OFF)
project(scm_abl)
enable_language(Fortran)
enable_language(C CXX)
set(CMAKE_CXX_STANDARD 14)
if(BUILD_DOC)
find_package(Doxygen)
if (DOXYGEN_FOUND)
set(DOXYGEN_IN ${CMAKE_CURRENT_SOURCE_DIR}/doxygen/config) # doxygen config file
# option ALL allows to build the docs together with the code
add_custom_target( pbl_doc ALL
COMMAND ${DOXYGEN_EXECUTABLE} ${DOXYGEN_IN}
WORKING_DIRECTORY ${CMAKE_CURRENT_SOURCE_DIR} ) # WORKING_DIRECTORY: Execute the command with the given current working directory. If it is a relative path it will be interpreted relative to the build tree directory corresponding to the current source directory.
else()
message("Doxygen need to be installed to generate the doxygen documentation")
endif (DOXYGEN_FOUND)
endif(BUILD_DOC)
set(lib_files )
\ No newline at end of file
README.md
+ 17
− 79
View file @ d76dea77
# scm-abl # scm-abl
## Atmospheric signle-column boundary layer routines and tests
this is a library (and a set of sets test cases) for parameterization of ABL vertical mixing, shallow and deep convection in ABL
## Getting started for use in INMCM6 and similar climate models
To make it easy for you to get started with GitLab, here's a list of recommended next steps.
Already a pro? Just edit this README.md and make it your own. Want to make it easy? [Use the template at the bottom](#editing-this-readme)!
## Add your files
- [ ] [Create](https://docs.gitlab.com/ee/user/project/repository/web_editor.html#create-a-file) or [upload](https://docs.gitlab.com/ee/user/project/repository/web_editor.html#upload-a-file) files
- [ ] [Add files using the command line](https://docs.gitlab.com/ee/gitlab-basics/add-file.html#add-a-file-using-the-command-line) or push an existing Git repository with the following command:
```
cd existing_repo
git remote add origin http://tesla.parallel.ru/inmcm60_pbl/scm-abl.git
git branch -M main
git push -uf origin main
```
## Integrate with your tools
- [ ] [Set up project integrations](http://tesla.parallel.ru/inmcm60_pbl/scm-abl/-/settings/integrations)
## Collaborate with your team
- [ ] [Invite team members and collaborators](https://docs.gitlab.com/ee/user/project/members/)
- [ ] [Create a new merge request](https://docs.gitlab.com/ee/user/project/merge_requests/creating_merge_requests.html)
- [ ] [Automatically close issues from merge requests](https://docs.gitlab.com/ee/user/project/issues/managing_issues.html#closing-issues-automatically)
- [ ] [Enable merge request approvals](https://docs.gitlab.com/ee/user/project/merge_requests/approvals/)
- [ ] [Set auto-merge](https://docs.gitlab.com/ee/user/project/merge_requests/merge_when_pipeline_succeeds.html)
## Test and Deploy
Use the built-in continuous integration in GitLab.
- [ ] [Get started with GitLab CI/CD](https://docs.gitlab.com/ee/ci/quick_start/index.html)
- [ ] [Analyze your code for known vulnerabilities with Static Application Security Testing (SAST)](https://docs.gitlab.com/ee/user/application_security/sast/)
- [ ] [Deploy to Kubernetes, Amazon EC2, or Amazon ECS using Auto Deploy](https://docs.gitlab.com/ee/topics/autodevops/requirements.html)
- [ ] [Use pull-based deployments for improved Kubernetes management](https://docs.gitlab.com/ee/user/clusters/agent/)
- [ ] [Set up protected environments](https://docs.gitlab.com/ee/ci/environments/protected_environments.html)
***
# Editing this README
When you're ready to make this README your own, just edit this file and use the handy template below (or feel free to structure it however you want - this is just a starting point!). Thanks to [makeareadme.com](https://www.makeareadme.com/) for this template.
## Suggestions for a good README
Every project is different, so consider which of these sections apply to yours. The sections used in the template are suggestions for most open source projects. Also keep in mind that while a README can be too long and detailed, too long is better than too short. If you think your README is too long, consider utilizing another form of documentation rather than cutting out information.
## Name
Choose a self-explaining name for your project.
## Description
Let people know what your project can do specifically. Provide context and add a link to any reference visitors might be unfamiliar with. A list of Features or a Background subsection can also be added here. If there are alternatives to your project, this is a good place to list differentiating factors.
## Badges
On some READMEs, you may see small images that convey metadata, such as whether or not all the tests are passing for the project. You can use Shields to add some to your README. Many services also have instructions for adding a badge.
## Visuals
Depending on what you are making, it can be a good idea to include screenshots or even a video (you'll frequently see GIFs rather than actual videos). Tools like ttygif can help, but check out Asciinema for a more sophisticated method.
## Installation ## Installation
Within a particular ecosystem, there may be a common way of installing things, such as using Yarn, NuGet, or Homebrew. However, consider the possibility that whoever is reading your README is a novice and would like more guidance. Listing specific steps helps remove ambiguity and gets people to using your project as quickly as possible. If it only runs in a specific context like a particular programming language version or operating system or has dependencies that have to be installed manually, also add a Requirements subsection. This is a Cmake based library, so standard cmake rules apply
list of option and their defaults:
```cmake
option(BUILD_DOC "Build documentation" OFF)
option(USE_CONFIG_PARSER "Build config parser" ON)
option(USE_NETCDF "Enable netcdf library" ON)
option(WITH_RRTMG "with(out) RRTMG radiation transfer" OFF)
```
## Usage ## Usage
Use examples liberally, and show the expected output if you can. It's helpful to have inline the smallest example of usage that you can demonstrate, while providing links to more sophisticated examples if they are too long to reasonably include in the README. see api examples
## Support ## Support
Tell people where they can go to for help. It can be any combination of an issue tracker, a chat room, an email address, etc.
## Roadmap
If you have ideas for releases in the future, it is a good idea to list them in the README.
## Contributing ## Contributing
State if you are open to contributions and what your requirements are for accepting them. Please contribute
For people who want to make changes to your project, it's helpful to have some documentation on how to get started. Perhaps there is a script that they should run or some environment variables that they need to set. Make these steps explicit. These instructions could also be useful to your future self.
You can also document commands to lint the code or run tests. These steps help to ensure high code quality and reduce the likelihood that the changes inadvertently break something. Having instructions for running tests is especially helpful if it requires external setup, such as starting a Selenium server for testing in a browser.
## Authors and acknowledgment ## Authors and acknowledgment
Show your appreciation to those who have contributed to the project.
## License ## License
For open source projects, say how it is licensed. GPL-3
## Project status ## Project status
If you have run out of energy or time for your project, put a note at the top of the README saying that development has slowed down or stopped completely. Someone may choose to fork your project or volunteer to step in as a maintainer or owner, allowing your project to keep going. You can also make an explicit request for maintainers. in development
\ No newline at end of file
doxygen/config 0 → 100644
+ 2855
− 0
View file @ d76dea77
This diff is collapsed.
radiation/CMakeLists.txt 0 → 100644
+ 166
− 0
View file @ d76dea77
cmake_minimum_required(VERSION 3.15)
project(rrtm)
set(CMAKE_CXX_STANDARD 14)
include_directories(modules)
include_directories(src)
enable_language(Fortran)
if ("${CMAKE_Fortran_COMPILER_ID}" STREQUAL "Intel")
set(dialect " -free -qoverride-limits")
set(bounds "-free -check-bounds -qoverride-limits")
elseif("${CMAKE_Fortran_COMPILER_ID}" STREQUAL "GNU")
##Adding compiler options for GNU
set(dialect " -ffree-form -finit-real=zero")
set(bounds "-ffree-form -fbounds-check -fcheck=all -fbacktrace --param=max-vartrack-size=1000000 -override-limits -fno-range-check -finit-real=zero")
elseif("${CMAKE_Fortran_COMPILER_ID}" STREQUAL "IntelLLVM")
##Adding compiler options for Intel OneAPI
set(dialect " -free -qoverride-limits")
set(bounds "-free -check-bounds -qoverride-limits")
else ()
message(SEND_ERROR "Couldn't parse Fortran compiler ID")
endif()
set(CMAKE_Fortran_FLAGS_DEBUG "-O0 ${CMAKE_Fortran_FLAGS_DEBUG} ${bounds}")
set(CMAKE_Fortran_FLAGS_RELEASE "-O1 -finline-functions ${dialect}" )
add_library(rrtmlw STATIC
modules/parkind.f90
modules/parrrtm.f90
modules/rrlw_cld.f90
modules/rrlw_con.f90
modules/rrlw_kg01.f90
modules/rrlw_kg02.f90
modules/rrlw_kg03.f90
modules/rrlw_kg04.f90
modules/rrlw_kg05.f90
modules/rrlw_kg06.f90
modules/rrlw_kg07.f90
modules/rrlw_kg08.f90
modules/rrlw_kg09.f90
modules/rrlw_kg10.f90
modules/rrlw_kg11.f90
modules/rrlw_kg12.f90
modules/rrlw_kg13.f90
modules/rrlw_kg14.f90
modules/rrlw_kg15.f90
modules/rrlw_kg16.f90
modules/rrlw_ref.f90
modules/rrlw_tbl.f90
modules/rrlw_vsn.f90
modules/rrlw_wvn.f90
src/rrtmg_lw_setcoef.f90
src/rrtmg_lw_cldprmc.f90
src/rrtmg_lw_cldprop.f90
src/rrtmg_lw_init.f90
src/rrtmg_lw_k_g.f90
src/rrtmg_lw_rtrn.f90
src/rrtmg_lw_rtrnmr.f90
src/rrtmg_lw_taumol.f90
src/rrtmg_lw_rad.nomcica.f90)
add_library(rrtmsw STATIC
modules/parkind.f90
modules/parrrsw.f90
modules/rrsw_aer.f90
modules/rrsw_cld.f90
modules/rrsw_con.f90
modules/rrsw_kg16.f90
modules/rrsw_kg17.f90
modules/rrsw_kg18.f90
modules/rrsw_kg19.f90
modules/rrsw_kg20.f90
modules/rrsw_kg21.f90
modules/rrsw_kg22.f90
modules/rrsw_kg23.f90
modules/rrsw_kg24.f90
modules/rrsw_kg25.f90
modules/rrsw_kg26.f90
modules/rrsw_kg27.f90
modules/rrsw_kg28.f90
modules/rrsw_kg29.f90
modules/rrsw_ref.f90
modules/rrsw_tbl.f90
modules/rrsw_vsn.f90
modules/rrsw_wvn.f90
src/rrtmg_sw_setcoef.f90
src/rrtmg_sw_cldprmc.f90
src/rrtmg_sw_cldprop.f90
src/rrtmg_sw_init.f90
src/rrtmg_sw_k_g.f90
src/rrtmg_sw_reftra.f90
src/rrtmg_sw_vrtqdr.f90
src/rrtmg_sw_spcvmc.f90
src/rrtmg_sw_spcvrt.f90
src/rrtmg_sw_taumol.f90
src/rrtmg_sw_rad.nomcica.f90)
add_library(rrtm STATIC
modules/parkind.f90
modules/parrrtm.f90
modules/rrlw_cld.f90
modules/rrlw_con.f90
modules/rrlw_kg01.f90
modules/rrlw_kg02.f90
modules/rrlw_kg03.f90
modules/rrlw_kg04.f90
modules/rrlw_kg05.f90
modules/rrlw_kg06.f90
modules/rrlw_kg07.f90
modules/rrlw_kg08.f90
modules/rrlw_kg09.f90
modules/rrlw_kg10.f90
modules/rrlw_kg11.f90
modules/rrlw_kg12.f90
modules/rrlw_kg13.f90
modules/rrlw_kg14.f90
modules/rrlw_kg15.f90
modules/rrlw_kg16.f90
modules/rrlw_ref.f90
modules/rrlw_tbl.f90
modules/rrlw_vsn.f90
modules/rrlw_wvn.f90
src/rrtmg_lw_setcoef.f90
src/rrtmg_lw_cldprmc.f90
src/rrtmg_lw_cldprop.f90
src/rrtmg_lw_init.f90
src/rrtmg_lw_k_g.f90
src/rrtmg_lw_rtrn.f90
src/rrtmg_lw_rtrnmr.f90
src/rrtmg_lw_taumol.f90
src/rrtmg_lw_rad.nomcica.f90
modules/parrrsw.f90
modules/rrsw_aer.f90
modules/rrsw_cld.f90
modules/rrsw_con.f90
modules/rrsw_kg16.f90
modules/rrsw_kg17.f90
modules/rrsw_kg18.f90
modules/rrsw_kg19.f90
modules/rrsw_kg20.f90
modules/rrsw_kg21.f90
modules/rrsw_kg22.f90
modules/rrsw_kg23.f90
modules/rrsw_kg24.f90
modules/rrsw_kg25.f90
modules/rrsw_kg26.f90
modules/rrsw_kg27.f90
modules/rrsw_kg28.f90
modules/rrsw_kg29.f90
modules/rrsw_ref.f90
modules/rrsw_tbl.f90
modules/rrsw_vsn.f90
modules/rrsw_wvn.f90
src/rrtmg_sw_setcoef.f90
src/rrtmg_sw_cldprmc.f90
src/rrtmg_sw_cldprop.f90
src/rrtmg_sw_init.f90
src/rrtmg_sw_k_g.f90
src/rrtmg_sw_reftra.f90
src/rrtmg_sw_vrtqdr.f90
src/rrtmg_sw_spcvmc.f90
src/rrtmg_sw_spcvrt.f90
src/rrtmg_sw_taumol.f90
src/rrtmg_sw_rad.nomcica.f90)
\ No newline at end of file
radiation/data/rrtmg_lw.nc 0 → 100644
+ 0
− 0
View file @ d76dea77
File added
radiation/modules/parkind.f90 0 → 100644
+ 20
− 0
View file @ d76dea77
module parkind
use, intrinsic :: iso_c_binding
implicit none
save
! integer kinds
! ------------
integer, parameter :: kind_ib = c_long ! 8 byte integer
integer, parameter :: kind_im = c_int ! 4 byte integer
!
! real kinds
! ----------
!
integer, parameter :: kind_rb = c_double ! 8 byte real
integer, parameter :: kind_rm = c_float ! 4 byte real
integer, parameter :: kind_bool = c_bool
end module parkind
radiation/modules/parrrsw.f90 0 → 100644
+ 119
− 0
View file @ d76dea77
module parrrsw
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!------------------------------------------------------------------
! rrtmg_sw main parameters
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!------------------------------------------------------------------
! name type purpose
! ----- : ---- : ----------------------------------------------
! mxlay : integer: maximum number of layers
! mg : integer: number of original g-intervals per spectral band
! nbndsw : integer: number of spectral bands
! naerec : integer: number of aerosols (iaer=6, ecmwf aerosol option)
! ngptsw : integer: total number of reduced g-intervals for rrtmg_lw
! ngNN : integer: number of reduced g-intervals per spectral band
! ngsNN : integer: cumulative number of g-intervals per band
!------------------------------------------------------------------
integer(kind=im), parameter :: mxlay = 603 !jplay, klev
integer(kind=im), parameter :: mg = 16 !jpg
integer(kind=im), parameter :: nbndsw = 14 !jpsw, ksw
integer(kind=im), parameter :: naerec = 6 !jpaer
integer(kind=im), parameter :: mxmol = 39
integer(kind=im), parameter :: nstr = 2
integer(kind=im), parameter :: nmol = 7
! Use for 112 g-point model
integer(kind=im), parameter :: ngptsw = 112 !jpgpt
! Use for 224 g-point model
! integer(kind=im), parameter :: ngptsw = 224 !jpgpt
! may need to rename these - from v2.6
integer(kind=im), parameter :: jpband = 29
integer(kind=im), parameter :: jpb1 = 16 !istart
integer(kind=im), parameter :: jpb2 = 29 !iend
integer(kind=im), parameter :: jmcmu = 32
integer(kind=im), parameter :: jmumu = 32
integer(kind=im), parameter :: jmphi = 3
integer(kind=im), parameter :: jmxang = 4
integer(kind=im), parameter :: jmxstr = 16
! ^
! Use for 112 g-point model
integer(kind=im), parameter :: ng16 = 6
integer(kind=im), parameter :: ng17 = 12
integer(kind=im), parameter :: ng18 = 8
integer(kind=im), parameter :: ng19 = 8
integer(kind=im), parameter :: ng20 = 10
integer(kind=im), parameter :: ng21 = 10
integer(kind=im), parameter :: ng22 = 2
integer(kind=im), parameter :: ng23 = 10
integer(kind=im), parameter :: ng24 = 8
integer(kind=im), parameter :: ng25 = 6
integer(kind=im), parameter :: ng26 = 6
integer(kind=im), parameter :: ng27 = 8
integer(kind=im), parameter :: ng28 = 6
integer(kind=im), parameter :: ng29 = 12
integer(kind=im), parameter :: ngs16 = 6
integer(kind=im), parameter :: ngs17 = 18
integer(kind=im), parameter :: ngs18 = 26
integer(kind=im), parameter :: ngs19 = 34
integer(kind=im), parameter :: ngs20 = 44
integer(kind=im), parameter :: ngs21 = 54
integer(kind=im), parameter :: ngs22 = 56
integer(kind=im), parameter :: ngs23 = 66
integer(kind=im), parameter :: ngs24 = 74
integer(kind=im), parameter :: ngs25 = 80
integer(kind=im), parameter :: ngs26 = 86
integer(kind=im), parameter :: ngs27 = 94
integer(kind=im), parameter :: ngs28 = 100
integer(kind=im), parameter :: ngs29 = 112
! Use for 224 g-point model
! integer(kind=im), parameter :: ng16 = 16
! integer(kind=im), parameter :: ng17 = 16
! integer(kind=im), parameter :: ng18 = 16
! integer(kind=im), parameter :: ng19 = 16
! integer(kind=im), parameter :: ng20 = 16
! integer(kind=im), parameter :: ng21 = 16
! integer(kind=im), parameter :: ng22 = 16
! integer(kind=im), parameter :: ng23 = 16
! integer(kind=im), parameter :: ng24 = 16
! integer(kind=im), parameter :: ng25 = 16
! integer(kind=im), parameter :: ng26 = 16
! integer(kind=im), parameter :: ng27 = 16
! integer(kind=im), parameter :: ng28 = 16
! integer(kind=im), parameter :: ng29 = 16
! integer(kind=im), parameter :: ngs16 = 16
! integer(kind=im), parameter :: ngs17 = 32
! integer(kind=im), parameter :: ngs18 = 48
! integer(kind=im), parameter :: ngs19 = 64
! integer(kind=im), parameter :: ngs20 = 80
! integer(kind=im), parameter :: ngs21 = 96
! integer(kind=im), parameter :: ngs22 = 112
! integer(kind=im), parameter :: ngs23 = 128
! integer(kind=im), parameter :: ngs24 = 144
! integer(kind=im), parameter :: ngs25 = 160
! integer(kind=im), parameter :: ngs26 = 176
! integer(kind=im), parameter :: ngs27 = 192
! integer(kind=im), parameter :: ngs28 = 208
! integer(kind=im), parameter :: ngs29 = 224
! Source function solar constant
real(kind=rb), parameter :: rrsw_scon = 1.36822e+03 ! W/m2
end module parrrsw
radiation/modules/parrrtm.f90 0 → 100644
+ 112
− 0
View file @ d76dea77
module parrrtm
use parkind ,only : im => kind_im
implicit none
save
!------------------------------------------------------------------
! rrtmg_lw main parameters
!
! Initial version: JJMorcrette, ECMWF, Jul 1998
! Revised: MJIacono, AER, Jun 2006
! Revised: MJIacono, AER, Aug 2007
! Revised: MJIacono, AER, Aug 2008
!------------------------------------------------------------------
! name type purpose
! ----- : ---- : ----------------------------------------------
! mxlay : integer: maximum number of layers
! mg : integer: number of original g-intervals per spectral band
! nbndlw : integer: number of spectral bands
! maxxsec: integer: maximum number of cross-section molecules
! (e.g. cfcs)
! maxinpx: integer:
! ngptlw : integer: total number of reduced g-intervals for rrtmg_lw
! ngNN : integer: number of reduced g-intervals per spectral band
! ngsNN : integer: cumulative number of g-intervals per band
!------------------------------------------------------------------
integer(kind=im), parameter :: mxlay = 603
integer(kind=im), parameter :: mg = 16
integer(kind=im), parameter :: nbndlw = 16
integer(kind=im), parameter :: maxxsec= 4
integer(kind=im), parameter :: mxmol = 39
integer(kind=im), parameter :: maxinpx= 39
integer(kind=im), parameter :: nmol = 7
! Use for 140 g-point model
integer(kind=im), parameter :: ngptlw = 140
! Use for 256 g-point model
! integer(kind=im), parameter :: ngptlw = 256
! Use for 140 g-point model
integer(kind=im), parameter :: ng1 = 10
integer(kind=im), parameter :: ng2 = 12
integer(kind=im), parameter :: ng3 = 16
integer(kind=im), parameter :: ng4 = 14
integer(kind=im), parameter :: ng5 = 16
integer(kind=im), parameter :: ng6 = 8
integer(kind=im), parameter :: ng7 = 12
integer(kind=im), parameter :: ng8 = 8
integer(kind=im), parameter :: ng9 = 12
integer(kind=im), parameter :: ng10 = 6
integer(kind=im), parameter :: ng11 = 8
integer(kind=im), parameter :: ng12 = 8
integer(kind=im), parameter :: ng13 = 4
integer(kind=im), parameter :: ng14 = 2
integer(kind=im), parameter :: ng15 = 2
integer(kind=im), parameter :: ng16 = 2
integer(kind=im), parameter :: ngs1 = 10
integer(kind=im), parameter :: ngs2 = 22
integer(kind=im), parameter :: ngs3 = 38
integer(kind=im), parameter :: ngs4 = 52
integer(kind=im), parameter :: ngs5 = 68
integer(kind=im), parameter :: ngs6 = 76
integer(kind=im), parameter :: ngs7 = 88
integer(kind=im), parameter :: ngs8 = 96
integer(kind=im), parameter :: ngs9 = 108
integer(kind=im), parameter :: ngs10 = 114
integer(kind=im), parameter :: ngs11 = 122
integer(kind=im), parameter :: ngs12 = 130
integer(kind=im), parameter :: ngs13 = 134
integer(kind=im), parameter :: ngs14 = 136
integer(kind=im), parameter :: ngs15 = 138
! Use for 256 g-point model
! integer(kind=im), parameter :: ng1 = 16
! integer(kind=im), parameter :: ng2 = 16
! integer(kind=im), parameter :: ng3 = 16
! integer(kind=im), parameter :: ng4 = 16
! integer(kind=im), parameter :: ng5 = 16
! integer(kind=im), parameter :: ng6 = 16
! integer(kind=im), parameter :: ng7 = 16
! integer(kind=im), parameter :: ng8 = 16
! integer(kind=im), parameter :: ng9 = 16
! integer(kind=im), parameter :: ng10 = 16
! integer(kind=im), parameter :: ng11 = 16
! integer(kind=im), parameter :: ng12 = 16
! integer(kind=im), parameter :: ng13 = 16
! integer(kind=im), parameter :: ng14 = 16
! integer(kind=im), parameter :: ng15 = 16
! integer(kind=im), parameter :: ng16 = 16
! integer(kind=im), parameter :: ngs1 = 16
! integer(kind=im), parameter :: ngs2 = 32
! integer(kind=im), parameter :: ngs3 = 48
! integer(kind=im), parameter :: ngs4 = 64
! integer(kind=im), parameter :: ngs5 = 80
! integer(kind=im), parameter :: ngs6 = 96
! integer(kind=im), parameter :: ngs7 = 112
! integer(kind=im), parameter :: ngs8 = 128
! integer(kind=im), parameter :: ngs9 = 144
! integer(kind=im), parameter :: ngs10 = 160
! integer(kind=im), parameter :: ngs11 = 176
! integer(kind=im), parameter :: ngs12 = 192
! integer(kind=im), parameter :: ngs13 = 208
! integer(kind=im), parameter :: ngs14 = 224
! integer(kind=im), parameter :: ngs15 = 240
! integer(kind=im), parameter :: ngs16 = 256
end module parrrtm
radiation/modules/rrlw_cld.f90 0 → 100644
+ 35
− 0
View file @ d76dea77
module rrlw_cld
use parkind, only : rb => kind_rb
implicit none
save
!------------------------------------------------------------------
! rrtmg_lw cloud property coefficients
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!------------------------------------------------------------------
! name type purpose
! ----- : ---- : ----------------------------------------------
! abscld1: real :
! absice0: real :
! absice1: real :
! absice2: real :
! absice3: real :
! absliq0: real :
! absliq1: real :
!------------------------------------------------------------------
real(kind=rb) :: abscld1
real(kind=rb) , dimension(2) :: absice0
real(kind=rb) , dimension(2,5) :: absice1
real(kind=rb) , dimension(43,16) :: absice2
real(kind=rb) , dimension(46,16) :: absice3
real(kind=rb) :: absliq0
real(kind=rb) , dimension(58,16) :: absliq1
end module rrlw_cld
radiation/modules/rrlw_con.f90 0 → 100644
+ 42
− 0
View file @ d76dea77
module rrlw_con
use parkind, only : rb => kind_rb
implicit none
save
!------------------------------------------------------------------
! rrtmg_lw constants
! Initial version: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!------------------------------------------------------------------
! name type purpose
! ----- : ---- : ----------------------------------------------
! fluxfac: real : radiance to flux conversion factor
! heatfac: real : flux to heating rate conversion factor
!oneminus: real : 1.-1.e-6
! pi : real : pi
! grav : real : acceleration of gravity
! planck : real : planck constant
! boltz : real : boltzmann constant
! clight : real : speed of light
! avogad : real : avogadro constant
! alosmt : real : loschmidt constant
! gascon : real : molar gas constant
! radcn1 : real : first radiation constant
! radcn2 : real : second radiation constant
! sbcnst : real : stefan-boltzmann constant
! secdy : real : seconds per day
!------------------------------------------------------------------
real(kind=rb) :: fluxfac, heatfac
real(kind=rb) :: oneminus, pi, grav
real(kind=rb) :: planck, boltz, clight
real(kind=rb) :: avogad, alosmt, gascon
real(kind=rb) :: radcn1, radcn2
real(kind=rb) :: sbcnst, secdy
end module rrlw_con
radiation/modules/rrlw_kg01.f90 0 → 100644
+ 74
− 0
View file @ d76dea77
module rrlw_kg01
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!-----------------------------------------------------------------
! rrtmg_lw ORIGINAL abs. coefficients for interval 1
! band 1: 10-250 cm-1 (low - h2o; high - h2o)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefao: real
!fracrefbo: real
! kao : real
! kbo : real
! kao_mn2 : real
! kbo_mn2 : real
! selfrefo: real
! forrefo : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: no1 = 16
real(kind=rb) :: fracrefao(no1) , fracrefbo(no1)
real(kind=rb) :: kao(5,13,no1)
real(kind=rb) :: kbo(5,13:59,no1)
real(kind=rb) :: kao_mn2(19,no1) , kbo_mn2(19,no1)
real(kind=rb) :: selfrefo(10,no1), forrefo(4,no1)
!-----------------------------------------------------------------
! rrtmg_lw COMBINED abs. coefficients for interval 1
! band 1: 10-250 cm-1 (low - h2o; high - h2o)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefa : real
!fracrefb : real
! ka : real
! kb : real
! absa : real
! absb : real
! ka_mn2 : real
! kb_mn2 : real
! selfref : real
! forref : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: ng1 = 10
real(kind=rb) :: fracrefa(ng1) , fracrefb(ng1)
real(kind=rb) :: ka(5,13,ng1) , absa(65,ng1)
real(kind=rb) :: kb(5,13:59,ng1), absb(235,ng1)
real(kind=rb) :: ka_mn2(19,ng1) , kb_mn2(19,ng1)
real(kind=rb) :: selfref(10,ng1), forref(4,ng1)
equivalence (ka(1,1,1),absa(1,1)), (kb(1,13,1),absb(1,1))
end module rrlw_kg01
radiation/modules/rrlw_kg02.f90 0 → 100644
+ 70
− 0
View file @ d76dea77
module rrlw_kg02
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!-----------------------------------------------------------------
! rrtmg_lw ORIGINAL abs. coefficients for interval 2
! band 2: 250-500 cm-1 (low - h2o; high - h2o)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefao: real
!fracrefbo: real
! kao : real
! kbo : real
! selfrefo: real
! forrefo : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: no2 = 16
real(kind=rb) :: fracrefao(no2) , fracrefbo(no2)
real(kind=rb) :: kao(5,13,no2)
real(kind=rb) :: kbo(5,13:59,no2)
real(kind=rb) :: selfrefo(10,no2) , forrefo(4,no2)
!-----------------------------------------------------------------
! rrtmg_lw COMBINED abs. coefficients for interval 2
! band 2: 250-500 cm-1 (low - h2o; high - h2o)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefa : real
!fracrefb : real
! ka : real
! kb : real
! absa : real
! absb : real
! selfref : real
! forref : real
!
! refparam: real
!-----------------------------------------------------------------
integer(kind=im), parameter :: ng2 = 12
real(kind=rb) :: fracrefa(ng2) , fracrefb(ng2)
real(kind=rb) :: ka(5,13,ng2) , absa(65,ng2)
real(kind=rb) :: kb(5,13:59,ng2), absb(235,ng2)
real(kind=rb) :: selfref(10,ng2), forref(4,ng2)
real(kind=rb) :: refparam(13)
equivalence (ka(1,1,1),absa(1,1)),(kb(1,13,1),absb(1,1))
end module rrlw_kg02
radiation/modules/rrlw_kg03.f90 0 → 100644
+ 75
− 0
View file @ d76dea77
module rrlw_kg03
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!-----------------------------------------------------------------
! rrtmg_lw ORIGINAL abs. coefficients for interval 3
! band 3: 500-630 cm-1 (low - h2o,co2; high - h2o,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefao: real
!fracrefbo: real
! kao : real
! kbo : real
! kao_mn2o: real
! kbo_mn2o: real
! selfrefo: real
! forrefo : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: no3 = 16
real(kind=rb) :: fracrefao(no3,9) ,fracrefbo(no3,5)
real(kind=rb) :: kao(9,5,13,no3)
real(kind=rb) :: kbo(5,5,13:59,no3)
real(kind=rb) :: kao_mn2o(9,19,no3), kbo_mn2o(5,19,no3)
real(kind=rb) :: selfrefo(10,no3)
real(kind=rb) :: forrefo(4,no3)
!-----------------------------------------------------------------
! rrtmg_lw COMBINED abs. coefficients for interval 3
! band 3: 500-630 cm-1 (low - h2o,co2; high - h2o,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefa : real
!fracrefb : real
! ka : real
! kb : real
! ka_mn2o : real
! kb_mn2o : real
! selfref : real
! forref : real
!
! absa : real
! absb : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: ng3 = 16
real(kind=rb) :: fracrefa(ng3,9) ,fracrefb(ng3,5)
real(kind=rb) :: ka(9,5,13,ng3) ,absa(585,ng3)
real(kind=rb) :: kb(5,5,13:59,ng3),absb(1175,ng3)
real(kind=rb) :: ka_mn2o(9,19,ng3), kb_mn2o(5,19,ng3)
real(kind=rb) :: selfref(10,ng3)
real(kind=rb) :: forref(4,ng3)
equivalence (ka(1,1,1,1),absa(1,1)),(kb(1,1,13,1),absb(1,1))
end module rrlw_kg03
radiation/modules/rrlw_kg04.f90 0 → 100644
+ 64
− 0
View file @ d76dea77
module rrlw_kg04
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!-----------------------------------------------------------------
! rrtmg_lw ORIGINAL abs. coefficients for interval 4
! band 4: 630-700 cm-1 (low - h2o,co2; high - o3,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefao: real
!fracrefbo: real
! kao : real
! kbo : real
! selfrefo: real
! forrefo : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: no4 = 16
real(kind=rb) :: fracrefao(no4,9) ,fracrefbo(no4,5)
real(kind=rb) :: kao(9,5,13,no4)
real(kind=rb) :: kbo(5,5,13:59,no4)
real(kind=rb) :: selfrefo(10,no4) ,forrefo(4,no4)
!-----------------------------------------------------------------
! rrtmg_lw COMBINED abs. coefficients for interval 4
! band 4: 630-700 cm-1 (low - h2o,co2; high - o3,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
! absa : real
! absb : real
!fracrefa : real
!fracrefb : real
! ka : real
! kb : real
! selfref : real
! forref : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: ng4 = 14
real(kind=rb) :: fracrefa(ng4,9) ,fracrefb(ng4,5)
real(kind=rb) :: ka(9,5,13,ng4) ,absa(585,ng4)
real(kind=rb) :: kb(5,5,13:59,ng4),absb(1175,ng4)
real(kind=rb) :: selfref(10,ng4) ,forref(4,ng4)
equivalence (ka(1,1,1,1),absa(1,1)),(kb(1,1,13,1),absb(1,1))
end module rrlw_kg04
radiation/modules/rrlw_kg05.f90 0 → 100644
+ 76
− 0
View file @ d76dea77
module rrlw_kg05
use parkind ,only : im => kind_im, rb => kind_rb
implicit none
save
!-----------------------------------------------------------------
! rrtmg_lw ORIGINAL abs. coefficients for interval 5
! band 5: 700-820 cm-1 (low - h2o,co2; high - o3,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefao: real
!fracrefbo: real
! kao : real
! kbo : real
! kao_mo3 : real
! selfrefo: real
! forrefo : real
! ccl4o : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: no5 = 16
real(kind=rb) :: fracrefao(no5,9) ,fracrefbo(no5,5)
real(kind=rb) :: kao(9,5,13,no5)
real(kind=rb) :: kbo(5,5,13:59,no5)
real(kind=rb) :: kao_mo3(9,19,no5)
real(kind=rb) :: selfrefo(10,no5)
real(kind=rb) :: forrefo(4,no5)
real(kind=rb) :: ccl4o(no5)
!-----------------------------------------------------------------
! rrtmg_lw COMBINED abs. coefficients for interval 5
! band 5: 700-820 cm-1 (low - h2o,co2; high - o3,co2)
!
! Initial version: JJMorcrette, ECMWF, jul1998
! Revised: MJIacono, AER, jun2006
! Revised: MJIacono, AER, aug2008
!-----------------------------------------------------------------
!
! name type purpose
! ---- : ---- : ---------------------------------------------
!fracrefa : real
!fracrefb : real
! ka : real
! kb : real
! ka_mo3 : real
! selfref : real
! forref : real
! ccl4 : real
!
! absa : real
! absb : real
!-----------------------------------------------------------------
integer(kind=im), parameter :: ng5 = 16
real(kind=rb) :: fracrefa(ng5,9) ,fracrefb(ng5,5)
real(kind=rb) :: ka(9,5,13,ng5) ,absa(585,ng5)
real(kind=rb) :: kb(5,5,13:59,ng5),absb(1175,ng5)
real(kind=rb) :: ka_mo3(9,19,ng5)
real(kind=rb) :: selfref(10,ng5)
real(kind=rb) :: forref(4,ng5)
real(kind=rb) :: ccl4(ng5)
equivalence (ka(1,1,1,1),absa(1,1)),(kb(1,1,13,1),absb(1,1))
end module rrlw_kg05
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment