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

Last change on this file since 253 was 253, checked in by emillour, 13 years ago

Generic GCM

  • Massive update to version 0.7

EM+RW
File size: 1.3 KB
RevLine 
[253]1      SUBROUTINE surface_nature(obliquit,qsurf,qsurfliquid
2     &   ,qsurfsnow,rnat,oceanarea)
3      IMPLICIT none
4
5!==================================================================
6!     
7!     Purpose
8!     -------
9!     Defines a few things
10!     
11!     Authors
12!     -------
13!     B. Charnay (2010)
14!     
15!     Called by
16!     ---------
17!     physiq.F
18!     
19!     Calls
20!     -----
21!     none
22!
23!     Notes
24!     -----
25!     rnat is terrain type: 0-ocean; 1-continent; 2-continental ice
26!     
27!==================================================================
28
29#include "dimensions.h"
30#include "dimphys.h"
31#include "comcstfi.h"
32#include "callkeys.h"
33#include "tracer.h"
34#include "fisice.h"
35#include "comgeomfi.h"
36#include "surfdat.h"
37#include "comsoil.h"
38
39        REAL qsurf(ngridmx,nqmx),ps(ngridmx)
40        REAL qsurfliquid(ngridmx)
41        REAL qsurfsnow(ngridmx)
42        INTEGER iq, ig
43        INTEGER rnat(ngridmx)
44        REAL oceanarea
45        REAL obliquit
46
47        do ig=1,ngridmx
48           rnat(ig)=1
49           dryness(ig)=1        !(coefficient for evaporation)
50           if (inertiedat(ig,1).gt.1E4) then
51              rnat(ig)=0
52           end if
53        end do
54
55! surface of all the oceans
56       
57        oceanarea=0.
58        do ig=1,ngridmx
59           if (rnat(ig).eq.0)then
60              oceanarea=oceanarea+area(ig)
61           end if
62        enddo
63
64        return
65        end
66
Note: See TracBrowser for help on using the repository browser.