1 | ! |
---|
2 | ! $Id: change_srf_frac_mod.F90 4004 2021-11-09 16:53:19Z lguez $ |
---|
3 | ! |
---|
4 | MODULE change_srf_frac_mod |
---|
5 | |
---|
6 | IMPLICIT NONE |
---|
7 | |
---|
8 | CONTAINS |
---|
9 | ! |
---|
10 | ! Change Surface Fractions |
---|
11 | ! Author J Ghattas 2008 |
---|
12 | |
---|
13 | SUBROUTINE change_srf_frac(itime, dtime, jour, & |
---|
14 | pctsrf, evap, z0m, z0h, agesno, & |
---|
15 | alb_dir, alb_dif, tsurf, ustar, u10m, v10m, pbl_tke & |
---|
16 | #ifdef ISO |
---|
17 | ,xtevap & |
---|
18 | #endif |
---|
19 | & ) |
---|
20 | |
---|
21 | |
---|
22 | |
---|
23 | ! |
---|
24 | ! This subroutine is called from physiq.F at each timestep. |
---|
25 | ! 1- For each type of ocean (force, slab, couple) receive new fractions only if |
---|
26 | ! it's time to modify (is_modified=true). Otherwise nothing is done (is_modified=false). |
---|
27 | ! If received new fraction : |
---|
28 | ! 2- Tests and ajustements are done on the fractions |
---|
29 | ! 3- Initialize variables where a new fraction(new or melted ice) has appered, |
---|
30 | ! |
---|
31 | |
---|
32 | USE dimphy |
---|
33 | USE surface_data, ONLY : type_ocean,version_ocean |
---|
34 | USE limit_read_mod |
---|
35 | USE pbl_surface_mod, ONLY : pbl_surface_newfrac |
---|
36 | USE cpl_mod, ONLY : cpl_receive_frac |
---|
37 | USE ocean_slab_mod, ONLY : fsic, ocean_slab_frac |
---|
38 | USE indice_sol_mod |
---|
39 | USE print_control_mod, ONLY: lunout |
---|
40 | #ifdef ISO |
---|
41 | USE infotrac_phy, ONLY: ntraciso |
---|
42 | #endif |
---|
43 | |
---|
44 | INCLUDE "YOMCST.h" |
---|
45 | !albedo SB >>> |
---|
46 | include "clesphys.h" |
---|
47 | !albedo SB <<< |
---|
48 | |
---|
49 | |
---|
50 | |
---|
51 | ! Input arguments |
---|
52 | !**************************************************************************************** |
---|
53 | INTEGER, INTENT(IN) :: itime ! current time step |
---|
54 | INTEGER, INTENT(IN) :: jour ! day of the year |
---|
55 | REAL, INTENT(IN) :: dtime ! length of time step (s) |
---|
56 | |
---|
57 | ! In-Output arguments |
---|
58 | !**************************************************************************************** |
---|
59 | |
---|
60 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: pctsrf ! sub-surface fraction |
---|
61 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: evap, agesno ! sub-surface fraction |
---|
62 | |
---|
63 | REAL, DIMENSION(klon,nbsrf+1), INTENT(INOUT) :: z0m,z0h ! sub-surface fraction |
---|
64 | !albedo SB >>> |
---|
65 | REAL, DIMENSION(klon,nsw,nbsrf), INTENT(INOUT) :: alb_dir,alb_dif |
---|
66 | !albedo SB <<< |
---|
67 | #ifdef ISO |
---|
68 | REAL, DIMENSION(ntraciso,klon,nbsrf), INTENT(INOUT) :: xtevap |
---|
69 | #endif |
---|
70 | |
---|
71 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: tsurf |
---|
72 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: ustar |
---|
73 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: u10m |
---|
74 | REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: v10m |
---|
75 | !jyg< |
---|
76 | !! REAL, DIMENSION(klon,klev+1,nbsrf), INTENT(INOUT) :: pbl_tke |
---|
77 | REAL, DIMENSION(klon,klev+1,nbsrf+1), INTENT(INOUT) :: pbl_tke |
---|
78 | !>jyg |
---|
79 | |
---|
80 | ! Loccal variables |
---|
81 | !**************************************************************************************** |
---|
82 | INTEGER :: i, nsrf |
---|
83 | LOGICAL :: is_modified ! true if pctsrf is modified at this time step |
---|
84 | LOGICAL :: test_sum=.FALSE. |
---|
85 | LOGICAL, DIMENSION(klon,nbsrf) :: new_surf |
---|
86 | REAL, DIMENSION(klon,nbsrf) :: pctsrf_old ! fraction from previous time-step |
---|
87 | REAL :: tmpsum |
---|
88 | |
---|
89 | pctsrf_old(:,:) = pctsrf(:,:) |
---|
90 | !**************************************************************************************** |
---|
91 | ! 1) |
---|
92 | ! For each type of ocean (force, slab, couple) receive new fractions only if it's time |
---|
93 | ! to modify (is_modified=true). Otherwise nothing is done (is_modified=false). |
---|
94 | !**************************************************************************************** |
---|
95 | SELECT CASE (type_ocean) |
---|
96 | CASE ('force') |
---|
97 | ! Read fraction from limit.nc |
---|
98 | CALL limit_read_frac(itime, dtime, jour, pctsrf, is_modified) |
---|
99 | CASE ('slab') |
---|
100 | ! Get fraction from slab module |
---|
101 | CALL ocean_slab_frac(itime, dtime, jour, pctsrf, is_modified) |
---|
102 | CASE ('couple') |
---|
103 | ! Get fraction from the coupler |
---|
104 | CALL cpl_receive_frac(itime, dtime, pctsrf, is_modified) |
---|
105 | END SELECT |
---|
106 | |
---|
107 | |
---|
108 | !**************************************************************************************** |
---|
109 | ! 2) |
---|
110 | ! Tests and ajustements on the new fractions : |
---|
111 | ! - Put to zero fractions that are too small |
---|
112 | ! - Test total fraction sum is one for each grid point |
---|
113 | ! |
---|
114 | !**************************************************************************************** |
---|
115 | IF (is_modified) THEN |
---|
116 | |
---|
117 | ! Test and exit if a fraction is negative |
---|
118 | IF (MINVAL(pctsrf(:,:)) < 0.) THEN |
---|
119 | WRITE(lunout,*)'Warning! One or several fractions are negative, itime=',itime |
---|
120 | WRITE(lunout,*)'at point = ',MINLOC(pctsrf(:,:)) |
---|
121 | WRITE(lunout,*)'value = ',MINVAL(pctsrf(:,:)) |
---|
122 | CALL abort_physic('change_srf_frac','Negative fraction',1) |
---|
123 | END IF |
---|
124 | |
---|
125 | ! Optional test on the incoming fraction |
---|
126 | IF (test_sum) THEN |
---|
127 | DO i= 1, klon |
---|
128 | tmpsum = SUM(pctsrf(i,:)) |
---|
129 | IF (ABS(1. - tmpsum) > 0.05) CALL abort_physic('change_srf_frac','Total fraction not equal 1.',1) |
---|
130 | END DO |
---|
131 | END IF |
---|
132 | |
---|
133 | ! Test for too small fractions of the sum land+landice and ocean+sea-ice |
---|
134 | WHERE ((pctsrf(:,is_ter) + pctsrf(:,is_lic)) < 2*EPSFRA) |
---|
135 | pctsrf(:,is_ter) = 0. |
---|
136 | pctsrf(:,is_lic) = 0. |
---|
137 | END WHERE |
---|
138 | |
---|
139 | WHERE ((pctsrf(:,is_oce) + pctsrf(:,is_sic)) < 2*EPSFRA) |
---|
140 | pctsrf(:,is_oce) = 0. |
---|
141 | pctsrf(:,is_sic) = 0. |
---|
142 | END WHERE |
---|
143 | |
---|
144 | ! Normalize to force total fraction to be equal one |
---|
145 | DO i= 1, klon |
---|
146 | tmpsum = SUM(pctsrf(i,:)) |
---|
147 | DO nsrf = 1, nbsrf |
---|
148 | pctsrf(i,nsrf) = pctsrf(i,nsrf) / tmpsum |
---|
149 | END DO |
---|
150 | END DO |
---|
151 | |
---|
152 | ! Test for too small fractions at each sub-surface |
---|
153 | WHERE (pctsrf(:,is_ter) < EPSFRA) |
---|
154 | pctsrf(:,is_lic) = pctsrf(:,is_ter) + pctsrf(:,is_lic) |
---|
155 | pctsrf(:,is_ter) = 0. |
---|
156 | END WHERE |
---|
157 | |
---|
158 | WHERE (pctsrf(:,is_lic) < EPSFRA) |
---|
159 | pctsrf(:,is_ter) = pctsrf(:,is_ter) + pctsrf(:,is_lic) |
---|
160 | pctsrf(:,is_lic) = 0. |
---|
161 | END WHERE |
---|
162 | |
---|
163 | WHERE (pctsrf(:,is_oce) < EPSFRA) |
---|
164 | pctsrf(:,is_sic) = pctsrf(:,is_oce) + pctsrf(:,is_sic) |
---|
165 | pctsrf(:,is_oce) = 0. |
---|
166 | END WHERE |
---|
167 | |
---|
168 | WHERE (pctsrf(:,is_sic) < EPSFRA) |
---|
169 | pctsrf(:,is_oce) = pctsrf(:,is_oce) + pctsrf(:,is_sic) |
---|
170 | pctsrf(:,is_sic) = 0. |
---|
171 | END WHERE |
---|
172 | ! Send fractions back to slab ocean if needed |
---|
173 | IF (type_ocean == 'slab'.AND. version_ocean.NE.'sicINT') THEN |
---|
174 | WHERE (1.-zmasq(:)>EPSFRA) |
---|
175 | fsic(:)=pctsrf(:,is_sic)/(1.-zmasq(:)) |
---|
176 | END WHERE |
---|
177 | END IF |
---|
178 | |
---|
179 | !**************************************************************************************** |
---|
180 | ! 3) |
---|
181 | ! Initialize variables where a new fraction has appered, |
---|
182 | ! i.e. where new sea ice has been formed |
---|
183 | ! or where ice free ocean has appread in a grid cell |
---|
184 | ! |
---|
185 | !**************************************************************************************** |
---|
186 | |
---|
187 | CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, & |
---|
188 | evap, z0m, z0h, agesno, & |
---|
189 | tsurf, alb_dir,alb_dif, ustar, u10m, v10m, pbl_tke & |
---|
190 | #ifdef ISO |
---|
191 | ,xtevap & |
---|
192 | #endif |
---|
193 | & ) |
---|
194 | |
---|
195 | |
---|
196 | ELSE |
---|
197 | ! No modifcation should be done |
---|
198 | pctsrf(:,:) = pctsrf_old(:,:) |
---|
199 | |
---|
200 | END IF ! is_modified |
---|
201 | |
---|
202 | END SUBROUTINE change_srf_frac |
---|
203 | |
---|
204 | |
---|
205 | END MODULE change_srf_frac_mod |
---|