source: trunk/LMDZ.TITAN/libf/phytitan/surfini.F @ 3709

Last change on this file since 3709 was 1647, checked in by jvatant, 9 years ago

+ Major clean of the new LMDZ.TITAN from too-generic options and routines (water, co2, ocean, surface type ...)
+ From this revision LMDZ.TITAN begins to be really separated from LMDZ.GENERIC
+ Partial desactivation of aerosols, only the dummy case is still enabled to keep the code running ( new aerosol routines to come in followings commits )

JVO
File size: 1.8 KB
Line 
1      SUBROUTINE surfini(ngrid,nq,qsurf,albedo,albedo_bareground)
2
3      USE surfdat_h, only: albedodat
4      use planetwide_mod, only: planetwide_maxval, planetwide_minval
5      use radinc_h, only : L_NSPECTV
6
7      IMPLICIT NONE
8     
9     
10ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
11cccccccccccccc                                                                 cccccccccccccc
12cccccccccccccc   Spectral Albedo Initialisation - Routine modified by MT2015.  cccccccccccccc
13cccccccccccccc                                                                 cccccccccccccc
14ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
15
16
17c--------------------
18c   Declarations:
19c--------------------
20
21      INTEGER,INTENT(IN) :: ngrid
22      INTEGER,INTENT(IN) :: nq
23      REAL,INTENT(OUT) :: albedo(ngrid,L_NSPECTV)
24      REAL,INTENT(OUT) :: albedo_bareground(ngrid)
25      REAL,INTENT(IN) :: qsurf(ngrid,nq) ! tracer on surface (kg/m2)
26
27      INTEGER :: ig,nw
28      REAL :: min_albedo,max_albedo
29
30c=======================================================================
31
32      ! We get the bare ground albedo from the start files.
33      DO ig=1,ngrid
34         albedo_bareground(ig)=albedodat(ig)
35         DO nw=1,L_NSPECTV
36            albedo(ig,nw)=albedo_bareground(ig)
37         ENDDO
38      ENDDO
39      call planetwide_minval(albedo_bareground,min_albedo)
40      call planetwide_maxval(albedo_bareground,max_albedo)
41      write(*,*) 'surfini: minimum bare ground albedo',min_albedo
42      write(*,*) 'surfini: maximum bare ground albedo',max_albedo
43
44   
45      call planetwide_minval(albedo,min_albedo)
46      call planetwide_maxval(albedo,max_albedo)
47      write(*,*) 'surfini: minimum corrected initial albedo',min_albedo
48      write(*,*) 'surfini: maximum corrected initial albedo',max_albedo
49
50
51      END
Note: See TracBrowser for help on using the repository browser.