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

Last change on this file since 1477 was 1397, checked in by milmd, 10 years ago

In LMDZ.GENERIC replacement of all phystd .h files by module files.
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 comgeomfi_h
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+area(ig)
59           end if
60        enddo
61
62        return
63        end
64
Note: See TracBrowser for help on using the repository browser.