Last change
on this file since 2916 was
2311,
checked in by Ehouarn Millour, 9 years ago
|
Further modifications to enforce physics/dynamics separation:
- moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
- created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to physics (should be used rather than iniprint.h)
- created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead when in physics.
- reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and print_control_mod flags.
EM
|
-
Property copyright set to
Name of program: LMDZ Creation date: 1984 Version: LMDZ5 License: CeCILL version 2 Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539 See the license file in the root directory
-
Property svn:eol-style set to
native
-
Property svn:keywords set to
Author Date Id Revision
|
File size:
1.6 KB
|
Rev | Line | |
---|
[1403] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
[972] | 4 | SUBROUTINE thermcell_init(ngrid,nlay,ztv,zlay,zlev, & |
---|
[878] | 5 | & lalim,lmin,alim_star,alim_star_tot,lev_out) |
---|
| 6 | |
---|
| 7 | !---------------------------------------------------------------------- |
---|
| 8 | !thermcell_init: calcul du profil d alimentation du thermique |
---|
| 9 | !---------------------------------------------------------------------- |
---|
[2311] | 10 | USE print_control_mod, ONLY: lunout |
---|
[878] | 11 | IMPLICIT NONE |
---|
[1026] | 12 | #include "thermcell.h" |
---|
[938] | 13 | |
---|
[878] | 14 | INTEGER l,ig |
---|
| 15 | !arguments d entree |
---|
| 16 | INTEGER ngrid,nlay |
---|
| 17 | REAL ztv(ngrid,nlay) |
---|
[972] | 18 | REAL zlay(ngrid,nlay) |
---|
| 19 | REAL zlev(ngrid,nlay+1) |
---|
[878] | 20 | !arguments de sortie |
---|
| 21 | INTEGER lalim(ngrid) |
---|
| 22 | INTEGER lmin(ngrid) |
---|
| 23 | REAL alim_star(ngrid,nlay) |
---|
| 24 | REAL alim_star_tot(ngrid) |
---|
| 25 | integer lev_out ! niveau pour les print |
---|
| 26 | |
---|
[972] | 27 | REAL zzalim(ngrid) |
---|
[878] | 28 | !CR: ponderation entrainement des couches instables |
---|
| 29 | !def des alim_star tels que alim=f*alim_star |
---|
| 30 | |
---|
| 31 | |
---|
[1403] | 32 | write(lunout,*)'THERM INIT V20C ' |
---|
[1026] | 33 | |
---|
[1403] | 34 | alim_star_tot(:)=0. |
---|
| 35 | alim_star(:,:)=0. |
---|
| 36 | lmin(:)=1 |
---|
| 37 | lalim(:)=1 |
---|
[878] | 38 | |
---|
[1403] | 39 | do l=1,nlay-1 |
---|
[878] | 40 | do ig=1,ngrid |
---|
[1403] | 41 | if (ztv(ig,l)> ztv(ig,l+1) .and. ztv(ig,1)>=ztv(ig,l) ) then |
---|
| 42 | alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.) & |
---|
[878] | 43 | & *sqrt(zlev(ig,l+1)) |
---|
[1403] | 44 | lalim(:)=l+1 |
---|
| 45 | alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l) |
---|
[878] | 46 | endif |
---|
| 47 | enddo |
---|
| 48 | enddo |
---|
| 49 | do l=1,nlay |
---|
| 50 | do ig=1,ngrid |
---|
[1403] | 51 | if (alim_star_tot(ig) > 1.e-10 ) then |
---|
[878] | 52 | alim_star(ig,l)=alim_star(ig,l)/alim_star_tot(ig) |
---|
| 53 | endif |
---|
| 54 | enddo |
---|
| 55 | enddo |
---|
[1403] | 56 | alim_star_tot(:)=1. |
---|
[1026] | 57 | |
---|
[1403] | 58 | return |
---|
[878] | 59 | end |
---|
Note: See
TracBrowser
for help on using the repository browser.