1 |
molod |
1.2 |
C $Header: /u/gcmpack/MITgcm/pkg/fizhi/update_earth_exports.F,v 1.1 2004/06/04 16:23:37 molod Exp $ |
2 |
molod |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
subroutine update_earth_exports (myTime, myIter, myThid) |
5 |
|
|
c---------------------------------------------------------------------- |
6 |
|
|
c Subroutine update_earth_exports - 'Wrapper' routine to update |
7 |
|
|
c the fields related to the earth's surface that are needed |
8 |
|
|
c by fizhi. |
9 |
|
|
c Also: Set up "bi, bj loop" and some timers and clocks here. |
10 |
|
|
c |
11 |
|
|
c Call: getalb (Set the 4 albedos based on veg type and time) |
12 |
|
|
c getemiss (Set the surface emissivity based on the veg type |
13 |
|
|
c and the snow depth) |
14 |
|
|
c getlgr (Set the leaf area index and surface greenness, |
15 |
|
|
c based on veg type and month) |
16 |
|
|
c----------------------------------------------------------------------- |
17 |
|
|
implicit none |
18 |
|
|
#include "CPP_OPTIONS.h" |
19 |
|
|
#include "SIZE.h" |
20 |
|
|
#include "GRID.h" |
21 |
molod |
1.2 |
#include "fizhi_land_SIZE.h" |
22 |
molod |
1.1 |
#include "fizhi_SIZE.h" |
23 |
|
|
#include "DYNVARS.h" |
24 |
|
|
#include "fizhi_coms.h" |
25 |
|
|
#include "gridalt_mapping.h" |
26 |
molod |
1.2 |
#include "fizhi_land_coms.h" |
27 |
molod |
1.1 |
#include "EEPARAMS.h" |
28 |
|
|
|
29 |
|
|
integer myTime, myIter, myThid |
30 |
|
|
|
31 |
|
|
integer i, j, L, bi, bj |
32 |
|
|
integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2 |
33 |
|
|
|
34 |
|
|
im1 = 1-OLx |
35 |
|
|
im2 = sNx+OLx |
36 |
|
|
jm1 = 1-OLy |
37 |
|
|
jm2 = sNy+OLy |
38 |
|
|
idim1 = 1 |
39 |
|
|
idim2 = sNx |
40 |
|
|
jdim1 = 1 |
41 |
|
|
jdim2 = sNy |
42 |
|
|
|
43 |
|
|
do bj = myByLo(myThid), myByHi(myThid) |
44 |
|
|
do bi = myBxLo(myThid), myBxHi(myThid) |
45 |
|
|
|
46 |
|
|
enddo |
47 |
|
|
enddo |
48 |
|
|
|
49 |
|
|
return |
50 |
|
|
end |