Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (3 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad
- Files:
-
- 109 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/drhook/yomhook_dummy.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation 8 8 ! nor does it submit to any jurisdiction. 9 ! 9 10 10 ! DR_HOOK is a profiling and debugging system for the IFS, and should 11 11 ! be called at the beginning and end of each subroutine. This is a -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver.F90
r5158 r5159 1 1 ! ecrad_driver.F90 - Driver for offline ECRAD radiation scheme 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! ECRAD is the radiation scheme used in the ECMWF Integrated 16 16 ! Forecasting System in cycle 43R3 and later. Several solvers are -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver_config.F90
r5158 r5159 1 1 ! ecrad_driver_config.F90 - Configure driver for offline ecRad radiation scheme 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_driver_read_input.F90
r5158 r5159 1 1 ! ecrad_driver_read_input.F90 - Read input structures from NetCDF file 2 ! 2 3 3 ! (C) Copyright 2018- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_ifs_driver.F90
r5158 r5159 1 1 ! ecrad_ifs_driver.F90 - Driver for offline ECRAD radiation scheme 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! ECRAD is the radiation scheme used in the ECMWF Integrated 16 16 ! Forecasting System in cycle 43R3 and later. Several solvers are … … 26 26 ! 2) Name of a NetCDF file containing one or more atmospheric profiles 27 27 ! 3) Name of output NetCDF file 28 ! 28 29 29 ! This version uses the infrastructure of the IFS, such as computing 30 30 ! effective radius and cloud overlap from latitude and other … … 36 36 ! satisfying this requirement is test/ifs/ecrad_meridian.nc in the 37 37 ! ecRad package. 38 ! 38 39 39 ! Note that the purpose of this file is simply to demonstrate the use 40 40 ! of the setup_radiation_scheme and radiation_scheme routines; all the -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ecrad_ifs_driver_blocked.F90
r5158 r5159 1 1 ! ecrad_ifs_driver_blocked.F90 - Driver for offline ECRAD radiation scheme 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! ECRAD is the radiation scheme used in the ECMWF Integrated 16 16 ! Forecasting System in cycle 43R3 and later. Several solvers are … … 26 26 ! 2) Name of a NetCDF file containing one or more atmospheric profiles 27 27 ! 3) Name of output NetCDF file 28 ! 28 29 29 ! This version uses the infrastructure of the IFS, such as computing 30 30 ! effective radius and cloud overlap from latitude and other … … 36 36 ! satisfying this requirement is test/ifs/ecrad_meridian.nc in the 37 37 ! ecRad package. 38 ! 38 39 39 ! Note that the purpose of this file is simply to demonstrate the use 40 40 ! of the setup_radiation_scheme and radiation_scheme routines as well -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/ifs_blocking.F90
r5158 r5159 1 1 ! ifs_blocking.F90 - Reshuffle ecRad data into an NPROMA-blocked data structure 2 ! 2 3 3 ! (C) Copyright 2022- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Balthasar Reuter 13 13 ! Email: balthasar.reuter@ecmwf.int … … 78 78 endif 79 79 80 ! 80 81 81 ! RADINTG 82 82 ! … … 344 344 345 345 ! ------------------------------------------------------- 346 ! 346 347 347 ! INPUT LOOP 348 ! 348 349 349 ! ------------------------------------------------------- 350 350 … … 524 524 525 525 ! ------------------------------------------------------- 526 ! 526 527 527 ! OUTPUT LOOP 528 ! 528 529 529 ! ------------------------------------------------------- 530 530 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_aerosol_optics_description.F90
r4773 r5159 1 1 ! (C) Copyright 2022- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation 8 8 ! nor does it submit to any jurisdiction. 9 ! 9 10 10 ! Author: Robin Hogan 11 11 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_cloud_generator.F90
r5158 r5159 1 1 ! (C) Copyright 2017- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/driver/test_spartacus_math.F90
r5158 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/cloud_overlap_decorr_len.F90
r4773 r5159 4 4 5 5 ! CLOUD_OVERLAP_DECORR_LEN 6 ! 6 7 7 ! (C) Copyright 2016- ECMWF. 8 ! 8 9 9 ! This software is licensed under the terms of the Apache Licence Version 2.0 10 10 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 11 ! 11 12 12 ! In applying this licence, ECMWF does not waive the privileges and immunities 13 13 ! granted to it by virtue of its status as an intergovernmental organisation 14 14 ! nor does it submit to any jurisdiction. 15 ! 15 16 16 ! PURPOSE 17 17 ! ------- 18 18 ! Calculate the cloud overlap decorrelation length as a function of 19 19 ! latitude for use in the radiation scheme 20 ! 20 21 21 ! INTERFACE 22 22 ! --------- 23 23 ! CLOUD_OVERLAP_DECORR_LEN is called from RADLSWR and RADIATION_SCHEME 24 ! 24 25 25 ! AUTHOR 26 26 ! ------ 27 27 ! Robin Hogan, ECMWF (using code extracted from radlswr.F90) 28 28 ! Original: 2016-02-16 29 ! 29 30 30 ! MODIFICATIONS 31 31 ! ------------- 32 ! 32 33 33 ! ------------------------------------------------------------------- 34 34 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/cos_sza.F90
r4773 r5159 2 2 3 3 !**** *COS_SZA* 4 ! 4 5 5 ! (C) Copyright 2015- ECMWF. 6 ! 6 7 7 ! This software is licensed under the terms of the Apache Licence Version 2.0 8 8 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 9 ! 9 10 10 ! In applying this licence, ECMWF does not waive the privileges and immunities 11 11 ! granted to it by virtue of its status as an intergovernmental organisation 12 12 ! nor does it submit to any jurisdiction. 13 ! 13 14 14 ! Purpose. 15 15 ! -------- … … 64 64 ! ---------- 65 65 ! ECMWF Research Department documentation of the IFS 66 ! 66 67 67 ! See also: Zhou, L., M. Zhang, Q. Bao, and Y. Liu (2015), On 68 68 ! the incident solar radiation in CMIP5 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/easy_netcdf_read_mpi.F90
r4773 r5159 1 1 ! easy_netcdf_read_mpi.f90 - Read netcdf file on one task and share with other tasks 2 ! 2 3 3 ! (C) Copyright 2017- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/fcttim.func.h
r4773 r5159 1 1 ! (C) Copyright 1989- ECMWF. 2 2 ! (C) Copyright 1989- Meteo-France. 3 ! 3 4 4 ! This software is licensed under the terms of the Apache Licence Version 2.0 5 5 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 6 ! 6 7 7 ! In applying this licence, ECMWF does not waive the privileges and immunities 8 8 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/fcttre.func.h
r5158 r5159 1 1 ! (C) Copyright 1988- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/ice_effective_radius.F90
r4773 r5159 5 5 6 6 ! ICE_EFFECTIVE_RADIUS 7 ! 7 8 8 ! (C) Copyright 2016- ECMWF. 9 ! 9 10 10 ! This software is licensed under the terms of the Apache Licence Version 2.0 11 11 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 12 ! 12 13 13 ! In applying this licence, ECMWF does not waive the privileges and immunities 14 14 ! granted to it by virtue of its status as an intergovernmental organisation 15 15 ! nor does it submit to any jurisdiction. 16 ! 16 17 17 ! PURPOSE 18 18 ! ------- 19 19 ! Calculate effective radius of ice clouds 20 ! 20 21 21 ! AUTHOR 22 22 ! ------ 23 23 ! Robin Hogan, ECMWF (using code extracted from radlswr.F90) 24 24 ! Original: 2016-02-24 25 ! 25 26 26 ! MODIFICATIONS 27 27 ! ------------- 28 ! 29 ! 28 29 30 30 ! ------------------------------------------------------------------- 31 31 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/liquid_effective_radius.F90
r4773 r5159 6 6 7 7 ! LIQUID_EFFECTIVE_RADIUS 8 ! 8 9 9 ! (C) Copyright 2015- ECMWF. 10 ! 10 11 11 ! This software is licensed under the terms of the Apache Licence Version 2.0 12 12 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 13 ! 13 14 14 ! In applying this licence, ECMWF does not waive the privileges and immunities 15 15 ! granted to it by virtue of its status as an intergovernmental organisation 16 16 ! nor does it submit to any jurisdiction. 17 ! 17 18 18 ! PURPOSE 19 19 ! ------- 20 20 ! Calculate effective radius of liquid clouds 21 ! 21 22 22 ! AUTHOR 23 23 ! ------ 24 24 ! Robin Hogan, ECMWF (using code extracted from radlswr.F90) 25 25 ! Original: 2015-09-24 26 ! 26 27 27 ! MODIFICATIONS 28 28 ! ------------- 29 ! 30 ! 29 30 31 31 ! ------------------------------------------------------------------- 32 32 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/radiation_scheme.F90
r4773 r5159 19 19 20 20 ! RADIATION_SCHEME - Interface to modular radiation scheme 21 ! 21 22 22 ! (C) Copyright 2015- ECMWF. 23 ! 23 24 24 ! This software is licensed under the terms of the Apache Licence Version 2.0 25 25 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 26 ! 26 27 27 ! In applying this licence, ECMWF does not waive the privileges and immunities 28 28 ! granted to it by virtue of its status as an intergovernmental organisation 29 29 ! nor does it submit to any jurisdiction. 30 ! 30 31 31 ! PURPOSE 32 32 ! ------- … … 36 36 ! it to the radiation scheme. It returns net fluxes and surface 37 37 ! flux components needed by the rest of the model. 38 ! 38 39 39 ! Lower case is used for variables and types taken from the 40 40 ! radiation library 41 ! 41 42 42 ! INTERFACE 43 43 ! --------- … … 45 45 ! SETUP_RADIATION_SCHEME routine (in the RADIATION_SETUP module) 46 46 ! populates the YRADIATION object, and should have been run first. 47 ! 47 48 48 ! AUTHOR 49 49 ! ------ 50 50 ! Robin Hogan, ECMWF 51 51 ! Original: 2015-09-16 52 ! 52 53 53 ! MODIFICATIONS 54 54 ! ------------- … … 62 62 ! 2019-02-04 R. Hogan Pass out surface longwave downwelling in each emissivity interval 63 63 ! 2019-02-07 R. Hogan SPARTACUS cloud size from PARAM_CLOUD_EFFECTIVE_SEPARATION_ETA 64 ! 64 65 65 !----------------------------------------------------------------------- 66 66 … … 293 293 ! differences between the effective temperature of this lowest layer 294 294 ! and the true value in the model. 295 ! 295 296 296 ! We may approximate the temperature profile in the lowest model level 297 297 ! as piecewise linear between the top of the layer T[k-1/2], the -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/radiation_setup.F90
r4773 r5159 2 2 3 3 ! RADIATION_SETUP - Setting up modular radiation scheme 4 ! 4 5 5 ! (C) Copyright 2015- ECMWF. 6 ! 6 7 7 ! This software is licensed under the terms of the Apache Licence Version 2.0 8 8 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 9 ! 9 10 10 ! In applying this licence, ECMWF does not waive the privileges and immunities 11 11 ! granted to it by virtue of its status as an intergovernmental organisation 12 12 ! nor does it submit to any jurisdiction. 13 ! 13 14 14 ! PURPOSE 15 15 ! ------- … … 19 19 ! radiation scheme, plus a small number of additional variables 20 20 ! needed for its implemenation in the IFS. 21 ! 21 22 22 ! INTERFACE 23 23 ! --------- … … 25 25 ! scheme is actually run using the RADIATION_SCHEME routine (not in 26 26 ! this module). 27 ! 27 28 28 ! AUTHOR 29 29 ! ------ 30 30 ! Robin Hogan, ECMWF 31 31 ! Original: 2015-09-16 32 ! 32 33 33 ! MODIFICATIONS 34 34 ! ------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/satur.F90
r4876 r5159 1 1 ! (C) Copyright 1996- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoe_spectral_planck.F90
r4773 r5159 1 1 ! (C) Copyright 2019- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation … … 11 11 12 12 ! YOE_SPECTRAL_PLANCK 13 ! 13 14 14 ! PURPOSE 15 15 ! ------- … … 19 19 ! spectral emissivity on the high-resolution model grid (rather than 20 20 ! the lower resolution grid seen by the radiation scheme). 21 ! 21 22 22 ! INTERFACE 23 23 ! --------- … … 26 26 ! temperatures at which the Planck function is required. FREE then 27 27 ! deallocates memory. 28 ! 28 29 29 ! AUTHOR 30 30 ! ------ 31 31 ! Robin Hogan, ECMWF 32 32 ! Original: 2019-02-04 33 ! 33 34 34 ! MODIFICATIONS 35 35 ! ------------- … … 214 214 ENDIF 215 215 ! Planck's law 216 ! 216 217 217 ! The exponential term is computed in double precision to avoid 218 218 ! overflow. The final result should still be in the range of a single -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoecld.F90
r4773 r5159 1 1 ! (C) Copyright 1989- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoephy.F90
r4853 r5159 1 1 ! (C) Copyright 1991- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoerad.F90
r4773 r5159 139 139 ! the following only available in newer modular radiation scheme: 140 140 ! 4 = SW/LW SOCRATES scheme 141 ! 141 142 142 ! LONEWSW: LOGICAL : .T. IF NEW SW CODE IS ACTIVE 143 143 ! LECSRAD: LOGICAL : .T. IF CLEAR-SKY RADIATION IS ARCHIVED AS PEXTR2 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoerdu.F90
r4773 r5159 1 1 ! (C) Copyright 1989- ECMWF. 2 2 ! (C) Copyright 1989- Meteo-France. 3 ! 3 4 4 ! This software is licensed under the terms of the Apache Licence Version 2.0 5 5 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 6 ! 6 7 7 ! In applying this licence, ECMWF does not waive the privileges and immunities 8 8 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yoethf.F90
r4876 r5159 1 1 ! (C) Copyright 1988- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifs/yomrip.F90
r4773 r5159 1 1 ! (C) Copyright 1989- ECMWF. 2 2 ! (C) Copyright 1989- Meteo-France. 3 ! 3 4 4 ! This software is licensed under the terms of the Apache Licence Version 2.0 5 5 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 6 ! 6 7 7 ! In applying this licence, ECMWF does not waive the privileges and immunities 8 8 ! granted to it by virtue of its status as an intergovernmental organisation … … 158 158 ! REAL(KIND=JPRB) :: RCOVSRLU 159 159 ! REAL(KIND=JPRB) :: RSIVSRLU 160 ! 160 161 161 ! !! things put here because their values are time-interpolated, i.e. non-constant 162 162 ! ! defined in YOEOZOC -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/abor1.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/mpl_module.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/parkind1.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation … … 9 9 10 10 MODULE PARKIND1 11 ! 11 12 12 ! *** Define usual kinds for strong typing *** 13 ! 13 14 14 IMPLICIT NONE 15 15 PUBLIC 16 16 SAVE 17 ! 17 18 18 ! Integer Kinds 19 19 ! ------------- 20 ! 20 21 21 INTEGER, PARAMETER :: JPIT = SELECTED_INT_KIND(2) 22 22 INTEGER, PARAMETER :: JPIS = SELECTED_INT_KIND(4) … … 32 32 #endif 33 33 34 ! 34 35 35 ! Real Kinds 36 36 ! ---------- 37 ! 37 38 38 INTEGER, PARAMETER :: JPRT = SELECTED_REAL_KIND(2,1) 39 39 INTEGER, PARAMETER :: JPRS = SELECTED_REAL_KIND(4,2) -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yomcst.F90
r4853 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yomdyncore.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yomlun.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yomlun_ifsaux.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yommp0.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsaux/yomtag.F90
r4773 r5159 1 1 ! (C) Copyright 2014- ECMWF. 2 ! 2 3 3 ! This software is licensed under the terms of the Apache Licence Version 2.0 4 4 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 5 ! 5 6 6 ! In applying this licence, ECMWF does not waive the privileges and immunities 7 7 ! granted to it by virtue of its status as an intergovernmental organisation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/modify_wv_continuum.F90
r4773 r5159 2 2 3 3 ! MODIFY_WV_CONTINUUM - Adjust the shortwave continuum coefficients 4 ! 4 5 5 ! PURPOSE 6 6 ! ------- … … 12 12 ! coefficients to match the CAVIAR continuum. Further details on the 13 13 ! impact were provided by Hogan et al. (2017, ECMWF Tech. Memo. 816). 14 ! 14 15 15 ! INTERFACE 16 16 ! --------- … … 19 19 ! argument is 1 then it implements the CAVIAR continuum by scaling 20 20 ! coefficients within the relevant SRTM modules. 21 ! 21 22 22 ! AUTHOR 23 23 ! ------ 24 24 ! Robin Hogan, ECMWF 25 25 ! Original: 2018-02-21 26 ! 26 27 27 ! MODIFICATIONS 28 28 ! ------------- 29 ! 29 30 30 ! ----------------------------------------------------------------------- 31 31 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/parsrtm.F90
r4773 r5159 70 70 71 71 !INTEGER(KIND=JPIM), PARAMETER :: JPGPT = 56 72 ! 72 73 73 !INTEGER(KIND=JPIM), PARAMETER :: NG16 = 3 74 74 !INTEGER(KIND=JPIM), PARAMETER :: NG17 = 6 … … 88 88 !-- configuration with 14 spectral intervals 89 89 ! and a total of 112 g-points (14xvariable number) 90 ! 90 91 91 !INTEGER(KIND=JPIM), PARAMETER :: JPGPT = 112 92 ! 92 93 93 !INTEGER(KIND=JPIM), PARAMETER :: NG16 = 6 94 94 !INTEGER(KIND=JPIM), PARAMETER :: NG17 = 12 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/rrtm_taumol1.F90
r4773 r5159 133 133 ! band 1: 10-350 cm-1 (low key - h2o; low minor - n2) 134 134 ! (high key - h2o; high minor - n2) 135 ! 135 136 136 ! note: previous versions of rrtm band 1: 137 137 ! 10-250 cm-1 (low - h2o; high - h2o) -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/rrtm_taumol2.F90
r4773 r5159 18 18 !********* 19 19 ! band 2: 350-500 cm-1 (low key - h2o; high key - h2o) 20 ! 20 21 21 ! note: previous version of rrtm band 2: 22 22 ! 250 - 500 cm-1 (low - h2o; high - h2o) 23 ! 23 24 24 ! --------------------------------------------------------------------------- 25 25 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/srtm_gas_optical_depth.F90
r4773 r5159 47 47 ! ADAPTED FROM SRTM_SPCVRT_MCICA (BY JEAN-JACQUES MORCRETTE) BY 48 48 ! ROBIN HOGAN 49 ! 49 50 50 ! MODIFICATIONS. 51 51 ! -------------- -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/surrtftr.F90
r4773 r5159 3 3 ! Adapted from Eli J. Mlawer, Atmospheric & Environmental Research. 4 4 ! by JJMorcrette, ECMWF 5 ! 5 6 6 ! JJMorcrette 20110613 flexible number of g-points 7 7 ! ABozzo 201306 updated to rrtmg v4.85 8 ! 8 9 9 ! ------------------------------------------------------------------ 10 10 … … 21 21 ! & NGS1, NGS2, NGS3, NGS4, NGS5, NGS6, NGS7, NGS8,& 22 22 ! & NGS9, NGS10,NGS11,NGS12,NGS13,NGS14,NGS15,NGS16 23 ! 23 24 24 ! ------------------------------------------------------------------ 25 25 … … 230 230 231 231 !------------------------------------------------------------------------------- 232 ! 232 233 233 WT( :) = (/& 234 234 & 0.1527534276_JPRB,0.1491729617_JPRB,0.1420961469_JPRB,0.1316886544_JPRB,& -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/yoerrtftr.F90
r4773 r5159 27 27 !INTEGER(KIND=JPIM) :: NGC(JPBAND) 28 28 !INTEGER(KIND=JPIM) :: NGS(JPBAND) 29 ! 29 30 30 !INTEGER(KIND=JPIM) :: NGN(JPGMAX) 31 31 !INTEGER(KIND=JPIM) :: NGB(JPGMAX) 32 ! 32 33 33 !INTEGER(KIND=JPIM) :: NGM(JPG*JPBAND) 34 34 !REAL(KIND=JPRB) :: WT(JPG) -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/ifsrrtm/yoesrtwn.F90
r4773 r5159 56 56 ! NGM : INTEGER : The index of each new g-point relative to the 57 57 ! original 16 g-points for each band. 58 ! 58 59 59 ! WT : REAL : RRTM weights for 16 g-points. 60 60 ! WTSUM : REAL : Sum of the weights -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/aeropt_5wv_ecrad.F90
r5158 r5159 1 1 !aeropt_5wv_ecrad.F90 2022-09-20 A. Idelkadi et O. Boucher 2 ! 2 3 3 SUBROUTINE AEROPT_5WV_ECRAD(istartcol,iendcol,istartlev,iendlev, & 4 4 config,thermodynamics,aerosol) -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/calcul_cloud_overlap_decorr_len.F90
r5133 r5159 18 18 ! Ld=mid_decorrelation_length si pres entre [440,680] 19 19 ! Ld=high_decorrelation_length si pres < 440hPa 20 ! 20 21 21 ! TO DO 22 22 ! Que faire de PDECORR_LEN_WATER_M ? -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/radiation_scheme_mod.F90
r5133 r5159 687 687 688 688 ! RADIATION_SCHEME - Interface to modular radiation scheme 689 ! 689 690 690 ! (C) Copyright 2015- ECMWF. 691 ! 691 692 692 ! This software is licensed under the terms of the Apache Licence Version 2.0 693 693 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 694 ! 694 695 695 ! In applying this licence, ECMWF does not waive the privileges and immunities 696 696 ! granted to it by virtue of its status as an intergovernmental organisation 697 697 ! nor does it submit to any jurisdiction. 698 ! 698 699 699 ! PURPOSE 700 700 ! ------- … … 704 704 ! it to the radiation scheme. It returns net fluxes and surface 705 705 ! flux components needed by the rest of the model. 706 ! 706 707 707 ! Lower case is used for variables and types taken from the 708 708 ! radiation library 709 ! 709 710 710 ! INTERFACE 711 711 ! --------- … … 713 713 ! SETUP_RADIATION_SCHEME routine (in the RADIATION_SETUP module) 714 714 ! should have been run first. 715 ! 715 716 716 ! AUTHOR 717 717 ! ------ 718 718 ! Robin Hogan, ECMWF 719 719 ! Original: 2015-09-16 720 ! 720 721 721 ! MODIFICATIONS 722 722 ! ------------- 723 ! 723 724 724 ! TO DO 725 725 ! ----- 726 ! 726 727 727 !----------------------------------------------------------------------- 728 728 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/radiation_setup.F90
r5133 r5159 2 2 3 3 ! RADIATION_SETUP - Setting up modular radiation scheme 4 ! 4 5 5 ! AUTHOR 6 6 ! ------ 7 7 ! Robin Hogan, ECMWF 8 8 ! Original: 2015-09-16 9 ! 9 10 10 ! MODIFICATIONS 11 11 ! ------------- 12 12 ! Abderrahmane Idelkadi LMD, juillet 2023 13 ! 13 14 14 !----------------------------------------------------------------------- 15 15 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/readaerosol_optic_ecrad.F90
r5154 r5159 1 1 ! $Id: readaerosol_optic_ecrad.F90 2 ! 2 3 3 SUBROUTINE readaerosol_optic_ecrad(debut, aerosol_couple, ok_alw, ok_volcan, & 4 4 flag_aerosol, flag_bc_internal_mixture, itap, rjourvrai, & … … 102 102 ! 103 103 !**************************************************************************************** 104 ! 105 ! 104 105 106 106 IF (aerosol_couple) THEN !--we get aerosols from tr_seri array from INCA 107 ! 107 108 108 !--copy fields from INCA tr_seri 109 109 !--convert to ug m-3 unit for consistency with offline fields 110 ! 110 111 111 itr = 0 112 112 DO iq = 1,nqtot … … 145 145 nitrcoarse(:,:) = tr_seri(:,:,id_CSNO3M) *zrho(:,:)*1.e9 ! CSNO3M 146 146 nitrinscoarse(:,:)= tr_seri(:,:,id_CINO3M) *zrho(:,:)*1.e9 ! CINO3M 147 ! 147 148 148 bcsol_pi(:,:) = 0.0 ! ASBCM pre-ind 149 149 pomsol_pi(:,:) = 0.0 ! ASPOMM pre-ind … … 159 159 nitrcoarse_pi(:,:) = 0.0 ! CSNO3M pre-ind 160 160 nitrinscoarse_pi(:,:)= 0.0 ! CINO3M 161 ! 161 162 162 ELSE !--not aerosol_couple 163 ! 163 164 164 ! Read and interpolate sulfate 165 165 IF ( flag_aerosol .EQ. 1 .OR. flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN … … 222 222 loaddust=0. 223 223 ENDIF 224 ! 224 225 225 ! Read and interpolate asno3m, csno3m, cino3m 226 226 IF (flag_aerosol .EQ. 6 .OR. flag_aerosol .EQ. 7 ) THEN … … 240 240 loadno3(:)=0.0 241 241 ENDIF 242 ! 242 243 243 ! CSSO4M is set to 0 as not reliable 244 244 sulfcoarse(:,:) = 0.0 ! CSSO4M (=SO4) + CSMSAM (=MSA) … … 247 247 ENDIF !--not aerosol_couple 248 248 249 ! 249 250 250 ! Store all aerosols mixing ratios in one variable for radiation scheme (unit kg/kg for ECRAD) 251 251 ! present-day values … … 281 281 m_allaer_pi(:,:,id_STRAT_phy) = 0.0 282 282 283 ! 283 284 284 ! Calculate the total mass of all soluble aersosols (in unit ug /m3) 285 285 ! to be revisited for AR6 … … 289 289 !**************************************************************************************** 290 290 ! 2) Calculate optical properties for the aerosols 291 ! 291 292 292 !**************************************************************************************** 293 293 DO k = 1, klev -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/lmdz/readaerosolstrato_ecrad.F90
r5154 r5159 1 ! 1 2 2 ! $Id: readaerosolstrato_ecrad.F90 tlurton $ 3 ! 3 4 4 SUBROUTINE readaerosolstrato_ecrad(config, debut, ok_volcan) 5 5 … … 438 438 439 439 IF (.NOT. ok_volcan) THEN 440 ! 440 441 441 !--this is the default case 442 442 !--stratospheric aerosols are added to both index 2 and 1 for double radiation calls … … 464 464 ENDWHERE 465 465 ENDDO 466 ! 466 467 467 ELSE 468 ! 468 469 469 !--this is the VOLMIP case 470 470 !--stratospheric aerosols are only added to index 2 in this case … … 523 523 ENDWHERE 524 524 ENDDO 525 ! 525 526 526 ELSE 527 ! 527 528 528 ! mod ThL 529 529 !--this is the VOLMIP case -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/ecrad_config.h
r4859 r5159 1 1 ! ecrad_config.h - Preprocessor definitions to configure compilation ecRad -*- f90 -*- 2 ! 2 3 3 ! (C) Copyright 2023- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! This file should be included in Fortran source files that require 16 16 ! different optimizations or settings for different architectures and -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_adding_ica_lw.F90
r5158 r5159 1 1 ! radiation_adding_ica_lw.F90 - Longwave adding method in independent column approximation 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive emission/albedo rather than planck/emissivity -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_adding_ica_sw.F90
r5158 r5159 1 1 ! radiation_adding_ica_sw.F90 - Shortwave adding method in independent column approximation 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-10-23 R. Hogan Renamed single-character variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_aerosol.F90
r4773 r5159 1 1 ! radiation_aerosol.F90 - Derived type describing aerosol 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2018-04-15 R. Hogan Add "direct" option -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_aerosol_optics.F90
r5158 r5159 1 1 ! radiation_aerosol_optics.F90 - Computing aerosol optical properties 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2018-04-15 R. Hogan Add "direct" option -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_aerosol_optics_data.F90
r5158 r5159 1 1 ! radiation_aerosol_optics_data.F90 - Type to store aerosol optical properties 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-10-23 R. Hogan Renamed single-character variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_aerosol_optics_description.F90
r5158 r5159 1 1 ! radiation_aerosol_optics_description.F90 - Type to store aerosol optics metadata 2 ! 2 3 3 ! (C) Copyright 2022- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_check.F90
r4773 r5159 1 1 ! radiation_check.F90 - Checking routines 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloud.F90
r5158 r5159 1 1 ! radiation_cloud.F90 - Derived type to store cloud/precip properties 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2019-01-14 R. Hogan Added inv_inhom_effective_size variable -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloud_cover.F90
r5158 r5159 1 1 ! radiation_cloud_cover.F90 - Compute cumulative cloud cover for McICA 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Generate profiles of the cumulative cloud cover as seen from TOA, 16 16 ! used in the McICA cloud generator. 17 ! 17 18 18 ! Modifications 19 19 ! 2020-10-07 R. Hogan Ensure iobj1 initialized in case of alpha_obj==0 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloud_generator.F90
r5158 r5159 1 1 ! radiation_cloud_generator.F90 - Generate water-content or optical-depth scalings for McICA 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Generate clouds for McICA using a method modified from Raisanen et 16 16 ! al. (2002) 17 ! 17 18 18 ! Modifications 19 19 ! 2018-02-22 R. Hogan Call masked version of PDF sampler for speed -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloud_optics.F90
r5158 r5159 1 1 ! radiation_cloud_optics.F90 - Computing cloud optical properties 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-07-22 R. Hogan Added Yi et al. ice optics model -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloud_optics_data.F90
r4853 r5159 1 1 ! radiation_cloud_optics_data.F90 - Type to store cloud optical properties 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloudless_lw.F90
r5158 r5159 1 1 ! radiation_cloudless_lw.F90 - Longwave homogeneous cloudless solver 2 ! 2 3 3 ! (C) Copyright 2019- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_cloudless_sw.F90
r5158 r5159 1 1 ! radiation_cloudless_sw.F90 - Shortwave homogeneous cloudless solver 2 ! 2 3 3 ! (C) Copyright 2019- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_config.F90
r5158 r5159 1 1 ! radiation_config.F90 - Derived type to configure the radiation scheme 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-07-22 R. Hogan Added Yi et al. ice optics model … … 28 28 ! 2020-05-18 R. Hogan Moved out_of_bounds_* to radiation_check.F90 29 29 ! 2021-07-04 R. Hogan Numerous changes for ecCKD and general cloud/aerosol optics 30 ! 30 31 31 ! Note: The aim is for ecRad in the IFS to be as similar as possible 32 32 ! to the offline version, so if you make any changes to this or any … … 354 354 ! SPARTACUS solver? This essentially means that in a situation 355 355 ! like this 356 ! 356 357 357 ! 000111 358 358 ! 222222 359 ! 359 360 360 ! Radiation downwelling from region 1 may be reflected back into 361 361 ! region 0 due to some degree of homogenization of the radiation -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_constants.F90
r4773 r5159 1 1 ! radiation_constants.F90 - Constants used in radiation calculations 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_delta_eddington.h
r5158 r5159 1 1 ! radiation_delta_eddington.h - Delta-Eddington scaling -*- f90 -*- 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! This file is intended to be included inside a module to ensure that 16 16 ! these simple functions may be inlined -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ecckd.F90
r5158 r5159 1 1 ! radiation_ecckd.F90 - ecCKD generalized gas optics model 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ecckd_gas.F90
r4853 r5159 1 1 ! radiation_ecckd_gas.F90 - type representing a single ecCKD gas 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ecckd_interface.F90
r5158 r5159 1 1 ! radiation_ecckd_interface.F90 - Interface to ecCKD gas optics model 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_flux.F90
r5158 r5159 1 1 ! radiation_flux.F90 - Derived type to store the output fluxes 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-09-08 R. Hogan Store g-point fluxes -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_gas.F90
r5158 r5159 1 1 ! radiation_gas.F90 - Derived type to store the gas mixing ratios 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2019-01-14 R. Hogan Added out_of_physical_bounds routine -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_gas_constants.F90
r4773 r5159 1 1 ! radiation_gas_constants.F90 - Molar mases and ID codes of the various gases 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_general_cloud_optics.F90
r5158 r5159 1 1 ! radiation_general_cloud_optics.F90 - Computing generalized cloud optical properties 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_general_cloud_optics_data.F90
r5158 r5159 1 1 ! radiation_general_cloud_optics_data.F90 - Type to store generalized cloud optical properties 2 ! 2 3 3 ! (C) Copyright 2019- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_homogeneous_lw.F90
r5158 r5159 1 1 ! radiation_homogeneous_lw.F90 - Longwave homogeneous-column (no cloud fraction) solver 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive emission/albedo rather than planck/emissivity -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_homogeneous_sw.F90
r5158 r5159 1 1 ! radiation_homogeneous_sw.F90 - Shortwave homogeneous-column (no cloud fraction) solver 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive albedos at g-points -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ice_optics_baran.F90
r4773 r5159 1 1 ! radiation_ice_optics_fu.F90 - Scheme for ice optical properties adapted from Baran's data 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ice_optics_baran2016.F90
r4773 r5159 1 1 ! radiation_ice_optics_baran2016.F90 - Baran et al. (2016) scheme for ice optical properties 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ice_optics_baran2017.F90
r4773 r5159 1 1 ! radiation_ice_optics_baran2017.F90 - 2017 parameterization of Baran's ice optical properties 2 ! 2 3 3 ! (C) Copyright 2017- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ice_optics_fu.F90
r5158 r5159 1 1 ! radiation_ice_optics_fu.F90 - Fu's scheme for ice optical properties 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2020-08-10 R. Hogan Bounded re to be <= 100um and g to be < 1.0 -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ice_optics_yi.F90
r4773 r5159 1 1 ! radiation_ice_optics_yi.F90 - Yi et al. (2013) ice optical properties 2 ! 2 3 3 ! (C) Copyright 2017- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Authors: Mark Fielding and Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! The reference for this ice optics parameterization is Yi, B., 16 16 ! P. Yang, B.A. Baum, T. L'Ecuyer, L. Oreopoulos, E.J. Mlawer, -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_ifs_rrtm.F90
r5158 r5159 1 1 ! radiation_ifs_rrtm.F90 - Interface to IFS implementation of RRTM-G 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive "surface" dummy argument -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_interface.F90
r4853 r5159 1 1 ! radiation_interface.F90 - Public interface to radiation scheme 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Changes to enable generalized surface description 17 17 ! 2017-09-08 R. Hogan Reverted some changes 18 ! 18 19 19 ! To use the radiation scheme, create a configuration_type object, 20 20 ! call "setup_radiation" on it once to load the look-up-tables and -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_liquid_optics_slingo.F90
r4773 r5159 1 1 ! radiation_liquid_optics_slingo.F90 - Slingo SW & Lindner-Li LW parameterization of liquid droplet optics 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_liquid_optics_socrates.F90
r5158 r5159 1 1 ! radiation_liquid_optics_socrates.F90 - SOCRATES method for parameterizing liquid droplet optics 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2020-08-10 R. Hogan Bounded re to be >=1.2um and <=50um -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_lw_derivatives.F90
r5158 r5159 1 1 ! radiation_lw_derivatives.F90 - Compute longwave derivatives for Hogan and Bozzo (2015) method 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! This module provides routines to compute the rate of change of 16 16 ! broadband upwelling longwave flux at each half level with respect to … … 22 22 ! skin temperature, following the method of Hogan and Bozzo (JAMES 23 23 ! 2015). Separate routines are provided for each solver. 24 ! 24 25 25 ! Note that currently a more approximate calculation is performed from 26 26 ! the exact one in Hogan and Bozzo (2015); here we assume that a … … 28 28 ! when in reality there is a change in shape of the Planck function in 29 29 ! addition to an overall increase in the total emission. 30 ! 30 31 31 ! Modifications 32 32 ! 2017-10-23 R. Hogan Renamed single-character variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_matrix.F90
r5158 r5159 1 1 ! radiation_matrix.F90 - SPARTACUS matrix operations 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2018-10-15 R. Hogan Added fast_expm_exchange_[23] 17 ! 17 18 18 ! This module provides the neccessary mathematical functions for the 19 19 ! SPARTACUS radiation scheme: matrix multiplication, matrix solvers -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_mcica_lw.F90
r5158 r5159 1 1 ! radiation_mcica_lw.F90 - Monte-Carlo Independent Column Approximation longtwave solver 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive emission/albedo rather than planck/emissivity -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_mcica_sw.F90
r5158 r5159 1 1 ! radiation_mcica_sw.F90 - Monte-Carlo Independent Column Approximation shortwave solver 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive albedos at g-points -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_monochromatic.F90
r5158 r5159 1 1 ! radiation_interface.F90 - Monochromatic gas/cloud optics for testing 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive "surface" dummy argument -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_optical_depth_scaling.h
r4773 r5159 1 1 ! radiation_optical_depth_scaling.h - Cloud optical-depth scaling for Tripleclouds 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 12 ! 11 12 13 13 ! Author: Robin Hogan 14 14 ! Email: r.j.hogan@ecmwf.int 15 ! 15 16 16 ! Modifications 17 17 ! 2017-07-14 R. Hogan Incorporate gamma distribution option 18 ! 18 19 19 ! This file is intended to be included inside a module to ensure that 20 20 ! this simple routine may be inlined -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_overlap.F90
r5158 r5159 1 1 ! radiation_overlap.F90 - Module to compute cloud overlap quantities 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-10-23 R. Hogan Renamed single-character variables -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_pdf_sampler.F90
r5158 r5159 1 1 ! radiation_pdf_sampler.F90 - Get samples from a PDF for McICA 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_random_numbers.F90
r5158 r5159 1 1 ! radiation_random_numbers.F90 - Generate random numbers for McICA solver 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! The derived type "rng_type" is a random number generator that uses 16 16 ! either (1) Fortran's built-in random_number function, or (2) a … … 25 25 ! mod(istate[i]*A, M). Uniform deviates are returned that normalize 26 26 ! the state elements to the range 0-1. 27 ! 27 28 28 ! The MINSTD generator was coded because the random_numbers_mix 29 29 ! generator in the IFS was found not to vectorize well on some … … 31 31 ! implementation should really be looked at and improved by someone 32 32 ! who knows what they are doing. 33 ! 33 34 34 ! Reference for MINSTD: Park, Stephen K.; Miller, Keith 35 35 ! W. (1988). "Random Number Generators: Good Ones Are Hard To Find" 36 36 ! (PDF). Communications of the ACM. 31 (10): 37 37 ! 1192-1201. doi:10.1145/63039.63042 38 ! 38 39 39 ! Modifications 40 40 ! 2022-12-01 R. Hogan Fixed zeroed state in single precision -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_regions.F90
r5158 r5159 1 1 ! radiation_regions.F90 -- Properties of horizontal regions in Tripleclouds & SPARTACUS 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-07-14 R. Hogan Incorporate gamma distribution option -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_save.F90
r5158 r5159 1 1 ! radiation_save.F90 - Save data to NetCDF files 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-22 R. Hogan Adapt for new way of describing longwave properties -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_single_level.F90
r5158 r5159 1 1 ! radiation_single_level.F90 - Derived type for single-level fields 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2019-01-14 R. Hogan Added out_of_physical_bounds routine -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spartacus_lw.F90
r5158 r5159 1 1 ! radiation_spartacus_lw.F90 - SPARTACUS longwave solver 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive emission/albedo rather than planck/emissivity … … 34 34 ! Sides (SPARTACUS), which can represent 3D effects using a matrix 35 35 ! form of the two-stream equations. 36 ! 36 37 37 ! Sections: 38 38 ! 1: Prepare general variables and arrays -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spartacus_sw.F90
r5158 r5159 1 1 ! radiation_spartacus_sw.F90 - SPARTACUS shortwave solver 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive albedos at g-points … … 46 46 ! Cloud Sides (SPARTACUS), which can represent 3D effects using a 47 47 ! matrix form of the two-stream equations. 48 ! 48 49 49 ! Sections: 50 50 ! 1: Prepare general variables and arrays -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_spectral_definition.F90
r5158 r5159 1 1 ! radiation_spectral_definition.F90 - Derived type to describe a spectral definition 2 ! 2 3 3 ! (C) Copyright 2020- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_thermodynamics.F90
r5158 r5159 1 1 ! radiation_thermodynamics.F90 - Derived type for pressure & temperature 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-05-11 R. Hogan Fix startcol/endcol for get_layer_mass -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_tripleclouds_lw.F90
r5158 r5159 1 1 ! radiation_tripleclouds_lw.F90 - Longwave "Tripleclouds" solver 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-28 R. Hogan Receive emission/albedo rather than planck/emissivity -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_tripleclouds_sw.F90
r5158 r5159 1 1 ! radiation_tripleclouds_sw.F90 - Shortwave "Tripleclouds" solver 2 ! 2 3 3 ! (C) Copyright 2016- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-11 R. Hogan Receive albedos at g-points -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_two_stream.F90
r5158 r5159 1 1 ! radiation_two_stream.F90 - Compute two-stream coefficients 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-05-04 P Dueben/R Hogan Use JPRD where double precision essential -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/utilities/easy_netcdf.F90
r5158 r5159 1 1 ! easy_netcdf.F90 - Module providing convenient NetCDF read/write capability 2 ! 2 3 3 ! (C) Copyright 2014- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! Modifications 16 16 ! 2017-04-28 R. Hogan Fix "reshape" when writing 3D array -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/utilities/radiation_io.F90
r4773 r5159 1 1 ! radiation_io.F90 - Provides logging and abort functionality 2 ! 2 3 3 ! (C) Copyright 2015- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation 10 10 ! nor does it submit to any jurisdiction. 11 ! 11 12 12 ! Author: Robin Hogan 13 13 ! Email: r.j.hogan@ecmwf.int 14 ! 14 15 15 ! This file provides an interface to the provision of file units used 16 16 ! for logging (nulout and nulerr) and for reading data files … … 18 18 ! appropriate for the environment in which the radiation scheme is 19 19 ! embedded. 20 ! 20 21 21 ! Rewrite this file as appropriate if the radiation scheme is to be 22 22 ! embedded into a model other than the ECMWF Integrated Forecasting -
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/utilities/random_numbers_mix.F90
r4773 r5159 2 2 3 3 ! (C) Copyright 2002- ECMWF. 4 ! 4 5 5 ! This software is licensed under the terms of the Apache Licence Version 2.0 6 6 ! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. 7 ! 7 8 8 ! In applying this licence, ECMWF does not waive the privileges and immunities 9 9 ! granted to it by virtue of its status as an intergovernmental organisation … … 85 85 ! Surface Radiative Transport at Large Scale via Monte Carlo. 86 86 ! Annual Review of Heat Transfer, Vol 9. 87 ! 87 88 88 ! Petersen W.P., 1994, Lagged Fibonacci Series Random Number Generator 89 89 ! for the NEC SX-3. International Journal of High Speed Computing 90 90 ! Vol. 6, No. 3, pp387-398. 91 ! 91 92 92 ! Mascagni M., Cuccaro S.A., Pryor D.V., Robinson M.L., 1995, 93 93 ! A Fast, High Quality and Reproducible Parallel Lagged-Fibonacci … … 157 157 ! efficient than the Fibonacci representation. The magic numbers 31 and 87 158 158 ! define the shift register primitive polynomial=(32,7,5,3,2,1,0). 159 ! 159 160 160 ! To ensure that different seeds produce distinct initial buffer states in 161 161 ! canonical form, bits 0...jpmm-2 of the initial seed (after XORing with jpmask … … 314 314 !-------------------------------------------------------------------------------- 315 315 ! Generate normally-distributed random numbers using the Box-Muller method. 316 ! 316 317 317 ! NB: this routine does not use buffering. This means that the following calls: 318 318 ! call gaussian_distribution (zx(1:k))
Note: See TracChangeset
for help on using the changeset viewer.