source: dynamico_lmdz/aquaplanet/LMDZ5/libf/phymar/PHY_SISVAT_OUT.f90 @ 4031

Last change on this file since 4031 was 3990, checked in by millour, 8 years ago

An updated LMDZ5 (based on LMDZ rev 2786) to play with.
EM
File size: 2.9 KB
Line 
1      subroutine PHY_SISVAT_OUT
2
3!------------------------------------------------------------------------------+
4!                                                         Wed 26-Jun-2013  MAR |
5!     SubRoutine PHY_SISVAT_OUT writes main variables from                     |
6!                    SISVAT    (Soil                                           |
7!                               Ice                                            |
8!                               Snow                                           |
9!                               Vegetation                                     |
10!                               Atmosphere                                     |
11!                               Transfer   Scheme)                             |
12!                                                                              |
13!                                                                              |
14!     version 3.p.4.1 created by H. Gallee,               Mon 24-Jun-2013      |
15!           Last Modification by H. Gallee,               Wed 26-Jun-2013      |
16!                                                                              |
17!------------------------------------------------------------------------------+
18
19
20! Global  Variables
21! =================
22
23      use Mod_Real
24      use Mod_PHY____grd
25      use Mod_SISVAT_grd
26      use Mod_SISVAT_dat
27      use Mod_SISVAT_kkl
28      use Mod_SISVAT_flx
29      use Mod_SISVAT_cdf
30
31
32
33
34! Local   Variables
35! =================
36
37! Dummy
38! -----
39
40      integer  ::  i     ,j     ,nm    ,ikp
41      integer  ::  isl
42
43
44
45
46! From 1D to 2D Horizontal Grid
47! =============================
48
49        DO ikp=1,kcolp
50           i   = ii__AP(ikp)
51           j   = jj__AP(ikp)
52        DO nm=1,mwp
53
54! Energy Fluxes                                             (OUTPUT/NetCDF)
55! ^^^^^^^^^^^^^                                              ^^^^^^^^^^^^^
56           SOsoNC_xyn(i,j,nm)      = SOsoKL(ikp,nm)                     ! Absorb.Sol.Rad.            [W/m2]
57           IRsoNC_xyn(i,j,nm)      = IRsoKL(ikp,nm)                     ! Absorb.IR  Rad.            [W/m2]
58           HSsoNC_xyn(i,j,nm)      = HSsoKL(ikp,nm)                     ! HS                         [W/m2]
59           HLsoNC_xyn(i,j,nm)      = HLsoKL(ikp,nm)                     ! HL                         [W/m2]
60           HLs_NC_xyn(i,j,nm)      = HLs_KL(ikp,nm)                     ! Evaporation           [mm w.e./s]
61           HLv_NC_xyn(i,j,nm)      = HLv_KL(ikp,nm)                     ! Transpiration              [W/m2]
62
63           eta_NC_xyn(i,j,nm)      = 0.                                 !
64         DO isl = -nsoil,0                                              !
65           eta_NC_xyn(i,j,nm)      = eta_NC_xyn(i,j,nm)                &! Soil Moisture
66     &    +eta_SV(ikp,nm,isl)      * dz_dSV(    isl)                    !
67         END DO                                                         !
68
69        END DO
70        END DO
71
72
73
74
75      end subroutine PHY_SISVAT_OUT
Note: See TracBrowser for help on using the repository browser.