|
Last change
on this file since 1477 was
1403,
checked in by emillour, 10 years ago
|
|
All models: Reorganizing the physics/dynamics interface.
- makelmdz and makelmdz_fcm scripts adapted to handle the new directory
settings
- misc: (replaces what was the "bibio" directory)
- Should only contain extremely generic (and non physics or
dynamics-specific) routines
- Therefore moved initdynav.F90, initfluxsto.F, inithist.F, writedynav.F90,
write_field.F90, writehist.F to "dyn3d_common"
- dynlonlat_phylonlat: (new interface directory)
- This directory contains routines relevent to physics/dynamics grid
interactions, e.g. routines gr_dyn_fi or gr_fi_dyn and calfis
- Moreover the dynlonlat_phylonlat contains directories "phy*" corresponding
to each physics package "phy*" to be used. These subdirectories should
only contain specific interfaces (e.g. iniphysiq) or main programs (e.g.
newstart)
- phy*/dyn1d: this subdirectory contains the 1D model using physics
from phy*
EM
|
|
File size:
1.2 KB
|
| Line | |
|---|
| 1 | SUBROUTINE wind_scal(pbaru,pbarv,us,vs) |
|---|
| 2 | c======================================================================= |
|---|
| 3 | c |
|---|
| 4 | c |
|---|
| 5 | c Subject: |
|---|
| 6 | c ------ |
|---|
| 7 | c On ramene les flux de masse /vents aux points scalaires. |
|---|
| 8 | c |
|---|
| 9 | c======================================================================= |
|---|
| 10 | IMPLICIT NONE |
|---|
| 11 | c----------------------------------------------------------------------- |
|---|
| 12 | c Declararations: |
|---|
| 13 | c --------------- |
|---|
| 14 | |
|---|
| 15 | #include "dimensions.h" |
|---|
| 16 | #include "paramet.h" |
|---|
| 17 | #include "comgeom.h" |
|---|
| 18 | |
|---|
| 19 | c Arguments: |
|---|
| 20 | c ---------- |
|---|
| 21 | |
|---|
| 22 | REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm) |
|---|
| 23 | REAL us(ip1jmp1,llm), vs (ip1jmp1,llm) |
|---|
| 24 | |
|---|
| 25 | c Local: |
|---|
| 26 | c ------ |
|---|
| 27 | |
|---|
| 28 | INTEGER ij,l |
|---|
| 29 | |
|---|
| 30 | c----------------------------------------------------------------------- |
|---|
| 31 | |
|---|
| 32 | c transport zonal: |
|---|
| 33 | c ---------------- |
|---|
| 34 | DO l=1,llm |
|---|
| 35 | DO ij=2,ip1jmp1 |
|---|
| 36 | us(ij,l)=.5*(pbaru(ij,l)+pbaru(ij-1,l)) |
|---|
| 37 | ENDDO |
|---|
| 38 | ENDDO |
|---|
| 39 | CALL SCOPY(jjp1*llm,us(iip1,1),iip1,us(1,1),iip1) |
|---|
| 40 | |
|---|
| 41 | |
|---|
| 42 | c Transport meridien: |
|---|
| 43 | c ------------------- |
|---|
| 44 | DO l=1,llm |
|---|
| 45 | DO ij=iip2,ip1jm |
|---|
| 46 | vs(ij,l)=.5*(pbarv(ij,l)+pbarv(ij-iip1,l)) |
|---|
| 47 | ENDDO |
|---|
| 48 | DO ij=1,iip1 |
|---|
| 49 | vs(ij,l)=0. |
|---|
| 50 | vs(ip1jm+ij,l)=0. |
|---|
| 51 | ENDDO |
|---|
| 52 | ENDDO |
|---|
| 53 | |
|---|
| 54 | RETURN |
|---|
| 55 | END |
|---|
Note: See
TracBrowser
for help on using the repository browser.
