source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/surface_nature.F @ 253

Last change on this file since 253 was 227, checked in by milmd, 10 years ago

Last LMDZ version (1315) with OpenMP directives and other stuff
File size: 1.3 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!#include "dimensions.h"
36!#include "dimphys.h"
37#include "comcstfi.h"
38#include "callkeys.h"
39
40        integer ngrid,nq
41
42        REAL qsurf(ngrid,nq),ps(ngrid)
43        REAL qsurfliquid(ngrid)
44        REAL qsurfsnow(ngrid)
45        INTEGER iq, ig
46        INTEGER rnat(ngrid)
47        REAL oceanarea
48        REAL obliquit
49
50        do ig=1,ngrid
51           rnat(ig)=1
52           dryness(ig)=1        !(coefficient for evaporation)
53           if (inertiedat(ig,1).gt.1E4) then
54              rnat(ig)=0
55           end if
56        end do
57
58! surface of all the oceans
59       
60        oceanarea=0.
61        do ig=1,ngrid
62           if (rnat(ig).eq.0)then
63              oceanarea=oceanarea+area(ig)
64           end if
65        enddo
66
67        return
68        end
69
Note: See TracBrowser for help on using the repository browser.