| Directory: | ./ |
|---|---|
| File: | phys/surf_land_bucket_mod.f90 |
| Date: | 2022-01-11 19:19:34 |
| Exec | Total | Coverage | |
|---|---|---|---|
| Lines: | 38 | 45 | 84.4% |
| Branches: | 38 | 56 | 67.9% |
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | ! | ||
| 2 | MODULE surf_land_bucket_mod | ||
| 3 | ! | ||
| 4 | ! Surface land bucket module | ||
| 5 | ! | ||
| 6 | ! This module is used when no external land model is choosen. | ||
| 7 | ! | ||
| 8 | IMPLICIT NONE | ||
| 9 | |||
| 10 | CONTAINS | ||
| 11 | |||
| 12 | 248640 | SUBROUTINE surf_land_bucket(itime, jour, knon, knindex, debut, dtime,& | |
| 13 | tsurf, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, & | ||
| 14 | spechum, petAcoef, peqAcoef, petBcoef, peqBcoef, pref, & | ||
| 15 | u1, v1, gustiness, rugoro, swnet, lwnet, & | ||
| 16 | snow, qsol, agesno, tsoil, & | ||
| 17 | 480 | qsurf, z0_new, alb1_new, alb2_new, evap, & | |
| 18 | fluxsens, fluxlat, tsurf_new, dflux_s, dflux_l) | ||
| 19 | |||
| 20 | USE limit_read_mod | ||
| 21 | USE surface_data | ||
| 22 | USE fonte_neige_mod | ||
| 23 | USE calcul_fluxs_mod | ||
| 24 | USE cpl_mod | ||
| 25 | USE dimphy | ||
| 26 | USE geometry_mod, ONLY: longitude,latitude | ||
| 27 | USE mod_grid_phy_lmdz | ||
| 28 | USE mod_phys_lmdz_para | ||
| 29 | USE indice_sol_mod | ||
| 30 | !**************************************************************************************** | ||
| 31 | ! Bucket calculations for surface. | ||
| 32 | ! | ||
| 33 | INCLUDE "clesphys.h" | ||
| 34 | INCLUDE "dimsoil.h" | ||
| 35 | INCLUDE "YOMCST.h" | ||
| 36 | |||
| 37 | ! Input variables | ||
| 38 | !**************************************************************************************** | ||
| 39 | INTEGER, INTENT(IN) :: itime, jour, knon | ||
| 40 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex | ||
| 41 | LOGICAL, INTENT(IN) :: debut | ||
| 42 | REAL, INTENT(IN) :: dtime | ||
| 43 | REAL, DIMENSION(klon), INTENT(IN) :: tsurf | ||
| 44 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay | ||
| 45 | REAL, DIMENSION(klon), INTENT(IN) :: tq_cdrag | ||
| 46 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow | ||
| 47 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum | ||
| 48 | REAL, DIMENSION(klon), INTENT(IN) :: petAcoef, peqAcoef | ||
| 49 | REAL, DIMENSION(klon), INTENT(IN) :: petBcoef, peqBcoef | ||
| 50 | REAL, DIMENSION(klon), INTENT(IN) :: pref | ||
| 51 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness | ||
| 52 | REAL, DIMENSION(klon), INTENT(IN) :: rugoro | ||
| 53 | REAL, DIMENSION(klon), INTENT(IN) :: swnet, lwnet | ||
| 54 | |||
| 55 | ! In/Output variables | ||
| 56 | !**************************************************************************************** | ||
| 57 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow, qsol | ||
| 58 | REAL, DIMENSION(klon), INTENT(INOUT) :: agesno | ||
| 59 | REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil | ||
| 60 | |||
| 61 | ! Output variables | ||
| 62 | !**************************************************************************************** | ||
| 63 | REAL, DIMENSION(klon), INTENT(OUT) :: qsurf | ||
| 64 | REAL, DIMENSION(klon), INTENT(OUT) :: z0_new | ||
| 65 | REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new, alb2_new | ||
| 66 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat | ||
| 67 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new | ||
| 68 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l | ||
| 69 | |||
| 70 | ! Local variables | ||
| 71 | !**************************************************************************************** | ||
| 72 | 960 | REAL, DIMENSION(klon) :: soilcap, soilflux | |
| 73 | 960 | REAL, DIMENSION(klon) :: cal, beta, dif_grnd, capsol | |
| 74 | 960 | REAL, DIMENSION(klon) :: alb_neig, alb_lim | |
| 75 | 960 | REAL, DIMENSION(klon) :: zfra | |
| 76 | 960 | REAL, DIMENSION(klon) :: radsol ! total net radiance at surface | |
| 77 | 960 | REAL, DIMENSION(klon) :: u0, v0, u1_lay, v1_lay | |
| 78 | 960 | REAL, DIMENSION(klon) :: dummy_riverflow,dummy_coastalflow | |
| 79 | INTEGER :: i | ||
| 80 | ! | ||
| 81 | !**************************************************************************************** | ||
| 82 | |||
| 83 | |||
| 84 | ! | ||
| 85 | !* Read from limit.nc : albedo(alb_lim), length of rugosity(z0_new) | ||
| 86 | ! | ||
| 87 | CALL limit_read_rug_alb(itime, dtime, jour,& | ||
| 88 | knon, knindex, & | ||
| 89 | 480 | z0_new, alb_lim) | |
| 90 | ! | ||
| 91 | !* Calcultaion of fluxes | ||
| 92 | ! | ||
| 93 | |||
| 94 | ! calculate total absorbed radiance at surface | ||
| 95 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | radsol(:) = 0.0 |
| 96 |
2/2✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
|
248160 | radsol(1:knon) = swnet(1:knon) + lwnet(1:knon) |
| 97 | |||
| 98 | ! calculate constants | ||
| 99 | 480 | CALL calbeta(dtime, is_ter, knon, snow, qsol, beta, capsol, dif_grnd) | |
| 100 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | if (type_veget=='betaclim') then |
| 101 | ✗ | CALL calbeta_clim(knon,jour,latitude(knindex(1:knon)),beta) | |
| 102 | endif | ||
| 103 | |||
| 104 | ! calculate temperature, heat capacity and conduction flux in soil | ||
| 105 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | IF (soil_model) THEN |
| 106 | CALL soil(dtime, is_ter, knon, snow, tsurf, qsol, & | ||
| 107 |
5/6✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 247680 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 247680 times.
✓ Branch 5 taken 480 times.
|
495840 | & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux) |
| 108 | |||
| 109 |
2/2✓ Branch 0 taken 480 times.
✓ Branch 1 taken 247680 times.
|
248160 | DO i=1, knon |
| 110 | 247680 | cal(i) = RCPD / soilcap(i) | |
| 111 | 248160 | radsol(i) = radsol(i) + soilflux(i) | |
| 112 | END DO | ||
| 113 | ELSE | ||
| 114 | ✗ | cal(:) = RCPD * capsol(:) | |
| 115 | ✗ | IF (klon_glo .EQ. 1) THEN | |
| 116 | ✗ | cal(:) = 0. | |
| 117 | ENDIF | ||
| 118 | ENDIF | ||
| 119 | |||
| 120 | ! Suppose zero surface speed | ||
| 121 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | u0(:)=0.0 |
| 122 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | v0(:)=0.0 |
| 123 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | u1_lay(:) = u1(:) - u0(:) |
| 124 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | v1_lay(:) = v1(:) - v0(:) |
| 125 | |||
| 126 | CALL calcul_fluxs(knon, is_ter, dtime, & | ||
| 127 | tsurf, p1lay, cal, beta, tq_cdrag, tq_cdrag, pref, & | ||
| 128 | precip_rain, precip_snow, snow, qsurf, & | ||
| 129 | radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, & | ||
| 130 | 1.,petAcoef, peqAcoef, petBcoef, peqBcoef, & | ||
| 131 | 480 | tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l) | |
| 132 | |||
| 133 | ! | ||
| 134 | !* Calculate snow height, run_off, age of snow | ||
| 135 | ! | ||
| 136 | CALL fonte_neige( knon, is_ter, knindex, dtime, & | ||
| 137 | tsurf, precip_rain, precip_snow, & | ||
| 138 | 480 | snow, qsol, tsurf_new, evap) | |
| 139 | ! | ||
| 140 | !* Calculate the age of snow | ||
| 141 | ! | ||
| 142 | 480 | CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:)) | |
| 143 | |||
| 144 |
4/4✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 183073 times.
✓ Branch 3 taken 64607 times.
|
248160 | WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0. |
| 145 | |||
| 146 |
2/2✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
|
248160 | DO i=1, knon |
| 147 | 247680 | zfra(i) = MAX(0.0,MIN(1.0, snow(i)/(snow(i)+10.0))) | |
| 148 | 248160 | alb_lim(i) = alb_neig(i) *zfra(i) + alb_lim(i)*(1.0-zfra(i)) | |
| 149 | END DO | ||
| 150 | |||
| 151 | ! | ||
| 152 | !* Return albedo : | ||
| 153 | ! alb1_new and alb2_new are here given the same values | ||
| 154 | ! | ||
| 155 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | alb1_new(:) = 0.0 |
| 156 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | alb2_new(:) = 0.0 |
| 157 |
2/2✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
|
248160 | alb1_new(1:knon) = alb_lim(1:knon) |
| 158 |
2/2✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
|
248160 | alb2_new(1:knon) = alb_lim(1:knon) |
| 159 | |||
| 160 | ! | ||
| 161 | !* Calculate the rugosity | ||
| 162 | ! | ||
| 163 |
2/2✓ Branch 0 taken 247680 times.
✓ Branch 1 taken 480 times.
|
248160 | DO i = 1, knon |
| 164 | 248160 | z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2)) | |
| 165 | END DO | ||
| 166 | |||
| 167 | !* Send to coupler | ||
| 168 | ! The run-off from river and coast are not calculated in the bucket modele. | ||
| 169 | ! For testing purpose of the coupled modele we put the run-off to zero. | ||
| 170 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | IF (type_ocean=='couple') THEN |
| 171 | ✗ | dummy_riverflow(:) = 0.0 | |
| 172 | ✗ | dummy_coastalflow(:) = 0.0 | |
| 173 | CALL cpl_send_land_fields(itime, knon, knindex, & | ||
| 174 | ✗ | dummy_riverflow, dummy_coastalflow) | |
| 175 | ENDIF | ||
| 176 | |||
| 177 | ! | ||
| 178 | !* End | ||
| 179 | ! | ||
| 180 | 480 | END SUBROUTINE surf_land_bucket | |
| 181 | ! | ||
| 182 | !**************************************************************************************** | ||
| 183 | ! | ||
| 184 | END MODULE surf_land_bucket_mod | ||
| 185 |