| [1647] | 1 | SUBROUTINE surfini(ngrid,nq,qsurf,albedo,albedo_bareground) |
|---|
| [787] | 2 | |
|---|
| [1482] | 3 | USE surfdat_h, only: albedodat |
|---|
| [1216] | 4 | use planetwide_mod, only: planetwide_maxval, planetwide_minval |
|---|
| [1482] | 5 | use radinc_h, only : L_NSPECTV |
|---|
| [787] | 6 | |
|---|
| [135] | 7 | IMPLICIT NONE |
|---|
| [1482] | 8 | |
|---|
| 9 | |
|---|
| 10 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 11 | cccccccccccccc cccccccccccccc |
|---|
| 12 | cccccccccccccc Spectral Albedo Initialisation - Routine modified by MT2015. cccccccccccccc |
|---|
| 13 | cccccccccccccc cccccccccccccc |
|---|
| 14 | ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc |
|---|
| 15 | |
|---|
| 16 | |
|---|
| 17 | c-------------------- |
|---|
| [135] | 18 | c Declarations: |
|---|
| [1482] | 19 | c-------------------- |
|---|
| 20 | |
|---|
| [1216] | 21 | INTEGER,INTENT(IN) :: ngrid |
|---|
| 22 | INTEGER,INTENT(IN) :: nq |
|---|
| [1482] | 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) |
|---|
| [135] | 26 | |
|---|
| [1482] | 27 | INTEGER :: ig,nw |
|---|
| [1216] | 28 | REAL :: min_albedo,max_albedo |
|---|
| [1482] | 29 | |
|---|
| [135] | 30 | c======================================================================= |
|---|
| 31 | |
|---|
| [1647] | 32 | ! We get the bare ground albedo from the start files. |
|---|
| [135] | 33 | DO ig=1,ngrid |
|---|
| [1482] | 34 | albedo_bareground(ig)=albedodat(ig) |
|---|
| 35 | DO nw=1,L_NSPECTV |
|---|
| 36 | albedo(ig,nw)=albedo_bareground(ig) |
|---|
| 37 | ENDDO |
|---|
| [135] | 38 | ENDDO |
|---|
| [1482] | 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 |
|---|
| [135] | 43 | |
|---|
| [1647] | 44 | |
|---|
| [1482] | 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 |
|---|
| [135] | 49 | |
|---|
| 50 | |
|---|
| 51 | END |
|---|
