Last change
on this file since 2537 was
1543,
checked in by emillour, 9 years ago
|
All models: Further adaptations to keep up with changes in LMDZ5 concerning
physics/dynamics separation:
- dyn3d:
- adapted gcm.F so that all physics initializations are now done in iniphysiq.
- dyn3dpar:
- adapted gcm.F so that all physics initializations are now done in iniphysiq.
- updated calfis_p.F to follow up with changes.
- copied over updated "bands.F90" from LMDZ5.
- dynphy_lonlat:
- calfis_p.F90, mod_interface_dyn_phys.F90,
follow up of changes in phy_common/mod_* routines
- phy_common:
- added "geometry_mod.F90" to store information about the grid (replaces
phy*/comgeomphy.F90) and give variables friendlier names: rlond =>
longitude , rlatd => latitude, airephy => cell_area,
cuphy => dx , cvphy => dy
- added "physics_distribution_mod.F90"
- updated "mod_grid_phy_lmdz.F90", "mod_phys_lmdz_mpi_data.F90",
"mod_phys_lmdz_para.F90", "mod_phys_lmdz_mpi_transfert.F90",
"mod_grid_phy_lmdz.F90", "mod_phys_lmdz_omp_data.F90",
"mod_phys_lmdz_omp_transfert.F90", "write_field_phy.F90"
and "ioipsl_getin_p_mod.F90" to LMDZ5 versions.
- phy[venus/titan/mars/std]:
- removed "init_phys_lmdz.F90", "comgeomphy.F90"; adapted routines to use
geometry_mod (longitude, latitude, cell_area, etc.)
EM
|
File size:
1.2 KB
|
Line | |
---|
1 | SUBROUTINE surface_nature(ngrid,nq,obliquit,qsurf,qsurfliquid |
---|
2 | & ,qsurfsnow,rnat,oceanarea) |
---|
3 | |
---|
4 | USE surfdat_h |
---|
5 | USE comsoil_h |
---|
6 | USE geometry_mod, ONLY: cell_area |
---|
7 | USE tracer_h |
---|
8 | |
---|
9 | IMPLICIT none |
---|
10 | |
---|
11 | !================================================================== |
---|
12 | ! |
---|
13 | ! Purpose |
---|
14 | ! ------- |
---|
15 | ! Defines a few things |
---|
16 | ! |
---|
17 | ! Authors |
---|
18 | ! ------- |
---|
19 | ! B. Charnay (2010) |
---|
20 | ! |
---|
21 | ! Called by |
---|
22 | ! --------- |
---|
23 | ! physiq.F |
---|
24 | ! |
---|
25 | ! Calls |
---|
26 | ! ----- |
---|
27 | ! none |
---|
28 | ! |
---|
29 | ! Notes |
---|
30 | ! ----- |
---|
31 | ! rnat is terrain type: 0-ocean; 1-continent; 2-continental ice |
---|
32 | ! |
---|
33 | !================================================================== |
---|
34 | |
---|
35 | integer ngrid,nq |
---|
36 | |
---|
37 | REAL qsurf(ngrid,nq),ps(ngrid) |
---|
38 | REAL qsurfliquid(ngrid) |
---|
39 | REAL qsurfsnow(ngrid) |
---|
40 | INTEGER iq, ig |
---|
41 | INTEGER rnat(ngrid) |
---|
42 | REAL oceanarea |
---|
43 | REAL obliquit |
---|
44 | |
---|
45 | do ig=1,ngrid |
---|
46 | rnat(ig)=1 |
---|
47 | dryness(ig)=1 !(coefficient for evaporation) |
---|
48 | if (inertiedat(ig,1).gt.1E4) then |
---|
49 | rnat(ig)=0 |
---|
50 | end if |
---|
51 | end do |
---|
52 | |
---|
53 | ! surface of all the oceans |
---|
54 | |
---|
55 | oceanarea=0. |
---|
56 | do ig=1,ngrid |
---|
57 | if (rnat(ig).eq.0)then |
---|
58 | oceanarea=oceanarea+cell_area(ig) |
---|
59 | end if |
---|
60 | enddo |
---|
61 | |
---|
62 | return |
---|
63 | end |
---|
64 | |
---|
Note: See
TracBrowser
for help on using the repository browser.