Parent Directory
| 
Revision Log
Removed procedure sortvarc0. Called sortvarc with an additional argument resetvarc instead. (Following LMDZ.) Moved current time computations and some printing statements from sortvarc to caldyn. Could then remove arguments itau and time_0 of sortvarc, and could remove "use dynetat0". Better to keep "dynetat0.f" as a gcm-only file. Moved some variables from module ener to module sortvarc. Split file "mathelp.f" into single-procedure files. Removed unused argument nadv of adaptdt. Removed dimension arguments of bernoui. Removed unused argument nisurf of interfoce_lim. Changed the size of argument lmt_sst of interfoce_lim from klon to knon. Removed case when newlmt is false. dynredem1 is called only once in each run, either ce0l or gcm. So variable nb in call to nf95_put_var was always 1. Removed variable nb. Removed dimension arguments of calcul_fluxs. Removed unused arguments precip_rain, precip_snow, snow of calcul_fluxs. Changed the size of all the arrays in calcul_fluxs from klon to knon. Removed dimension arguments of fonte_neige. Changed the size of all the arrays in fonte_neige from klon to knon. Changed the size of arguments tsurf and tsurf_new of interfsurf_hq from klon to knon. Changed the size of argument ptsrf of soil from klon to knon.
| 1 | module dteta1_m |
| 2 | |
| 3 | IMPLICIT NONE |
| 4 | |
| 5 | contains |
| 6 | |
| 7 | SUBROUTINE dteta1(teta, pbaru, pbarv, dteta) |
| 8 | |
| 9 | ! From LMDZ4/libf/dyn3d/dteta1.F, version 1.1.1.1, 2004/05/19 12:53:06 |
| 10 | ! Authors: P. Le Van, F. Forget |
| 11 | |
| 12 | ! Calcul du terme de convergence horizontale du flux d'enthalpie |
| 13 | ! potentielle. |
| 14 | |
| 15 | USE dimens_m, ONLY: iim, jjm, llm |
| 16 | USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1, jjp1 |
| 17 | USE filtreg_m, ONLY: filtreg |
| 18 | |
| 19 | REAL, intent(in):: teta(ip1jmp1, llm) |
| 20 | REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm) |
| 21 | REAL, intent(out):: dteta(iim + 1, jjm + 1, llm) |
| 22 | |
| 23 | ! Local: |
| 24 | INTEGER l, ij |
| 25 | REAL hbyv(ip1jm, llm), hbxu(ip1jmp1, llm) |
| 26 | |
| 27 | !---------------------------------------------------------------- |
| 28 | |
| 29 | DO l = 1, llm |
| 30 | DO ij = iip2, ip1jm - 1 |
| 31 | hbxu(ij, l) = pbaru(ij, l) * 0.5 * (teta(ij, l) + teta(ij + 1, l)) |
| 32 | end DO |
| 33 | |
| 34 | DO ij = iip1+ iip1, ip1jm, iip1 |
| 35 | hbxu(ij, l) = hbxu(ij - iim, l) |
| 36 | end DO |
| 37 | |
| 38 | DO ij = 1, ip1jm |
| 39 | hbyv(ij, l)= pbarv(ij, l) * 0.5 * (teta(ij, l) + teta(ij + iip1, l)) |
| 40 | end DO |
| 41 | end DO |
| 42 | |
| 43 | CALL convflu(hbxu, hbyv, llm, dteta) |
| 44 | |
| 45 | ! stockage dans dh de la convergence horizontale filtrée du flux |
| 46 | ! d'enthalpie potentielle |
| 47 | CALL filtreg(dteta, jjp1, llm, 2, 2, .true.) |
| 48 | |
| 49 | END SUBROUTINE dteta1 |
| 50 | |
| 51 | end module dteta1_m |
| ViewVC Help | |
| Powered by ViewVC 1.1.21 |