Parent Directory
| 
Revision Log
Just encapsulated SUBROUTINE vlsplt in a module and cleaned it. In procedure vlx, local variables dxqu and adxqu only need indices iip2:ip1jm. Otherwise, just cleaned vlx. Procedures dynredem0 and dynredem1 no longer have argument fichnom, they just operate on a file named "restart.nc". The programming guideline here is that gcm should not be more complex than it needs by itself, other programs (ce0l etc.) just have to adapt to gcm. So ce0l now creates files "restart.nc" and "restartphy.nc". In order to facilitate decentralizing the writing of "restartphy.nc", created a procedure phyredem0 out of phyredem. phyredem0 creates the NetCDF header of "restartphy.nc" while phyredem writes the NetCDF variables. As the global attribute itau_phy needs to be filled in phyredem0, at the beginnig of the run, we must compute its value instead of just using itap. So we have a dummy argument lmt_pas of phyredem0. Also, the ncid of "startphy.nc" is upgraded from local variable of phyetat0 to dummy argument. phyetat0 no longer closes "startphy.nc". Following the same decentralizing objective, the ncid of "restart.nc" is upgraded from local variable of dynredem0 to module variable of dynredem0_m. "restart.nc" is not closed at the end of dynredem0 nor opened at the beginning of dynredem1. In procedure etat0, instead of creating many vectors of size klon which will be filled with zeroes, just create one array null_array. In procedure phytrac, instead of writing trs(: 1) to a text file, write it to "restartphy.nc" (following LMDZ). This is better because now trs(: 1) is next to its coordinates. We can write to "restartphy.nc" from phytrac directly, and not add trs(: 1) to the long list of variables in physiq, thanks to the decentralizing of "restartphy.nc". In procedure phyetat0, we no longer write to standard output the minimum and maximum values of read arrays. It is ok to check input and abort on invalid values but just printing statistics on input seems too much useless computation and out of place clutter.
| 1 | guez | 3 | |
| 2 | guez | 81 | ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/limz.F,v 1.1.1.1 2004/05/19 |
| 3 | ! 12:53:07 lmdzadmin Exp $ | ||
| 4 | guez | 3 | |
| 5 | guez | 81 | SUBROUTINE limz(s0, sz, sm, pente_max) |
| 6 | guez | 3 | |
| 7 | guez | 81 | ! Auteurs: P.Le Van, F.Hourdin, F.Forget |
| 8 | guez | 3 | |
| 9 | guez | 81 | ! ******************************************************************** |
| 10 | ! Shema d'advection " pseudo amont " . | ||
| 11 | ! ******************************************************************** | ||
| 12 | ! nq,iq,q,pbaru,pbarv,w sont des arguments d'entree pour le s-pg .... | ||
| 13 | guez | 3 | |
| 14 | |||
| 15 | guez | 81 | ! -------------------------------------------------------------------- |
| 16 | USE dimens_m | ||
| 17 | USE paramet_m | ||
| 18 | USE comconst | ||
| 19 | USE disvert_m | ||
| 20 | USE conf_gcm_m | ||
| 21 | USE comgeom | ||
| 22 | IMPLICIT NONE | ||
| 23 | guez | 3 | |
| 24 | |||
| 25 | |||
| 26 | guez | 81 | ! Arguments: |
| 27 | ! ---------- | ||
| 28 | REAL pente_max | ||
| 29 | REAL s0(ip1jmp1, llm), sm(ip1jmp1, llm) | ||
| 30 | REAL sz(ip1jmp1, llm) | ||
| 31 | guez | 3 | |
| 32 | guez | 81 | ! Local |
| 33 | ! --------- | ||
| 34 | guez | 3 | |
| 35 | guez | 155 | INTEGER ij, l |
| 36 | guez | 3 | |
| 37 | guez | 81 | REAL q(ip1jmp1, llm) |
| 38 | REAL dzq(ip1jmp1, llm) | ||
| 39 | guez | 3 | |
| 40 | guez | 81 | REAL dzqw(ip1jmp1) |
| 41 | REAL adzqw(ip1jmp1), dzqmax(ip1jmp1) | ||
| 42 | guez | 3 | |
| 43 | guez | 155 | REAL ssum |
| 44 | guez | 81 | INTEGER ismax, ismin |
| 45 | EXTERNAL ssum, convflu, ismin, ismax | ||
| 46 | |||
| 47 | DO l = 1, llm | ||
| 48 | DO ij = 1, ip1jmp1 | ||
| 49 | q(ij, l) = s0(ij, l)/sm(ij, l) | ||
| 50 | dzq(ij, l) = sz(ij, l)/sm(ij, l) | ||
| 51 | END DO | ||
| 52 | END DO | ||
| 53 | |||
| 54 | ! calcul de la pente en haut et en bas de la maille | ||
| 55 | DO ij = 1, ip1jmp1 | ||
| 56 | DO l = 1, llm - 1 | ||
| 57 | dzqw(l) = q(ij, l+1) - q(ij, l) | ||
| 58 | END DO | ||
| 59 | dzqw(llm) = 0. | ||
| 60 | |||
| 61 | DO l = 1, llm | ||
| 62 | adzqw(l) = abs(dzqw(l)) | ||
| 63 | END DO | ||
| 64 | |||
| 65 | ! calcul de la pente maximum dans la maille en valeur absolue | ||
| 66 | |||
| 67 | DO l = 2, llm - 1 | ||
| 68 | dzqmax(l) = pente_max*min(adzqw(l-1), adzqw(l)) | ||
| 69 | END DO | ||
| 70 | |||
| 71 | ! calcul de la pente avec limitation | ||
| 72 | |||
| 73 | DO l = 2, llm - 1 | ||
| 74 | IF (dzqw(l-1)*dzqw(l)>0. .AND. dzq(ij,l)*dzqw(l)>0.) THEN | ||
| 75 | dzq(ij, l) = sign(min(abs(dzq(ij,l)),dzqmax(l)), dzq(ij,l)) | ||
| 76 | ELSE | ||
| 77 | ! extremum local | ||
| 78 | dzq(ij, l) = 0. | ||
| 79 | END IF | ||
| 80 | END DO | ||
| 81 | |||
| 82 | DO l = 1, llm | ||
| 83 | sz(ij, l) = dzq(ij, l)*sm(ij, l) | ||
| 84 | END DO | ||
| 85 | |||
| 86 | END DO | ||
| 87 | |||
| 88 | RETURN | ||
| 89 | END SUBROUTINE limz |
| ViewVC Help | |
| Powered by ViewVC 1.1.21 |