source: trunk/LMDZ.GENERIC/libf/phystd/surface_nature.F @ 2537

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.