Last change
on this file since 4013 was
3908,
checked in by idelkadi, 4 years ago
|
Online implementation of the radiative transfer code ECRAD in the LMDZ model.
- Inclusion of the ecrad directory containing the sources of the ECRAD code
- interface routine : radiation_scheme.F90
- Adaptation of compilation scripts :
- compilation under CPP key CPP_ECRAD
- compilation with option "-rad ecard" or "-ecard true"
- The "-rad old/rtm/ecran" build option will need to replace the "-rrtm true" and "-ecrad true" options in the future.
- Runing LMDZ simulations with ecrad, you need :
- logical key iflag_rrtm = 2 in physiq.def
- namelist_ecrad (DefLists?)
- the directory "data" containing the configuration files is temporarily placed in ../libfphylmd/ecrad/
- Compilation and execution are tested in the 1D case. The repository under svn would allow to continue the implementation work: tests, verification of the results, ...
|
File size:
1.3 KB
|
Rev | Line | |
---|
[3908] | 1 | /** |
---|
| 2 | * (C) Copyright 2014- ECMWF. |
---|
| 3 | * |
---|
| 4 | * This software is licensed under the terms of the Apache Licence Version 2.0 |
---|
| 5 | * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. |
---|
| 6 | * |
---|
| 7 | * In applying this licence, ECMWF does not waive the privileges and immunities |
---|
| 8 | * granted to it by virtue of its status as an intergovernmental organisation |
---|
| 9 | * nor does it submit to any jurisdiction. |
---|
| 10 | */ |
---|
| 11 | |
---|
| 12 | // cas.h |
---|
| 13 | // |
---|
| 14 | // compare_and_swap -based locks |
---|
| 15 | // |
---|
| 16 | // Thanks to https://github.com/majek/dump/blob/master/msqueue/queue_lock_myspinlock1.c |
---|
| 17 | // |
---|
| 18 | |
---|
| 19 | #include <signal.h> |
---|
| 20 | #include <sched.h> |
---|
| 21 | |
---|
| 22 | #ifndef INLINE |
---|
| 23 | #define INLINE __inline__ |
---|
| 24 | #endif |
---|
| 25 | |
---|
| 26 | #if defined(__GNUC__) && !defined(__NEC__) |
---|
| 27 | |
---|
| 28 | #define CAS(lock,oldval,newval) __sync_bool_compare_and_swap(lock,oldval,newval) |
---|
| 29 | |
---|
| 30 | #else |
---|
| 31 | |
---|
| 32 | #warning *** CAS-locks self-implemented *** |
---|
| 33 | |
---|
| 34 | static INLINE int CAS(volatile sig_atomic_t *lock, int oldval, int newval) |
---|
| 35 | { |
---|
| 36 | int tmp = *lock; |
---|
| 37 | if (tmp == oldval) *lock = newval; |
---|
| 38 | return tmp; |
---|
| 39 | } |
---|
| 40 | |
---|
| 41 | #endif |
---|
| 42 | |
---|
| 43 | static INLINE void cas_init(volatile sig_atomic_t *lock) |
---|
| 44 | { |
---|
| 45 | if (lock) *lock = 0; |
---|
| 46 | } |
---|
| 47 | |
---|
| 48 | static INLINE void cas_lock(volatile sig_atomic_t *lock) |
---|
| 49 | { |
---|
| 50 | while (1) { |
---|
| 51 | int i; |
---|
| 52 | for (i=0; i < 10000; ++i) { |
---|
| 53 | if (CAS(lock, 0, 1)) return; |
---|
| 54 | } |
---|
| 55 | sched_yield(); |
---|
| 56 | } |
---|
| 57 | } |
---|
| 58 | |
---|
| 59 | static INLINE void cas_unlock(volatile sig_atomic_t *lock) |
---|
| 60 | { |
---|
| 61 | CAS(lock, 1, 0); |
---|
| 62 | } |
---|
Note: See
TracBrowser
for help on using the repository browser.